功-率-谱-估-计程序.pptVIP

m=k时，有： 算法的程序编写流程图参见讲义P79或丁P154 图 3.3.1 利用列文森递推法计算功率谱的流程图 表明： 系统函数所有极点均在单位圆内 讨论： fs=1; N=10000; t=0:1/fs:N; A=20^0.5; B=A; C=1;　 x1=A*sin(2*0.15*pi*t); x2=B*sin(2*0.2*pi*t); noise=C*randn(size(t)); x=x1+x2+noise; %被分析信号 练习：自编程序，请理解、消化 M=15; %阶数 M1=20; %阶数 %求自相关函数 for m=0:M-1, % 对r赋初值 temp=0; for n=1:N-M, temp=x(n)*x(n+m)+temp; end r(m+1)=(1/N)*temp; end a(m+1,m+1)=0; a(1,1)=-r(2)/r(1); xgm(1)=(1-(a(1,1))^2)*r(1); % p=2; while pM p=p+1; temp=0; for i=1:p-2 temp=a(p-2,i)*r(p-i)+temp; end a(p-1,p-1)=-(r(p)+temp)/xgm(p-2); for k=1