LU法求解線性方程組,matlab程式設計

2023-01-23 08:30:40 字數 3194 閱讀 8846

1樓:

a=[1 4 0 1;1 5 1 0;-3 2 0 3;-4 0 1 4];

b=[11 12 7 5]'; %

[l,u]=lu(a); %計算上三角和下三角x=u\(l\b) %求解x

matlab中怎樣利用lu分解求解方程組

2樓:匿名使用者

如果a,b是引數,那麼x,y是未知數一個方程兩個未知數,有無限多的解啊換句話說是二維座標系中符合x^2=y^3-2^a*7^b的曲線上的每一點的(x,y)都是解

如何用matlab求解線性方程組

3樓:匿名使用者

給你舉個例子,已知兩個方程表示式,可以利用solve函式求解

[x,y] = solve('-0.43*x -0.67*y -90=0','-0.32*x-0.38*y-280=0')

用matlab程式設計求解方程組 100

4樓:匿名使用者

求解方程組的執行結果找不到解的主要問題有:

1、用fsolve函式求解的初值問題,解決的方法是通過調整獲得。

2、求和的問題,在呼叫myfun(x)時,始終在原地求解當i等於某值時的y值,而沒有累加求和。

按上述方法,修改後執行可求解得到

x1=-0.78894,x2=0.1253,f1=1.7053e-13,f2=1.6342e-13

matlab程式設計解方程組。。寫出matlab程式

5樓:工科天蠍男

a=[3,2,5;-1,4,3;1,-1,3];

b=[6;5;1];

x=inv(a)*b

結果:x =

0.5000

1.0000

0.5000

即用矩陣的乘除法(這裡用的是逆矩陣)來計算。

或者用matlab自己的solve函式來計算syms x1 x2 x3;

[x1,x2,x3]=solve('3*x1+2*x2+5*x3=6','-x1+4*x2+3*x3=5','x1-x2+3*x3=1')

結果為:

x1 =

1/2x2 =

1x3 =1/2

用matalb編寫程式,用追趕法求解三對角線性方程組:

6樓:青山客_李文楊

function x=trid(a,b,c,d)% 追趕法求解三對角的線性方程組 ax=d% b為主對角線元素,a,c分別為次對角線元素,d為右端項% a=[ b1 c1

% a2 b2 c2

% ......

% a_(n-1) b_(n-1) c_(n-1)

% a_(n) b_(n) ]

% b=[b1...b_(n)]

% a=[0 a2...a_(n)]

% c=[c1...c_(n-1)]

n=length(b);

u(1)=b(1);

for i=2:n

l(i)=a(i)/u(i-1);

u(i)=b(i)-l(i)*c(i-1);

endy(1)=d(1);

for i=2:n

y(i)=d(i)-l(i)*y(i-1);

endx(n)=y(n)/u(n);

for i=n-1:-1:1

x(i)=(y(i)-c(i)*x(i+1))/u(i);end

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

7樓:刺友互

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)求得的線性方程組解是一樣的。

8樓:匿名使用者

%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=rbendend

求解線性方程組

首先說明這個方程組的解不是唯一的,因為四個未知數卻只有三個方程組,一般情況下是肯定沒有唯一解的。用線性代數的方法求解如下,主要就是用到矩陣的初等行變換,樓主想必應該知道 否則也不會做到這道題吧 可以參考 大學數學基礎教程 三 線性代數與空間解析幾何 1,將該方程組的增廣矩陣為 a b 0 1 3 1...

數學求解線性方程組的通解,求線性方程組的通解請寫下過程謝謝

增廣矩陣抄經襲行變換化成bai 字數限制du 1 0 3 7 13 7 13 7 0 1 2 7 4 7 4 7 0 0 0 0 0 通解zhi為 13 7,4 7,0,0 c1 3,2,7,0 c2 13,4,0,7 c1,c2 為任意dao常數 求線性方程組的通解 請寫下過程謝謝 方程來組的通解...

關於線性方程組XA O

首先你要明確一個問題x是嚴格來說是一個列向量,一個列向量左乘以一個與之同型號的向量是啥。你自己琢磨琢磨,全國人民都知道。求基礎解析的方法 首先你要明確一個問題,基礎解析是一個極大無關組,而且可以作為ax 0的部分解 由於是極大無關的,基礎解析是線性無關的,既然是線性無關的,ax 0的任意解都可以由基...