%----------------------------------------------------------------------------
% exa050402 for example 5.4.2 and fig5.4.3
% to explain the minimum,maxmum-phase system;
%--------------------------------------------------------------------------
clear

% 構(gòu)造了關(guān)于單位圓對(duì)稱的零點(diǎn)
N=20;
a=0.8;
b=0.6;
b1=[1 -b];
a1=[1 -a];
b2=[b -1];

% 比較二者幅頻響應(yīng)和相頻響應(yīng)
[H1 w]=freqz(b1,a1,256,1);      % 最小相位系統(tǒng)
[H2 w]=freqz(b2,a1,256,1);      % 最大相位系統(tǒng)
fupin1=abs(H1);
fupin2=abs(H2);
subplot(221)
plot(w,fupin1.*fupin1,'b.',w,fupin2.*fupin2,'b.');grid on;
subplot(222)
plot(w,unwrap(angle(H1)),w,unwrap(angle(H2)));grid on;

% 計(jì)算兩個(gè)系統(tǒng)的沖激響應(yīng)
h1=impz(b1,a1,N);
h2=impz(b2,a1,N);
t=0:N-1;
subplot(223)
stem(t,h1,'.')
hold on;
plot(t,zeros(size(h1)))
subplot(224)
stem(t,h2,'.')
hold on;
plot(t,zeros(size(h2)))