如何用matlab解方程組的解析解

時間 2022-01-13 05:55:06

1樓:匿名使用者

用solve函式

先舉一例,解方程"x^2+100*x+99=0"

在matlab 」command window"中輸入如下命令:

x=solve('x^2+100*x+99=0','x')首先來求一個二元一次方程組

9x+8y=10 式1

13x+14y=12 式2

[x,y]=solve('9*x+8*y=10','13*x+14*y=12','x','y')

如何用matlab解方程組的解析解

2樓:年智茂賦

用solve函式

先舉一例,解方程"x^2+100*x+99=0"

在matlab

」command

window"中輸入如下命令:

x=solve('x^2+100*x+99=0','x')首先來求一個二元一次方程組

9x+8y=10

式113x+14y=12

式2[x,y]=solve('9*x+8*y=10','13*x+14*y=12','x','y')

如何用matlab解方程組?

3樓:程順寬

s=solve('2*x*y=1,x+2=y+z,x+y-z=4','x,y,z'); %前面的參復

數是方程制組列表,後面是未知變數列bai表s.x %輸出未知數

dux的值

zhis.y %輸出未知數y的值

s.z %輸出未知數z的值

f=@(x)2*x; %定義一個匿名函式y=2x,其中@(x)表示x是匿名函式的

dao自變數

fplot(f,[-10,10]) %畫圖函式,第一個引數是函式名,第二引數要畫函式的區間[-10,10]

4樓:郯璇俟玟玉

沒有解抄析解,可bai

嘗試du數值解zhi

法dao

>>[a,b,c,d,e,f,g,h,k]=solve('g+4*k+8.314*log(a)=-19.3','2*k+h+8.

314*log(b)=192.6','g+h+8.314*log(c)=200.

6','g+2*h+8.314*log(d)=395.9','2*k+8.

314*log(e)=0','a+c+d-2*f^(-1)=0','4*a+2*b+2*e-14*f^(-1)=0','b+c+2*d-3*f^(-1)=0','a+b+c+d+e=1')

warning:

explicit

solution

could

notbe

found.

>insolve

at98a=

[empty

sym]b=

c=d=

e=f=

g=h=k=

5樓:老語開悅遠

>>symsa%

定義a為符號變數

f=4*cos(a)-2*sin(a)-6*exp(-2*a);

%方程內關係式

容solve(f)

%解方程

ans=

.86889441616212189601300889832946+0.*i

如何用matlab求下面方程組的解析解

6樓:匿名使用者

clear

clcsyms iq ir p1 q1 u1 r xf1=((p1-(p1^2+q1^2)*r/u1^2)/sqrt((u1+(p1*r+q1*x)/u1)^2+(p1*x-q1*r)^2/u1^2))-iq;

f2=((q1-(p1^2+q1^2)*x/u1^2)/sqrt((u1+(p1*r+q1*x)/u1)^2+(p1*x-q1*r)^2/u1^2))-ir;

[p1,q1]= solve(f1,f2,p1,q1)你要的解沒有意義。

兩個解的解析表示式在word裡佔了a4紙100多頁。

如果令iq=ir=0,解就簡單的多。

求高手解答用matlab 解複雜方程組的解

7樓:叢葉

syms x y z a b c

[x,y,z]=solve('72.8*(1+cos(a))=2*(sqrt(x*21.8)+sqrt(25.

5*y)+sqrt(25.5*z))','50.8*(1+cos(b))=2*sqrt(x*50.

8)','48*(1+cos(c))=2*(sqrt(x*29)+sqrt(47*y)+sqrt(1.92*z))','x,y,z')

x =12.700000000000000000000000000000+25.400000000000000000000000000000*cos(b)+12.

700000000000000000000000000000*cos(b)^2

y =-.98416931170728906540241988580218*cos(c)+9.7643439805890752790240685157559*cos(a)*cos(b)-23.

462398366292790782350185322144*cos(c)*cos(b)+7.1488996766620738069752563394849*cos(b)^2+3.3341639162497969191291895638814*cos(a)^2+.

40958164398297044880227538559881*cos(a)+.59974496762043211062439587258181*cos(b)-16.023090169105698668480172257920*cos(c)*cos(a)+19.

250659611845600192713968847131*cos(c)^2+.12578649948056747012125686189223e-1

z =.15194565438596826916006412565125e-193*(.58477278352492431981436322974990e98+.

58477278352492431981436322974990e98*cos(a)-.26731041191058757560106827855831e98*((cos(b)+1.)^2)^(1/2)-1.

*(-.64771139107889760381210034820223e194*cos(c)+.64262081202967154457022801098613e195*cos(a)*cos(b)-.

15441309237244940736312726425476e196*cos(c)*cos(b)+.47049056490306928474554525841306e195*cos(b)^2+.21943134403702282352778341593461e195*cos(a)^2+.

82783874266682206962299704487537e192+.26955798481907364090410143645319e194*cos(a)+.39471018111316040430045885264655e194*cos(b)-.

10545277016218098275353846992100e196*cos(c)*cos(a)+.12669437431192070703927236534687e196*cos(c)^2)^(1/2))^2

______________________________

就是代入求值就行了,例如:(接著上面的求x)

x=subs(x,,)

x =127/10+127/5*cos(333/10)+127/10*cos(333/10)^2

__________________________________

其實在一開始就給定a,b,c的值也是可以的,如果不要求一定要求解析解的話,樓下說的fsolve倒也是一種選擇

8樓:

用fsolve函式可以解出數值解

matlab 解析解轉換數值解

9樓:匿名使用者

subs 替換變數為數值

double 將解析解符號轉變為雙精度數值

10樓:

你可能誤解了解析解與數值解的**。

solve給你的是解析解,你得到解析解以後,給引數賦值,然後再使用eval()函式,即可得到結果(但不是你說的數值解)。

如果需要數值解,用積分來說吧,我們有辛普森(simposon),trapzoidal等方法,解方程還有gradient based等。

這些matlab也有相關的函式。

建議你看一本書,叫做:computaitonal methods for engineers,裡面都是用matlab實現數值解的方法。

11樓:

你可以自己編個演算法,比如計算方法裡講過的二分法和切線法都是數值解的解法

12樓:匿名使用者

vpa(解析解)=數值解

13樓:匿名使用者

都有解析解了你還要數值解幹嘛?把變數帶進去不就是數了麼,不理解你是啥意思。

用matlab解方程組,無解析解,求助!

14樓:我行我素

你的算式有問題,三個方程式

vm1+y0=a0*exp(-b*(vt1+t0)),vm2+y0=a0*exp(-b*(vt2+t0)),vm3+y0=a0*exp(-b*(vt3+t0))

把已知代入後,0.0149+y0=a0*exp(-0.125*(0.

0044+t0)),0.0136+y0=a0*exp(-0.125*(0.

0103+t0)),0.0130+y0=a0*exp(-0.125*(0.

0162+t0))

1式減2式,2式減3式得:0.0013=a0*[exp(-0.

125*(0.0044+t0))-exp(-0.125*(0.

0103+t0))],0.0006=a0*[exp(-0.125*(0.

0103+t0))-exp(-0.125*(0.0162+t0))]

兩式相除,左邊=2.1667,

右邊=[exp(-0.125*(0.0044+t0))-exp(-0.

125*(0.0103+t0))]/[exp(-0.125*(0.

0103+t0))-exp(-0.125*(0.0162+t0))]含t0項約掉

=[exp(-0.125*0.0044)-exp(-0.

125*0.0103)]/[exp(-0.125*0.

0103)-exp(-0.125*0.0162)]=1.

007左邊≠右邊

matlab解方程,數值解或解析解

15樓:嚕嚕晗寶

用solve()函式可以求得解析解。使用方法如下:

syms n0 n2 n1 r phi

solve(r-((n0-n2)^2*(cos(phi/2))^2+(n0*n2/n1-n1)^2*(sin(phi/2))^2)/((n0+n2)^2*(cos(phi/2))^2+(n0*n2/n1+n1)^2*(sin(phi/2))^2),'n1')

ans =

1/(-2*tan(1/2*phi)^2+2*r*tan(1/2*phi)^2)*tan(1/2*phi)*(-(-2+2*r)*(r*n0^2+r*n2^2+2*r*tan(1/2*phi)^2*n0*n2+2*r*n0*n2-n0^2+2*n0*n2-n2^2+2*tan(1/2*phi)^2*n0*n2-(6*n0^2*n2^2+6*r^2*n0^2*n2^2+8*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^2*n0^2*n2^2+r^2*n0^4-2*r*n0^4+r^2*n2^4-2*r*n2^4-4*n0^3*n2-4*n0*n2^3+4*r^2*n0^3*tan(1/2*phi)^2*n2+4*r^2*n2^3*tan(1/2*phi)^2*n0+8*r^2*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^4*n0^2*n2^2+4*r^2*n0^3*n2+4*r*n0^2*n2^2+4*r^2*n2^3*n0-4*n0^3*tan(1/2*phi)^2*n2-4*n2^3*tan(1/2*phi)^2*n0+n0^4+n2^4)^(1/2)))^(1/2)

-1/(-2*tan(1/2*phi)^2+2*r*tan(1/2*phi)^2)*tan(1/2*phi)*(-(-2+2*r)*(r*n0^2+r*n2^2+2*r*tan(1/2*phi)^2*n0*n2+2*r*n0*n2-n0^2+2*n0*n2-n2^2+2*tan(1/2*phi)^2*n0*n2-(6*n0^2*n2^2+6*r^2*n0^2*n2^2+8*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^2*n0^2*n2^2+r^2*n0^4-2*r*n0^4+r^2*n2^4-2*r*n2^4-4*n0^3*n2-4*n0*n2^3+4*r^2*n0^3*tan(1/2*phi)^2*n2+4*r^2*n2^3*tan(1/2*phi)^2*n0+8*r^2*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^4*n0^2*n2^2+4*r^2*n0^3*n2+4*r*n0^2*n2^2+4*r^2*n2^3*n0-4*n0^3*tan(1/2*phi)^2*n2-4*n2^3*tan(1/2*phi)^2*n0+n0^4+n2^4)^(1/2)))^(1/2)

1/(-2*tan(1/2*phi)^2+2*r*tan(1/2*phi)^2)*tan(1/2*phi)*(-(-2+2*r)*(r*n0^2+r*n2^2+2*r*tan(1/2*phi)^2*n0*n2+2*r*n0*n2-n0^2+2*n0*n2-n2^2+2*tan(1/2*phi)^2*n0*n2+(6*n0^2*n2^2+6*r^2*n0^2*n2^2+8*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^2*n0^2*n2^2+r^2*n0^4-2*r*n0^4+r^2*n2^4-2*r*n2^4-4*n0^3*n2-4*n0*n2^3+4*r^2*n0^3*tan(1/2*phi)^2*n2+4*r^2*n2^3*tan(1/2*phi)^2*n0+8*r^2*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^4*n0^2*n2^2+4*r^2*n0^3*n2+4*r*n0^2*n2^2+4*r^2*n2^3*n0-4*n0^3*tan(1/2*phi)^2*n2-4*n2^3*tan(1/2*phi)^2*n0+n0^4+n2^4)^(1/2)))^(1/2)

-1/(-2*tan(1/2*phi)^2+2*r*tan(1/2*phi)^2)*tan(1/2*phi)*(-(-2+2*r)*(r*n0^2+r*n2^2+2*r*tan(1/2*phi)^2*n0*n2+2*r*n0*n2-n0^2+2*n0*n2-n2^2+2*tan(1/2*phi)^2*n0*n2+(6*n0^2*n2^2+6*r^2*n0^2*n2^2+8*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^2*n0^2*n2^2+r^2*n0^4-2*r*n0^4+r^2*n2^4-2*r*n2^4-4*n0^3*n2-4*n0*n2^3+4*r^2*n0^3*tan(1/2*phi)^2*n2+4*r^2*n2^3*tan(1/2*phi)^2*n0+8*r^2*tan(1/2*phi)^2*n0^2*n2^2+16*r*tan(1/2*phi)^4*n0^2*n2^2+4*r^2*n0^3*n2+4*r*n0^2*n2^2+4*r^2*n2^3*n0-4*n0^3*tan(1/2*phi)^2*n2-4*n2^3*tan(1/2*phi)^2*n0+n0^4+n2^4)^(1/2)))^(1/2)

如何用ecel解方程組,如何用Excel解方程組?

以解二元一次方程組為例 a1 x b1 y c1 a2 x b2 y c2 在excel中,a1 c1分別為a1 b1 c1,a2 c2分別為a2 b2 c2 計算a4 mdeterm a1 b2 若不為0,則方程組有唯一解 若為0,則方程組有多個解 若方程組有唯一解,則 在a6中輸入 minver...

解方程組解不等式組

1 x 3y 1 3x 2y 8 從 得x 1 3y 移項,即同時減去3y可得 將x代入 y 1 合併同類項,並同時除以11可得 x 2 將x代入y可得 2 6x 15 2 4x 3 2x 1 3 x 2 2 3 從 得6x 15 8x 6 x 2 9 從 得2 2x 1 3x 4 同時乘以6 x ...

用matlab解方程組的問題,方程求解未知數

把係數寫成矩陣,然後直接用等號右邊的矩陣除以係數矩陣,或者直接用solve函式,看看幫助,裡面寫的很詳細 k打醬油 方程數大於未知量個數,屬於超定方程組,一般沒有精確解,只能求最小二乘解。matlab中可以利用左除或偽逆法求解。clear c1 10496.108 c2 10500.275 c3 1...