| 查看: 2124 | 回复: 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]}] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求mathematica编程
已经有6人回复
- 【求助】开平方问题【已完结】
已经有3人回复
- 【求助】求解氢原子电子密度的时候遇到的4个问题
已经有4人回复
- 【求助】如何利用matlab计算这个无穷级数的和
已经有12人回复
- 【求助】100金币求助高人一个程序解方程
已经有15人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
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
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
20
