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...求齊次線性方程組的解,要具體過程
求齊次線性方程組的基礎解系,並求方程組的通解
非齊次線性方程組的解向量個數的問題