Question

跪求下面2题如何用MATLAB编程！！！小女子不胜感激！！！

1.已知H_a (s)=(s+1)/(s^2+5s+6),分别用脉冲响应不变法,双线性变换法,利用impinvar及bilinear函数求系统函数H(z),选择T=1s.
2.设计一个切比雪夫I型数字带通滤波器，要求：带通范围100~250Hz，阻带上限300Hz,下限50Hz，通带呢内波纹小于3dB，阻带为-30dB。（设计采样频率为1000Hz）

Answer 1

b=[1 1];
a=[1 5 6];
T=1;
[bz,az]=impinvar(b,a,T) %脉冲响应不变法
>>bz =

1.0000 -0.2209

az =

1.0000 -0.1851 0.0067
即H(z)=[1.0000-0.2209z^(-1)]/[1.0000-0.1851z^(-1)+0.0067z^(-2)]
双线性法：
b=[1 1];
a=[1 5 6];
T=1;Fs=1/T;
[bz,az]=bilinear(b,a,Fs) %双线性

切比雪夫：
%数字带通滤波器程序
fp= [100,250];fs=[50,300]; %模拟通带、阻带频率
rp=3;rs=30; %通带波动、阻带衰减
fo=1000; %采样频率

%频率指标变换
wp=2*pi*fp/fo; %数字通带频率
ws=2*pi*fs/fo; %数字阻带频率

%切比雪夫1型数字带通滤波器直接设计
[n,wn]=cheb1ord(wp/pi,ws/pi,rp,rs);
[b,a]=cheby1(n,rp,wn); %数字带通滤波器系数

%切比雪夫1型数字带通滤波器频率响应
[h,w]=freqz(b,a);
mag=abs(h);pha=angle(h);
subplot(1,2,1);plot(w,mag); %幅频曲线
subplot(1,2,2);plot(w,pha); %相频曲线