文档介绍：--------------------------校验:_____________-----------------------日期:_____________信号处理实验三报告实验三抽样抽样引起的混叠【实验内容】f0=500HZ,fs分别取100HZ,1kHZ,10kHZ,绘出x[n]及其DTFT。【程序】%---文件%问题实验内容closeall;clearall;f0=500;%初相n=0:20;fs1=100;%采样频率100HZt1=n/fs1;x1=sin(2*pi*f0*t1+B);[X1,W1]=dtft(x1,1000);%对采样信号做DTFTfs2=1000;%采样频率1kHZt2=n/fs2;x2=sin(2*pi*f0*t2+B);[X2,W2]=dtft(x2,1000);%对采样信号做DTFTfs3=10000;%采样频率10kHZt3=n/fs3;x3=sin(2*pi*f0*t3+B);[X3,W3]=dtft(x3,1000);%对采样信号做DTFThf=figure;%绘制x1时域采样结果subplot(2,3,1);stem(n,x1);title('fs=100HZ时域采样结果');xlabel('n');ylabel('x1[n]');subplot(2,3,4);%绘制x1频域解图plot(W1,X1);gridon;title('fs=100HZ采样信号的DTFT');xlabel('W1');ylabel('X1');%绘制x2时域采样结果subplot(2,3,2);stem(n,x2);title('fs=1kHZ时域采样结果');xlabel('n');ylabel('x2[n]');subplot(2,3,5);%绘制x2频域解图plot(W2,X2);gridon;title('fs=1kHZ采样信号的DTFT');xlabel('W2');ylabel('X2');%绘制x3时域采样结果subplot(2,3,3);stem(n,x3);title('fs=10kHZ时域采样结果');xlabel('n');ylabel('x3[n]');%绘制x3频域解图subplot(2,3,6);plot(W3,X3);gridon;title('fs=10kHZ采样信号的DTFT');xlabel('W3');ylabel('X3');saveas(【结果及分析】如图所示,前三幅图分别为欠采样、临界采样和完全采样的结果,后三幅图为他们分别对应的DTFT。图正弦信号的不同频率采样及DTFT【实验内容】求出并绘制其傅里叶变换xa(jΩ)(n),x2(n);画出他们的DTFT并与xa(jΩ)。【程序】%---文件%问题a%傅里叶变换结果X(W)=2000/(1000000+W^2)clearall;%模拟信号;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)];%频率介于-Wmax和Wamx之间Xa=[fliplr(Xa),Xa(2:501)];%Xa介于-Wmax和Wamx之间hf=figure;subplot(2,2,1);plot(t*1000,xa);xlabel('t(ms)');ylabel('xa(t)');gridon;title('模拟信号');subplot(2,2,2);plot(W/(2*pi*1000),Xa*1000);xlabel('f(kHZ)');ylabel('Xa(jw)');gridon;title('连续时间傅里叶变换');;n=-25:25;x1=exp(-1000*abs(n*Ts1));%离散信号K=500;k=0:K;w=pi*k/K;X1=x1*exp(-j*n'*w);X1=real(X1);%离散时间傅里叶变换w=[-fliplr(w),w(2:K+1)];X1=[fliplr(X1),X1(2:K+1)];subplot(2,2,3);plot(w/pi,X1);xlabel('w(rad/s)');ylabel('X1(w)');gridon;title('Ts=5000HZ');;n=-25:25;x2=exp(-1000*abs(n*Ts2));%离散信号K=500;k=0:K;w=pi*k/K;X2=x2*exp(-j*n'*w);X2=real(X2);%离散时间傅里叶变换w=[-fliplr(w),w(2:K+1)];X2=[fliplr(X2