📄 25t.m
字号:
Rp=1;Rs=25;wp1=0.2*pi;ws1=0.4*pi;T=0.001;
wp=(2/T)*tan(wp1/2);ws=(2/T)*tan(ws1/2);
[n,wn]=cheb2ord(wp,ws,Rp,Rs,'s') %从此处可以计算阶数n
[b,a]=cheby2(n,Rs,wn,'low','s') %由b,a的值可以得到系统函数
[bz,az]=bilinear(b,a,1/T);
[b0,B,A]=dir2cas(bz,az)
[db,mag,pha,grd,w]=freqz_m(bz,az);
subplot(2,1,1);plot(w/pi,db);
axis([0,1,-65,10]);
xlabel('');
ylabel('');
title('幅频相应')
subplot(2,1,2);
plot(w/pi,pha);
xlabel('频率(单位:pi)');
ylabel('相位');title('相频相应');
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -