hc825.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第八章例8.2.5程序hc825
% 椭圆型模拟滤波器设计
% 电子工业出版社出版  陈怀琛编著 2004年8月
%
OmegaP=2*pi*5000; OmegaS=2*pi*12000;     % 输入设计指标
Rp=1; As=30;
[N, OmegaC] = ellipord(OmegaP, OmegaS, Rp, As, 's')   % 调用椭圆设计函数
[z0,p0,k0] = ellipap(N,Rp,As)          % 调用椭圆原型函数,得到左半平面零极点
b0 = real(poly(z0))			       % 由零点计算分子系数向量
a0 = real(poly(p0))               % 由极点计算分母系数向量
p = p0*OmegaC; z = z0*OmegaC;   		% 将零极点乘以Omegac，得到非归一化零极点
Nm=length(a0)-length(b0);		% 求极点个数与零点个数之差Nm
[Ha,w0]=freqs(k0*b0,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(p,'x'),hold on,plot(z,'o')     % 画出极点分布图
axis square, axis equal, grid on        % 使x,y等比例，以便看出极点的椭圆分布
xlabel('\sigma'),ylabel('j\Omega')
[sos,G]=tf2sos([zeros(1,Nm),b0],a0)  % 求级联形式的归一化传递函数