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...