matlab信号程序,信号matlab程序

中南大学赵亚湘信号处理课程设计matlab

⑴请设计实验对采样定理进行验证。

① 请画出其傅立叶变换图；.

② 请画出信号重建图(可以采用不同的内插方法)

思路：采样定理

对连续信号进行等间隔采样形成采样信号，采样信号的频率是原连续信号的频率以采样频率为周期进行周期延拓

设连续信号xa(t)属于带限信号，最高截止频率为Ω

发生频谱混叠。若Ωc，如果采样角频率Ωs>=Ωc则不会s<=Ωc，则发生频谱混叠。

用内插函数恢复出原模拟信号

取信号xa(t)=e 1000|t|，并求出其傅利叶变换，原信号最高频率2KHZ，令采样频率为5000次每秒

%原模拟信号

Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t));

%连续时间傅利叶变换

Wmax=2*pi*2000;K=500;k=0:1:K;W=k*Wmax/K;

Xa=xa*exp(-j*t'*W)*Dt;Xa=real(Xa);

W=[-fliplr(W),W(2:501)];

Xa=[fliplr(Xa),Xa(2:501)];

%离散时间信号

Ts=0.0002;Fs=1/Ts;n=-25:1:25;nTs=n*Ts;x=exp(-1000*abs(n*Ts));

%离散时间信号DTFT

w=k*pi/K;

X=x*exp(-j*n'*w);X=real(X);

w=[-fliplr(w),w(2:K+1)];

X=[fliplr(X),X(2:K+1)];

%信号重构

x1=x*sinc(Fs*(ones(length(n),1)*t-nTs'*ones(1,length(t))));

%画图

subplot(3,2,1);plot(t*1000,xa);xlabel('t毫秒');ylabel('xa(t)');title('模拟信号');

subplot(3,2,2);plot(W/(2*pi*1000),Xa*1000);

xlabel('频率(kHZ)');ylabel('Xa(jw)*1000');title('原信号的傅里叶变换');

subplot(3,2,3);stem(n*Ts*1000,x);

xlabel('n');ylabel('x(n)');title('采样信号');

subplot(3,2,4);plot(w/pi,X);

xlabel('以pi为单位的频率');ylabel('X(e^jw)');title('离散时间傅里叶变换');

subplot(3,1,3);plot(t*1000,x1);title('重建出来的信号');