24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1713  |  回复: 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]}]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhw2211

金虫 (正式写手)

我按你的方法改了程序,但是还是没有结果,感觉是不是拉普拉斯变换后解出的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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

jjdg

版主 (知名作家)

优秀版主

你是用z,t作图的吧?
但你的z值是从哪里来?
努力学习!以正当途径!获得需要的知识!
2楼2011-07-25 17:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhw2211

金虫 (正式写手)

最后一句应该是lot[m, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}]
3楼2011-07-25 18:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhw2211

金虫 (正式写手)

晕,居然出现一个笑脸。最后一句话z应该是m,我写错了。
4楼2011-07-25 18:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见