fe52.m

不用太多解释了吧

function dE52=fE52(E1,E3,E4,E5,E6,E7,Rm,Rn,Plat,Ph,Ve,Vu,Vn,Dn,Fe,Fu)
Wie=15*pi/(180*3600);
dE52=Fu*E1-Fe*E3-2*(Wie*sin(Plat)+Ve*tan(Plat)/(Rn+Ph))*E4-(Vu/(Rm+Ph))*E5-(Vn/(Rm+Ph))*E6-(2*Wie*cos(Plat)+Ve/((Rn+Ph)*(cos(Plat))^2))*Ve*E7+Dn;