?? samp8_6.m
字號:
%Samp8_6
[b,a]=butter(4,0.4); %產生4階Butterworth濾波器
[h,w]=freqz(b,a,64); %計算頻率響應
[bb,aa]=invfreqz(h,w,3,3); %采用invfreqz求解濾波器的系數
wt=ones(size(w)); %產生權向量
niter=30; %迭代次數
[bbb,aaa]=invfreqz(h,w,3,3,wt,niter); %采用權向量和迭代次數求濾波器系數
[h1,w1]=freqz(bb,aa,w); %計算模型1的頻率響應
[h2,w2]=freqz(bbb,aaa,w); %計算模型2的頻率響應
plot(w/pi,abs(h),'k',w1/pi,abs(h1),'k:',w2/pi,abs(h2),'b.'); %繪制頻率響應
xlabel('歸一化頻率');ylabel('振幅')
legend('理想模型','未采用權向量和迭代次數','采用權向量和迭代次數');
grid on
sumerr1=sum(abs(h-h1).^2) %輸出模型1誤差
sumerr2=sum(abs(h-h2).^2) %輸出模型2誤差
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -