hc824.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第八章例8.2.4程序hc824
% 切贝雪夫II模拟滤波器设计
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
OmegaP=2*pi*5000; OmegaS=2*pi*12000;     % 输入设计指标
Rp=1; As=30;
[N, OmegaC] = cheb2ord(OmegaP, OmegaS, Rp, As, 's')   % 调用切贝雪夫II设计函数
[z0,p0,k0] = cheb2ap(N, As)          % 调用切贝雪夫II原型函数,得到左半平面零极点
b0 = k0*real(poly(z0))			       % 由零点计算分子系数向量
a0 = real(poly(p0))               % 由极点计算分母系数向量
[Ha,w0]=freqs(b0*k0,a0);      % 求归一化频率响应
subplot(1,2,1), 
set(gcf,'color','w')              % 置图形背景色为白
semilogx(w0*OmegaC,20*log10(abs(Ha)/max(abs(Ha)))) % 画出归一化分贝频率特性
grid on                     % 画出非归一化对数频率特性
xlabel('\Omega'),ylabel('dBHa(\Omega)')
subplot(1,2,2), 
plot(p0*OmegaC,'x'),hold on,plot(z0*OmegaC,'o')     % 画出极点分布图
xlabel('\sigma'),ylabel('j\Omega')
axis square, axis equal, grid on        % 使x,y等比例，以便看出极点的椭圆分布
Nm=length(a0)-length(b0);
[sos,G]=tf2sos([zeros(1,Nm),b0],a0)  % 求级联形式的归一化传递函数