| 查看: 2104 | 回复: 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]}] |
回复此楼» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有8人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有9人回复
-
放榜前的不淡定
已经有18人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求mathematica编程
已经有6人回复
- 【求助】开平方问题【已完结】
已经有3人回复
- 【求助】求解氢原子电子密度的时候遇到的4个问题
已经有4人回复
- 【求助】如何利用matlab计算这个无穷级数的和
已经有12人回复
- 【求助】100金币求助高人一个程序解方程
已经有15人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
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
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
5楼2011-07-25 21:47:56
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
|
我按你的方法改了程序，但是还是没有结果，感觉是不是拉普拉斯变换后解出的Z表达式太复杂，mathematica无法进行拉普拉斯反变换？？？下面是程序运行后得到的m m={1/2 Sqrt[\[Pi]/5] InverseLaplaceTransform[( E^(5 + s + s^2/20) (1 - 840/((7 + s) (120 + s))) Erfc[(10 + s)/( 2 Sqrt[5])])/( 1 - 840/((7 + s) (120 + s)) - ( 110 (1 - 1/2 E^(1/20 (10 + s)^2) Sqrt[\[Pi]/5] s Erfc[(10 + s)/(2 Sqrt[5])]))/(110 + s)), s, t]} |
6楼2011-07-26 09:25:25
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
7楼2011-07-26 09:26:19
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
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
★
jjdg(金币+1): 感谢参与 2011-07-30 14:50:08
jjdg(金币+1): 感谢参与 2011-07-30 14:50: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[m, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] |
9楼2011-07-26 09:33:09
xhw2211
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 9240.1
- 红花: 1
- 帖子: 301
- 在线: 56.3小时
- 虫号: 570836
- 注册: 2008-06-08
- 专业: 机械结构强度学
10楼2011-07-26 09:46:47
