sigmacalca.m

高频地波雷达完整仿真matlab程序,可以实现地波传播模块、目标回波功率模块等6个模块的仿真

function sigemaf=SigmaCalcA(f0,u,q,w)

%主函数
g=9.8;
c=3*10^8;
% u=24*0.5144;              %由用户给出
% f0=10*10^6;               %由用户给出
% q=5*pi/6;                 %由用户给出
% w=4;                      %由用户给出
f0=f0*10^6;
k0=2*pi*f0/c;
wB=sqrt(2*g*k0);
eta=-4:0.01:4;
N=length(eta);
sigema=zeros(1,N);
for i=1:N
    if eta(i)==-1
        sigema(i)=sigema1L(q,u,f0)/16/2;
    elseif eta(i)==1
        sigema(i)=sigema1R(q,u,f0)/16/2;
    else
        [sitaL,sitaR]=limitfunc(eta(i));
        sigema(i)=Rombergfunc(sitaL,sitaR,eta(i),u,q,f0)/16/2;
    end
end 

 x=ifft(sigema);
 ww=gausswin(801,1);
 y=x.*ww.';
 s2=fft(y);

 for i=1:N
    if abs(s2(i))<0.0000000001
       s2(i)=0.0000000001;
    end
 end
 
 
 nn = 20;
 n1 = 400-nn;
 n2 = 400+nn;
 NWild = 0;
 for i=n1:n2
     if i>n1
         ss(1:i-n1) = abs(s2(n1:i-1));
         ss(i-n1+1:n2-n1) = abs(s2(i+1:n2));
         mm = mean(ss);
         if abs(s2(i))>4*mm
             NWild = NWild+1;
             wildVal(NWild) = i;
         end
         if abs(s2(i))<mm/4
             NWild = NWild+1;
             wildVal(NWild) = i;
         end
     end
 end
 
 sumS = sum(abs(s2(n1:n2)));
 for i=1:NWild
     sumS = sumS-abs(s2(wildVal(i)));
 end
 mm = sumS/(n2-n1+1-NWild);
 for i=1:NWild
     s2(wildVal(i)) = mm;
 end
 
 nn = length(s2);
 for i=1:nn
     if i<nn-2
         s2(i) = mean(abs(s2(i:i+2)));
     else
         s2(i) = mean(abs(s2(i-2:i)));
     end
 end
 
 sigema=s2;

% for i=1:N
%     if sigema(i)<0.00000001
%        sigema(i)=0.00000001;
%     end
% end
% sigemadB=10*log10(sigema);
% sigemadB=sigemadB';
% figure;
% plot(eta,sigemadB);
% grid;

%%bragg频率与目标doppler频率同为角频率或Hz时
% if (w<-4*wB)|(w>4*wB)
%     sigemaf=0;
% else 
%     tem=floor(100*w/wB)+400;
%     sigemaf=sigema(tem);

%%bragg频率为角频率，而目标doppler频率为Hz时   
if (w*2*pi<-4*wB)|(w*2*pi>4*wB)
    sigemaf=10^(-10);
else 
    tem=floor(100*2*pi*w/wB)+400;
    sigemaf=sigema(tem);
end