[討論] 畢旅
這次統算作業比較簡單,相信大家都可以自己寫出來
但以防萬一 (碰!! 阿!看錯,以為是一萬)
還是PO個想法給大家參考參考
先講第二題好了,
要用到R內建的函式simplex,可能大部分的人會不知道後面變數要怎麼打,
我就來教教大家
要先裝一個包裝,簡稱package,這邊要用的包裝是"boot",靴子
(指令為library(boot))
然後joke以用"?simplex"來查詢函式內容,可以看到變數有a、A1、b1、A2、b2
、A3、b3、maxi...
a 是指要maximize的函數的係數向量
A1 是限制式為"<="的係數矩陣
b1 是限制式"<="右邊的那個東東矩陣
A2 則跟A1相反是指">="的限制式係數矩陣
b2 是限制式">="右邊的那個西西矩陣
A3 應該是"="吧
b3 就是右邊那個吧,嗯!
maxi 如果是TURE就是找最大,如果是FALSE就是找最小
其他的我看不懂,
所以這次作業的話就類似以下這種寫法:
library(boot)
s=c(2,2,3)
e1=c(-2,1,1)
e2=c(4,-1,3)
e=rbind(e1,e2);d=(1,3)
simplex(s,e,d,maxi=T)
這樣應該就會跑出來了!!YA
再回來說第一題,
可以試試basic variables取1跟2 剩下就是nonbasic variables
你就會發現做出來的w_{3}、w_{4}、w_{5}都會大於0喔^.<
就 STOP!!!!!!!!!!!!!!
回去解出x,y帶進原函式得最大值,就降!
這次作業大概就是這麼一回事,
如果有任何疑問可以到410找Min-Han,隨時可以為您解答
又到時候講掰掰~12點以後不要打給我
因為我睡啦~~~~~~~~~~~~~~~~~~~~~~~~~~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 1.34.178.19
→
05/07 22:35, , 1F
05/07 22:35, 1F
討論串 (同標題文章)