數(shù)值計(jì)算實(shí)驗(yàn)word版(已修改)
2025-01-19 16:08 本頁面
【正文】 數(shù) 值 實(shí) 驗(yàn) 數(shù)值實(shí)驗(yàn) 1 線性方程組求解 一、方法與程序 1． LUPA? ：帶選主元的分解法 （ MATLAB 程序） Function x =lufact(A,b) % Input A is an N? N matrix % b is N? 1 matrix %Output x is an N? 1 matrix containing the solution to Ax=b %Initialize x, y, the temporary storage matrix C, and the row %permutation information matrix R [N,N]=size(A)。 x=zeros(N,1)。 y=zeros(N,1)。 c=zeros(N,1)。 R=1:1。 For k=1:N1 % Find the pivot row for colum p [max1,j]=max (abs(A(p:N,p)))。 % Interchange row p and j c=A(p,:)。 A(p,:)=A(j+p1,:)。 A(j+p1,:)=c。 D=R(p)。 R(p)= R(j+p1)。 R(j+p1)=d。 if a(p,p)==0 ‘A is singular. No unique solution’ break end %Calculate multiplier and place in sundiagonal portion of A For k=p+1:N mult=A(k,p)/A(p,p)。 A(p,:)=A(j+p1,:)。 A(k,p)=mult。 A(k,p+1:N)= A(k,p+1:N)mult* A(p,p+1:N)。 end end %solve for y y(1)=b(R(1))。 For k=2:N Y(k)=b(R(k))A(k,1:k1)*y(1:k1) end %solve for x x(N)=y(N)/A(N,N)。 For k=N1:1:1 x(k)=(y(k)A(k,k+1:N)*x(k+1:N))/A(k,k)。 end 二、數(shù)值試驗(yàn)內(nèi)容 1） 用 帶選主元的分解法 求解線性方程組 bAx? ,其中 ???????????????????1362520053127531A 和 ???????????????4321b 使用 MATLAB 中的 [L， U， P]=lu(A)命令檢查得到的答案． 2） 使用帶選主元的分解法 求解線性方程組 bAx? ,其中 NNijaA ?? ][ ，1?? jij ia ， ,][ 111 Nba Nij ?? 當(dāng) 2?i 時 )1/()1(1 ??? iib Ni ．對于 11,7,3?N 的情況分別求解． 精確解為 Tx )1111( ?? ．對得到的結(jié)果與精確解的差異進(jìn)行解釋． 數(shù)值實(shí)驗(yàn) 2 Lagrange 插值 數(shù)值實(shí)驗(yàn)內(nèi)容：對一組數(shù)據(jù)做 Lagrange 插值，根據(jù)插值多項(xiàng)式估計(jì)函數(shù)值 . 調(diào)
公司管理相關(guān)推薦
鄂ICP備17016276號-1