顯示廣告
隱藏 ✕
回列表(←)
Disp BBS guest 註冊 登入(i) 線上人數: 859
首頁(home) 上頁(↑) 下頁(↓) 末頁(end)
※ 本文為 MindOcean 轉寄自 ptt.cc 更新時間: 2021-05-22 23:39:12
看板 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
噓 MotleyCrue: 文組會來噓妳1F 05/22 20:37
→ guitar0225: 你為什麼要讓我想起來！！2F 05/22 20:37
→ htc10: https://i.imgur.com/DvUsbWT.jpg3F 05/22 20:37
推 appleball200: 推XD4F 05/22 20:38
→ kumori: 跨謀拉 講重點好嗎5F 05/22 20:38
推 sandiato: 這個可能要找t2d來解答6F 05/22 20:38
推 sam90205: 太鬼了吧7F 05/22 20:38
推 hipmyhop: 笑死8F 05/22 20:39
推 alcpeon911: 笑死 不要這麼認真好嗎9F 05/22 20:39
→ ridecule: 文組應該是推吧 因為看不懂 好神奇10F 05/22 20:40
推 lightaurora: 你是鬼吧11F 05/22 20:40
噓 weekman: 用拉式轉換前提要是線性系統12F 05/22 20:40
→ Ceferino: 直接說你覺得今天實際數多少13F 05/22 20:40
噓 xiaoyao: 老闆：講重點14F 05/22 20:40
推 Neisseria: 問教授可不可以當碩論啊15F 05/22 20:40
推 plusonezero: 嗯嗯，跟我想的一樣16F 05/22 20:40
推 Cwren: 笑死17F 05/22 20:40
→ rdmtb: 太認真了。18F 05/22 20:40
推 billliu: ID 好像女生要不要校正回歸一下~~19F 05/22 20:40
推 JackTheRippe: 快推 不然會被說看不懂20F 05/22 20:41
推 aivxg: 跟我想的一樣 笑死21F 05/22 20:41
噓 lucifiel1618: 沒料的就愛貼一堆圖講不出東西22F 05/22 20:41
推 strongrain: 你太認真了吧XD23F 05/22 20:42
噓 f17690815: 講中文24F 05/22 20:47
推 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
推 AlianF: 公三小29F 05/22 20:55
→ NCUking: https://i.imgur.com/UtH2g7v.jpg30F 05/22 20:57
推 taylor0607: 樣本太少 這樣泛化能力應該挺差的31F 05/22 20:58
推 hshs305: 推...神人快跪32F 05/22 21:00
推 aletheia: 笑死33F 05/22 21:05
推 clsied: 我前天剛考完你不要這樣34F 05/22 21:06
→ hw1: 幹...看不董辣35F 05/22 21:06
噓 thelittleone: 你老師沒教你要寫結論嗎36F 05/22 21:07
推 duriamon: 電梯向上37F 05/22 21:13
噓 SeaForest: 濫用爆卦38F 05/22 21:14
→ FreedomTrail: XD目標波形校正教程XD39F 05/22 21:16
推 Tianmu: 假設錯誤。 因為有外力改變絕對沒辦法用公式推算。 中間40F 05/22 21:17
→ Tianmu: 已經有增加防疫為三級。
推 MrCool5566: 霸氣趨緩42F 05/22 21:19
推 parkeryang: 推z transform43F 05/22 21:38
→ kennetharies: 講中文好嗎！XD44F 05/22 21:54
推 haniah: 痾 看不懂45F 05/22 21:59
推 champagnechu: 模型結果推測明天確診數會是多少？46F 05/22 22:32
推 Sam0453: 求幫預測明天數字47F 05/22 22:35
推 waterydan: 吃飽太閒喔=  =48F 05/22 23:17
推 hodj: 浪費才能系列開始49F 05/22 23:24

--
※ 看板: Gossiping　文章推薦值: -1 目前人氣: 0 累積人氣: 192　
分享網址: 複製 已複製
(￣︿￣)ｐ JosephC0227 說瞎！
1樓 時間: 2021-05-23 00:32:28 (台灣) 
→
JosephC0227
　 05-23 00:32 TW 
這是白癡解釋吧！
r)回覆 e)編輯 d)刪除 M)收藏 ^x)轉錄 同主題: =)首篇 [)上篇 ])下篇