| 查看: 2118 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xhw2211金虫 (正式写手)
|
[求助]
mathematica 编程求助
|
||
|
头一次用mathtematica编的程序，总是没有画出图形，不知道错在什么地方了，大牛看一下。 z1 = LaplaceTransform[Exp[-10*t - 5*t^2], t, s]; z2 = LaplaceTransform[(10 + 10*t)*Exp[-10*t - 5*t^2], t, s]; z3 = Simplify[z2]; equ1 = z1 + z3*y - z; y1 = LaplaceTransform[7*Exp[-7*t], t, s]; y2 = LaplaceTransform[110*Exp[-110*t], t, s]; equ2 = y1*x + y2*z - y; x1 = LaplaceTransform[120*Exp[-120*t], t, s]; equ3 = x1*y - x; Solve[{equ1 == 0, equ2 == 0, equ3 == 0}, {x, y, z}] m = InverseLaplaceTransform[z, s, t] Plot[z, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求mathematica编程
已经有6人回复
- 【求助】开平方问题【已完结】
已经有3人回复
- 【求助】求解氢原子电子密度的时候遇到的4个问题
已经有4人回复
- 【求助】如何利用matlab计算这个无穷级数的和
已经有12人回复
- 【求助】100金币求助高人一个程序解方程
已经有15人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
|
把修改后的程序贴出来： z1 = LaplaceTransform[Exp[-10*t - 5*t^2], t, s]; z2 = LaplaceTransform[(10 + 10*t)*Exp[-10*t - 5*t^2], t, s]; z3 = Simplify[z2]; equ1 = z1 + z3*y - z; y1 = LaplaceTransform[7*Exp[-7*t], t, s]; y2 = LaplaceTransform[110*Exp[-110*t], t, s]; equ2 = y1*x + y2*z - y; x1 = LaplaceTransform[120*Exp[-120*t], t, s]; equ3 = x1*y - x; rr = Solve[{equ1 == 0, equ2 == 0, equ3 == 0}, {x, y, z}]; z = z /. rr; m = InverseLaplaceTransform[z, s, t] Plot[z, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] |
8楼2011-07-26 09:27:33
2楼2011-07-25 17:51:26
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
lot[m, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] 3楼2011-07-25 18:33:04
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
4楼2011-07-25 18:35:01
30