SOLVE函数使用几点注意

jinghelen

用Matlab编程遇到问题，来到振动论坛的Matlab版，版主和其它版友的热心帮助让我十分的感动，现将自己问的问题，以及通过实践得到的一些心得与大家分享。

solve关于解方程的一些技巧。

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

donkeyxu

好！

花如月

总结得不错，对初学者比较有用。
---------------------------------------------------------
本版面一向鼓励原创和经验分享
---------------------------------------------------------