【正文】 s(L(int(L)*63))*3600。//經(jīng)度L同當(dāng)前子午線的差值單位為秒 double s=p2*p2*p2*p2。 double v=l2*l2*l2*l2。 double u=cos(B)*cos(B)*cos(B)。 x=CompexS(x)+N1/(2*(p2)*(p2))*sin(B)*cos(B)*(l2)*(l2)+N/(24*(p2)*(p2))*sin(B)*u*(5t*t+9*n1*n1+4*n1*n1*n1*n1)*v+N1/(720*p2*p2*s)*sin(B)*cos(B)*cos(B)*u*(6158*t*t+t*t*t*t)*l2*l2*v。 y=N1/(p2)*cos(B)*(l2)+N*p2/(6*s)*u*(1t*t+n*n)*v/l2+N1/(120*p2*s)*cos(B)*cos(B)*u*(518*t*t+t*t*t*t58*n1*n1*t*t)*l2*v。 }void Transformation::fangsuan(double x,double y,int type,double n,double amp。b,double amp。l){ parameter(type)。 double Bf=CompexBf(x)。 double Y=y50000。 double X=x。 double Tf=tan(Bf)。 double nf=e1*cos(Bf)。 double Nf=a/sqrt(1e*e*sin(Bf)*sin(Bf))。 double Mf=Nf/(1+e1*e1*cos(Bf)*cos(Bf))。 double g=Nf*Nf*Nf*Nf。 double h=Y*Y*Y*Y。 double c=Tf*Tf*Tf*Tf。 double B=BfTf/(2*Mf*Nf)*Y*Y+Tf/(4*Mf*pow(Nf,3))*(5+3*Tf*Tf+nf*nf9*nf*nf*Tf*Tf)*hTf/(720*Mf*Nf*g)*(61+90*Tf*Tf+45*c)*h*Y*Y。 double L=1/(Nf*cos(Bf))*YNf/(6*g*cos(Bf))*(1+2*Tf*Tf+nf*nf)*h/Y+1/(120*Nf*g*cos(Bf))*(5+28*Tf*Tf+24*c+6*nf*nf+8*nf*nf*Tf*Tf)*h*Y。 b=dgreetodms(B)。 l=dgreetodms(L)。}void Transformation::clear(double x,double y,int type,double n,double b,double l){ x=0。 y=0。 type=0。 n=0。 b=0。 l=0。}include include 0145110615 ymh include 0145110615 ymh include include iostreaminclude void CMy0145110615ymh7Dlg::OnCbnSelchangeCombo1(){ // TODO: 在此添加控件通知處理程序代碼 type=()。}void CMy0145110615ymh7Dlg::OnBnClickedOk(){ // TODO: 在此添加控件通知處理程序代碼 //OnOK()。 UpdateData(true)。 Transformation sb。 (x, y, type, n, x, y)。 UpdateData(false)。}void CMy0145110615ymh7Dlg::OnBnClickedButton1()//正算{ // TODO: 在此添加控件通知處理程序代碼 UpdateData(true)。 Transformation sb。 ( b, l, type, n, x,y)。 UpdateData(false)。}void CMy0145110615ymh7Dlg::OnBnClickedButton2()//反算{ UpdateData(true)。 Transformation sb。 ( b, l, type, n, x,y)。 UpdateData(false)。 // TODO: 在此添加控件通知處理程序代碼} 運(yùn)行結(jié)果：53總 結(jié) 這個(gè)實(shí)驗(yàn)雖然運(yùn)行出來了可是結(jié)果卻是不正確的，正算和反算好像都出了問題，清除按鈕下面我也調(diào)用了類當(dāng)中清除的函數(shù)可是點(diǎn)清除按鈕的時(shí)候也沒反應(yīng)。我也仔細(xì)檢查過可是最終還是沒能解決問題。 在編程的過程中我發(fā)現(xiàn)一個(gè)問題就是程序沒運(yùn)行出來的時(shí)候大家都會(huì)想盡一切辦法去解決掉他們，當(dāng)程序運(yùn)行出來了程序內(nèi)部實(shí)現(xiàn)要是存在一些問題的話則很少有人會(huì)去改正這些錯(cuò)誤。 實(shí)驗(yàn)八 平差程序設(shè)計(jì)基礎(chǔ)一、 實(shí)驗(yàn)?zāi)康? 鞏固過程的定義與調(diào)用? 鞏固類的創(chuàng)建與使用? 鞏固間接平差模型及平差計(jì)算? 掌握平差程序設(shè)計(jì)的基本技巧與步驟二、 實(shí)驗(yàn)內(nèi)容水準(zhǔn)網(wǎng)平差程序設(shè)計(jì)。設(shè)計(jì)一個(gè)水準(zhǔn)網(wǎng)平差的程序，要求數(shù)據(jù)從文件中讀取，計(jì)算部分與界面無關(guān)。1. 水準(zhǔn)網(wǎng)間接平差模型：2. 計(jì)算示例：近似高程計(jì)算：3. 水準(zhǔn)網(wǎng)平差計(jì)算一般步驟(1)讀取觀測(cè)數(shù)據(jù)和已知數(shù)據(jù)；(2)計(jì)算未知點(diǎn)高程近似值；(3)列高差觀測(cè)值誤差方程；(4)根據(jù)水準(zhǔn)路線長(zhǎng)度計(jì)算高差觀測(cè)值的權(quán)；(5)組成法方程；(6)解法方程，求得未知點(diǎn)高程改正數(shù)及平差后高程值；(7)求高差觀測(cè)值殘差及平差后高差觀測(cè)值；(8)精度評(píng)定；(9)輸出平差結(jié)果。4. 水準(zhǔn)網(wǎng)高程近似值計(jì)算算法 實(shí)驗(yàn)代碼：pragma onceclass LevelControlPoint{public: LevelControlPoint(void)。 ~LevelControlPoint(void)。public: CString strName。//點(diǎn)名 CString strID。//點(diǎn)號(hào) float H。 bool flag。//標(biāo)記是否已經(jīng)計(jì)算出近似高程值，若計(jì)算出則為，否則為}。class CDhObs{public: CDhObs(void)。 ~CDhObs(void)。public: LevelControlPoint* cpBackObj。//后視點(diǎn) LevelControlPoint* cpFrontObj。//前視點(diǎn) double ObsValue。//高差值 double Dist。//測(cè)站的距離}。include include LevelControlPoint::LevelControlPoint(void){ strName=_T()。 strID=_T()。 H=0。 flag=0。}LevelControlPoint::~LevelControlPoint(void){}CDhObs::CDhObs(void){}CDhObs::~CDhObs(void){}pragma onceincludeinclude class AdjustLevel{public: AdjustLevel(void)。 ~AdjustLevel(void)。public: LevelControlPoint* m_pKnownPoint。//已知點(diǎn)數(shù)組 int m_iKnownPointCount。//已知點(diǎn)個(gè)數(shù) LevelControlPoint* m_pUnknownPoint。//未知點(diǎn)數(shù)組 int m_iUnknownPointCount。//未知點(diǎn)個(gè)數(shù) CDhObs* m_pDhObs。//高差觀測(cè)值數(shù)組 int m_iDhObsCount。//高差觀測(cè)值個(gè)數(shù)public: void SetKnownPointSize(int size)。//創(chuàng)建大小為size的已知點(diǎn)數(shù)組 void SetUnkonwnPointSize(int size)。//創(chuàng)建大小為size的未知點(diǎn)數(shù)組 void SetDhObsSize(int size)。//創(chuàng)建大小為size的觀測(cè)值數(shù)組 bool LoadObsData(const CStringamp。 strFile)。//讀入觀測(cè)文件 CString* SplitString(CString str, char split, intamp。 iSubStrs)。 void ApproHeignt(void)。//計(jì)算近似值private: LevelControlPoint* SearchKnownPointUsingID(CString ID)。 LevelControlPoint* SearchUnknownPointUsingID(CString ID)。 LevelControlPoint* SearchPointUsingID(CString ID)。 CMatrix LevleWeight(void)。//計(jì)算權(quán)矩陣 public: void FormErrorEquation(CMatrix amp。B, CMatrix amp。L)。//組成誤差方程 void EquationCompute(CMatrix amp。x)。//計(jì)算法方程 void Accuracy_Assessment(double amp。r0,CMatrix amp。Qxx)。//精度評(píng)定 void CompAdjust(double amp。r0,CMatrix Qx[])。}。include include include include include AdjustLevel::AdjustLevel(void){ m_pKnownPoint=NULL。//已知點(diǎn)數(shù)組 m_iKnownPointCount=0。//已知點(diǎn)個(gè)數(shù) m_pUnknownPoint=NULL。//未知點(diǎn)數(shù)組 m_iUnknownPointCount=0。//未知點(diǎn)個(gè)數(shù) m_pDhObs=NULL。//高差觀測(cè)值數(shù)組 m_iDhObsCount=0。//高差觀測(cè)值個(gè)數(shù)}AdjustLevel::~AdjustLevel(void){ if(m_pKnownPoint!=NULL) { delete[] m_pKnownPoint。 m_pKnownPoint=NULL。 } if(m_pUnknownPoint!=NULL) { delete[] m_pUnknownPoint。 m_pUnknownPoint=NULL。 } if(m_pDhObs!=NULL) { delete[] m_pDhObs。 m_pDhObs=NULL。 }}void AdjustLevel::SetKnownPointSize(int size){ m_pKnownPoint=new LevelControlPoint[size]。//創(chuàng)建動(dòng)態(tài)指針 m_iKnownPointCount=size。}void AdjustLevel::SetUnkonwnPointSize(int size){ m_pUnknownPoint=new LevelControlPoint[size]。 m_iUnknownPointCount=size。}void AdjustLevel::SetDhObsSize(int size){ m_pDhObs=new CDhObs[size]。 m_iDhObsCount=size。//高差觀測(cè)值個(gè)數(shù)}bool AdjustLevel::LoadObsData(const CStringamp。 strFile){ CStdioFile sf。 if(!(strFile,CFile::modeRead)) return false。//創(chuàng)建并打開文件對(duì)象 CString strLine。 bool bEOF=(strLine)。//讀取第一行，即已知點(diǎn)的數(shù)目 SetKnownPointSize(_ttoi(strLine))。//根據(jù)已知點(diǎn)的數(shù)目，創(chuàng)建已知點(diǎn)數(shù)組； int n=0。 for(int i=0。im_iKnownPointCount。i++)//讀取已知點(diǎn)的點(diǎn)名和高程值 { (strLine)。 CString *pstrData=SplitString(strLine,39。,39。,n)。 m_pKnownPoint[i].strName=pstrData[0]。 m_pKnownPoint[i].strID=pstrData[0]。 m_pKnownPoint[i].H=_tstof(pstrData[1])。 m_pKnownPoint[i].flag=1。//已知點(diǎn)不用平差，故將其的flag設(shè)置為 delete[] pstrData。 pstrData=NULL。 } (strLine)。//讀取未知點(diǎn)的個(gè)數(shù) SetUnkonwnPointSize(_ttoi(strLine))。//根據(jù)未知點(diǎn)的個(gè)數(shù)創(chuàng)建未知點(diǎn)數(shù)組 (str