图3-2-7预畸变在双线性变换中的作用.txt

二阶数字谐振器的频率特性和滤波器的幅频特性和相位特性和预畸变在双线性变换中的作用

%图3-2-4双线性变换模拟滤波器和数字滤波器的幅频特性曲线
clear;
fs=1500;
bs=2000;
as=[1 2000];
%freqs(bs,as);

figure(1);
omeg=linspace(0,fs*pi,1000);  %FREQ中频率需输入弧度OMEG=0~0.5fs
h=freqs(bs,as,omeg);
habs=10*log(abs(h));
f=omeg/2/pi;  %归一化频率1=0.5FS
%plot(f,habs,'k-');
plot(f,abs(h),'k-');
hold on;


hold on;
%figure(2);
%冲激不变法
[bd,ad]=bilinear(bs,as,fs); 				%fs修正幅度变化
%[bd,ad]=impinvar(bs,as,fs); 				%fs修正幅度变化
[h,w]=freqz(bd,ad,1000);			%数字滤波器频响 0~0.5FS
habs2=10*log(abs(h));
w=w'/pi/2*fs;
w=w(1:1000);
habs2=habs2(1:1000);
plot(w,abs(h),'k:');

hold on;
bs=2360.53;
as=[1,2360.53];
[bd,ad]=bilinear(bs,as,fs); 				%fs修正幅度变化
%[bd,ad]=impinvar(bs,as,fs); 				%fs修正幅度变化
[h,w]=freqz(bd,ad,1000);			%数字滤波器频响 0~0.5FS
habs2=10*log(abs(h));
w=w'/pi/2*fs;
w=w(1:1000);
habs2=habs2(1:1000);
plot(w,abs(h),'k-.');