数学建模社区-数学中国

标题: 大家好我想问如何在解集中保留实数解剔除复数解 [打印本页]
作者: lhrman    时间: 2015-5-17 22:44
标题: 大家好我想问如何在解集中保留实数解剔除复数解
[img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\289623144\QQ\WinTemp\RichOle\N`H1F`X@4FAYE0A41M_[O~3.jpg[/img][img=0,1]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\289623144\QQ\WinTemp\RichOle\N`H1F`X@4FAYE0A41M_[O~3.jpg[/img]quen := {s0*s2 = 3/8, s0+s1+s2 = 1, 2*s0*s1+2*s1*s2 = -9/2}；计算三元方程组的解
, T, I3 S7 a4 j, G) K7 {
) G0 A' Q. m0 f* C
: d! y! |0 o2 m* Q: |# a9 B# \5 wquen2 := fsolve(quen);  用fsolve只能得到一个解# f+ s% j# K" F2 U. d
/ m- w" i* l7 ]/ L) l
quen3 := allvalues(solve(quen));可以将所有的解求出，但是其中涵盖太多的复数解，我们只需要实数解（实数解不只一个）  Q; D, U- X: l  i" f( D
: S3 s* f  }# y2 i: m4 f
map(subs, [quen3], quen)；验证方程的解，quen3中所有的解都是方程的解。
7 j4 M' `/ |; v# Z& x+ o# u, V7 d5 `  m) |* Z( u
我想求出所有解之后把复数解都删掉4 m* ]. ?' M7 p( H- K
! h( E. i2 Q9 E3 C* V
作者: haoyongle    时间: 2015-5-18 15:36
quen3(imag(quen3)~=0)=[];
* x2 C* D& }3 j



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5