程序設(shè)計(jì)復(fù)習(xí)題 (2).docx

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