s1func.m

高频地波雷达完整仿真matlab程序 网上找到的

function [s1]=S1(y,sita,m1,k0,u,q)
  k1=y^2;
  g=9.8;
  kc=(g/u^2)/(2*k0);
  if m1==1
      S=(0.0081/k1^4)*exp(-0.74*kc^2/k1^2);
      G=(4/(3*pi))*(cos(sita/2-q/2))^4;
      s1=S*G; 
  elseif m1==-1
      S=(0.0081/k1^4)*exp(-0.74*kc^2/k1^2);
      G=(4/(3*pi))*(cos((sita+pi-q)/2))^4;
      s1=S*G;
  end