用高斯消元法解線性方程組 的matlab程式

時間 2021-08-11 18:15:20

1樓:刺友互

1、下圖是需要求解的線性方程組。

2、開啟matlab,利用左除法(\)求解上述線性方程組。輸入如下**:close all; clear all; clc% matlab左除法(\)求解線性方程組,a = [1 2 3;-1 3 7;9 0 3];b = [1 4 7]';x = a\b。

3、儲存和執行上述**,利用左除法(\)得到線性方程組的解。

4、用求逆法(inv)求解線性方程組,輸入如下**:close all; clear all; clc,% matlab求逆法(inv)求解線性方程組,% a是線性方程組等號左邊係數構成的矩陣。

5、儲存和執行上述**,利用求逆法(inv)得到線性方程組的解如下。

6、最後,可以看到左除法(\)和求逆法(inv)求得的線性方程組解是一樣的。

2樓:匿名使用者

%function [ra,rb,n,x]=gaus(a,b)

a=[1 2 3;

1 4 9;

1 8 27];

b=[123

];b=[a b];

n=length(b);

ra=rank(a); %a的秩

rb=rank(b); %b的秩

zhica=rb-ra;

if zhica>0

disp('請注意:因為ra~=rb,所以此方程組無解.');

return;

endif ra==rb %方程有唯一解

if ra==n

disp('請注意:因為ra=rb=n,所以此方程組有唯一解.')

x=zeros(n,1);

%生成上三角矩陣

for p= 1:n-1

for k=p+1:n

m= b(k,p)/ b(p,p); %要將b(k,p)化為零的係數

b(k,p:n+1)= b(k,p:n+1)-m* b(p,p:n+1); %整行乘以這個係數,最

endend

b=b(1:n,n+1);

a=b(1:n,1:n);

%求解上三角矩陣

x(n)=b(n)/a(n,n);

for q=n-1:-1:1

x(q)=(b(q)-sum(a(q,q+1:n)*x(q+1:n)))/a(q,q);

endelse %多解

disp('請注意:因為ra=rb

endend

找列主元高斯消去法來求解線性代數方程組解的matlab程式

3樓:匿名使用者

function x=gauss_pivot(a,b)% 用gauss列主主元消去法解線性方程組ax=b%x是未知向量

n=length(b);

x=zeros(n,1);

c=zeros(1,n);

d1=0

for i=1:n-1

max=abs(a(i,i));

m=i;

for j=i+1:n

if max

max=abs(a(j,i));

m=j;

endend

if(m~=i)

for k=i:n

c(k)=a(i,k);

a(i,k)=a(m,k);

a(m,k)=c(k);

endd1=b(i);

b(i)=b(m);

b(m)=d1;

endfor k=i+1:n

for j=i+1:n

a(k,j)=a(k,j)-a(i,j)*a(k,i)/a(i,i);

endb(k)=b(k)-b(i)*a(k,i)/a(i,i);

a(k,i)=0;

endend

%回代求解

x(n)=b(n)/a(n,n);

for i=n-1:-1:1

sum=0;

for j=i+1:n

sum=sum+a(i,j)*x(j);

endx(i)=(b(i)-sum)/a(i,i);end

4樓:eleven紅毛

function x=gauss_pivot(a,b)n=length(b);

x=zeros(n,1);

c=zeros(1,n);

d1=0

for i=1:n-1

max=abs(a(i,i));

m=i;

for j=i+1:n

if max~

62616964757a686964616fe59b9ee7ad9431333264663032=i)

for k=i:n

c(k)=a(i,k);

a(i,k)=a(m,k);

a(m,k)=c(k);

endd1=b(i);

b(i)=b(m);

b(m)=d1;

endfor k=i+1:n

for j=i+1:n

a(k,j)=a(k,j)-a(i,j)*a(k,i)/a(i,i);

endb(k)=b(k)-b(i)*a(k,i)/a(i,i);

a(k,i)=0;

endend

x(n)=b(n)/a(n,n);

for i=n-1:-1:1

sum=0;

for j=i+1:n

sum=sum+a(i,j)*x(j);

endx(i)=(b(i)-sum)/a(i,i);end

線性代數—用高斯消元法解線性方程組 15

5樓:匿名使用者

a =1 2 3 20 -2 -5 -3

3 4 3 -2

a =1 2 3 20 -2 -5 -3

0 -2 -6 -8

a =1 2 3 20 -2 -5 -3

0 0 -1 -5

a =1 2 3 20 -2 -5 -3

0 0 1 5

a =1 2 0 -130 -2 0 22

0 0 1 5

a =1 0 0 90 -2 0 22

0 0 1 5

a =1 0 0 90 1 0 -11

0 0 1 5

x1 = 9

x2 = -11

x3 = 5

6樓:匿名使用者

進行行變化用a的列減去r(a)得基礎解系

求齊次線性方程組的解,要具體過程

設解向量為x x1,x2,x3 初等變換之後 1,1,2 因為x是3維向量,x的方程組係數矩陣的秩為1,所以基礎解系含解個數為3 1 2。同解方程組是 x1 x2 2 x3 0 通解為x1 1 k1 2 k2 x2 1 k1 x3 1 k2 k1,k2是任意常數 於是基礎解系就是n1 1,1,0 t...

求齊次線性方程組的基礎解系,並求方程組的通解

墨汁諾 係數矩陣 a 2 3 1 5 3 1 2 4 1 2 3 1 初等行變換為 1 2 3 1 2 3 1 5 3 1 2 4 初等行變換為 1 2 3 1 0 7 7 7 0 7 7 7 初等行變換為 1 0 1 1 0 1 1 1 0 0 0 0 方程組同解變形為 x1 x3 x4,x2 x...

非齊次線性方程組的解向量個數的問題

條件沒有問題.非齊次方程的解與對應的齊次方程的基礎解系是線性無關的,也就是說非齊次方程ax b的解向量組成的向量組的秩 n 秩 a 1,n是未知數個數.記得同濟版線性代數課後有相關的習題.對於本題來說,秩 a 1時,ax b就可以找到四個線性無關的解.例如,a 1 0 0 0 0 0 0 0 0 0...