次方系數(shù)對(duì)方程根的影響課程設(shè)計(jì)報(bào)告書(shū)-wenkub
2023-04-07 09:47:24 本頁(yè)面
【正文】 。 return。 end end L Uend4. 結(jié)果 LU_de(A)m = 4n = 4L = 1 0 0 0 2 1 0 0 1 2 1 0 3 3 2 1U = 2 4 2 6 0 1 2 3 0 0 3 6 0 0 0 1在編寫(xiě)程序過(guò)程中由于角標(biāo)較多因此在運(yùn)行過(guò)程中出現(xiàn)了不少角標(biāo)不對(duì)的錯(cuò)誤題目；給出函數(shù)f(x)=1/(1+25x^2)，求f（x在[1,1]上取5個(gè)和9個(gè)等距節(jié)點(diǎn)，做最小二乘擬合，得出均方誤差方誤差。enda=polyfit(x,y,3)syms xf1=a(1)*x^3+a(2)*x^2+a(3)*x+a(4) x=[1 0 1]。enda=polyfit(x,y,3)syms xf2=a(1)*x^3+a(2)*x^2+a(3)*x+a(4) x=[1 0 1]。h=(ba)/2。 h=(ba)/n。 x1=x+h。h=(xtx0)/(pointnum1)。 f1=f1(:)39。 f3=f3(:)39。endfunction[x,y]=euler(fun,x0,xt,y0,pointnum)if nargin4 y0=0。y(1,:)=y0(:)39。endRunge_kutta方法 x0=0。,39。)。（1）分別將五個(gè)點(diǎn)的數(shù)據(jù)代入橢圓一般方程中，寫(xiě)出五個(gè)待定系數(shù)滿(mǎn)足的等式，整理后寫(xiě)出線(xiàn)性方程組AX=R以及方程組的系數(shù)矩陣和右端項(xiàng)b；（2）用MARLAB求低階方程的指令A(yù)＼b求出待定系數(shù)a1,a2,a3,a4,a5.答：1)2)一．源程序：A=[53605^2 2*53605*6026 6026^2 2*53605 2*6026。 68894^2 2*68894*68894 68894^2 2*68894 2*68894]。首先運(yùn)用Romberg積分方法對(duì)給出定積分進(jìn)行積分，然后對(duì)得到的結(jié)果用插值方法，分別求出Lagrange插值多項(xiàng)式和Newton插值多項(xiàng)式，再運(yùn)用最小二乘法的思想求出擬合多項(xiàng)式，最后對(duì)這些不同類(lèi)型多項(xiàng)式進(jìn)行比較，找出它們各自的優(yōu)劣。二．問(wèn)題解決： 依據(jù)問(wèn)題出現(xiàn)先后順序，對(duì)三個(gè)小問(wèn)進(jìn)行討論。Romberg算法的計(jì)算過(guò)程列出如下表：kh(步長(zhǎng))T0(k)T1(k)T2(k)T3(k)0hT0(0)①1h/2T0(1)②T1(0)③2h/22T0(2)④T1(1)⑤T2(0)⑥3h/23T0(3)⑦T1(2)⑧T2(1)⑨T3(0)⑩………………表中①~⑩表示計(jì)算順序，k表示二分次數(shù)。h=ba。R(1,1)=h*(feval(39。,p,b))/2。 。ik)|i4 i=i+1。,p,a)+feval(39。i=0。a=0。 由復(fù)合梯形公式的展開(kāi)定理，得到如下關(guān)系式：T1(h)I=aah2+a2h4+a3h6+…+amh2m+…其中，I=,T1(h)=Tn.利用Richardson外推定理對(duì)T1(h)進(jìn)行加速，注意這里取m=1,q=,有利用T0()和T0()可以得到實(shí)際上T1(h)就是復(fù)化拋物線(xiàn)求積公式，一般的
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1