| 查看: 537 | 回复: 3 | |||
[交流]
求助mathematica的问题 已有2人参与
|
|
现在碰上一个问题，要解多个方程组，最后要求出clenda的表达式，表达式里含有fi。折腾了一天，没搞定，没办法，不知出什么问题，求助于版上高人。程序如下： Clear[eq1, eq2, eq3, eq4，c, fp, nos, nop, lendafp, fs, fi，deltak]; c = 3*10^8;(*c为光速*) fp = 555*10^12;(*fp为pump光的频率*) eq1 = Solve[{fp == c/(nop*lendafp), (lendafp*10^3)*(-0.8) + 2.74 == nop}, {nop, lendafp}](*nop为pump光的折射率，lendafp为pump光的波长*) eq2 = Solve[{fp == fs + fi, fi == c/(nei*lendafi), (lendafi*10^3)*(-0.03) + 2.185 == nei}, {nei, lendafi, fs}](*nei为idle光的折射率，lendafi为idle光的波长，fs为signal光的频率*) eq3 = Solve[{fp == fs + fi, fs == c/(nos*lendafs), (lendafs*10^3)*(-0.04) + 2.735 == nos}, {fs, nos, lendafs}] (*nos为signal光的折射率，lendafs为signal光的波长，fs为signal光的频率*) eq4 = Solve[{0 == (2 Pi/c)*(2.739* fp - (0.0025 (-3.03585*10^17 + 547 fi + Sqrt[ 9.2162786625*10^34 - 3.3212007*10^20 fi + 299209 fi^2]))/(-5.55*10^14 + fi)*(5.55*10^14 - fi) - ( 0.0025 (437 fi + Sqrt[fi] Sqrt[-1.44*10^15 + 190969 fi]))/fi* fi) - (2 Pi/clenda)}, fi] |
回复此楼» 猜你喜欢
一志愿同济大学323分（080500）求调剂
已经有3人回复
-
330分求调剂
已经有7人回复
-
336求调剂
已经有3人回复
-
285求调剂
已经有5人回复
-
调剂推荐
已经有9人回复
-
材料调剂
已经有6人回复
-
材料专硕322分
已经有3人回复
-
321求调剂
已经有5人回复
-
301求调剂
已经有6人回复
-
085600，321分求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个Mathematica积分的问题（重金求解）
已经有6人回复
- 【求助】MATLAB、Mathematica和Maple这三大数学软件的优缺点
已经有41人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】请教一个关于Mathematica中计算小量的问题
已经有4人回复
- 【求助】请问怎么求函数方程的函数表达式？
已经有11人回复
- 【求助】mathematica化简问题【已完结】
已经有15人回复
- 【求助】Mathematica怎么对不等式解的最小值集合作图？
已经有8人回复
- 【求助】关于Mathematica绘图问题【已完结】
已经有3人回复
- 【求助】Mathematica画图问题
已经有7人回复
- 【求助】先解常微分方程 然后画图【已完结】
已经有11人回复
- 【求助】求解非线性方程
已经有8人回复
- 【求助】Mathematica7.0
已经有12人回复
- 【求助】origin和mathematica的问题！
已经有6人回复
- 【求助】紧急求助有关MATHEMATICA的问题
已经有9人回复
- 【求助】如何提高Mathematica运行速度？
已经有6人回复
walk1997
主管区长
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
把程序里面的小数全部改写成分数，应该就可以了 如果是实在非常伏复杂的情况，只能用小数的话，就用函数定义来求解。你这个情况，用解析的速度也够快。。 我改写下，你检查下方程里面的数字有没有改错？ 不知道是不是有问题，fi的结果对clenda的取值是个复数，实部很大，虚部很小......... －----------------------------------------------------- Clear["Global`*"]; c = 3*10^8; fp = 555*10^12; eq1 = Solve[{fp == c/(nop*lendafp), (lendafp*10^3)*(-8/10) + 274/10 == nop}, {nop, lendafp}]; eq2 = Solve[{fp == fs + fi, fi == c/(nei*lendafi), (lendafi*10^3)*(-3/100) + 2185/1000 == nei}, {nei, lendafi, fs}]; eq3 = Solve[{fp == fs + fi, fs == c/(nos*lendafs), (lendafs*10^3)*(-4/100) + 2735/1000 == nos}, {fs, nos, lendafs}] ; eq4 = Solve[{0 == (2 Pi/c)*(2739/1000* fp - (25/ 10000 (-303585*10^12 + 547 fi + Sqrt[92162786625*10^24 - 33212007*10^13 fi + 299209 fi^2]))/(-555*10^12 + fi)*(555*10^12 - fi) - (25 / 10000 (437 fi + Sqrt[fi] Sqrt[-144*10^13 + 190969 fi]))/ fi*fi) - (2 Pi/clenda)}, fi]; eq1 = eq1 // Simplify; eq2 = eq2 // Simplify; eq3 = eq3 // Simplify; TimeUsed[] eq4 = eq4 // Simplify; TimeUsed[] fit[x_] := fi /. eq4 /. {clenda -> x}; fit[1.0] Plot[Re[fit[x]][[1]], {x, 1, 20}] Plot[Re[fit[x]][[2]], {x, 1, 20}] Plot[Re[fit[x]][[3]], {x, 1, 20}] －-----------------------------------------－－－－ |
2楼2012-10-27 01:33:20
walk1997
超级版主
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
3楼2012-10-27 01:34:37
4楼2012-11-17 11:57:18
