你的式子很复杂,我这里简单的以第一个式子计算PL为例子
由于其他参数已知,所以可以简化为
P=A*e*f/(1+f)^2;
其中e就是ε,而f就是φ,进一步简化,设这里A=1(具体数值自己计算代入)
[e f]=meshgrid(0:0.01:0.2,0:0.1:2);
P=e.*f./(1+f).^2;
surf(f,e,P);
这里e的取样是从0~0.2 以0.01为间隔
f的取样是0~2,以0.1为取样间隔
图像如下