作者TriQ ()
看板MATLAB
标题[心得] solve函数使用注意
时间Tue Oct 30 20:49:02 2012
今天在用solve函数解方程式，遇到一些问题，找到下面资料解决，po给大家参考。
原文：
http://www.chinavib.com/thread-53383-1-1.html
1.函数为符号函数
(当a,b,c为符号时,运算同样成立) 例题:
>> solve('a*x^2+b*x+c=0')
ans =
1/2/a*(-b+(b^2-4*a*c)^(1/2))
1/2/a*(-b-(b^2-4*a*c)^(1/2))
2.一般默认x为指定未知量；也可特指,如：特指a为未知量
>>solve('a*x^2+b*x+c=0','a')
ans =
-(b*x+c)/x^2
3.如果a,b,c有赋值，加入subs或eval语句
例题1
>> a=1;b=3;c=2;subs(solve('a*x^2+b*x+c=0'))
ans =
-1
-2
例题2
>> a=[1;1];b=[3;2];c=[2;1];
>> for m=1:2
>> eval(solve('a(m)*x^2+b(m)*x+c(m)=0'))
>> end
ans =
-1
-2
ans =
-1
-1
4.运算中出现复数 一方程如下，其中除xc外都已经赋值，xc为未知量
>>xc=subs(solve('As*fy=fc*b*xc*(xc*ey/(e0*(h0-xc))-(xc*ey/(h0-xc))^2/(3*e0^2))
+Es*As*(xc-as)/(h0-xc)*ey'));
得到的结果:
xc =
1.0e+002 *
1.7517 - 0.0000i
-1.0593 - 0.0000i
0.6221 + 0.0000i
如果复数值非常小的话，可以忽略，根据实际需要取所需实数。
例如：取第三个值的实数部分
>>real(xc(3))
ans=62.21
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.41.178.71