a=6378245
b=6356863
e=(a^2-b^2)/a^2
N1=a/sqrt(1-e^2*(sind(25))^2)
r1=N1*cosd(25)
N2=a/sqrt(1-e^2*(sind(47))^2)
r2=N2*cosd(47)
s1=asind(e*sind(25))/2
U1=tand(45+25/2)/(tand(45+s1))^e
s2=asind(e*sind(47))/2
U2=tand(45+47/2)/(tand(45+s2))^e
c=(log10(r2)-log10(r1))/(log10(U1)-log10(U2))
K=r1*(U1)^c/c
s3=asind(e*sind(15))/2
Us=tand(45+15/2)/(tand(45+s3))^e
Ps=K/(Us)^c
Q1=[15:5:55]';
Q2=70:5:145;
V=c*Q2;
x=(Ps-K./(tand(45+Q1./2)./(tand(45+asind(e*sind(Q1))./2)).^e).^c*cosd(V))/10000000;
y=K./(tand(45+Q1./2)./(tand(45+asind(e*sind(Q1))./2)).^e).^c*sind(V)/10000000;
z=0*Q1*Q2/10000000;
mesh(x,y,z)
view([0,90])
版权所有:编程辅导网 2021 All Rights Reserved 联系方式:QQ:99515681 微信:codinghelp 电子信箱:99515681@qq.com
免责声明:本站部分内容从网络整理而来,只供参考!如有版权问题可联系本站删除。