fe12.m

不用太多解释了吧

function dE12=fE12(E2,E3,E5,E9,Rm,Rn,Plat,Ph,Ve,Vn,Ee)
Wie=15*pi/(180*3600);
dE12=(Wie*sin(Plat)+Ve*tan(Plat)/(Rn+Ph))*E2-(Wie*cos(Plat)+Ve/(Rn+Ph))*E3+(1/(Rm+Ph))*E5-(Vn/((Rm+Ph)^2))*E9+Ee;