FIR滤波器的设计作业.docVIP

生医《数字信号处理》期末大作业 一．设计FIR带阻滤波器 1.1数据的产生 1.用randn()函数产生均值为零，功率为0.1（即方差为0.1）的高斯分布的白噪声， clear; u1=randn(1,500000); u1=u1*sqrt(0.1); power_u1=var(u1); subplot(2,1,1); plot(u1(1:200)); grid on; xlabel(n); ylabel(x(n)); u2=randn(1,500000); u2=u2*sqrt(0.1); power_u2=var(u2); subplot(2,1,2); plot(u2(1:200)); grid on; xlabel(n); ylabel(x(n)); 和通过由5个FIR子系统（）级联而成的FIR系统，得到输出，令，在上加四个已知幅度（）和归一化频率（）的复正弦信号，得到已知的试验信号。 分析可得：由已知的5个FIR子系统的，可知系统的，从而得到总的FIR系统的系统函数，最终的试验信号，求幅频响应则调用freqz（）函数。 clear; b1=[1,1.98,0.9801]; b2=[1,-1.98,0.9801]; h1=conv(b1,b2); b3=[1,-1.8418,0.9801]; h2=conv(h1,b3); b4=[1,-1.5,0.9801]; h3=co