fseries.m

SIRP法相干相关K分布雷达杂波的建模与仿真

format long;
sigma_f=20;%功率谱3dB带宽
Fs=512;%采样频率
T0=1/Fs;%采样间隔
f=0:T0:Fs/2;
abs_H_f=0;
for n=0:13;
    A(n+1)=2*sigma_f*T0*sqrt(pi)*exp(-4*sigma_f^2*pi^2*n^2*T0^2);
    if n==0;
        abs_H_f=abs_H_f+A(n+1);
    else
        abs_H_f=abs_H_f+2*A(n+1)*cos(2*pi*n*f*T0);
        if n==3;
            temp=abs_H_f;
        end
    end
end
powerf_theory=exp(-f.^2/(2*sigma_f^2));
powerf_theory=powerf_theory/max(powerf_theory);%归一化
powerf_fseries=abs_H_f.^2;
powerf_fseries=powerf_fseries/max(powerf_fseries);%归一化
temp=temp/max(temp);%归一化
plot(f,powerf_theory,f,powerf_fseries,'r',f,temp.^2,'g');
xlabel('f');
ylabel('S(f)')
legend('理论值','N=9时的拟合值','N=3时的拟合值');
% grid on;
text(12,0.6,'1');
text(20,0.3,'2');
text(50,0.4,'3');