本帖最后由 lijinfeng042 于 2010-11-3 19:26 编辑
使用Mupad可以非常方便的求解非线性方程的数值解,简化输入的步骤,同时添加各种限制~只要有一下几个函数numeric::linsolve, numeric::polyroots, numeric::polysysroots, numeric::realroot, numeric::realroots, numeric::solve
注意:以下例子均在Mupad环境进行,如需在Command Window 运行,请作如下处理
M='代码'
result = evalin(symengine,M)
%这样调用也是有缺点的,就是少数命令不支持...
主要有以下几个方面内容,呵呵 做个目录
float 对结果求近似值
numeric::solve
numeric::fsolve
numeric::linsolve
numeric::polyroots
numeric::polysysroots
numeric::realroot
numeric::realroots
注意:在数值解法里,将会忽略assume及assuming对未知数的条件限制,例如x>0将不可以用assume(x>0)表示 可以在区间里限定x=0..infinite
float 对结果求近似值
这个很基本的,类似double, 用法 float(object, n)
float(solve([x^3 + x^2 + 2*x = y, y^2 = x^2], [x, y]))复制代码