博士畢業 Day7

老實說從第一天之後就沒有進度阿阿阿~
會空七天是因為我耍廢兩天在生理期上，然後前往新竹跟家人露營
然後回台北各種跳舞痾阿阿阿~
但是逝者已矣，來者猶可追!!!

所以今天來個進度~

------------------------------
總而言之，Lucky me ，我找到另一個matlab code : fbm1d (網頁請點)
一個one dimensional 的 fractional Brownian motion
然後發現依然無論我模擬多少次(次數上升到1000000都沒收斂) ，

最後我找到他的numerical code來源~(LUCKY ME)

% In Stochastic Geometry, Spatial Statistics and Random Fields(pp. 369-404)

% Springer International Publishing, DOI: 10.1007/978-3-319-10064-7_12

想到了不收斂得原因(完全沒從比較有系統的數值方法下去搜尋阿我! )
時間不夠長?!
這個意思比較像是步數不夠多! 而我回去看他source code的解釋也是一個明白了!
他用一個正規化normalize的方法對時間進行處理! 所以t0 = 0 < 1/n < 2/n <...< tn=1 (n : grid points)
當n不夠大，就是切的不夠細，感覺就像你模擬一個很長時間的process但是中間只有三個點，導致很難收斂，太粗! 目前是讓 n > = 10^4 ，然後讓realizations > 100,000次
下面是他怎麼做(p413~414; Schmidt, V., 2014)

Schmidt, V. (2014). Stochastic geometry, spatial statistics and random fields. Springer.

#話說因為找到source code，所以我就沒有用matlab那個function，自己貼了做成一個function，叫做FBMfromBook，其實只是因為fbm1d的n不知道為什麼規定一定要是2的次方倍數~?!

#話說matlab用矩陣一個限制很多阿阿阿阿!!! 是否該學個平行運算~

🔻🔻🔻
所以目前我本人的想法是這樣，讓n很大! 意思是切很細，每一個事件的時間都標準化到0~1再乘開!!!最長的事件，乘開之後的1/n*T<=0.05就好! 這樣比較短的事件，只要照著dt 跳著取值就好!
🔺🔺🔺

-----------------
題外話，我今天一個看了這首歌的MV然後眼淚就掉下來~
閃靈、元千歲-暮沉武德殿(民謠版) CHTHONIC-Defenders of Bú-Tik Palace(acoustic)

今天正好是言論自由日呢!
我想，如果真的被大陸控制，我會毫不猶豫去赴死~看是要怎麼樣能夠造成影響，即使造成一點點漣漪也可以! 雖然對親友很抱歉，但是沒有一條生命會白白送死，即使是我們一輩子都不會知道名子的人都一樣!更何況是在這個資訊爆炸的時代，誰都能肉搜別人的時代，不會有一條生命不留下一些地走，但是必須要懷抱志氣的離開，必須要，而且不會後悔!
講這麼壯烈的話是沒有用的，所以現在必須要造成一點影響才行，所以要多傳一點訊息給家裡的LINE群組!!