该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
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