[爆卦] Step by Step 求解校正回歸
問題:
校正回歸圖
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 (韓國)
※ 文章網址: https://www.ptt.cc/bbs/Gossiping/M.1621687024.A.C2F.html
噓
05/22 20:37,
4年前
, 1F
05/22 20:37, 1F
→
05/22 20:37,
4年前
, 2F
05/22 20:37, 2F
→
05/22 20:37,
4年前
, 3F
05/22 20:37, 3F

推
05/22 20:38,
4年前
, 4F
05/22 20:38, 4F
→
05/22 20:38,
4年前
, 5F
05/22 20:38, 5F
推
05/22 20:38,
4年前
, 6F
05/22 20:38, 6F
推
05/22 20:38,
4年前
, 7F
05/22 20:38, 7F
推
05/22 20:39,
4年前
, 8F
05/22 20:39, 8F
推
05/22 20:39,
4年前
, 9F
05/22 20:39, 9F
→
05/22 20:40,
4年前
, 10F
05/22 20:40, 10F
推
05/22 20:40,
4年前
, 11F
05/22 20:40, 11F
噓
05/22 20:40,
4年前
, 12F
05/22 20:40, 12F
→
05/22 20:40,
4年前
, 13F
05/22 20:40, 13F
噓
05/22 20:40,
4年前
, 14F
05/22 20:40, 14F
推
05/22 20:40,
4年前
, 15F
05/22 20:40, 15F
推
05/22 20:40,
4年前
, 16F
05/22 20:40, 16F
推
05/22 20:40,
4年前
, 17F
05/22 20:40, 17F
→
05/22 20:40,
4年前
, 18F
05/22 20:40, 18F
推
05/22 20:40,
4年前
, 19F
05/22 20:40, 19F
推
05/22 20:41,
4年前
, 20F
05/22 20:41, 20F
推
05/22 20:41,
4年前
, 21F
05/22 20:41, 21F
噓
05/22 20:41,
4年前
, 22F
05/22 20:41, 22F
推
05/22 20:42,
4年前
, 23F
05/22 20:42, 23F
噓
05/22 20:47,
4年前
, 24F
05/22 20:47, 24F
推
05/22 20:48,
4年前
, 25F
05/22 20:48, 25F
推
05/22 20:50,
4年前
, 26F
05/22 20:50, 26F
推
05/22 20:51,
4年前
, 27F
05/22 20:51, 27F
推
05/22 20:54,
4年前
, 28F
05/22 20:54, 28F
推
05/22 20:55,
4年前
, 29F
05/22 20:55, 29F
→
05/22 20:57,
4年前
, 30F
05/22 20:57, 30F

推
05/22 20:58,
4年前
, 31F
05/22 20:58, 31F
推
05/22 21:00,
4年前
, 32F
05/22 21:00, 32F
推
05/22 21:05,
4年前
, 33F
05/22 21:05, 33F
推
05/22 21:06,
4年前
, 34F
05/22 21:06, 34F
→
05/22 21:06,
4年前
, 35F
05/22 21:06, 35F
噓
05/22 21:07,
4年前
, 36F
05/22 21:07, 36F
推
05/22 21:13,
4年前
, 37F
05/22 21:13, 37F
噓
05/22 21:14,
4年前
, 38F
05/22 21:14, 38F
→
05/22 21:16,
4年前
, 39F
05/22 21:16, 39F
推
05/22 21:17,
4年前
, 40F
05/22 21:17, 40F
→
05/22 21:18,
4年前
, 41F
05/22 21:18, 41F
推
05/22 21:19,
4年前
, 42F
05/22 21:19, 42F
推
05/22 21:38,
4年前
, 43F
05/22 21:38, 43F
→
05/22 21:54,
4年前
, 44F
05/22 21:54, 44F
推
05/22 21:59,
4年前
, 45F
05/22 21:59, 45F
推
05/22 22:32,
4年前
, 46F
05/22 22:32, 46F
推
05/22 22:35,
4年前
, 47F
05/22 22:35, 47F
推
05/22 23:17,
4年前
, 48F
05/22 23:17, 48F
推
05/22 23:24,
4年前
, 49F
05/22 23:24, 49F
推
05/24 15:20,
4年前
, 50F
05/24 15:20, 50F