PTT推薦

[爆卦] Step by Step 求解校正回歸

angelcherry 發表於 2021/5/22 下午8:37:02
看板Gossiping標題[爆卦] Step by Step 求解校正回歸作者
angelcherry
 (cherry)時間推噓25 推:32 噓:7 →:11

問題:
校正回歸圖
https://i.imgur.com/hzFZ8VP.png

圖https://i.imgur.com/hzFZ8VP.jpg, [爆卦] Step by Step 求解校正回歸

目標：
我們想要求解這個黑盒子(Black box)，這個Function是什麼？
https://i.imgur.com/AcL7p7P.png
圖https://i.imgur.com/AcL7p7P.jpg, [爆卦] Step by Step 求解校正回歸

準備軟體：
Excel，Matlab

首先，必須要先有個假設，這個黑盒子會隨時間改動，符合確診數校正不一致。
那麼每天校正一次，即一天等於86400秒間隔。

第一步：
準備資料集，將校正回歸之後跟之前的數據分別設定為輸出和輸入。
開啟一個新的Excel，第一列輸入input和output，輸入完畢後存成blackbox.xlsx。
https://i.imgur.com/cyD92bF.png

圖https://i.imgur.com/cyD92bF.jpg, [爆卦] Step by Step 求解校正回歸

第二步：
執行Matlab，選擇Import Data，選擇剛剛的blackbox.xlsx。
https://i.imgur.com/rmZicxh.png
圖https://i.imgur.com/rmZicxh.jpg, [爆卦] Step by Step 求解校正回歸

第三步：
資料集會顯示出來，預設會抓取所有Data，選擇Import Selection。
https://i.imgur.com/2qwwjcI.png
圖https://i.imgur.com/2qwwjcI.jpg, [爆卦] Step by Step 求解校正回歸

第四步：
產生一個6x2的table變數，這樣可以開始建構所需的資料物件(Data object)。
https://i.imgur.com/46rfJzU.png
圖https://i.imgur.com/46rfJzU.jpg, [爆卦] Step by Step 求解校正回歸

第五步:
建構一個iddata object，採樣時間為86400秒，輸入為blackbox.input這一行，輸出同理。
會產生y1和u1代表輸出和輸入。
https://i.imgur.com/GoJktCw.png
圖https://i.imgur.com/GoJktCw.jpg, [爆卦] Step by Step 求解校正回歸

第六步：
開啟systemIdentification工具來嘗試求解這個黑盒子模型。
https://i.imgur.com/PbxOmZf.png
圖https://i.imgur.com/PbxOmZf.jpg, [爆卦] Step by Step 求解校正回歸

第七步：
下拉Import data，選擇Data objcet...。
https://i.imgur.com/2uC49cF.png
圖https://i.imgur.com/2uC49cF.jpg, [爆卦] Step by Step 求解校正回歸

第八步：
Object文字框中填入data，即剛剛創建的資料物件，按下Import按鈕。
https://i.imgur.com/a6bi53R.png
圖https://i.imgur.com/a6bi53R.jpg, [爆卦] Step by Step 求解校正回歸
圖上顯示資料集的波形圖。
https://i.imgur.com/R3zXs5o.png
圖https://i.imgur.com/R3zXs5o.jpg, [爆卦] Step by Step 求解校正回歸

第九步:
勾選Time plot可以查看輸入輸出資料波形。
https://i.imgur.com/G76UrIV.png
圖https://i.imgur.com/G76UrIV.jpg, [爆卦] Step by Step 求解校正回歸

第十步：
下拉Estimate-->，選擇Transfer Function Models...(轉移函數)。
https://i.imgur.com/mHZhMSq.png
圖https://i.imgur.com/mHZhMSq.jpg, [爆卦] Step by Step 求解校正回歸

第十一步：
更改0個zeros(零點)，選擇Discrete-time(離散時間)，按下Estimate。
https://i.imgur.com/9orFQuY.png
圖https://i.imgur.com/9orFQuY.jpg, [爆卦] Step by Step 求解校正回歸
系統就開始估計。
https://i.imgur.com/ZWKCJGQ.png
圖https://i.imgur.com/ZWKCJGQ.jpg, [爆卦] Step by Step 求解校正回歸
之後會產生估計後的模型。
https://i.imgur.com/HRC2mNy.png
圖https://i.imgur.com/HRC2mNy.jpg, [爆卦] Step by Step 求解校正回歸

第十二步：
點選估計後的模型，會顯示近似解和Fit的程度85.11%。
1 - 0.5842 z^-1 - 0.4157 z^-2
https://i.imgur.com/PW0xldN.png
圖https://i.imgur.com/PW0xldN.jpg, [爆卦] Step by Step 求解校正回歸

第十三步：
勾選Model output，會顯示和目標波形Fit程度。
https://i.imgur.com/4mJtOIw.png
圖https://i.imgur.com/4mJtOIw.jpg, [爆卦] Step by Step 求解校正回歸
可嘗試不同的估計模型，但目前這個是最佳解。
https://i.imgur.com/vYr0CS6.png
圖https://i.imgur.com/vYr0CS6.jpg, [爆卦] Step by Step 求解校正回歸

下台一鞠躬，康撒哈米達:)

--

※ PTT留言評論
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 222.100.19.84 (韓國)
※ 文章網址: https://www.ptt.cc/Gossiping/M.1621687024.A.C2F
噓

MotleyCrue05/22 20:37文組會來噓妳

→

guitar022505/22 20:37你為什麼要讓我想起來！！

→

htc1005/22 20:37https://i.imgur.com/DvUsbWT.jpg

圖https://i.imgur.com/DvUsbWT.jpg, [爆卦] Step by Step 求解校正回歸
推

appleball20005/22 20:38推XD

→

kumori05/22 20:38跨謀拉 講重點好嗎

推

sandiato05/22 20:38這個可能要找t2d來解答

推

sam9020505/22 20:38太鬼了吧

推

hipmyhop05/22 20:39笑死

推

alcpeon91105/22 20:39笑死 不要這麼認真好嗎

→

ridecule05/22 20:40文組應該是推吧 因為看不懂 好神奇

推

lightaurora05/22 20:40你是鬼吧

噓

weekman05/22 20:40用拉式轉換前提要是線性系統

→

Ceferino05/22 20:40直接說你覺得今天實際數多少

噓

xiaoyao05/22 20:40老闆：講重點

推

Neisseria05/22 20:40問教授可不可以當碩論啊

推

plusonezero05/22 20:40嗯嗯，跟我想的一樣

推

Cwren05/22 20:40笑死

→

rdmtb05/22 20:40太認真了。

推

billliu05/22 20:40ID 好像女生要不要校正回歸一下~~

推

JackTheRippe05/22 20:41快推 不然會被說看不懂

推

aivxg05/22 20:41跟我想的一樣 笑死

噓

lucifiel161805/22 20:41沒料的就愛貼一堆圖講不出東西

推

strongrain05/22 20:42你太認真了吧XD

噓

f1769081505/22 20:47講中文

推

yling551805/22 20:48假設疫情減緩，找模型滿足不是難事。

推

yl321805/22 20:50下星期再來驗證這個公式是不是對的

推

azeroth05/22 20:51XDDDDDDDDD八卦閒人真的很多

推

omoideomoi05/22 20:54可以不要到這邊浪費才能嗎 認吱作戰單位不會信的啦

推

AlianF05/22 20:55公三小

→

NCUking05/22 20:57https://i.imgur.com/UtH2g7v.jpg

圖https://i.imgur.com/UtH2g7v.jpg, [爆卦] Step by Step 求解校正回歸
推

taylor060705/22 20:58樣本太少 這樣泛化能力應該挺差的

推

hshs30505/22 21:00推...神人快跪

推

aletheia05/22 21:05笑死

推

clsied05/22 21:06我前天剛考完你不要這樣

→

hw105/22 21:06幹...看不董辣

噓

thelittleone05/22 21:07你老師沒教你要寫結論嗎

推

duriamon05/22 21:13電梯向上

噓

SeaForest05/22 21:14濫用爆卦

→

FreedomTrail05/22 21:16XD目標波形校正教程XD

推

Tianmu05/22 21:17假設錯誤。 因為有外力改變絕對沒辦法用公式推算。 中間

→

Tianmu05/22 21:18已經有增加防疫為三級。

推

MrCool556605/22 21:19霸氣趨緩

推

parkeryang05/22 21:38推z transform

→

kennetharies05/22 21:54講中文好嗎！XD

推

haniah05/22 21:59痾 看不懂

推

champagnechu05/22 22:32模型結果推測明天確診數會是多少？

推

Sam045305/22 22:35求幫預測明天數字

推

waterydan05/22 23:17吃飽太閒喔= =

推

hodj05/22 23:24浪費才能系列開始

推

iam30719 05/24 15:20太神啦 不寫篇論文嗎 讓國際看看台灣的

Gossiping最新熱門推薦
🔥🔥🔥