数学建模社区-数学中国
标题:
大家好我想问如何在解集中保留实数解剔除复数解
[打印本页]
作者:
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}；计算三元方程组的解
9 q2 [' Q0 F e7 ?/ U4 e. ?
) y0 p- G" C. C# Q% t
6 a' E, i9 ?& w Y2 j, a$ e s/ b8 V
quen2 := fsolve(quen); 用fsolve只能得到
一个解
+ R7 y n. h* G( E5 x# w
B z5 q+ x! \( ` g
quen3 := allvalues(solve(quen));可以将所有的解求出，但是其中涵盖太
多的复数解
，我们只需要实数解（
实数解不只一个
）
& X$ L3 @7 ~9 H- [" N
) N: ~* v6 O' D3 R& z5 [$ I
map(subs, [quen3], quen)；验证方程的解，quen3中所有的解都是方程的解。
1 p' W: z! [7 k: O1 _
# R# Q: U* }4 d8 W3 g$ l, ~
我想求出所有解之后把复数解都删掉
+ Z4 G- `, b6 ^) J0 n
: E% G8 m; l Y+ e
作者:
haoyongle
时间:
2015-5-18 15:36
quen3(imag(quen3)~=0)=[];
; W. u/ S- [( {
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5