| 查看: 1199 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
Mathematica求助 泰勒展开和求导 已有1人参与
|
|||
|
求教Mathematica： 多元函数展开后(Series[])，比如x展开到m次，y展开到n次，怎么只保留m+n<4的项？ 多谢！！ |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个关于 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
- 专业: 粒子物理学和场论
|
对这种多项式处理不是很熟 另外对模式的理解个人还不是很深入 只能用比较笨的方式 下面的代码是比较笨的方法解决你这个问题 例子是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
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
- 专业: 粒子物理学和场论
|
上面包含了i] 出错 改成k 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[#[[k]]^ToExpression[("n" <> ToString[k] <> "_." ], {k, 1, Length[#]}]) /; Evaluate[ Total[Table[ ToExpression[("n" <> ToString[k])], {k, 1, Length[#]}]] >= n]) -> 0 &, jihe] Do[ test = ExpandAll[Normal[Series[test, {x[[k]], 0, n - 1}]]]; test = test /. rule; , {k, 1, Length[x]}] test |
5楼2015-05-05 22:07:11
25