Maple 公式推导
2014-12-24 16:05
417 查看
Maple 公式推导
2012-03-24 13:07:35| 分类:学习 | 标签:maple
|举报
|字号 订阅
即以此功德,庄严佛净土。上报四重恩,下救三道苦。惟愿见闻者,悉发菩提心。在世富贵全,往生极乐国。
阿英按:个人认为公式推导最主要的就是“代入”和“化简”两步,大家抓住这两步其实就够用了,其他的用到了再查啊,阿弥陀佛!
① 选择哪个函数进行变量替换?
Maple中是subs函数,algsubs函数和simplify函数。
原来我不知道要把一个表达式带入另外一个表达式要用algsubs函数,subs函数仅仅只能替换一个变量而已,比如z=a+b, 如果x+y是一个通项中的operand,则有x + y = a + b + c ;algsubs一次能替换一个表达式中的某些中间项,比如algsubs( a+b=d, 1+a+b+c );中a+b就是表达式1+a+b+c中的中间项。
② 公式推导中要使用哪些函数?
factor, combine, expand, subs,algsubs,simplify, convert,alias
③ 变量的约束条件如何加?
变量的约束条件由assume给出,如assume(x>0); 经过assume后的x后边跟了一个flag,显示为"x~"; about(x)是查询x的约束条件是什么;在assume的基础上附加一个assumption,通过additionally函数,例如additionally(x<2); 去除所有约束采用x := 'x';
例1: subs 单变量替换
p := (x+y)^2 + 1/(x+y)^2;
normal(p); (* normal(f) The expression f is converted to factored normal form. This is the form numerator/denominator, where the numerator and denominator are relatively prime polynomials with integer coefficients. *)
我们希望p的分子上包含(x+y)^4的项,我们可以令 z = (x + y),然后代入p的表达式。
pz := subs(x+y=z,p); (*替代规则是将 x+y = z 等式左手边的项x+y用右手边的项z替代,然后代入多项式p中*)
npz := normal(pz);
q := subs(z = x + y, npz);
(*下面做一个简要的说明,关于op函数的用法*)
op([2,1], q); (*op是取出表达式q中的operand,前面[2,1]是个列表,表示以q的第二个操作数为新的子表达式,取出该子表达式中的第一个操作数,具体分解见下面执行结果*)
(* 言归正传:subsop - substitute for specified operands in an expression *)
subsop([2,1] = x-y,q); (* 将q中嵌套的操作数 [2,1] = (x+y) 用(x-y)替换 ,嵌套的operand用[2,1]表示我以前不知道 *)
例2: subs 嵌套替换和多个变量同时替换
p := a + 2*b + 3*c; (*定义一个新的多项式*)
subs(a=b, b=c, c=a, p); (*我以前不知道原来是从左到右依次嵌套替换,此句的效果与 subs(c=a,subs(b=c,subs(a=b,p))); 等价*)
(* The substitutions within a set or list are performed simultaneously,如果要同时替换p中的多个变量,需要将各个替换等式用花括号或者方括号括起来 *)
subs({a=b, b=c, c=a}, p);
subs([a=b, b=c, c=a], p);
例3:如何用一个子表达式替换表达式中的一部分,将一个表达式带入另外一个表达式不能用subs函数,而是要用algsubs函数。
p := a + 2*b + 3*c;
subs(a+b=c,p);
algsubs(a+b=c,p);
例4: expand 命令,matlab中也由此命令expand命令,matlab中 expand((x+1)^3) returns x^3+3*x^2+3*x+1
p := (a+b+c)*(x^3 + 9*x + 8);
expand(p); (* expand命令的作用是将多项式乘开 *)
(* 下面的内容演示如何在expand过程中保持第二个operand 不变 *)
op(2,p); (* p的第二个因子 ,第二个operand,这里打印出来看看是啥东东 *)
q := subsop(2=z,p); (* 将p的第二个operand用z代替 *)
eq := expand(q);
eq := subs(z=op(2,p),eq);
p2 := op(2,p);
irreduc(p2); (* p2是不可化简的 *)
alias(alpha=RootOf(p2)); % alias就是变量替换,这个在化简的过程中也有用
ap := (alpha + 1)^7;
expand(ap);
evala(%);
例5. factor的数域扩展,factor仅仅分解到有理数就停止了,下面的例子演示如何通过数域扩展,在无理数域进行 因式分解
p := x^2 + 2;
factor(p);
factor(p,{sqrt(2),I});
fp := polytools[split](p,x); % polytools[split] - complete Factorization of a Polynomial
rfp := convert(fp,radical);
map(evalc,rfp); % 我以前不知道要map一下,注意与evalc(rfp);的结果有差异
evalc(rfp);
相关文章推荐
- 利用Maple推导向量方程的微分公式
- LSTM公式推导csdn
- 信噪比(一些概念,公式推导,实验分析)
- 坐标旋转公式(推导)
- OPENGL NEHE Lesson11 11课的计算公式推导
- 机器学习公式推导
- 相机标定简介与MatLab相机标定工具箱的使用(未涉及原理公式推导)
- Lucene学习总结之六:Lucene打分公式的数学推导
- 深度学习算法之CNN、RNN、LSTM公式推导
- 反向传播算法(过程及公式推导)
- CF #278 (Div. 2) B.(暴力枚举+推导公式+数学构造)
- python 字典推导公式
- Stanford机器学习课程笔记——LR的公式推导和过拟合问题解决方案
- HDU 2298 Toxophily(公式推导)
- Lucene打分公式的数学推导
- 453. Minimum Moves to Equal Array Elements 数学公式推导
- 逻辑回归-参数迭代公式推导
- <深度学习系列>反向传播算法的公式推导
- 四旋翼飞行器旋转矩阵公式推导!
- 【快速幂 && 公式推导】LightOJ - 1282 Leading and Trailing