matlab解線性方程組的直接方法-在線瀏覽
2024-10-01 12:40本頁面
【正文】 :n)。請注意：因為A的r階主子式等于零，：39。RA,returnendend if h(1,i)~=0 disp(39。)hl。1 12 7。hl=pdLUfj(A)運行后輸出結(jié)果為請注意：因為A的各階主子式都不等于零，：RA = 3， hl = 1 10 48（2）在MATLAB工作窗口輸入程序 A=[1 2 3。4 5 6]。1 2 3。hl=pdLUfj(A)運行后輸出結(jié)果為請注意：因為A的n階行列式hl等于零，RA = 2， hl = 0 直接LU分解法及其MATLAB程序?qū)⒕仃囘M行直接LU分解的MATLAB程序function hl=zhjLU(A)[n n] =size(A)。 if RA~=ndisp(39。), RA,hl=det(A)。endhl=h(1:n)。請注意：因為A的r階主子式等于零，：39。RAreturnendend if h(1,i)~=0 disp(39。)for j=1:nU(1,j)=A(1,j)。L(i,i)=1。L(2,1)=A(2,1)/U(1,1)。L(i,k)=(A(i,k) L(i,1:k1)*U(1:k1,k))/U(k,k)。endendendendhl。0 1 0 1。0 1 0 3]。for p=1:nh(p)=det(A(1:p, 1:p))。zA=A39。請注意：因為A的各階順序主子式hl不全大于零，：39。zA,returnendendif h(1,i)0disp(39。) hl。1 2 3 4。5 7 8 9]。1 3 0 3。1 3 6 19],hl=zddc(A)運行后輸出結(jié)果A = 1 1 2 1 1 3 0 3 2 0 9 6 1 3 6 19 請注意： A是對稱矩陣請注意：因為A的各階順序主子式hl都大于零，：zA = 1 1 2 1 1 3 0 3 2 0 9 6 1 3 6 19 hl = 1 2 6 24 （3）在MATLAB工作窗口輸入程序 A=[1/sqrt(2) 1/sqrt(2) 0 0。 0 0 1/sqrt(2) 1/sqrt(2)。1 6 0。hl=zddc(A)運行后輸出結(jié)果 A = 2 1 1 1 6 0 1 0 4請注意： A是對稱矩陣請注意：因為A的各階順序主子式hl不全大于零，：zA = 2 1 1 hl = 2 11 38 1 6 0 1 0 4可見不是正定矩陣，是負定矩陣；因為 = T 因此，是對稱矩陣. 求解線性方程組的LU方法及其MATLAB程序 解線性方程組的楚列斯基（Cholesky）分解法及其MATLAB程序 先將矩陣進行楚列斯基分解,然后解矩陣方程，并用其他方法驗證..解 在工作窗口輸入A=[1 1 2 1。 2 0 9 6。b1=1:2:7。 R=chol(A)。*R,R1=inv(R)。 x=R1*R2*b,Rx=A\bx運行后輸出方程組的解和驗證結(jié)果x = Rx = * C = * 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 解線性方程組的直接LU分解法及其MATLAB程序 首先將矩陣直接進行LU分解，然后解矩陣方程，.解 (1) A=[1 0 2 0。1 2 4 3。b=[1。1。 hl=zhjLU(A),AL*U運行后輸出LU分解請注意：因為A的各階主子式都不等于零，：L = 1 0 0 0 0 1 0 0 1 2 1 00 1 0 1hl = 1 1 2 4RA = 4U = 1 0 2 0 0 1 0 1 0 0 2 1 0 0 0 2分解為一個單位下三角形矩陣和一個上三角形矩陣的積 .(2)在工作窗口輸入 U=[1 0 2 0。0 0 2 1。 L=[1 0 0 0。1 2 1 0。b=[1。1。U1=inv(U)。 X=U1*L1*b,x=A\b運行后輸出方程組的解X = 解線性方程組的選主元的LU方法及其MATLAB程序 先將矩陣進行LU分解，然后解矩陣方程 其中，.解 方法1 根據(jù)（）式編寫MATLAB程序，然后在工作窗口輸入 A=[ 2 3 4。11 21 13 41。 b=[1。1。 [L U P]=LU(A), U1=inv(U)。 X=U1* L1*P*bP = 0 0 1 0 0 1 0 0
黨政相關相關推薦
鄂ICP備17016276號-1