ar频率估计 matlab,求matlab高手帮小弟看下这个程序关于频率估计的..

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

clear all

clc

A=5;

f0=0.2;

w0=f0;

t=0:49;

N=length(t);

alfa=pi/4;

zao=randn(1,N)+randn(1,N)*j;

ww=zeros(1,3);

for i=1:2

y=A*exp(j*((w0-ww(1,i))*t)+alfa)+zao;

wt=3*N/(2*(N*N-1))*(1-((2*t-N-2)/N).^2);

for m=1:N-1;

if(real(y(m+1)>=0))

a(m)=atan(imag(y(m+1))/real(y(m+1)));

else

a(m)=atan(imag(y(m+1))/real(y(m+1)))+pi;

end

if(real(y(m)>=0))

b(m)=atan(imag(y(m))/real(y(m)));

else

b(m)=atan(imag(y(m))/real(y(m)))+pi;

end

if(y(m)<0&y(m+1)>0)

x1(m)=a(m)-b(m)+2*pi;

else

x1(m)=a(m)-b(m);

end

end

sum=0;

for m=1:N-1;

sum=sum+wt(m)*x1(m);

end

ww(1,i+1)=sum;

end

w=ww(1,i)+ww(1,i+1)

w =

0.1998

clear all

clc

A=1;

f0=10;

w0=f0;

N=100

t=0:N-1;

alfa=pi/4;

y=A*exp(j*(w0*t+alfa));

wt=3*N/(2*(N*N-1))*(1-((2*t-N-2)/N).^2);

for m=1:N-1;

x1(m)=((atan(imag(y(m+1))/real(y(m+1))))-(atan(imag(y(m))/real(y(m)))));

end

sum=0;

for m=1:N-1;

sum=sum+wt(m)*x1(m);

N =

100