请大家看看我的三维建模代码哪里出了问题
mathematica吧
全部回复
仅看楼主
level 4
这是我的文档,我是真的不知道为什么,我咋调数值这个图像就是只在y=-1000那里凹下去,明明很多地方都应该下沉,这个图像就只凹那一块,这个函数明明在x=0,y=500处Wxy值为-2000,可是它在图像上该处Wxy值为0,我真的想不明白
2023年12月08日 16点12分 1
level 4
r[i_] := h + (2 i - 1)/(2 n)*L2 Sin[\[Alpha]] Cot[\[Theta]];
\[Alpha] =
10 Degree; \[CurlyPhi]1 = 0; \[CurlyPhi]2 = 0; L2 = 65; L1 = 65; \
Wmax = -1.6; h = 900; n = 10; \[Theta] = 85 Degree;
W[x_, y_] :=
Sum[Wmax/(2 (1 +
Tan[\[Alpha]] Cot[\[Theta]]) (Erf[
Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - ((i - 1)/
n) L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))] -
Erf[Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - (i/n)*
L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))]) ) ((Erf[
Sqrt[Pi] (x + (i - 1)/n*L2 Sin[\[CurlyPhi]1])/r[i]] -
Erf[Sqrt[
Pi] (x + (i - 1)/n*
L2 Sin[\[CurlyPhi]1] - (L1 +
L2 (2 i - 1)/(2 n)) (Sin[\[CurlyPhi]1] +
Cos[\[CurlyPhi]1] + Tan[\[CurlyPhi]2]))/r[i]]))/2, {i,
1, 100}]; Plot3D[W[x, y], {x, -1000, 2000}, {y, -1000, 1000},
PlotRange -> All, AxesLabel -> {"x", "y", "W[x, y]"},
PlotLabel -> "W[x, y]的三维图"]
Print["在 x = 0 且 y = 500 处的具体值为: ", result]
2023年12月08日 16点12分 2
level 4
求助,为了建这个模型0我MMA学了快一个月了,失败了我也想知道我是怎么死的,函数式我检查过了没有问题
2023年12月08日 16点12分 3
信誓旦旦说式子没问题到头来还是式子错了的情况我见得多了。Plot3D绘图没有多少可错的点,十有八九是表达式本身不对。你按下Ctrl+Shift+N切成二维表达式再好好检查下。
2024年01月06日 05点01分
level 4
按道理应该呈现这种下沉三维图形的
我真的要长脑子了
2023年12月08日 16点12分 4
level 8
应该是Z轴单位太大了,你看你的W[x,y] 轴的单位是10^7,-2000多的值在上面也就10^-4级别,有但和零差别不大。重新调整Z轴范围效果如下。(作了一些改动,为了加快执行速度)
W[x_, y_] =
Total@N@Table[
Wmax/(2 (1 +
Tan[\[Alpha]] Cot[\[Theta]]) (Erf[
Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - ((i - 1)/
n) L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))] -
Erf[Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - (i/n)*
L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))])) ((Erf[
Sqrt[Pi] (x + (i - 1)/n*L2 Sin[\[CurlyPhi]1])/r[i]] -
Erf[Sqrt[
Pi] (x + (i - 1)/n*
L2 Sin[\[CurlyPhi]1] - (L1 +
L2 (2 i - 1)/(2 n)) (Sin[\[CurlyPhi]1] +
Cos[\[CurlyPhi]1] + Tan[\[CurlyPhi]2]))/r[i]]))/
2, {i, 1, 100}];
Plot3D[W[x, y], {x, -1000, 2000}, {y, -1000, 1000}, PlotPoints -> 10,
PlotRange -> {-50000, 0}, AxesLabel -> {"x", "y", "W[x, y]"},
PlotLabel -> "W[x, y]的三维图"]
2023年12月09日 03点12分 5
level 9
绝对是函数式或者参数错了
2024年01月06日 09点01分 6
1