| 查看: 1202 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
Mathematica求助 泰勒展开和求导 已有1人参与
|
|||
|
求教Mathematica： 多元函数展开后(Series[])，比如x展开到m次，y展开到n次，怎么只保留m+n<4的项？ 多谢！！ |
回复此楼» 猜你喜欢
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
感觉是下周放榜了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个关于 quadgk 函数求积分的问题
已经有7人回复
- Matlab 7.0 基础教程 清华大学
已经有29人回复
- 如何求解下列方程：xlnx=a(x-b)
已经有14人回复
- [求助]用Mathematica软件进行taylor展开
已经有3人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
运行结果 a - a^3/6 + b - (a^2 b)/2 - (a b^2)/2 - b^3/6 + c - (a^2 c)/2 - a b c - (b^2 c)/2 - (a c^2)/2 - (b c^2)/2 - c^3/6 + d1 - (a^2 d1)/2 - a b d1 - (b^2 d1)/2 - a c d1 - b c d1 - (c^2 d1)/2 - (a d1^2)/2 - ( b d1^2)/2 - (c d1^2)/2 - d1^3/6 + d2 - (a^2 d2)/2 - a b d2 - ( b^2 d2)/2 - a c d2 - b c d2 - (c^2 d2)/2 - a d1 d2 - b d1 d2 - c d1 d2 - (d1^2 d2)/2 - (a d2^2)/2 - (b d2^2)/2 - (c d2^2)/2 - ( d1 d2^2)/2 - d2^3/6 + d3 - (a^2 d3)/2 - a b d3 - (b^2 d3)/2 - a c d3 - b c d3 - (c^2 d3)/2 - a d1 d3 - b d1 d3 - c d1 d3 - ( d1^2 d3)/2 - a d2 d3 - b d2 d3 - c d2 d3 - d1 d2 d3 - (d2^2 d3)/2 - ( a d3^2)/2 - (b d3^2)/2 - (c d3^2)/2 - (d1 d3^2)/2 - ( d2 d3^2)/2 - d3^3/6 + d4 - (a^2 d4)/2 - a b d4 - (b^2 d4)/2 - a c d4 - b c d4 - (c^2 d4)/2 - a d1 d4 - b d1 d4 - c d1 d4 - ( d1^2 d4)/2 - a d2 d4 - b d2 d4 - c d2 d4 - d1 d2 d4 - (d2^2 d4)/2 - a d3 d4 - b d3 d4 - c d3 d4 - d1 d3 d4 - d2 d3 d4 - (d3^2 d4)/2 - ( a d4^2)/2 - (b d4^2)/2 - (c d4^2)/2 - (d1 d4^2)/2 - (d2 d4^2)/2 - ( d3 d4^2)/2 - d4^3/6 没好好看 应该是没问题 |
7楼2015-05-05 22:10:28
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
2楼2015-04-23 06:46:30
3楼2015-04-27 17:45:06
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
对这种多项式处理不是很熟 另外对模式的理解个人还不是很深入 只能用比较笨的方式 下面的代码是比较笨的方法解决你这个问题 例子是7个变量泰勒展开，只取总的次数小于5的项。 个人感觉mathematic应该有更简单的方式只需要构造单个的rule 挑选出所需的项，不过我没找到这样的函数，只好用字符的方式和多个rule来实现这种笨挑选。 f[x__] := Sin[Total[x]]; x = {a, b, c, d1, d2, d3, d4}; n = 5; test = f[x] jihe = Delete[Subsets[x], 1]; rule = Map[((Times @@ Table[#[]^ToExpression[("n" <> ToString <> "_."  ], {i, 1, Length[#]}]) /; Evaluate[ Total[Table[ ToExpression[("n" <> ToString)], {i, 1, Length[#]}]] >= n]) -> 0 &, jihe] Do[ test = ExpandAll[Normal[Series[test, {x[], 0, n - 1}]]]; test = test /. rule; , {i, 1, Length[x]}] test |
4楼2015-05-05 22:05:30
