求解三對角方程組的追趕法(Matlab程式)

2021-06-05 10:37:06 字數 695 閱讀 2947

clear all;clc;

fprintf('輸入n:(10,20,30)\n');

n=input('');

na=zeros(1,n);b=zeros(1,n);c=zeros(1,n-1);

a=zeros(n,n);

for i=1:n

for j=1:n

if j==i

a(i,j)=2;

elseif (j-i==-1)||(j-i==1)

a(i,j)=1;

else

a(i,j)=0;

endend

endfprintf('係數矩陣:');a

fprintf('主對角線元素:');b=diag(a)'

fprintf('第1條對角線元素:');c=diag(a,1)'

fprintf('第-1條對角線元素:');a(2:n)=c

d=zeros(1,n);d(1)=1;d(n)=((-1)^(n+1));

fprintf('給定的係數矩陣(轉置後):');d

y=zeros(n,1);x=y;

u=zeros(1,n);l=u;%u矩陣u的主對角線元素,l矩陣l的第-1條對角線元素

%---------追過程--------

u(1)=b(1);y(1)=d(1);

for i=2:n

l(i

三對角矩陣壓縮儲存 注意對角元素的下標

對角矩陣的壓縮儲存 對角矩陣是指所有非零元素全部集中在中心幾條對角線上的矩陣。下面以三對角矩陣 所有非零元素集中在中心三條對角線上 為例描述對角矩陣的壓縮儲存方法。圖2 8是乙個三對角矩陣,使用一維陣列a m 來壓縮儲存矩陣資訊,則陣列中的元素依次為a11,a12,a21,a22,a23,a32,a...

三對角陣的一維陣列表示

資料結構基礎 c語言版 第2版,ellis horowitz著,朱仲濤譯 2.9節,page78,習題5 三對角陣是除主對角線及其相鄰上下兩條對角線之外都是零元素的方陣,三條對角線中的元素可以按行儲存在一維陣列b中,a 0 0 對應b 0 程式設計實現在b中訪問a i j 0 i,j 如下 incl...

高斯消去法求解方程組的實現 C

高斯消去法,其實就是利用初等變換把矩陣轉換成上三角,便於求行列式,求逆,求線性方程組的解 只要動筆演算一下相關矩陣元素的索引,就能使得核心 特別短,初等變換如下 for int i 0 i column 1 i 方程求解 如下 double result new double column 1 fo...