main.m

来自「基于超外差的瞬时测频方法的研究。请多多指教」· M 代码 · 共 41 行

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

pulseLen=1000;
freq=20;
snr=-3;
p_signal=10^(snr/10);
amp=sqrt(p_signal*2);
k=1;
dt=5e-9;
signal=creatlfm(pulseLen,freq,amp,k,dt);
noise1=randn(1,pulseLen);
noise2=randn(1,pulseLen);
noise=noise1+j*noise2;
x=signal+noise;
[startFreq,LFMk]=delfm(x,dt)

load number;
num1(jj)=startFreq;
num2(jj)=LFMk;
jj=jj+1;
save number jj num1 num2;
clear;
end

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