信号与系统11 - 大学课件.docxVIP

设计脉冲响应不变法、双线性变换法设计IIR数字巴特沃斯低通数字滤波器。例如：通带截止频率0.2*pi

设计脉冲响应不变法、双线性变换法设计IIR数字巴特沃斯低通数字滤波器。例如：通带截止频率0.2*pi

1(N,rp,wp,high);freqz(Bz,Az)六、※FIR滤波器的设计、※FIR滤

(wp,ws,rp,rs,s);[b,a]=butter(n,ws,s);[bz,az]

2*pi/2);ws=2/t*tan(0.3*pi/2);rp=1;rs=10;[N,wc]

信号处理与系统分析课程设计内容

一、※三个系统的系统函数分别为h1(

h2(s)=(4s+5)/[s(5s+1)];h3(s)=(4s+5)/[5s2+5s+5)]用simulink来仿真三个系统

的阶跃响应，并分析系统的稳定性。

1、a=[5];

b=[51];

[zpk]=tf2zp(a,b)

z=Emptymatrix:0-by-1

p=-0.2000

k=1

极点p=-0.2，在s域的左半开平面，所以系统稳定！

2、a=[45];

b=[510];

通数字滤波器，要求通带截止频率ωp=0.8πrad,

通数字滤波器，要求通带截止频率ωp=0.8πrad,通带最大衰减αp=3db,阻带截止频率ωs=0.

k=0.8000极点p=-0.2在s域的左半开平面，但p=0不在左半开平面，所以系统不稳定！a=

1,2)plot(w,ang

];m=[1,0];dev=[0.02,0.01];[N,f0,m0,w]=firpmord

[zpk]=tf2zp(a,b)

z=-1.2500

p=0

-.2000

k=0.8000

极点p=-0.2在s域的左半开平面，但p=0不在左半开平面，所以系统不稳定！

3、a=[45];

b=[555];

[zpk]=tf2zp(a,b)

z=-1.2500

p=-0.5000+0.8660i

-0.5000-0.8660i

k=0.8000

由极点p可知，-0.5000+0.8660i和0.5000-0.8660i并不都在s域的左半开平面，所以该系统不稳定！

];m=[1,0];dev=[0.02,0.01];[N,f0,m0,w]=firpmordpi);ylabel(相位(rad));Wp=3/8;dev1=0.02;Ws=0.5;(2

];m=[1,0];dev=[0.02,0.01];[N,f0,m0,w]=firpmord

pi);ylabel(相位(rad));Wp=3/8;dev1=0.02;Ws=0.5;

(2πfn/f)+sin(2πfn/f)，分析其频谱。f1=2;f2=2.05;fs=10;

如图所示，要求设计一数字滤波器，将干扰滤除掉，指标是允许|s(f)|在0≤f≤15Khz频率范围中

1s2s

三、频谱分析

设一序列中含有两种频率成分，f1=2hz，f2=2.05hz，采样频率为fs=10Hz，即x(n)=sin(2

πfn/f)+sin(2πfn/f)，分析其频谱。

f1=2;f2=2.05;

fs=10;Ts=1/fs;Tp=20;N=fs*Tp;

n=[0:N-1];

xn=sin(2*pi*f1*Ts*n)+sin(2*pi*f2*Ts*n);

Xk=fft(xn,N);

stem(n,abs(Xk),*)

xlabel(k);ylabel(|X|)

信号：x(t)=s(t)+u(t),其中u(t)

信号：x(t)=s(t)+u(t),其中u(t)是干扰信号，s(t)+u(t)的频谱不混叠，其幅度谱

5πrad，阻带最小衰减αs=18dbwp=0.8*pi;rp=3;ws=0.5*pi;rs=1

khzclear;fc=15000;fs=20000;Fs=80000;f=[fc,fs

bs(Xk),*)xlabel(k);ylabel(|X|)四、※IIR滤波器的

四、※IIR滤波器的设计

脉冲响应不变法、双线性变换法设计IIR数字巴特沃斯低通数字滤波器。例如：通带截止频率0.2*pi,阻带截止频率0.3*pi,通带波动1dB；在频率0.3π到π之间的阻带衰减大于

10dB。

t=1;

wp=0.2*pi/t;ws=0.3*pi/t;rp=1;rs=10;

[n,wc]=buttord(wp,ws,rp,rs,s);

[b,a]=butter(n,ws,s);

[bz,az]=imp