作者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/m.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