作者kittor ()
看板Python
标题[问题] 求指正
时间Tue Apr 17 09:42:26 2018
import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
def model(y, u ,t):
2dydt=-y+u
return dydt
y0=1
t=np.linspace(0,10)
u=0
y1=odeint(model, y0, t, args=(u, ))
u=0.2
y2=odeint(model, y0, t, args=(u, ))
.
.
.
u=2.0
y11=odeint(model, y0, t, args=(u, ))
plt.plot(t, y1, 'g:', linewidth=5, label='u=0')
plt.plot(t, y2, 'g:', linewidth=5, label='u=0.2')
.
.
.
plt.plot(t, y11, 'g:', linewidth=5, label='u=2.0')
plt. xlabel('time')
plt. ylabel('y(t)')
plt.show()
不好意思，自学程式中的我，只能贴上来求指正
如果有错的地方，希望可以指点我一下
...的部分代表条件递移的相同程式码，直接省略掉了
--
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 118.168.165.35
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Python/M.1523929349.A.65C.html
1F：→ HybridSC: 怎麽你一开始的function看起来就不太对 04/17 16:42
2F：→ HybridSC: 2dydt ? 04/17 16:42
3F：→ kittor: 没关系了，我後来有找到解答了 04/17 16:44
4F：嘘 Blankfein: 邱个屁 04/18 06:25
是要邱什麽? B大看不见我写的"不好意思"和"希望"两个用词是吗?
5F：推 OnePiecePR: 加油 04/18 08:35
※ 编辑: kittor (118.166.153.138), 04/18/2018 09:11:07
6F：嘘 leo00134: 自己语气看不出有问题 04/18 15:28
所以是我的语气，你就看得出有问题?
※ 编辑: kittor (111.250.165.102), 04/18/2018 15:36:58
7F：→ ooooSEX: 祝你好运 04/18 16:40
8F：嘘 hadoop: 颇ㄏ 04/18 19:01
所以这位又有什麽高见了?
※ 编辑: kittor (111.248.168.148), 04/18/2018 19:49:26
9F：嘘 ooooSEX: 看来你很喜欢辩论 04/18 22:07
11F：→ kittor: 然後呢 04/18 22:09
12F：→ kittor: 你信不信，我曾有不是你认为爱辩论的我的时期 04/18 22:11
13F：嘘 laputaflutin: 04/18 23:27
14F：→ kittor: 所以又想说什麽? 04/18 23:33
15F：嘘 AlaRduTP: 什麽一二楼没给指正 那时我是在回N大的文 赞同他的内容 04/18 23:36
16F：→ AlaRduTP: 又不是在回答你的问题 现在才是特地回来嘘你 04/18 23:36
17F：→ kittor: 推 AlaRduTP: 真的呢~ 发问写的好也能提高大家回答意愿！ 04/18 23:38
18F：→ kittor: 所以到底这篇是n大的文，还是我的文? 04/18 23:38
19F：嘘 pups003: 奇文共赏 04/19 00:20
20F：嘘 ooooSEX: 亏你还是台大化工所毕业的 04/19 01:39
21F：→ Blankfein: 难怪30几岁还在搭讪7-11店员 04/19 06:34
22F：嘘 ACMANIAC: 讲一下学习态度差就气成这样，EQ要再加点油 04/19 08:16
23F：→ kittor: 我的学习态度是呈现在有人说了的是知识内容相关的部分 04/19 09:24
24F：→ kittor: 给b大，你确定是我在搭讪店员吗? 04/19 09:25
25F：→ coeric: 难道是搭讪店长？XDDDDDD 04/19 10:26
26F：→ kittor: 能够被看得这麽平凡，我应该要高兴吼 04/19 10:30
27F：嘘 Fanchien: 嘘 04/19 21:36
28F：嘘 leo00134: 真是奇文共赏 像如此不知觉的人还是第一次看到 04/19 21:37
29F：→ kittor: 是怎样地不自觉了，我都能想到没有给出资讯而能帮助的提问 04/19 21:40
30F：→ kittor: 了，结果一堆人只纠结在"资讯不足" 04/19 21:40
31F：嘘 leo00134: 我只纠结在你回文的语气 请你自己再看看下面那篇回文是 04/19 21:42
32F：→ leo00134: 否有问题 04/19 21:42
33F：→ kittor: 你又知道我真的不高兴的时候的用字和下面那篇回文的用字了 04/19 21:43
34F：→ kittor: ? 04/19 21:43
35F：→ kittor: 给你点出我哪几句你有意见的，来来来 04/19 21:44
36F：嘘 ringlwe: 态度不好 嘘 04/20 08:13
37F：嘘 hadoop: 2 04/20 12:58
38F：嘘 pxycho: 这样子问真的呵呵 04/20 18:44
39F：嘘 cpper: 嘘 04/21 09:23
40F：→ kittor: 这麽爱嘘，也没看见楼上几位指教别人，整个电脑操作是只会 04/21 09:58
41F：→ kittor: 按2推文而已? 04/21 09:58
42F：嘘 cpper: 态度这麽差 嘘 04/21 11:58
43F：→ kittor: 这样就差，是没看过不删码的问题是吗 04/21 14:34
44F：嘘 jn8029: 路过 04/21 19:38
45F：→ kittor: 这麽懂得给嘘，我看给的推也不会怎样的价值而已 04/21 20:28
46F：→ cpper: 呵呵我刚才在java板看到你问的问题，也是令人喷血 04/21 20:30
47F：→ cpper: 我给你一大段推文意见了 够意思了吧 多学学怎麽问问题吧 04/21 20:30
48F：→ cpper: 别再浪费大家时间了 04/21 20:30
49F：→ cpper: 另外学程式如果你自己连IDE环境都不装，凭空想像程式会怎麽 04/21 20:31
50F：→ cpper: 跑，那是永远不可能学好的 04/21 20:31
51F：→ kittor: 我敢手写就好了，反正程式语言不见得是我的主力，爽 04/21 20:37
52F：→ kittor: 间隔两个字看不见? 04/21 20:43
53F：→ kittor: 说真的，这个微方背後的数学模式也不是算难的，能被你这样 04/21 20:50
54F：→ kittor: 评，我还庆幸自己会工程数学勒 04/21 20:51
55F：→ kittor: 另外再提醒你一点，原来你也是要跑过才知道错在哪的? 04/21 20:52
56F：→ kittor: 也建议你，有些人的说法根本已经在打自己的脸了，不盖你 04/21 20:52
57F：嘘 Blankfein: 呵 04/21 21:02
58F：→ kittor: 继续啊，我每天都来动脑看你们嘘文的价值是什麽 04/21 21:02
59F：→ kittor: 只要还有一丝别人语法好不好不关己的事，就不知道为什麽2 04/21 21:03
60F：→ kittor: 按得心安理得了 04/21 21:03
61F：嘘 EGsux: 难怪三十几岁没楼没车没女友 04/22 01:02
62F：→ kittor: → kittor: 能够被看得这麽平凡，我应该要高兴吼 04/22 08:45
63F：→ kittor: 送给还在看这篇的人，量你们没人有胆敢让可能是人生胜利组 04/22 08:46
64F：→ kittor: 的自己成为人生失败组的 04/22 08:46
65F：嘘 leo00134: 人生到底多边缘才会这样刷存在感? 04/22 09:36
66F：→ kittor: 所以到底要多好骗才会被文字迷惑? 04/22 09:47
67F：嘘 Blankfein: 路过 04/22 19:19
69F：→ ACMANIAC: 有点像庄园大师，哈哈哈哈 04/22 20:25
70F：嘘 paul0407: 您很强~您很棒~任何指正都是多余! 04/23 11:53