%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誤差