数字信号处理作业.doc

数字信号处理作业 1 DFT各种现象的图形 2 混叠现象参数设置的疑问 0800030061 闫 严 模拟原始信号求N点DFT的幅度谱和相位谱： 有混叠现象的信号求N点DFT的幅度谱和相位谱： 有频谱泄漏现象的信号求N点DFT的时域图和频域图 n=[0:1:99]; %采样点数即矩形窗长度 xn=cos(3*pi/20*n); %原信号时域表达式 %生成信号序列 Xk=fft(xn); %对加矩形窗后的信号做N点的傅利叶变换 kx=[0:1:length(Xk)-1]; %求出频域图的X轴 subplot(2,1,1); %在两行一列的第一行图中画信号时域图 stem(n,xn); %绘制源曲线信号时域图 title(time domain signal); %做标题 %绘制源曲线信号 subplot(2,1,2); %在两行一列的第二行图中画信号频域图 stem(kx,abs(Xk)); %绘制源曲线信号相位谱图 title(frequence domain signal(n=100)); %做标题 加长矩形长度改善频漏现象的信号求N点DFT的时域图和频域图： n=[0:1:149]; %采样点数即矩形窗长度(比上面的长) xn=cos(3*pi/20*n); %原信号时域表达式 %生成信号序列 Xk=fft(xn); %对加矩形窗后的信号做N点的傅利叶变换 kx=[0:1:length(Xk)-1]; %求出频域图的X轴 subplot(2,1,1); %在两行一列的第一行图中画信号时域图 stem(n,xn); %绘制源曲线信号时域图 title(time domain signal); %做标题 %绘制源曲线信号 subplot(2,1,2); %在两行一列的第二行图中画信号频域图 stem(kx,abs(Xk)); %绘制源曲线信号相位谱图 title(frequence domain signal(n=150)); %做标题 栅栏现象和信号末尾补零改善的信号求N点DFT图： n=[0:1:29]; %补零之前的采样矢量 xn=cos(3*pi/20*n); %原信号时域表达式 %生成信号序列 gn=[xn,zeros(1,30)]; %对信号进行末尾补零防止栅栏现象 N=length(n); %补零之前的采样点数 M=N+30;