24小时热门版块排行榜    

查看: 2112  |  回复: 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

金虫 (正式写手)

把修改后的程序贴出来:
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
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 00:28 by winsaint
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见