【文章內(nèi)容簡介】 我們在原有的解題思路中加入C語言程序，并讓它來幫助我們解決矩陣的復(fù)雜運算。既然用到了程序，我們就必須保證其運算的簡潔性、正確性，尤其是在編寫過程中要認(rèn)真檢查，為程序順利運行打下基礎(chǔ)。另外在各個子程序調(diào)用過程中，我們要充分考慮其順序性并反復(fù)調(diào)試，以便得到理想結(jié)果。盡管在這次課程設(shè)計中遇到了很多困難，但我卻得到了不少收獲，并培養(yǎng)了自己正確應(yīng)用公式、綜合分析和解決問題的能力，同時也為今后步入社會打下了一定的基礎(chǔ)。另外，我們還要學(xué)會綜合利用自身所學(xué)的知識，并將它們聯(lián)系起來幫助自己有效地解決實際中的問題?？傊谶@次課程設(shè)計中我不但過了比較充實的一周，還收獲了不少知識。includeiostreamincludefstreamincludeincludeincludeiomanipincludeusing namespace std。 class SZWPC{ private: int gcz_zs。 //高差總數(shù) int szd_zs。 //總點數(shù) int yz_szd_zs。 //已知點數(shù) double s0。//單位權(quán)水準(zhǔn)路線長度 double m_pvv。 //[pvv] int *qsd_dh。 //高差起點號 int *zd_dh。 //高差終點號 char **dm。 //點名地址數(shù)組 double *gcz。 //觀測值數(shù)組 double *szd_gc。 //高程值數(shù)組 double *P。 //觀測值的權(quán) double *BTPB,*BTPL。 //法方程系數(shù)矩陣與自由項 double *dX。 //高程改正數(shù)、平差值 double *V。 //殘差 double zwc。 //單位權(quán)中誤差 public: SZWPC()。 ~SZWPC()。 int ij(int i,int j)。//對稱矩陣下標(biāo)計算函數(shù) bool inverse(double a[],int n)。//對稱正定矩陣求逆(僅存下三角元素)(參考他人） void inputdata(char *datafile)。//輸入原始數(shù)據(jù)函數(shù) int dm_dh(char *name)。 //點名轉(zhuǎn)點號 void ca_H0()。 //近似高程計算函數(shù) void ca_BTPB()。 //法方程組成函數(shù) void ca_dX()。 //高程平差值計算函數(shù) void printresult(char *resultfile)。 //精度估計與平差值輸出函數(shù) double ca_V()。 //殘差計算函數(shù) void zxecpc(char *resultfile)。//最小二乘平差函數(shù)}。//////// // 構(gòu)造函數(shù)SZWPC::SZWPC(){ gcz_zs=0。 szd_zs=0。 yz_szd_zs=0。}//// // 析構(gòu)函數(shù)SZWPC::~SZWPC(){ if(gcz_zs0) { delete []qsd_dh。 delete []zd_dh。 delete []gcz。 delete []P。 delete []V。 } if(szd_zs0) { delete []szd_gc。 delete []BTPB。 delete []BTPL。 delete []dX。 for(int i=0。 iszd_zs。i++) if(dm[i]!=NULL)delete[](dm[i])。 delete []dm。 }}////// // 對稱矩陣下標(biāo)計算函數(shù)int SZWPC::ij(int i,int j){ return (i=j)? i*(i+1)/2+j :j*(j+1)/2+i。} // 對稱正定矩陣求逆(僅存下三角元素)bool SZWPC::inverse(double a[],int n){ double *a0=new double[n]。 for(int k=0。kn。k++) { double a00=a[0]。 if(a00+==) { delete []a0。 return false。 } for(int i=1。in。i++) { double ai0 = a[i*(i+1)/2]。 if(i=nk1)a0[i]=ai0/a00。 else a0[i]=ai0/a00。 for(int j=1。j=i。j++) { a[(i1)*i/2+j1]=a[i*(i+1)/2+j]+ai0*a0[j]。 } } for(int i=1。in。i++) { a[(n1)*n/2+i1]=a0[i]。 } a[n*(n+1)/21]=。 } delete []a0。 return true。}////////// 原始數(shù)據(jù)輸入函數(shù)void SZWPC::inputdata(char *datafile){ ifstream infile(datafile,ios::in)。//聲明輸入句柄infile打開地址為datafile的文件并 if(!infile) { cerr Open error!endl。 } infilegcz_zsszd_zsyz_szd_zs。 infiles0。 szd_gc=new double [szd_zs]。 dX=new double [szd_zs]。 BTPB=new double [szd_zs*(szd_zs+1)/2]。 BTPL=new double [szd_zs]。 qsd_dh=new int [gcz_zs]。 zd_dh=new int [gcz_zs]。 gcz=new double [gcz_zs]。 V=new double [gcz_zs]。 P=new double [gcz_zs]。 dm=new char* [szd_zs]。 for(int i=0。iszd_zs