数学建模社区-数学中国

标题: 大家好我想问如何在解集中保留实数解剔除复数解 [打印本页]

作者: lhrman    时间: 2015-5-17 23:36
标题: 大家好我想问如何在解集中保留实数解剔除复数解
[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};计算三元方程组的解
1 n# t+ ^, [  R2 n
; K( ?$ k/ m4 I! M
& I+ i7 u8 P' U- Y1 K4 G1 G8 Rquen2 := fsolve(quen);  用fsolve只能得到一个解) O0 }; M; x6 Q1 @' L

$ p7 M$ f2 N. `) Nquen3 := allvalues(solve(quen));可以将所有的解求出,但是其中涵盖太多的复数解,我们只需要实数解(实数解不只一个)
: K# k$ w$ i9 z2 b$ M0 N
  V- e5 n& ^9 }! ^3 F$ p8 imap(subs, [quen3], quen);验证方程的解,quen3中所有的解都是方程的解。
5 n4 W/ z# s; }0 n, j1 u; d0 [  @6 i2 T/ Y. x: Y% G
我想求出所有解之后把复数解. K7 \9 k- z% ?8 _  V

& x) O9 G9 F/ Z) o) n1 C( [




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