| 查看: 1847 | 回复: 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]}] |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求mathematica编程
已经有6人回复
- 【求助】开平方问题【已完结】
已经有3人回复
- 【求助】求解氢原子电子密度的时候遇到的4个问题
已经有4人回复
- 【求助】如何利用matlab计算这个无穷级数的和
已经有12人回复
- 【求助】100金币求助高人一个程序解方程
已经有15人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51005.2
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
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
