看板 Gossiping作者 angelcherry (cherry)標題 [爆卦] Step by Step 求解校正回歸時間 Sat May 22 20:37:02 2021
問題:
校正回歸圖
https://i.imgur.com/hzFZ8VP.png
目標:
我們想要求解這個黑盒子(Black box),這個Function是什麼?
https://i.imgur.com/AcL7p7P.png
準備軟體:
Excel,Matlab
首先,必須要先有個假設,這個黑盒子會隨時間改動,符合確診數校正不一致。
那麼每天校正一次,即一天等於86400秒間隔。
第一步:
準備資料集,將校正回歸之後跟之前的數據分別設定為輸出和輸入。
開啟一個新的Excel,第一列輸入input和output,輸入完畢後存成blackbox.xlsx。
https://i.imgur.com/cyD92bF.png
第二步:
執行Matlab,選擇Import Data,選擇剛剛的blackbox.xlsx。
https://i.imgur.com/rmZicxh.png
第三步:
資料集會顯示出來,預設會抓取所有Data,選擇Import Selection。
https://i.imgur.com/2qwwjcI.png
第四步:
產生一個6x2的table變數,這樣可以開始建構所需的資料物件(Data object)。
https://i.imgur.com/46rfJzU.png
第五步:
建構一個iddata object,採樣時間為86400秒,輸入為blackbox.input這一行,輸出同理
。
會產生y1和u1代表輸出和輸入。
https://i.imgur.com/GoJktCw.png
第六步:
開啟systemIdentification工具來嘗試求解這個黑盒子模型。
https://i.imgur.com/PbxOmZf.png
第七步:
下拉Import data,選擇Data objcet...。
https://i.imgur.com/2uC49cF.png
第八步:
Object文字框中填入data,即剛剛創建的資料物件,按下Import按鈕。
https://i.imgur.com/a6bi53R.png
圖上顯示資料集的波形圖。
https://i.imgur.com/R3zXs5o.png
第九步:
勾選Time plot可以查看輸入輸出資料波形。
https://i.imgur.com/G76UrIV.png
第十步:
下拉Estimate-->,選擇Transfer Function Models...(轉移函數)。
https://i.imgur.com/mHZhMSq.png
第十一步:
更改0個zeros(零點),選擇Discrete-time(離散時間),按下Estimate。
https://i.imgur.com/9orFQuY.png
系統就開始估計。
https://i.imgur.com/ZWKCJGQ.png
之後會產生估計後的模型。
https://i.imgur.com/HRC2mNy.png
第十二步:
點選估計後的模型,會顯示近似解和Fit的程度85.11%。
1 - 0.5842 z^-1 - 0.4157 z^-2
https://i.imgur.com/PW0xldN.png
第十三步:
勾選Model output,會顯示和目標波形Fit程度。
https://i.imgur.com/4mJtOIw.png
可嘗試不同的估計模型,但目前這個是最佳解。
https://i.imgur.com/vYr0CS6.png
下台一鞠躬,康撒哈米達:)
--
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 222.100.19.84 (韓國)
※ 文章代碼(AID): #1WgFhmml (Gossiping)
※ 文章網址: https://www.ptt.cc/bbs/Gossiping/M.1621687024.A.C2F.html
→ kumori: 跨謀拉 講重點好嗎5F 05/22 20:38
→ ridecule: 文組應該是推吧 因為看不懂 好神奇10F 05/22 20:40
噓 weekman: 用拉式轉換前提要是線性系統12F 05/22 20:40
→ Ceferino: 直接說你覺得今天實際數多少13F 05/22 20:40
推 Cwren: 笑死17F 05/22 20:40
→ rdmtb: 太認真了。18F 05/22 20:40
推 billliu: ID 好像女生要不要校正回歸一下~~19F 05/22 20:40
推 aivxg: 跟我想的一樣 笑死21F 05/22 20:41
推 yling5518: 假設疫情減緩,找模型滿足不是難事。25F 05/22 20:48
推 yl3218: 下星期再來驗證這個公式是不是對的26F 05/22 20:50
推 azeroth: XDDDDDDDDD八卦閒人真的很多27F 05/22 20:51
推 omoideomoi: 可以不要到這邊浪費才能嗎 認吱作戰單位不會信的啦28F 05/22 20:54
推 clsied: 我前天剛考完你不要這樣34F 05/22 21:06
→ hw1: 幹...看不董辣35F 05/22 21:06
推 Tianmu: 假設錯誤。 因為有外力改變絕對沒辦法用公式推算。 中間已經有增加防疫為三級。40F 05/22 21:17
推 haniah: 痾 看不懂45F 05/22 21:59
推 hodj: 浪費才能系列開始49F 05/22 23:24
--