作業(yè)4空間后方交會(huì)(參考版)
2025-03-27 07:03本頁(yè)面
【正文】 17 。雖然有些地方是效仿的，不過也是自己要弄懂以后，再試著自己一步步理出來。不斷調(diào)試，不斷改進(jìn)，這個(gè)過程還是有些難度的，最后還實(shí)現(xiàn)了用文本導(dǎo)入數(shù)據(jù)計(jì)算及文本輸出結(jié)果。其實(shí)主要總結(jié)起來就一個(gè)比較難的點(diǎn)：編程實(shí)現(xiàn)矩陣的運(yùn)算尤其是矩陣求逆。而編程方面所面臨的問題就比較多了，上一學(xué)期就有過一個(gè)類似的作業(yè)，當(dāng)時(shí)沒有攻破編程這一題，便用的Excel解決，但事后也沒有繼續(xù)研究，實(shí)為遺憾。此次作業(yè)，便是又一次的挑戰(zhàn)。return 0。delete []js。 pt++。j++) { *ptemp=*pt。i++){ for (j=0。for (i=0。 } }}ptemp=matrix。irow。j++) { swap(*(m+r*row+j),*(m+js[r]*row+j))。r){ if (js[r]!=r) { for (j=0。//恢復(fù)行列for (r=row1。i++) { if(i!=k) { *(m+i*row+k)*=(*(m+k*row+k))。 } } } } for (i=0。jrow。irow。j++) { if (j!=k) { *(m+k*row+j)*=*((m+k*row+k))。 for (j=0。i++) { swap(*(m+i*row+k),*(m+i*row+js[k]))。 } } //列交換 if (js[k]!=k) { for (i=0。irow。 } } } if (0 == max) { return 1。 is[k]=i。jrow。irow。k++){ double max=0。for (k=0。 }}int k。 ptemp++。jrow。irow。ptemp=matrix。double *ptemp,*pt=m。)**********************************************/int InverseMatrix(double *matrix,const int amp。}/***********************************************函數(shù)名：InverseMatrix *函數(shù)介紹：求矩陣的逆（高斯約當(dāng)法） *輸入?yún)?shù)：(in/out)matrix（矩陣首地址），*（in）row（矩陣階數(shù)）*輸出參數(shù)：matrix（原矩陣的逆矩陣）*返回值：int ，0成功，1失敗*調(diào)用函數(shù)：swap(doubleamp。a=b。a,double amp。return 0。delete []mm。delete []ATL。delete []AT。delete []L。}OutPut(Q,mm,m0,Xs,Ys,Zs,Phi,Omega,Kappa,amp。i6。double *mm=new double[6]。 VTV+=V[i]*V[i]。i2*Num。double VTV=0。MatrixMul(amp。//注：協(xié)因數(shù)陣，旋轉(zhuǎn)矩陣等計(jì)算本應(yīng)該使用最后外方位元素值，//由于變換很小忽略double *Q=ATA。 Kappa+=Xg[5]。 Phi+=Xg[3]。 Ys+=Xg[1]。 MatrixMul(ATA,6,ATL,1,6,Xg)。 if(InverseMatrix(ATA,6)) return 1。 } } //求ATA MatrixMul(AT,6,amp。j6。i2*Num。 A[2*i+1][5]=(xxx)。 A[2*i+1][3]=(xxx)*sin(Omega)(((yyy)/f)* ((xxx)*cos(Kappa)(yyy)*sin(Kappa)) f*sin(Kappa))*cos(Omega)。 A[2*i+1][1]=(R[1][1]*f+R[1][2]*(yyy))/Z。 A[2*i][5]=(yyy)。 A[2*i][3]=(yyy)*sin(Omega)(((xxx)/f)* ((xxx)*cos(Kappa)(yyy)*sin(Kappa))+ f*cos(Kappa))*cos(Omega)。 A[2*i][1]=(R[1][0]*f+R[1][2]*(xxx))/Z。 L[2*i+1]=Data[i][1](f*Y/Z)。 yyy=f*Y/Z。 double xxx,yyy。 double Y=R[0][1]*(Data[i][2]Xs)+R[1][1]*(Data[i][3]Ys)+ R[2][1]*(Data[i][4]Zs)。iNum。 R[2][2]=cos(Phi)*cos(Omega)。 R[2][0]=sin(Phi)*cos(Kappa)+cos(Phi)*sin(Omega)*sin(Kappa)。 R[1][1]=cos(Omega)*cos(Kappa)。 R[0][2]=sin(Phi)*cos(Omega)。//迭代計(jì)算do { //旋轉(zhuǎn)矩陣 R[0][0]=cos(Phi)*cos(Kappa)sin(Phi)*sin(Omega)*sin(Kappa)。double *ATL=new double[6]。double *AT=new double[2*Num*6]。typedef
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1