fe32.m

不用太多解释了吧

function dE32=fE32(E1,E2,E4,E7,E9,Rm,Rn,Plat,Ph,Ve,Vn,Eu)
Wie=15*pi/(180*3600);
dE32=(Wie*cos(Plat)+Ve/(Rn+Ph))*E1+(Vn/(Rm+Ph))*E2-(tan(Plat)/(Rn+Ph))*E4-(Wie*cos(Plat)+Ve/((Rn+Ph)*(cos(Plat))^2))*E7+(Ve*tan(Plat)/(Rn+Ph)^2)*E9+Eu;