數(shù)值分析課程課程設(shè)計(jì)-免費(fèi)閱讀
【正文】 50。)運(yùn)行結(jié)果：sol = x: [0 1] y: [2x5 double] yp: [2x5 double] solver: 39。)legend(39。ya39。t39。解析：ode**函數(shù)無(wú)能為力Matlab中，提供了 bvp解算器。x39。f39。endans=[T39。Y=zeros(1,N+1)。)。39。x(1,:)=x039。x=diff(cos(t))。h=ba。for i=1:length(x) for j=1:2 A(i,j)=x(i)^(j1)。y39。y=[1,3,4,5,6]。x39。 y2。a2=[1 ,x2,x2^2,x2^3]。Taylor approximation39。xx^3/639。xx.^3/6+x.^5/12039。y=x。taylor(xx^3/6+x^5/120)fplot(39。]。X=[1。U=[20 2 3。1。由三個(gè)列向量拼裝可得逆矩陣。b=[0 0 0]39。L(2,1)=A(2,1)/U(1,1)。)。if RA~=n disp(39。 %迭代后處于4k+1位置上的點(diǎn)的坐標(biāo)為迭代前的相應(yīng)坐標(biāo) p(2:4:m,:)=q+d。10 0]?，F(xiàn)在需要在直線段的中間依次插入三個(gè)點(diǎn)產(chǎn)生第一次迭代的圖形（圖14）。 if i=amp。 for k=1:5p=5*p/4+1。指 導(dǎo) 老 師 評(píng) 語(yǔ)建議：從學(xué)生的工作態(tài)度、工作量、設(shè)計(jì)(論文的)創(chuàng)造性、學(xué)術(shù)性、使用性及書(shū)面表達(dá)能力等方面給出評(píng)價(jià)。題6．4 常微分方程的解，dsolve()函數(shù)使用題8．2 差分法解常微分方程邊值問(wèn)題，ode函數(shù)無(wú)能為力，Matlab中提供bvp解算器。個(gè)別題目設(shè)計(jì)驗(yàn)證，加深對(duì)數(shù)值分析的理解。查閱資料時(shí)了解到一些MATLAB論壇。n= 39。i=39。1．3 繪制Koch分形曲線問(wèn)題描述：從一條直線段開(kāi)始，將線段中間的三分之一部分用一個(gè)等邊三角形的另兩條邊代替，形成具有5個(gè)結(jié)點(diǎn)的新的圖形（圖14）；在新的圖形中，又將圖中每一直線段中間的三分之一部分都用一個(gè)等邊三角形的另兩條邊代替，再次形成新的圖形（圖15），這時(shí)，圖形中共有17個(gè)結(jié)點(diǎn)。矩陣的第一列元素分別為5個(gè)結(jié)點(diǎn)的x坐標(biāo) ，第二列元素分別為5個(gè)結(jié)點(diǎn)的y坐標(biāo)。 %diff計(jì)算相鄰兩個(gè)點(diǎn)的坐標(biāo)之差,得到相鄰兩點(diǎn)確定的向量 %則d就計(jì)算出每個(gè)向量長(zhǎng)度的三分之一,與題中將線段三等分對(duì)應(yīng) m=4*n3。設(shè)消元過(guò)程可將矩陣A化為上三角矩陣U，試求出消元過(guò)程所用的乘數(shù)、并以如下格式構(gòu)造下三角矩陣L和上三角矩陣U驗(yàn)證：矩陣A可以分解為L(zhǎng)和U的乘積，即A=LU。 end hl=h(1:n)。 end for k=2:n for i=2:n for j=2:n L(1,1)=1。 endend以上上代碼保存為M文件，并在命令窗口輸入A=[20 2 3。2．2 用矩陣分解方法求上題中A的逆矩陣。1。L=[1,0,0。解：(1)H=[1 1/2 1/3。1/2 1/3 1/4。taylor(xx^3/6)fplot(39。sk39。,[0,pi/2,0,2],2e3,39。,39。)ylabel(39。根據(jù)經(jīng)驗(yàn)，一架水平飛行的飛機(jī)其降落曲線是一條三次曲線。a3。y=s(3)*x.^2s(4)*x.^3plot(x,y,39。下面數(shù)據(jù)中，第二行數(shù)據(jù)為晶片上晶體數(shù)目在不同年代與1959年時(shí)數(shù)目比較的倍數(shù)。)xlabel(39。4。（1）牛頓萊布尼茨公式；（2）梯形公式；（3）辛卜生公式；（4）復(fù)合梯形公式。 i4=i4+4/(1+x^2)。i=vpa(y)運(yùn)行結(jié)果：i = 6．1 用歐拉公式和四階龍格庫(kù)塔法分別求解下列初值問(wèn)題；解：（1）1歐拉公式：function [t,x]=Euler(fun,t0,tt,x0,N)h=(ttt0)/N。39。Dx=*x/(1+2*t)39。保存輸入：main_Euler運(yùn)行結(jié)果：ans = 0 2四階龍格庫(kù)塔法function R=rk4(f,a,b,ya,N)%y39。 k3=h*feval(f,T(j)+h/2,Y(j)+k2/2)。以39。,39。解： p=[1 55 1320 18150 157773 902055 3416930 8409500 12753576 10628640 6328800]。[y(2)。[ya(1)。y=deval(sol,t)。x39。1 7 1。y0=(x(2)/2)。設(shè)已知圓通過(guò)三點(diǎn)P1（1,3）,P2(1,7),P3(6,2),其方程為(xx0)2+(yy0)2=R2按下列步驟求圓的半徑R和圓心的坐標(biāo)（x0,y0）(1)設(shè)圓的一般方程為x2+y2+ax+by+c=0將圓上三個(gè)點(diǎn)的坐標(biāo)數(shù)據(jù)分別代入上面方程；（2）整理后寫關(guān)于未知數(shù)a,b,c的三元線性方程組AX=b（3）用MATLAB解方程組的指令A(yù)＼b求解；（4）通過(guò)配方求出圓的半徑R和圓心坐標(biāo)。解點(diǎn)39。)。)。sinit=bvpinit(0:1,[0,0])。f39。解v=dsolve(39。3
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1