24小时热门版块排行榜    

查看: 2077  |  回复: 4

小木虫:)

荣誉版主 (著名写手)

[交流] 【求助】请教一个关于Mathematica中计算小量的问题 已有3人参与

小虫自从用上Mathematica之后发现非常不错,于是平时手算的内容很多都移给Mathematica了,最近小虫正在计算一些变分和微扰问题,这里面涉及一些0阶、一阶之类的问题,但是小虫发现用Mathematica很难处理这样的问题,小虫就想问问有没有什么命令或者思路可以求解一个方程的N阶近似呢?具体而言也就是

比如,如下一个问题



其中




这里的 都是小量

也就是



我希望求出F的零阶近似


和一阶近似


这个问题可能比较过分,但是小虫还是想问,电脑能不能做这样的工作?

[ Last edited by 小木虫 on 2011-1-20 at 22:51 ]
回复此楼
霸道做事,厚道做人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
余泽成(金币+3): 谢谢分享经验! 2011-01-21 13:29:05
类似于微扰问题,mathematica应该能做,不过我也在找路的困惑中...
以前关于方程的零级解、一级解都是手推的,那个累就不说了,还怕错
现在的方法是我提供求解的方法,mathematica去推
说实话,这方法很笨,几乎没有通用性,但是起码不怕推错了,聊胜于无吧
2楼2011-01-20 14:19:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gsview

木虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
余泽成(金币+3): 谢谢参与应助! 2011-01-21 13:29:23
对这个问题我可以给出一点小建议,通用性上还需要加强!聊胜于无!

A0 = Abar + x*dA;
B0 = Bbar + x*dB;
C0 = Cbar + x*dC;
F = A0*B0*C0;
F = Expand[F];
F = F /. {x^3 -> 0, x^2 -> 0};
F = F /. {x -> 1};
Print[F]
3楼2011-01-20 20:41:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫:)

荣誉版主 (著名写手)

谢谢两位大侠的回帖,我觉得gsview的方法很巧,可以实用,非常感谢,这下再算微扰问题可就容易了

程序可以这样编

A[x_]=A0+x*A1;
B[x_]=B0+x*B1;
H[x_]=H0+x*H1;
F=Expand[A[x]*B[x]*H[x]];
F0=F/.{x^3 -> 0, x^2 -> 0, x -> 0}
F1=F /. {x^3 -> 0, x^2 -> 0}

结果F0就是0级近似,F1就是1级近似。(这里我发现C[x]好象是个特殊字符,不宜使用,我用H替代了)

这个小程序应该很实用,可以减轻许多虫子计算的痛苦,谢谢咯
霸道做事,厚道做人
4楼2011-01-20 23:11:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫:)

荣誉版主 (著名写手)

★ ★
余泽成(金币+2): 鼓励交流,欢迎常来程序语言版! 2011-06-04 19:22:04
找到一个更好的方法了

A[x_]=A0+x*A1;
B[x_]=B0+x*B1;
H[x_]=H0+x*H1;
Fn=Normal[Series[A[x]*B[x]*H[x],{x,0,n}]];

Fn就是n级近似
霸道做事,厚道做人
5楼2011-06-03 11:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小木虫:) 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 6/300 2026-03-17 17:13 by 生物工程调剂
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 070300化学319求调剂 +3 锦鲤0909 2026-03-17 3/150 2026-03-17 15:01 by 我的船我的海
[考研] 289求调剂 +6 步川酷紫123 2026-03-11 6/300 2026-03-17 10:23 by Sammy2
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
信息提示
请填处理意见