p1func.m

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

function [p1]=P1(y,sita,m1,m2)
 temp=1+m1*m2*y*(y^2+cos(sita))/(1+2*y^2*cos(sita)+y^4)^(3/4);
 if temp==0
     temp=0.00000000001;
 end
 p1=1/abs(temp);