1、考試課程 數學實驗 2004.6.24下午班級 姓名 學號 得分 說明（1）第一、二、三題的答案直接填在試題紙上；（2）第四題將數學模型、簡要解題過程和結果寫在試題紙上；卷面空間不夠時，請寫在背面；（3）除非特別說明，所有計算結果小數點后保留4位數字。（4）考試時間為120分鐘。一、（10分）已知一組數據x00.10.20.30.40.5y05.31.82.64.96.3用數值積分求x在（0, 0.5）內y的積分為 1.7750 。根據表中數據，試建立2次多項式回歸模型，其回歸系數為。在顯著性水平為0.05情況下有無異常數據（若有，請指出第幾個點為異常數據）？第二個。數值積分：梯形公式：x=0

2、:0.1:0.5;y=0 5.3 1.8 2.6 4.9 6.3;S1=trapz(x,y)輸出結果：S1 =1.7750多項式回歸：y=0 5.3 1.8 2.6 4.9 6.3;x1=0:0.1:0.5; %切削時間x2=x1.2; %刀具厚度n=6; %已知的數據容量T=ones(n,1),x1', x2' %1與自變量組成的輸入矩陣b,bint,r,rint,s=regress(y',T); %回歸分析程序(=0.05)b,bint,s,rcoplot(r,rint)輸出結果：b = 1.4821 5.58216.6071二、(10分) 已知常微分方程組初值問題

3、，試用數值方法求_0.897893_(保留小數點后6位數字)。你用的方法是_龍格庫塔方法 _，其精度為_四階_。%待解常微分方程組函數M文件源程序：function dy=ff (x,y) dy=y(2);-y(1)*sin(x)- exp(x);%應用歐拉方法和龍格-庫塔方法求解該常微分方程：ts=0:0.1:0.4;y0=1,0;x,y=ode45(ff, ts,y0); %龍格-庫塔方法求數值解x, y(:,1)輸出結果：三、 (10分) 已知線性代數方程組Ax=b, 其中, , ,若方程組右端項有小擾動，試根據誤差估計式估計_0.06416_ (保留小數點后5位數字)；若取初值，則用高

4、斯-賽德爾迭代法求解Ax=b時，_(1.07103,0.37900,-0.00072,-0.06084)_，此迭代是否收斂_是_。線性代數方程組解的誤差分析：故其誤差上限為：A=10 -7 0 1 ;-3 22 6 2 ;5 -1 51 -1 ;2 1 0 25; b=8 5 5 1;db=0 0 0 0.1;d=cond(A)*norm(db)/norm(b)輸出結果：d = 0.061460749095949A=10 -7 0 1 ;-3 22 6 2 ;5 -1 51 -1 ;2 1 0 25;D=diag(diag(A); %從稀疏矩陣A中提取DL=-tril(A,-1); %從稀疏矩

5、陣A中提取LU=-triu(A,1); %從稀疏矩陣A中提取U b=8 5 5 1' %設定方程組右端項向量x= zeros(4,1); %設定方程組初始向量 m= inv(D-L)*U;n= inv(D-L)*b; %高斯-賽德爾迭代法for j2=1:5 y=m*(x(:,j2); for i=1:4 x(i,j2+1)=y(i,:)+n(i,:); end endt2=x(:,end) %輸出迭代法最終結果j2 輸出結果： t2 = 1.071032233651353 0.378998741589479 -0.000722946127134 -0.060842528355687

6、高斯-賽德爾迭代判斂法：1) 若求收斂域,需要計算系數矩陣A的正定與對稱性質；2) 若僅僅判斷收斂，直接計算迭代矩陣L的譜半徑與義的大小關系；lamda=eig(inv(D-L)*U)pubanjing=max(abs(lamda)輸出結果：pubanjing =0.218937948256270四、（20分）飛機在飛行過程中，能夠收到地面上各個監控臺發來的關于飛機當前位置的信息，根據這些信息可以比較精確地確定飛機的位置。如下圖所示，高頻多向導航設備（VOR）能夠得到飛機與該設備連線的角度信息；距離測量裝置（DME）能夠得到飛機與該設備的距離信息。圖中飛機接收到來自3個VOR給出的角度和1個D

7、ME給出的距離（括號內是相應設備測量的精度，即絕對誤差限），并已知這4種設備的x，y坐標（假設飛機和這些設備在同一平面上）。請你根據這些信息確定當前飛機的位置，要求建立相應的數學模型并給出解答。提示：對角度信息進行處理時，可以考慮使用MATLAB的 atan2 函數。提示：對角度信息進行處理時，可以考慮使用MATLAB的 atan2 函數。0yxVOR2x=629, y=375309.00 (1.30)864.3(2.0)飛機x=?, y=?VOR1x=764, y=1393161.20 (0.80)DMEx=155, y=987VOR3x=1571, y=25945.10 (0.60)北DM

8、Ex=155, y=987第四題圖：飛機與監控臺（圖中坐標和測量距離的單位是“公里”）考試課程 數學實驗 2004.6.24下午班級 學號 姓名 得分 說明（1）第一、二、三題的答案直接填在試題紙上；（2）第四題將數學模型、簡要解題過程和結果寫在試題紙上；卷面空間不夠時，請寫在背面；（3）除非特別說明，所有計算結果小數點后保留4位數字。（4）考試時間為120分鐘。一、（10分）已知一組數據x00.10.20.30.40.5y00.31.88.64.96.3用數值積分求x在（0, 0.5）內y的積分為 。根據表中數據，試建立2次多項式回歸模型，其回歸系數為。在顯著性水平為0.05情況下有無異常數

9、據（若有，請指出第幾個點為異常數據）？ 。二、(10分) 已知常微分方程組初值問題，試用數值方法求_(保留小數點后6位數字)。你用的方法是_ _，其精度為_ _。三、 (10分) 已知線性代數方程組Ax=b, 其中, , , 若方程組右端項有小擾動，試根據誤差估計式估計_ (保留小數點后5位數字)；若取初值，則用高斯-賽德爾迭代法求解Ax=b時，_，此迭代是否收斂_。四、（20分）飛機在飛行過程中，能夠收到地面上各個監控臺發來的關于飛機當前位置的信息，根據這些信息可以比較精確地確定飛機的位置。如下圖所示，高頻多向導航設備（VOR）能夠得到飛機與該設備連線的角度信息；距離測量裝置（DME）能夠得

10、到飛機與該設備的距離信息。圖中飛機接收到來自3個VOR給出的角度和1個DME給出的距離（括號內是相應設備測量的精度，即絕對誤差限），并已知這4種設備的x，y坐標（假設飛機和這些設備在同一平面上）。請你根據這些信息確定當前飛機的位置，要求建立相應的數學模型并給出解答。提示：對角度信息進行處理時，可以考慮使用MATLAB的 atan2 函數。0yxVOR2x=629, y=375309.00 (1.30)861.3(2.0)飛機x=?, y=?VOR1x=764, y=1393161.20 (0.80)DMEx=155, y=987VOR3x=1571, y=25945.10 (0.60)北DME

11、x=155, y=987第四題圖：飛機與監控臺（圖中坐標和測量距離的單位是“公里”）考試課程 數學實驗 參考答案與評分標準 2004.6.24A卷（班級-姓名-學號）一、1.7750; 1.4821, 5.5821, 6.6071; 第2個點為異常數據二、0.897892 (或0.897893), 龍格-庫塔方法，3級2階(ode23)(或5級4階(ode45)三、0.06146, 1.0710, 0.3790, -0.0007, -0.0608, 收斂四、模型為 得到飛機的坐標為（978.3070，723.9838），誤差的平方和為0.6685。附主要程序示例：（fun為上述函數，程序略）X

12、=746 629 1571 155;Y=1393 375 259 987;theta=161.2,45.1,309.0-360*2*pi/360; % 角度轉換sigma=0.8,0.6,1.3*2*pi/360;d4=864.3; %B卷為861.3sigma4=2;x0=900,700; % 初值x,norm,res,exit,out=lsqnonlin(fun,x0,X,Y,theta,sigma,d4,sigma4)B卷（班級-學號-姓名）一、1.8750; -0.9821, 28.5464, -27.3214; 第4個點為異常數據二、0.831371 (0.831373), 龍格-庫

13、塔方法，3級2階(ode23)( 5級4階(ode45)三、0.07726, 1.0610, 0.3776, 0.0026, -0.0600, 收斂四、模型同A卷。計算結果：飛機的坐標為（975.2474，723.1857），誤差的平方和為0.9132。評分標準：一、第一空4分，中間三空每空1分，最后一空3分。二、第一空6分，后兩空每空2分。三、第一空3分，中間一空5分，最后一空2分。四、 模型15分（基本擬合模型5分；角度轉換正確5分；無量綱化處理正確5分）；計算結果5分（若因模型錯誤導致結果錯誤，可依計算方法或程序是否正確、完整給分，最多給3分；若無計算方法或程序的說明，則不給分）。考試課

14、程 數學實驗 參考答案與評分標準 2004.6.24A卷（班級-姓名-學號）一、1.7750; 1.4821, 5.5821, 6.6071; 第2個點為異常數據二、0.897892 (或0.897893), 龍格-庫塔方法，3級2階(ode23)(或5級4階(ode45)三、0.06146, 1.0710, 0.3790, -0.0007, -0.0608, 收斂四、模型為 得到飛機的坐標為（978.3070，723.9838），誤差的平方和為0.6685。附主要程序示例：（fun為上述函數，程序略）X=746 629 1571 155;Y=1393 375 259 987;theta=16

15、1.2,45.1,309.0-360*2*pi/360; % 角度轉換sigma=0.8,0.6,1.3*2*pi/360;d4=864.3; %B卷為861.3sigma4=2;x0=900,700; % 初值x,norm,res,exit,out=lsqnonlin(fun,x0,X,Y,theta,sigma,d4,sigma4)B卷（班級-學號-姓名）一、1.8750; -0.9821, 28.5464, -27.3214; 第4個點為異常數據二、0.831371 (0.831373), 龍格-庫塔方法，3級2階(ode23)( 5級4階(ode45)三、0.07726, 1.0610, 0.3776, 0.002