信号处理与系统分析课程设计..docVIP

信号处理与系统分析课程设计 一 二 %三 %频谱分析 f1=2;f2=2.05;%模拟信号的频率分量 fs=10;Ts=1/fs;%取样频率和周期 Tp=100;N=fs*Tp;%信号观测时间和取样点数N n=[0:N-1]; xn=sin(2*pi*f1/fs*n)+sin(2*pi*f2/fs*n);%数字信号 Xk=fft(xn,N);%快速FFT变换 stem(n,abs(Xk),.) title(频谱分析) xlabel(k);ylabel % 四 IIR滤波器的设计 %脉冲响应不变法设计IIR数字巴特沃斯低通滤波器 T=1;%T=1s wp=0.2*pi;ws=0.3*pi;ap=1;as=10;%数字滤波器指标 Wp=wp/T;Ws=ws/T;ap=1;as=10;%T=1时的模拟滤波器指标 [N,Wc]=buttord(Wp,Ws,ap,as,s);%计算相应的模拟滤波器阶数N和3dB截止频率 [B,A]=butter(N,Wc,s);%计算响应的模拟滤波器系统函数 [Bz,Az]=impinvar(B,A,1/T);%用脉冲响应不变法将模拟滤波器转换成数字滤波器 [H,w]=freqz(Bz,Az); plot(w/pi,abs(H)) title(脉冲响应不变法设计IIR数字巴特沃斯低通滤波器) xlabel(w/pi)