hc932.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第九章例9.3.2程序hc932
% 系数量化误差对IIR滤波器极点分布的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear,close all
b=[1,0.4,-0.03,0.232];              % 分子系数向量
a=[1,-2.4760,2.9541,-1.9631,0.4900];% 分母系数向量
aq=brqtize(a,7), bq=brqtize(b,7);   % 量化后的分母、分子系数向量 
ra=(aq-a)./abs(a),                  % 量化后的分母系数的相对误差
p=roots(a),pq=roots(aq),	        % 分母系数量化前、后的极点
z=roots(b);zq=roots(bq);	        % 分子系数量化前、后的零点
rp=abs(pq-p)./abs(p)                % 量化前、后极点的相对误差
zplane(z,p), hold on,               % 画出量化前零、极点位置(o和x)
plot(pq,'*'),plot(zq,'+'),          % 画出量化后零、极点位置(+和*)
set(gcf,'color','w')                % 置图形背景色为白
pause,
[sos,G]=tf2sos(b,a);                % 化为二阶级联形式
a1=sos(1,[4:6]),a2=sos(2,[4:6])     % 取出sos中的二阶分母系数向量       
a1q=bqtize(a1,7);a2q=bqtize(a2,7);  % 量化后的二阶分母系数向量 
p1q=roots(a1q),p2q=roots(a2q),      % 二阶分母系数量化后的极点
abs(p1q),abs(p2q)                   % 极点的模
[H,w]=freqz(b,a);                   % 量化前的频率特性
[Hq,w]=freqz(bq,aq);                % 量化后的频率特性
figure(2),
subplot(2,1,1),plot(w,abs(H))
subplot(2,1,2),plot(w,abs(Hq))
set(gcf,'color','w')                 % 置图形背景色为白