大地测量学实习报告 - 图文

double eep=0.006738525414683;//第二偏心率的平方 double ee=0.006693421622966; double a=6378245.0000000000; double P=180*3600/PI;

double N1, A1 ,RA, S, D, H1, H2, h1, h2, D12 ,Bm,ym , Mm, Nm, Rm, x1, x2,y2 ,y1,f1 ,f2,f;//h1h2为正常高,D12高斯平面边长,Bm为平均纬度，Rm平均曲率半径

UpdateData(); h1=m_h1,h2=m_h2;

x1=m_x1,y1=m_y1-500000,x2=m_x2,y2=m_y2-500000; B1=(m_B1du+m_B1fen/60+m_B1miao/3600)*PI/180; B2=(m_B2du+m_B2fen/60+m_B2miao/3600)*PI/180; D=m_D;//实测长度 ym=0.5*(y1+y2);

A1=(m_A12du+m_A12fen/60+m_A12miao/3600)*PI/180; N1=a/sqrt(1-ee*sin(B1)*sin(B1));

RA=N1/(1+eep*cos(B1)*cos(B1)*cos(A1)*cos(A1)); H1=m_h1+m_h0; H2=m_h2+m_h0;

f=1-(H2-H1)/D*(H2-H1)/D; f1=1+H1/RA; f2=1+H2/RA;

S=D*sqrt(f/(f1*f2))+D*D*D/(24*RA*RA); Bm=0.5*(B1+B2);

Mm=a*(1-ee)/sqrt((1-ee*sin(Bm)*sin(Bm))*(1-ee*sin(Bm)*sin(Bm))*(1-ee*sin(Bm)*sin(Bm)));

Nm=a/sqrt(1-ee*sin(Bm)*sin(Bm)); Rm=sqrt(Mm*Nm);

D12=(1+ym*ym/(2*Rm*Rm)+(y2-y1)*(y2-y1)/(24*Rm*Rm)+ym*ym*ym*ym/(24*Rm*Rm*Rm*Rm))*S;

m_S=S; m_D12=D12;

UpdateData(FALSE); }

大地主题正反算主要代码：

void CMyDlg::OnButton1() {

// TODO: Add your control notification handler code here

double ee=0.006693421622966;//第一偏心率的平方 double eep=0.006738525414683;//第二偏心率的平方 double a=6378245.0000000000;

double PI=3.14159265358979323846264338327950288; double P=180*3600/PI; double

B1,L1,A12,S12,M1,N1,t,n,dertB0,dertL0,dertA0,Bm,Lm,Am,dertB,dertL,dertA,Mm,Nm,tm,nm,B2,L2,A21,Vm,f1,f2,f3,f4; UpdateData();

B1=(m_B1du+m_B1fen/60+m_B1miao/3600)*PI/180;//化弧度 L1=(m_L1du+m_L1fen/60+m_L1miao/3600)*PI/180;

A12=(m_A12du+m_A12fen/60+m_A12miao/3600)*PI/180; S12=m_S12;

f1=1-ee*sin(B1)*sin(B1); M1=a*(1-ee)/sqrt(f1*f1*f1);

N1=a/sqrt(1-ee*sin(B1)*sin(B1)); t=tan(B1);

n=sqrt(eep*cos(B1)*cos(B1)); dertB0=S12*cos(A12)/M1;

dertL0=S12*sin(A12)/(N1*cos(B1)); dertA0=dertL0*sin(B1);//单位弧度 Bm=B1+0.5*dertB0; Lm=L1+0.5*dertL0; Am=A12+0.5*dertA0;

f2=1-ee*sin(Bm)*sin(Bm); Mm=a*(1-ee)/sqrt(f2*f2*f2);

Nm=a/sqrt(1-ee*sin(Bm)*sin(Bm)); tm=tan(Bm);

Vm=sqrt(1+eep*cos(Bm)*cos(Bm)); nm=sqrt(eep*cos(Bm)*cos(Bm));

f4=3*cos(Am)*cos(Am)*nm*nm*(tm*tm-1-nm*nm-4*nm*nm*tm*tm);

f3=1+S12*S12/(24*Nm*Nm)*(sin(Am)*sin(Am)*(2+3*tm*tm+2*nm*nm)+f4); dertB=S12*cos(Am)/Mm*f3;

f4=tm*tm*sin(Am)*sin(Am)-cos(Am)*cos(Am)*(1+nm*nm-9*nm*nm*tm*tm+nm*nm*nm*nm);

f3=1+S12*S12/(24*Nm*Nm)*f4;

dertL=S12*sin(Am)/Nm/cos(Bm)*f3;

f4=cos(Am)*cos(Am)*(2+7*nm*nm+9*nm*nm*tm*tm+5*nm*nm*nm*nm)+sin(Am)*sin(Am)*(2+tm*tm+2*nm*nm); f3=1+S12*S12/(24*Nm*Nm)*f4; dertA=S12*sin(Am)*tm/Nm*f3; do {

dertB0=dertB;

dertL0=dertL; dertA0=dertA;

f4=3*cos(Am)*cos(Am)*nm*nm*(tm*tm-1-nm*nm-4*nm*nm*tm*tm);

f3=1+S12*S12/(24*Nm*Nm)*(sin(Am)*sin(Am)*(2+3*tm*tm+2*nm*nm)+f4); dertB=S12*cos(Am)/Mm*(f3);

f4=tm*tm*sin(Am)*sin(Am)-cos(Am)*cos(Am)*(1+nm*nm-9*nm*nm*tm*tm+nm*nm*nm*nm);

f3=1+S12*S12/(24*Nm*Nm)*f4;

dertL=S12*sin(Am)/Nm/cos(Bm)*f3;

f4=cos(Am)*cos(Am)*(2+7*nm*nm+9*nm*nm*tm*tm+5*nm*nm*nm*nm)+sin(Am)*sin(Am)*(2+tm*tm+2*nm*nm); f3=1+S12*S12/(24*Nm*Nm)*f4; dertA=S12*sin(Am)*tm/Nm*f3;

Bm=B1+0.5*dertB; Lm=L1+0.5*dertL; Am=A12+0.5*dertA;

Mm=a*(1-ee)/sqrt((1-ee*sin(Bm)*sin(Bm))*(1-ee*sin(Bm)*sin(Bm))*(1-ee*sin(Bm)*sin(Bm)));

Nm=a/sqrt(1-ee*sin(Bm)*sin(Bm)); tm=tan(Bm);

Vm=sqrt(1+eep*cos(Bm)*cos(Bm)); nm=sqrt(eep*cos(Bm)*cos(Bm)); }

while(((dertB-dertB0)>1.0e-10)||((dertL-dertL0)>1.0e-10)||((dertA-dertA0)>1.0e-10)); B2=B1+dertB; B2=B2*180/PI; L2=L1+dertL; L2=L2*180/PI; A21=A12+dertA+PI; A21=A21*180/PI; m_B2du=(int)B2;

m_B2fen=(int)((B2-(int)B2)*60);

m_B2miao=B2*3600-m_B2du*3600-m_B2fen*60; m_L2du=(int)L2;

m_L2fen=(int)((L2-(int)L2)*60);

m_L2miao=L2*3600-m_L2du*3600-m_L2fen*60; m_A21du=(int)A21;

m_A21fen=(int)((A21-(int)A21)*60);

m_A21miao=A21*3600-m_A21du*3600-m_A21fen*60; UpdateData(FALSE); }

void CMyDlg::OnButton2() {

// TODO: Add your control notification handler code here double ee=0.006693421622966;//第一偏心率的平方 double eep=0.006738525414683;//第二偏心率的平方 double a=6378245.0000000000;

double PI=3.14159265358979323846264338327950288; double P=180*3600/PI;

double B1, L1 , B2 ,L2 ,Bm, dertB,dertA, dertL, Mm, Nm ,r01, r21, r03,s10,s12, s30,t01,t21,t03,Am,U,Vm, V,t m, nm ,T,C,S,A1,A2,f1; UpdateData();

B1=(m_b1du+m_b1fen/60+m_b1miao/3600)*PI/180;//化弧度 L1=(m_l1du+m_l1fen/60+m_l1miao/3600)*PI/180;

B2=(m_b2du+m_b2fen/60+m_b2miao/3600)*PI/180;//化弧度 L2=(m_l2du+m_l2fen/60+m_l2miao/3600)*PI/180;

Bm=0.5*(B1+B2); dertB=B2-B1; dertL=L2-L1;

f1=1-ee*sin(Bm)*sin(Bm); Mm=a*(1-ee)/sqrt(f1*f1*f1);

Nm=a/sqrt(1-ee*sin(Bm)*sin(Bm)); nm=sqrt(eep*cos(Bm)*cos(Bm)); tm=tan(Bm);

Vm=sqrt(1+eep*cos(Bm)*cos(Bm)); r01=Nm*cos(Bm);

r21=Nm*cos(Bm)/(24*Vm*Vm*Vm*Vm)*(1+nm*nm-9*nm*nm*tm*tm+nm*nm*nm*nm); r03=-Nm*cos(Bm)*cos(Bm)*cos(Bm)*tm*tm/24; s10=Nm/(Vm*Vm);

s12=-Nm*cos(Bm)*cos(Bm)*(2+3*tm*tm+2*nm*nm)/(24*Vm*Vm); s30=Nm*(nm*nm-tm*tm*nm*nm)/(8*Vm*Vm*Vm*Vm*Vm*Vm); t01=tm*cos(Bm);

t21=cos(Bm)*tm*(2+7*nm*nm+9*tm*tm*nm*nm)/(24*Vm*Vm*Vm*Vm); t03=cos(Bm)*cos(Bm)*cos(Bm)*tm*(2+tm*tm+2*nm*nm)/24; U=r01*dertL+r21*dertB*dertB*dertL+r03*dertL*dertL*dertL; V=s10*dertB+s12*dertB*dertL*dertL+s30*dertB*dertB*dertB; dertA=t01*dertL+t21*dertB*dertB*dertL+t03*dertL*dertL*dertL; C=fabs(V/U);

if(fabs(dertB)>=fabs(dertL)) T=atan(U/V); else

T=PI/4+atan(fabs((1-C)/(1+C))); if(dertB>0&&dertL>=0)