24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xhw2211

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 9240.1
红花: 1
帖子: 301
在线: 56.3小时
虫号: 570836
注册: 2008-06-08
专业: 机械结构强度学

[求助] 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人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2011-07-25 13:38:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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]}