关于Mathematica软件的命令使用的讨论

chyanog

mathematica 发表于 2013-11-4 14:33
第二个版本的代码
(*利用递归解决问题!*)
Clear["Global`*"];(*Clear all variables*)

这个实现不错，改写一下：

1. Clear["`*"];
   
2. fun[{x_, xs___}] := Prepend[fun /@ {xs}, fun@x]
   
3. fun[x_] := FromCharacterCode@x
   
4. si = {{65, {66, {70}, 69}}, 67};
   
5. fun /@ si
   

复制代码

chyanog

chyanog 发表于 2013-11-4 14:02

@Mathematica，似乎没什么解释的，主要就模式匹配和RuleBased，不熟悉的符号或函数，按F1看看，然后试几个例子就明白了，判断整数是或小数的试试IntegerQ，NumericQ，Element，NumberQ等

mathematica

chyanog 发表于 2013-11-4 16:08
@Mathematica，似乎没什么解释的，主要就模式匹配和RuleBased，不熟悉的符号或函数，按F1看看，然后试几 ...

@chyanog
代码1:
If[ListQ@#,#,FromCharacterCode@#]&//@{{65, {66, {70}, 69}}, 67}
代码2:
If[Not@IntegerQ@#,#,FromCharacterCode@#]&//@{{65, {66, {70}, 69}}, 67}

你觉得这个如何?

mathematica

chyanog 发表于 2013-11-4 16:08
@Mathematica，似乎没什么解释的，主要就模式匹配和RuleBased，不熟悉的符号或函数，按F1看看，然后试几 ...

不可否认,这个也马马虎虎是一种新的代码,虽然是"抄袭"你的

wayne

mathematica 发表于 2013-11-4 16:29
@chyanog
代码1:
If&//@{{65, {66, {70}, 69}}, 67}

有点大动干戈了吧.
本质上就是让FromCharacterCode 具有 Listable 的属性.所以:

1. SetAttributes[FromCharacterCode, Listable];
   
2. FromCharacterCode[si]

复制代码

wayne

zeroieme 发表于 2013-11-4 13:25
借楼提问：
怎么让Mathematica求多元多项式的指数
比如a x y + x 以x、y为未知数。Exponent[a x y + x, { ...

既然Mathematica的 Exponent不支持多元函数,我们不妨自己写一个.
基本思路 就是获取所有的项,然后将所有自变量赋成一个变量,再用Exponent.

1. f = Max[Exponent[MonomialList[#] /. Thread[{x, y} -> t], t]] &

复制代码

MonomialList有一个用来排序的选项,可以按照 total degree来排序, 排序完第一个就是最高次项了.

1. MonomialList[#, {x, y}, "DegreeLexicographic"]

复制代码

mathematica

1. Map[If[ListQ@#, #, FromCharacterCode@#] &, {{65, {66, {70}, 69}},
   
2.   67}, Infinity]

复制代码

这也是一种办法!

zeroieme

wayne 发表于 2013-11-4 22:21
既然Mathematica的 Exponent不支持多元函数,我们不妨自己写一个.
基本思路 就是获取所有的项,然后将所 ...

但MonomialList容易溢出

wayne

zeroieme 发表于 2013-11-6 11:48
但MonomialList容易溢出

我测试了一下,发现3^2000项,前面提供的代码,MonomialList 运行需要55秒钟,
而下面的代码需要41秒钟

1. eqs = (a x y + y + x)^2000;
   
2. Total[First[CoefficientRules[eqs, {x, y}, "DegreeLexicographic"][[All, 1]]]]

复制代码

wayne

其实我们完全没必要展开所有的项.
只需要计算 每一个 参与乘积的因子的最高次幂就行了.