程序設計復習題 (2).docx

一、 解釋以下語句含義 A3： 數乘矩陣，矩陣A乘以2 A-B： 矩陣A與矩陣B之差 ab：比較數組a是否等于數組b cos(x): x的余弦函數Fix(x)：對x結尾去整b=ones(3,3)：生成3行3列全1元素矩陣AB： A乘以B或矩陣A與矩陣B之積A+B： 數字A與數字B相加或A，B兩個同階矩陣相加ab：比較數組a是否大于數組bsin(x): x的正弦函數sqrt(x)：x的平方根A=zeros(2,3)：生成2行3列全0元素矩陣二、 程序設計如右圖所示，導線控制測量，已知A，B點，觀測了角度和邊長，給出程序設計。1. 角度平差 ：1）先給出A、B點的起始數據2）定義一個4行3列的全0矩陣angle輸入已知角度，定義一個3行3列的全0矩陣l輸入已知距離。3）計算角度閉合差，然后平均分配到四個角當中求出角度平差值并放入到矩陣angle當中。2.計算坐標方位角：1）先定義一個4行4列的全0矩陣alfa。2）根據已知A.B點的坐標以及公式acos(Dx/sqrt(Dx*Dx+Dy*Dy)推算出起始方位角。3）根絕DY是否小于0正確推算出起始坐標方位角,根據公式alfa(i+1)=alfa(i)+pi+temp推算出各點的坐標方位角。并將左右的坐標方位角輸入到矩陣alfa當中。3.坐標平差：1）定義一個4行4列的全0矩陣dx, 定義一個4行4列的全0矩陣dy.2）根據公式dx(i)=dist(i)*cos(alfa(i+1)計算出各x軸增量，根據公式dy(i)=dist(i)*sin(alfa(i+1)計算出個y軸增量。3）計算出x軸y軸的累積增量來定權。4）計算出x軸y軸的數學增量來算距離平差值。5）最后根據距離平差值以及角度平差值來計算各點的坐標平差值。4.精度評定如右圖所示，水準控制測量，已知A點高程，采用四等水準觀測，給出程序設計。1.初始數據準備，計算閉合差：2. 定權，計算高差改正數：3. 計算高程：4.給出精度評定：三、程序編寫，要有注解1. 坐標正算：ax=get(handles.edit1,string)；%讀取文本框1中的信息ax=str2num(ax)； %把文本框1轉為數據格式，可調用ay=get(handles.edit2,string)；%讀取文本框2中的信息ay=str2num(ay)； %把文本框2轉為數據格式，可調用s=get(handles.edit3,string)； %讀取文本框3中的信息s=str2num(s)； %把文本框3轉為數據格式，可調用a=get(handles.edit4,string)； %讀取文本框4中的信息a=str2num(a)； %把文本框4轉為數據格式，可調用bx=ax+s*cos(pi*a/180)；by=ay+s*cos(pi*a/180)；轉換公式，得出數據set(handles.edit5,string,bx)；%把數據格式變?yōu)槲谋靖袷?，并展現(xiàn)出來set(handles.edit6,string,by) %把數據格式變?yōu)槲谋靖袷?，并展現(xiàn)出來2.坐標反算：ax=get(handles.edit1,string)；ax=str2num(ax)；ay=get(handles.edit2,string)；ay=str2num(ay)；bx=get(handles.edit3,string)；bx=str2num(bx)；by=get(handles.edit4,string)；by=str2num(by)；a=atan(abs(by-ay)/abs(bx-ax)*180/pi；if a=0&a=90&a=180&a=270&a=360 a=360-a；endset(handles.edit5,string,a)；I=sqrt(by-ay)2+(bx-ax)2)； %計算AB的水平距離set(handles.edit6,string,I) %把數據格式變?yōu)槲谋靖袷剑⒄宫F(xiàn)出來3.角度轉換為弧度：a=input(輸入a的數值) %輸入已知的角度值b=input(輸入b的數值) %輸入已知的分度值c=input(輸入c的數值) %輸入已知的秒度值y=(a+b/60+c/3600)*pi)/180 %利用轉換公式求出弧度值4.弧度轉化為角度：a=input(輸入a的數值) %輸入要轉化的弧度值b=a*180/pi %利用公式求解c=fix(b) %求出要求的角度值d=fix(b-c)*60) %求出要求的分度值f=(b-c)*3600-d*60 %求出要求的秒度值四、程序閱讀，說明程序功能，優(yōu)缺點x1=1000.000 %起始坐標數據y1=1000.000x4=1200.000y4=1000.000 angle=ones(4,3) %定義一個4行3列全0矩陣for i=1:4 %給出各角度數據 angle(i,1)=90 angle(i,2)=00 angle(i,3)=02enddist=zeros(3) %定義一個3行3列全0矩陣dist(1)=300.002 %給出距離值dist(2)=200.001dist(3)=300.000%*角度平差* deta=0.0for i=1:4 deta=deta+angle(i,1)*3600+angle(i,2)*60+angle(i,3)enddeta=deta-360*3600 %計算角度閉合差for i=1:4 %計算各角度平差值 angle(i,3)=angle(i,3)-deta/4end%*計算方位角* alfa=zeros(4); %定義一個4行4列全0矩陣Dx=x1-x4;Dy=y1-y4temp=acos(Dx/sqrt(Dx*Dx+Dy*Dy) %計算方位角if Dy0 %判斷坐標方位角 temp=2*pi-tempendalfa(1)=temp %給出起始坐標方位角值%*計算各點坐標方位角*for i=1:3 temp=angle(i,1)+angle(i,2)/60+angle(i,3)/3600 temp=temp*pi/180 %角度換算弧度 alfa(i+1)=alfa(i)+pi+tempenddisp(alfa) %顯示出alfa%*平面坐標平差* dx=zeros(4) %定義dx為一個4行4列全0矩陣dy=zeros(4) %定義dy為一個4行4列全0矩陣ddx=0ddy=0px=0py=0for i=1:3 dx(i)=dist(i)*cos(alfa(i+1) %x坐標增量 dy(i)=dist(i)*sin(alfa(i+1) %y坐標增量 px=px+abs(dx(i) %x坐標增量絕對值之和，定權用 py=py+abs(dy(i) %y坐標增量絕對值之和，定權用 ddx=ddx+dx(i) %x坐標增量累計(理論上=dx，實際有誤差） ddy=ddy+dy(i) %y坐標增量累計(理論上=dy，實際有誤差）enddetax=ddx+Dx %計算x坐標增量閉合差detay=ddy+Dy %計算y坐標增量閉合差for i=1:3 %計算各坐標增量改正值 dx(i)=dx(i)-abs(dx(i)/px)*detax dy(i)=dy(i)-abs(dy(i)/py)*detayendxc=zeros(4) %定義xc為一個4行4列全0矩陣 yc=zeros(4) %定義yc為一個4行4列全0矩陣xc(1)=x1 yc(1)=y1for i=1:3 %得出最終的坐標改正值 xc(i+1)=xc(i)+dx(i) yc(i+1)=yc(i)+dy(i)end- 第二個程序 -% *% * 高程數據：起算數據，高差數據，距離數據。 *% *H1 =100H12 =12.364D12 =647.3H23 =2.324D23 =842.1H34 =-3.683D34 =587.4H41 =-11.003D41 =732.3% *% * 計算高差閉合差deta，計算距離和sigma_D,計算高差平差值，計算高程平差值 *% *deta =H12+H23+H34+H41sigma_D=D12+D23+D34+D41H12 =H12-deta*(D12/sigma_D)H23 =H23-deta*(D23/sigma_D)H34 =H34-deta*(D34/sigma_D)H41 =H41-deta*(D41/sigma_D)H2 =H1+H12;H3 =H2+H23;H4 =H3+H34;H1 =H4+H41;% *% * 判斷是否合格，若合格輸出計算成果，否則查找錯誤 *% *fprintf(deta=%10.4fn,deta*1000)fprintf(STD=%10.4fn,4*sqrt(sigma_D/1000)if abs(deta*1000)4*sqrt(sigma_D/1000) fprintf(Good!n); fprin