hc821.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第八章例8.2.1程序hc821
% 离散信号的基本序列
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=input('N= ');                 % 输入滤波器参数N和OmegaC
OmegaC=input('OmegaC= ');
[z0,p0,k0] = buttap(N)          % 调用buttap函数,得到系统的左半平面零极点
p = p0*OmegaC; z = z0*OmegaC;   % 将零极点乘以OmegaC，得到非归一化零极点
k = k0*OmegaC^N;                % 将k0乘以OmegaC^N,得到非归一化k
b = real(poly(z));b = k*b       % 由零点计算分子系数向量
a = real(poly(p))               % 由极点计算分母系数向量
b=[zeros(1,length(a)-length(b)),b]; 	% 将b左端补零，使与a同长
[sos,G]=tf2sos(b,a)             	% 求分解为级联二阶环节的参数