2024年西南交通大学数字信号处理实验2

数字信号处理 MATLAB 上机实验第三章3-23 序列 x(n)={1, 2, 3, 3, 2, 1} 1〕 求出 x(n)的傅里叶变换 X(ejω)， 画出幅频特性和相频特性曲线(提示： 用 1024 点 FFT 近似 X(ejω))；  2〕 计算 x(n)的 N〔N≥6〕点离散傅里叶变换 X(k)， 画出幅频特性和相频特性曲线； 3〕 将 X(ejω)和 X(k)的幅频特性和相频特性曲线分别画在同一幅图中， 验证 X(k)是 X(ejω)的等间隔采样， 采样间隔为 2π/N； 4〕 计算 X(k)的 N 点 IDFT， 验证 DFT 和 IDFT 的惟一性。实验分析 〔1〕 题用 1024 点 DFT 近似 x(n)的傅里叶变换。 〔2〕 题用 36 点 DFT。 〔4〕 题求傅里叶反变换验证 IDFT 的惟一性。实验代码与截图1 到 3 问xn=[1 2 3 3 2 1];Xen=fft(xn,1024);n1=0:length(Xen)-1;amp = abs(Xen);phi = angle(Xen);Xkn=fft(xn,36);n2=0:length(Xkn)-1;amp2 = abs(Xkn);phi2 = angle(Xkn);subplot(221);plot(n1,amp)title('Xejw 幅频特性');xlabel('n');ylabel('Amp')subplot(222);plot(n1,phi)title('Xejw 相频特性');xlabel('n');ylabel('Phi')subplot(223);stem(n2,amp2,'.')title('Xk 幅频特性');xlabel('n');ylabel('Amp')subplot(224);stem(n2,phi2,'.')title('Xk 相频特性');xlabel('n');ylabel('Phi')截图如下第 4 问xn=[1 2 3 3 2 1];Xkn2=fft(xn,6);x6n=ifft(Xkn2);n2=0:length(x6n)-1;subplot(2,1,2);stem(n2,x6n,'.');title('X6k 傅里叶逆变换');xlabel('n');ylabel('x6n');Xkn1=fft(xn,16);x16n=ifft(Xkn1);n1=0:length(x16n)-1;subplot(2,1,1);stem(n1,x16n,'.');title('X16k 傅里叶逆变换');xlabel('n');ylabel('x16n')截图为3-25 序列 h(n)=R6(n), x(n)=nR8(n)。  1〕 计算 yc(n)=h(n) 8 x(n)；  2〕 计算 yc(n)=h(n) 16 x(n)和 y(n)=h(n)*x(n)； 3〕 画出 h(n)、 x(n)、 yc(n)和 y(n)的波形图， 观察总结循环卷积与线性卷积的关系。实验分析循环卷积为线性卷积的周期延拓序列的主值序列； 当循环卷积区间长度大于等于线性卷积序列长度时，二者相等。实验代码与截图hn=[1 1 1 1 1 1];xn=[0 1 2 3 4 5 6 7]; %用 DFT 计算 8 点循环卷积 yc8n:H8k=fft(hn,8); %计算 h(n)的 8 点 DFTX8k=fft(xn,8); %计算 x(n)的 8 点 DFTYc8k=H8k.*X8k;yc8n=ifft(Yc8k,8);%用 DFT 计算 16 点循...