matlab數值分析

線性代數的直接接法??追趕法求解線性方程組Ax=b，其中A是三對角方陣

%追趕法求解線性方程組Ax=b，其中A是三對角方陣 function x=tridiagsolver(A,b) [n,n]=size(A); for i=1:nif(i==1)l(i)=A(i,i);y(i)=b(i)/l(i);else i<nl(i)=A(i,i)-A(i,i-1)*u(i-1); y(i)=(b(i)-A(i,i-1)*y(i-1))/l(i);endif(i<n)u(i)=A(i,i+1)/l(i);end end x(n)=y(n); for j=n-1:-1:1x(j)=y(j)-x(j+1)*u(j); end

調用程序

clc clear A=[2,-1,0,0;-1,3,-2,0;0,-2,4,-3;0,0,-3,5]; b=[6;1;-2;1]; x= tridiagsolver(A,b)

總結

以上是生活随笔為你收集整理的2021-01-07 matlab数值分析 线性代数的直接接法  追赶法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。