% 最小相位序列復倒譜分析N=100;x=[1,1,2,3,4,5,6,7];n=1:8;k=1:N;w=exp(j*2*pi/N)*k;Xk=x*(exp(-j*2*pi/N).^(n'*(k-1)));FXk=log(Xk);FXrk=log(abs(Xk));COR=zeros(1,length(n));for i=2:N    if angle(Xk(i))-angle(Xk(i-1))>pi       COR(i)=COR(i-1)-2*pi;    elseif angle(Xk(i-1))-angle(Xk(i))>pi       COR(i)=COR(i-1)+2*pi;     else        COR(i)=COR(i-1);    endend   r=(1/pi)*Xk(N/2+1);sign=Xk(1)/abs(Xk(1));argX=angle(Xk)+COR-(2*pi/N)*r.*ones(1,N);fxn=sign*(1/N)*(FXrk+j*argX)*(exp(j*2*pi/N).^(k'*n));subplot(2,2,1);stem(n,x);title('x(n)');subplot(2,2,2);stem(n,fxn);title('fx(n)');subplot(2,2,3);stem(k,FXrk);title('復對數實部');subplot(2,2,4);stem(k,argX);title('復對數虛部');