appla 发表于 2006-10-4 11:30

请教关于符号函数不等式的问题

我现在做一个程序,用routh判剧,得到了一些值.b()是一个6*3的矩阵
现在要求说b(4,1)*b(4,2)>0;
b(4,1)=((5/3+(-4460109110491993/422212465065984+1011723334372798245593504385912315/162259276829213363391578010288128*pi)*T^2-5488135394582167/316659348799488*T)/(5/3*T+1)*(-707527188533663/140737488355328*T+760303746666911/70368744177664+(1/12+2111066468319334115083237618823685/162259276829213363391578010288128*pi)*T+52004897100201825/4503599627370496*pi)+41430476124406897618513672539851085/162259276829213363391578010288128*pi*T^2-(-46672912408830829681098729181745997/162259276829213363391578010288128*pi+10788829688449105/1266637395197952)*T-2157765937689821/422212465065984-2111066468319334115083237618823685/162259276829213363391578010288128*pi)/(5/3+(-4460109110491993/422212465065984+1011723334372798245593504385912315/162259276829213363391578010288128*pi)*T^2-5488135394582167/316659348799488*T)*(5/3*T+1)
b(4,2)也跟它差不多
其中只有一个符号变量T
设 syms T
我想用 maple('solve(b(4,1)*b(4,2)>=0,T)')来做
但是运行的结果是' '
不知道为什么会这样?
应该怎么多呢?

appla 发表于 2006-10-4 11:31

望高手指教一下,做不出来就毕不了业了

buaalqq 发表于 2006-10-5 11:30

试一下把语句maple('solve(b(4,1)*b(4,2)>=0,T)')改成maple('solve(b(4,1)*b(4,2)>=0,T);'),这是maple语法格式规定的,不这样偶尔会引起运行停滞,另外我觉得把b(4,1)和b(4,2)替换成它们代表的表达式吧。
我也没试过,不知道行不行。希望有帮助。

appla 发表于 2006-10-6 18:38

呵呵,多谢多谢,改为b(4,1)和b(4,2)的表达式后果然能出结果了.真是太谢谢了,帮了大忙了!
页: [1]
查看完整版本: 请教关于符号函数不等式的问题