數(shù)值分析課程課程設(shè)計-文庫吧資料
2025-01-24 15:57本頁面
【正文】 x4,y4)format longa1=[1 ,x1,x1^2,x1^3]。根據(jù)經(jīng)驗，一架水平飛行的飛機其降落曲線是一條三次曲線。3．2 繪制飛機的降落曲線一架飛機飛臨北京國際機場上空時，其水平速度為540km/h，飛行高度為1 000m。)結(jié)果分析：圖中紅色點線為正弦曲線，藍(lán)色的星線為一階泰勒逼近，綠色叉線為二階泰勒逼近，黑色正方形線為三階泰勒逼近。)title(39。)ylabel(39。,2)xlabel(39。,39。,39。,39。)hold offlegend(39。,[0,pi/2,0,2],2e3,39。)hold onfplot(39。,[0,pi/2,0,2],2e3,39。)hold onfplot(39。plot(x,y,39。hold onx=0::pi/2。sk39。y=sin(x)。xx^3/6+x^5/12039。,[0 pi/2]) (4) syms x。taylor(xx^3/6)fplot(39。taylor(sin(x))x=0:*pi:piplot(x,sin(x)) (2) syms x。X=H\bX = 與相差較大，矩陣為病態(tài)矩陣3．1 用泰勒級數(shù)的有限項逼近正弦函數(shù)用計算機繪出上面四個函數(shù)的圖形。b=[。1/2 1/3 1/4。b=H*Xb = 與題中相同(2) 先求出解X，與數(shù)據(jù)作比較 。1。1/3 1/4 1/5]。解：(1)H=[1 1/2 1/3。Y1=L\b1 X1=U\Y1 Y2=L\b2X2=U\Y2Y3=L\b3X3=U\Y3Y1 = X1 = Y2 =實驗驗證： [X1 X2 X3]ans = 而：inv（A）=[X1 X2 X3] 得證 0 X2 = Y3 = 0 X3 = 19 2．3 驗證希爾伯特矩陣的病態(tài)性：對于三階矩陣取右端向量，驗證：（1）向量是方程組的準(zhǔn)確解；（2）取右端向量b的三位有效數(shù)字得，求方程組的準(zhǔn)確解，并與X的數(shù)據(jù)作比較 。0 。,1]。L=[1,0,0。1,8,1。1]。 b3=[0。1。0]。解:MATLAB代碼如下：b1=[1。對第一個方程組，由于A=LU，所以先求解下三角方程組,再求解上三角方程組，則可得逆矩陣的第一列列向量；類似可解第二、第三方程組，得逆矩陣的第二列列向量的第三列列向量。2．2 用矩陣分解方法求上題中A的逆矩陣。 2 3 15]。h=zhjLU(A)程序運行結(jié)果：L = U= 0 0 0 0 0 0 實驗驗證： 可以直接使用MATLAB內(nèi)置LU分解A=[20 2 3。 2 3 15]。 endend以上上代碼保存為M文件，并在命令窗口輸入A=[20 2 3。 else U(k,j)=A(k,j)L(k,1:k1)*U(1:k1,j)。L(i,1)=A(i,1)/U(1,1)。 if ij L(1,1)=1。 end for k=2:n for i=2:n for j=2:n L(1,1)=1。)。 RA,hl return end end if h(1,i)~=0 disp(39。因為A的各階主子式等不等于零，：39。 end hl=h(1:n)。 RA,hl=det(A)。因為A的n階行列式hl等于零，：39。RA=rank(A)。設(shè)消元過程可將矩陣A化為上三角矩陣U，試求出消元過程所用的乘數(shù)、并以如下格式構(gòu)造下三角矩陣L和上三角矩陣U驗證：矩陣A可以分解為L和U的乘積，即A=LU。 %用向量方法計算迭代后處于4k位置上的點的坐標(biāo) n=m。 %用向量方法計算迭代后處于4k+2位置上的點的坐標(biāo) p(3:4:m,:)=q+d+d*A39。 %以原點為起點,前n1個點的坐標(biāo)為終點形成向量 p(5:4:m,:)=p(2:n,:)。 %diff計算相鄰兩個點的坐標(biāo)之差,得到相鄰兩點確定的向量 %則d就計算出每個向量長度的三分
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1