大概看了一下你的原始碼,你的程式本來就是用來時變系統的,似乎是一維連續釋放模式,因為C(1)程式執行執行過程沒有變動過(是CC(1),筆誤sorry)不過連續釋放模式CC(1)應該要是濃度最高的,你跑出你跑出來的模型很奇怪,可能你有變動程式碼過的關係
https://i.imgur.com/F2EvkVn.jpg 連續釋放模型的污染峰應該會像圖1一樣隨時間變化,中間確實會降不來你模型尺度和釋放時間沒配合好的話,OUTPUT數據繪製起來會感覺怪怪的,如圖2,你如果用較小的模型去看已經發展到T3的污染峰的話,看起來像是一條直線,但如果只看發展到T2的污染峰,看起來就很親切了。同理,如果你用如圖3的較大尺度觀察T3污染峰的話,圖形似乎也正常多了。建議以未修改的原始程式,先隨意調整END_X、TT,甚至是速度項的UU(1)、UU(2),來試試看是否可產生類似的圖形。但調END_X時要注意3個Xstop的配合。另外你把CC(N)設為0的話,那點會類似1個很強的sink,所以你跑的圖形可看到CC(N)前幾點濃度被硬拉下來又因為是1維模型,質量繞不過去那點,他表現得會像是一個污染物抽出量一直在變化的井,目的是讓那點的污染物質量(濃度)抽為0。