24小时热门版块排行榜    

查看: 2082  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 小木虫:) 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂(0805) +6 烟汐忆海 2026-03-16 15/750 2026-03-18 21:12 by Gaodh_82
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 0703化学调剂,求各位老师收留 +9 秋有木北 2026-03-14 9/450 2026-03-18 19:49 by macy2011
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 0703化学调剂 ,六级已过,有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见