欢迎光临程序员联合开发网(www.pudn.com)

				% 最小相位序列复倒谱分析
				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);
				    end
				end   
				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('复对数虚部');