main1.m

基于超外差的瞬时测频方法的研究。请多多指教

clear
clc
jj=1;
save number jj;
num=100;
for i=1:num
    

pulseLen=5000;
freq=20;
snr=6;
p_signal=10^(snr/10);
amp=sqrt(p_signal*2);
dt=5e-9;
a1=10;
a2=1;
a3=10;
signal=Creatnlfm(pulseLen,a1,a2,a3,amp,dt);
noise1=randn(1,pulseLen);
noise2=randn(1,pulseLen);
noise=noise1+j*noise2;
x=signal+noise;
len=length(x);
delay=round(0.25*len);
x1=x(delay+1:len).*conj(x(1:len-delay));
[startFreq,LFMk]=delfm(x1,dt);
a4=LFMk/(6*delay*dt);
a3=a4
tin1=1:5000;
y1=x.*exp(-j*2*pi*(a3*tin1.^3*dt*dt*dt));
[startFreq,LFMk]=delfm(y1,dt)
a1=startFreq;
a2=LFMk/2;


load number;
num1(jj)=a1;
num2(jj)=a2;
num3(jj)=a3;

jj=jj+1;
save number jj num1 num2 num3;
clear;
end

num=100;
load number;
num1
num2
num3
mean1=mean(num1)
mean2=mean(num2)
mean3=mean(num3)
rms1=num1-mean1;
rms1=rms1.^2;
rms1=sqrt(sum(rms1)/num)
rms2=num2-mean2;
rms2=rms2.^2;
rms2=sqrt(sum(rms2)/num)
rms3=num3-mean3;
rms3=rms3.^2;
rms3=sqrt(sum(rms3)/num)