作业帮 > 综合 > 作业

Mathematica求数值解,错在哪?

来源：学生作业帮 编辑：神马作文网作业帮 分类：综合作业 时间：2024/11/13 22:33:20
Mathematica求数值解,错在哪?
h = 4.13566743*(10^-15)
m = 0.1*9.109*(10^-31）
Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*
Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*
Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] +
Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*
Cos[Sqrt[8*Pi^2*m*y/h^2]*x/2] - Cos[x] == 0,{x,y}],{x,0,10,
1}]
Mathematica求数值解,错在哪?
1 你的m的赋值式里面用了中文标点.
2 Do是不会主动显示循环中的计算结果的,要使用Print语句,当然最方便的改法还是把Do改成Table.
3 最关键的一点,这个超越方程恐怕不是Solve和NSolve能解的了的.可以考虑用FindRoot,不过这就需要你对解的范围有大概的了解,要了解解的范围,一个可能的方法,是使用ContourPlot观察解的位置：
h = 4.13566743*(10^-15);
m = 0.1*9.109*(10^-31);
ContourPlot[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2*m*y/h^2]*x/2] - Cos[x] == 0, {x, 0, 11}, {y, 0, 20}, PlotPoints -> 30]


然后就根据图上显示的解的范围去寻根吧.
再问： 那如果要把==0换成
Mathematica求数值解,错在哪? mathematica里面求方程的数值解 问个Mathematica求方程数值解的问题 求数值微分的mathematica方法 mathematica如何对分段函数的积分求数值积分 使用mathematica解常微分方程（数值解） mathematica求解微分方程的数值解,做出图像, mathematica与matlab计算比较?各自优点在哪里?仅纯数值计算而言,谁更适合? 在mathematica中求指数为矩阵得计算 【Mathematica】这个方程哪错了? mathematica问题急求! 用mathematica求参数.