基于MATLAB的数字信号处理的课程设计报告.docVIP

课程设计掌握《数字信号处理基础》课程的基本理论；掌握应用MATLAB进行信号处理的程序设计： 低通数字滤波器，，， 要求： 选择合适的窗函数，数字包括个子图，滤波器的。用法，设计数字低通滤波器包括个子图，滤波器的。模拟信号 其中，，，，取采样频率。 要求： 以对进行取样，。用一个图形窗口，包括两个子图，以及()的波形； 用FFT对进行谱分析，要求频率分辨率不超过5Hz。求出一个记录长度中的最少点数，并用一个图形窗口，包括两个子图，以及的幅值数字进行滤波，求出滤波器的输出，并用FFT对进行谱分析，要求频率分辨率不超过5Hz。求出一个记录长度中的最少点数，并用一个图形窗口，包括个子图，（）、、（）和的幅值； 用要求2中设计的Chebyshev低通数字进行滤波，求出滤波器的输出，并用FFT对进行谱分析，要求频率分辨率不超过5Hz。求出一个记录长度中的最少点数，并用一个图形窗口，包括个子图，（）、、（）和的幅值。 备注： (1) 要求编写一个主程序，完成上述问题。 (2) 要求各结果图均标出图题以及横纵轴的名称。 (3) 要求给主要语句加上标注。 设计思路： 1． 在任务一中，最小阻带衰减As=50dB，则就此可以确定所要选择的窗函数为哈明窗。 2. 对于任务二，根据通带频率wp，阻带频率ws,通带衰减Rp，阻带衰减As，计算出所设计的切比雪夫滤波器阶数，极点等重要参数。 对连续时间信号进行采样，得到，根据要求绘制出子图。 频率分辨率等于采样频率除以采样点数，据此可以计算出一个记录中的最少点数，执行函数fft，可以将x（n）经过傅里叶变换求出X(k)。 用线性相位低通数字进行滤波，在设计低通数字，经过fft变换后，可以对进行频谱分析。 用Chebyshev低通数字进行滤波时，由于设计时能够将切比雪夫滤波器系统函数（z域，双线性变换变换的表达式）的零极点均已求出，那么就可以调用filter函数求出，最后将得到滤波器的输出进行fft变换，对进行频谱分析，按要求画出相应的子图。 实验结果： 数字 数字低通滤波器 以对进行取样进行谱分析 线性相位低通数字进行滤波 Chebyshev低通数字进行滤波 程序代码： 主程序： clear; %任务一： %窗口函数的设计 wp=0.26*pi; ws=0.41*pi; tr_width=ws-wp;%带宽 N1=ceil(6.6*pi/tr_width)+1;%阶数 n=0:N1-1; wc=(ws+wp)/2;%理想低通的截止频率 hd=ideal_lp(wc,N1);%理想低通的冲击响应 w_ham=(hamming(N1)); h=hd.*w_ham;%FIR滤波器的冲击响应 [db,mag,pha,grd,w]=freqz_m(h,[1]); delta_w=2*pi/1000; %作图 figure(1); subplot(2,2,1); stem(n,h,r); axis([0 46 -0.1 0.4]); title(单位冲击响应); subplot(2,2,2); plot(w/pi,pha,r); axis([0 1 -4 4]); title(相频响应); subplot(2,2,3); plot(w/pi,mag,r); axis([0 0.8 0 1.1]); title(幅频响应); subplot(2,2,4); plot(w/pi,db,r); axis([0 0.8 -100 10]); title(幅频响应(dB));xlabel(以\pi为单位的频率);ylabel(对数幅度/dB); %任务二; %切比雪夫滤波器的设计 Rp=0.75; As=50; T=1; Fs=1/T; omegap=(2/T)*tan(wp/2); omegas=(2/T)*tan(ws/2); [cs,ds]=cheb(omegap,omegas,Rp,As);%调用函数求解拉普拉斯变换的分子与分母多项式的系数数组 [b,a]=bilinear(cs,ds,Fs);%做双线性变换 [db,mag,pha,grd,w]=freqz_m(b,a);%频域上采样并记录幅频和相频响应 %作图 figure(2); subplot(3,1,1); plot(w/pi,mag); title(幅度响应);grid; axis([0 0.6 0 1.1]);ylabel(幅度);xlabel(以\pi为单位的频率); subplot(3,1,2); plot(w/pi,db);title(幅频响应(dB));grid; axis([0 0.6 -60 10]);xlabel(以\pi为单位的频率);ylabel(对数幅度/dB); subplot(3,1,3); plot(w/pi,pha);ti