采样与模拟信号重建MATLAB实例%关于连续函数求傅⾥叶变换
%⽤有限长序列近似原函数（利⽤e^-5 约为 0）
%从⽽确定出序列间隔T的范围，接下来要确定T的步进量
%要求：步进T<<；采样间隔
%先求出傅⾥叶变换换后幅值在什么（设为f）频率下趋向0
%步进T取⼀个值<<1/f\
%analog signal
dt = 0.00005;                                                      %时间步进量
t = -0.005:dt:0.005;                                              %时间范围
xa = exp(-1000*abs(t));                                        %求函数值
%contunites_time fourier transform
Wmax = 2*pi* 2000;                                            %观察的最⾼频率
K = 500;                                                              %500份频率值
k = 0: 1: K;
W = k*Wmax/K;
Xa = xa * exp(-1i * t'*W)*dt;
Xa = real(Xa);
W = [-fliplr(W),W(2:501)];              %Flip matrix left to right  倒置左右逐个交换
Xa = [fliplr(Xa),Xa(2:501)];            %xa over -Xa to Xa 合并矩阵
subplot(2,1,1);
plot(t*1000,xa);grid
title('analog signal');
xlabel('t  (ms)');
ylabel('xa(t)');
subplot(2,1,2);
plot(W/(2*pi*1000),Xa*1000);grid
title('continues_time fourier transform');
xlabel('f  (Khz)');
ylabel('Xa(jw) * 1000');
%analog signal
dt = 0.00005;
t = -0.005: dt: 0.005;
xa = exp(-1000*abs(t));
%discrete_time signal
ts = 0.0002;
n = -25:1:25;
x = exp(-1000*abs(n*ts));
%discrete-fourier transform
K = 500;
k = 0:1:K;
w = 2*pi*k/K;
X = x*exp(-1i *n'*w);
X= real(X);
w = [-fliplr(w),w(2:K+1)];
X = [fliplr(X),X(2:K+1)];                                      %要对应-w所求的值
subplot(2,1,1);
%hold on retains the current plot and certain axes properties so that
%subsequent graphing commands add to the existing graph. If no
%current axes exist before you call hold on, MATLAB creates new
%axes and retains the default properties. However, some axes
%properties change to accommodate additional graphics objects.
%For example, the axes' limits increase when the data requires them
%to do so. hold on sets the NextPlot property of the current figure and axes to add.
plot(t*1000,xa);
title('discrete signal');
xlabel('t in msec');
ylabel('x(n)');
hold on
%stem(X,Y)    stem  Plot discrete sequence data
stem(n*ts*1000,x);
%gtext  Mouse placement of text in 2-D view
%gtext displays a text string in the current figure window after you select
%a location with the mouse.
%gtext('string')matlab求傅里叶变换
%hold      Retain current graph in figure
%hold off resets axes properties to their defaults before drawing new
%plots. hold off is the default. hold off sets the NextPlot property of the current axes to replace. gtext('ts = 0.2msec'); hold off
subplot(2,1,2);
plot(w/pi, X);
title('discrete_time fourier transform');
xlabel('frequence in pi unit');
ylabel('X(w)');

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。