| 查看: 1310 | 回复: 1 | |||
[交流]
用mathematica求解微分方程
|
|
如何用mathematica解这两个方程，并画出xi与t的曲线图，（其中yi=i*h,m = 2.4/10^16; \[Beta] = 5/10^5; G = 82*10^9; v = 0.29; b = 2.5/10^10; D= G/2 \[Pi]/(1 - v); h = 10 b; w = 0.1; n = 3; d = (n + 1)*h; ff= 0.5*10^9;(ff替代了方程里的一个字符) 求高手帮忙  我的解法是：方程里的f1与f2没有赋值成功，求高手帮忙 m = 2.4/10^16; \[Beta] = 5/10^5; G = 82*10^9; v = 0.29; b = 2.5/10^10; dd = G/2 \[Pi]/(1 - v); h = 10 b; w = 0.1; n = 3; d = (n + 1)*h; ff = 0.5*10^9; ClearAll[y, x, f2, f1]; y = Table[i*h, {i, n}]; f2 = Table[1, {n}]; f1 = Table[1, {n}]; For[i = 1, i <= n, i++, For[k = 1, k <= n, k++, If[k == i, f2[[k]] = 0, f2[[k]] = (x[t] - x[k][t]) ((x[t] - x[k][t])^2 - (y[] - y[[k]])^2)/((x[t] - x[k][t])^2 + (y[] - y[[k]])^2)^2]]; f1[] = Apply[Plus, f2] ]; eqns = {Table[{m*x''[t] + \[Beta]*x'[t] - b*ff - dd*b^2*f1[] + dd*w*b*(x[t] y[]/(x[t]^2 + y[]^2) - x[t]*(y[] - d)/(x[t]^2 + (y[] - d)^2)) == 0}, {i, n}], Table[{x[0] == 0, x'[0] == 0}, {i, n}]}; NDSolve[eqns, Table[x, {i, n}], {t, 0, 0.1/10^9}]  360反馈意见截图1167302256711587.png  360反馈意见截图1699120795121125.png |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
范进中举一文的中心思想
已经有4人回复
-
跳槽后在研项目怎么办？
已经有7人回复
2楼2015-08-16 09:07:45
40