数，据此可以计算出一个记录中的最少点数，执行函数fft，可以将x（
）经过傅里叶变换求出Xk。5．用线性相位低通数字滤波器对x
进行滤波，在设计低通数字滤波器时已经求出了单位冲击响应序列h（
），那么将采样得到的序列x（
）与h（
）进行卷积，就可以得到滤波器的输出y1
，经过fft变换后，可以对y1
进行频谱分析。6．用Chebyshev低通数字滤波器对x
进行滤波时，由于设计时能够将切比雪夫滤波器系统函数（z域，双线性变换变换的表达式）的零极点均已求出，那么就可以调用filter函数求出y2
，最后将得到滤波器的输出进行fft变换，对y2
进行频谱分析，按要求画出相应的子图。实验结果：实验结果：数字线性相位FIR低通滤波器
f数字ChebyshevI型低通滤波器
f以fs10kHz对xt进行取样
用FFT对x
进行谱分析
f线性相位低通数字滤波器对x
进行滤波
Chebyshev低通数字滤波器对x
进行滤波
程序代码：程序代码
主程序：主程序：clear任务一：窗口函数的设计wp026piws041pitr_widthwswp带宽
fN1ceil66pitr_width1阶数
0N11wcwswp2理想低通的截止频率hdideal_lpwcN1理想低通的冲击响应w_hamhammi
gN1hhdw_hamFIR滤波器的冲击响应dbmagphagrdwfreqz_mh1delta_w2pi1000作图figure1subplot221stem
hraxis0460104title单位冲击响应subplot222plotwpipharaxis0144title相频响应subplot223plotwpimagraxis008011title幅频响应subplot224plotwpidbraxis00810010title幅频响应dBxlabel以pi为单位的频率ylabel对数幅度dB任务二切比雪夫滤波器的设计Rp075As50T1Fs1Tomegap2Tta
wp2omegas2Tta
ws2
fcsdschebomegapomegasRpAs调用函数求解拉普拉斯变换的分子与分母多项式的系数数组babili
earcsdsFs做双线性变换dbmagphagrdwfreqz_mba频域上采样并记录幅频和相频响应作图figure2subplot311plotwpimagtitle幅度响应gridaxis006011ylabel幅度xlabel以pi为单位的频率subplot312plotwpidbtitle幅频响应dBgridaxis0066010xlabel以pi为单位的频率ylabel对数幅度dBsubplot313plotwpiphatitle相频响应gridaxis00644xlabel以pi为单位的频率ylabel相位任务三：fs10000f1120f24980f33250f41150Nq512
q0Nq1xg2cos2pif1
qfs5si
2pif2
qfs8cos2pif3
qfs75cos2pif4
qfs模拟信号采样的表达式作图figure3subplot211fplotxat0511ctitr