數(shù)值計(jì)算實(shí)驗(yàn)word版-在線瀏覽
2025-02-24 16:08本頁(yè)面
【正文】 4321b 使用 MATLAB 中的 [L， U， P]=lu(A)命令檢查得到的答案． 2） 使用帶選主元的分解法 求解線性方程組 bAx? ,其中 NNijaA ?? ][ ，1?? jij ia ， ,][ 111 Nba Nij ?? 當(dāng) 2?i 時(shí) )1/()1(1 ??? iib Ni ．對(duì)于 11,7,3?N 的情況分別求解． 精確解為 Tx )1111( ?? ．對(duì)得到的結(jié)果與精確解的差異進(jìn)行解釋． 數(shù)值實(shí)驗(yàn) 2 Lagrange 插值 數(shù)值實(shí)驗(yàn)內(nèi)容：對(duì)一組數(shù)據(jù)做 Lagrange 插值，根據(jù)插值多項(xiàng)式估計(jì)函數(shù)值 . 調(diào)用格式： yi=Lagran_(x,y,xi) x,y: 數(shù)組形式的數(shù)據(jù)表 xi: 待計(jì)算函數(shù)值的橫坐標(biāo)數(shù)組 yi: 用 Lagrange 插值多項(xiàng)式算出的 y 值數(shù)組 Function fi=Lagran_(x,f,xi) fi=zeros(size(xi)) npl=length(f) for i=1:npl z=ones(size(xi)) for j=1:npl if I~=j,z=z.*(xix(j))/(x(i)x(j))。 T=zeros(1,N)。 復(fù)化 Simpson 求積公式 利用 )(xf 在 12 ?N 個(gè)等步長(zhǎng)采樣點(diǎn)：khaxk ?? Nkxfxfbfafhdxxf Nk kba Nk k 2,1,0),)(4)(2)()((3)( 1 1211 2 ?????? ?? ? ? ??? 的逼近積分．注意： bxax N ?? 20 , 復(fù)化 Simpson 求積算法（ MATLAB 程序 ） Function s =simpson(f,a,b,N) %Input f is the integrand input as a string ‘f’ % a and b are upper and lower limits of integration % N is the number of subintervals %Output s is the simpson rule sum h= (ba)/(2*N)。 s2=0。 s1=s1+feval(f,x) end For k=1:(N1) x=a+h*2*k。 二、數(shù)值試驗(yàn)內(nèi)容 1）使用 6 點(diǎn) GaussLengder 求積公式逼近積分 : ? ??? 30 22 )()()( dttvtxxxv 即用求積公式去離散積分方程中的積分項(xiàng) ,并求積以上積分方程 ,得到其解函數(shù)的近似表達(dá)式 . 2 ）用的 復(fù)化 Simpson 求積公式求下列積分，要求絕對(duì)誤差限為71021 ???? ． 1． ?? ??11 12 )1( dxx 2． ? ?40 2 dxex x 3． ? ??0 )2sin( dxex x 數(shù)值實(shí)驗(yàn) 4 微分方程數(shù)值解法 一 、方法與程序 4 階 Rungekutta 法 )22(6 43211 kkkkhuu kk ?????? 4 階 Rungekutta 算法（ MATLAB 程序） Function R =rk4(f,a,b,ya,N) %Input f is the
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1