文档介绍：【实验内容】f0=500HZ,fs分别取100HZ,1kHZ,10kHZ,绘出x[n]及其DTFT。【程序】%------%;clearall;f0=500;B=*pi;%初相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(hf,'');【结果及分析】,前三幅图分别为欠采样、临界采样和完全采样的结果,后三幅图为他们分别对应的DTFT。【实验内容】求出并绘制其傅里叶变换xa(jΩ)(n),x2(n);画出他们的DTFT并与xa(jΩ)。【程序】%------%问题a%傅里叶变换结果X(W)=2000/(1000000+W^2)clearall;Dt=;%模拟信号t=-:Dt:;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=[fl