error in my code

3 次查看（过去 30 天）

common fernando 2020-12-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/676488-error-in-my-code

评论： common fernando 2020-12-3

clear all; close all; clc
fs=200;     %sampling freq.
dt =1/fs;
N0=fs/3;	%number of samples/cycle
m=3;        %no. of cycles
t = linspace(0,200,1+N0*m); %data window
fi=50; %Frequency test
ww=wgn(201,1,-40);
size(transpose(ww))
x= sin(2*pi*fi*t + 0.3);
v = bsxfun( @plus, x , ww  );
%v=@(t) (sin(2*pi*fi.*t + 0.3)+ transpose(wgn(1+N0*m,1,-40)));
tmax=1;
n=N0-1:-1:0;
f0=50;
f=50.88;
Hc=2/N0*cos(2*pi*n/N0+pi/N0);
Hs=-2/N0*sin(2*pi*n/N0+pi/N0);
t_est=[];
f_est=[];
j_max=tmax*fs;
for j=1:j_max+1
	x=v((j-1:j+N0-2)*dt);
    c(j)=x*Hc';
    s(j)=x*Hs';    
    if(j>N0)
        Ac(j-N0)=sqrt(sum(c(end-N0+1:end).^2)/N0);   
        As(j-N0)=sqrt(sum(s(end-N0+1:end).^2)/N0);    
        cc(j-N0)=c(end-N0+1:end)*Hc';
        ss(j-N0)=c(end-N0+1:end)*Hs';
        if(j>2*N0)
            Acc(j-2*N0)=sqrt(sum(cc(end-N0+1:end).^2)/N0);   
            Ass(j-2*N0)=sqrt(sum(ss(end-N0+1:end).^2)/N0);   
            ccc(j-2*N0)=cc(end-N0+1:end)*Hc';
            ccs(j-2*N0)=cc(end-N0+1:end)*Hs';
            ssc(j-2*N0)=ss(end-N0+1:end)*Hc';
            sss(j-2*N0)=ss(end-N0+1:end)*Hs';
            ff=f0*N0/pi*atan(tan(pi/N0)*((ccc(j-2*N0).^2+ccs(j-2*N0).^2)./(ssc(j-2*N0).^2+sss(j-2*N0).^2)).^.25);
            t_est=[t_est;(j-1)*dt];
            f_est=[f_est;ff];
        end
    end
end
t_est;
f_est
     RMSE = sqrt(mean((f_est-fi).^2))
plot(t_est,f_est,'red') 
hold on
xlabel('time')
ylabel('frequency')
title('DE white noise')
plot (t,fi)
plot (t,fi)
hold off