| 查看: 1266 | 回复: 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 |
回复此楼» 猜你喜欢
化工学硕294分，求导师收留
已经有29人回复
-
求调剂
已经有10人回复
-
考研求调剂
已经有13人回复
-
求调剂
已经有3人回复
-
人工智能320调剂08工类还有机会吗
已经有17人回复
-
考研英一数一338分
已经有10人回复
-
求助调剂，跨调
已经有15人回复
-
085600材料与化工329分求调剂
已经有20人回复
-
085600材料与化工349分求调剂
已经有15人回复
-
求调剂
已经有13人回复
2楼2015-08-16 09:07:45
