[討論]關於fminsearch 求多變數的最小

作者: cglineage (超級阿嘴)   2014-12-26 16:42:18
小弟目前遇到的問題是
要做用mle來估計參數的的問題
目前在真實資料輸入後是一個6x6的矩陣
而這6x6的矩陣在帶入分配之後,會形成了一個也是6x6的符號矩陣
每一個元素中都有 A B C三個參數
再將矩陣中的每個元素用sum總和去算則變成一個 三個變數的多項式
目前 是用很拙劣的方法(記事本) 將矩陣總合的式子中的ABC
替換成x(1) x(2) x(3)
將整串複製入 fun.m
再用fminsearch 跑
[x,fval] = fminsearch(@fun, x0)
藉此得到一組參數和min
但在模擬資料的時候就出現問題了,模擬一千筆資料就會顯得我很智障
但查書之後並無找到fminserch有我理想中的用法
是否能直接使用命令視窗中所算出來的式子(符號變數是ABC)
用fminsearch
直接求出本來A,B,C的Min
作者: H45 (!H45)   2014-12-26 17:11:00
技術上可以,再怎麼難辦,寫個.m檔產生一千個fun.m檔總可以吧再說還有feval這個函數可以用阿抱歉,跟feval沒關係
作者: cglineage (超級阿嘴)   2014-12-26 17:20:00
可是每次都要複製到記事本改,再複製過去..
作者: H45 (!H45)   2014-12-26 20:47:00
不用吧,你的符號是什麼型式的? 總可以轉為double吧?再怎麼沒辦法,也可以用MATLAB讀取記事本資料才是阿抱歉,如果可以直接轉double,那複製到記事本這步就沒必要
作者: mp19990920 (郝渴連)   2014-12-27 08:54:00
您所述的 ABC "參數" 意思是用 syms 符號型態嗎?
作者: cglineage (超級阿嘴)   2014-12-27 16:56:00
沒錯在前面就是用syms定的樓樓上我試試看算完之後再用另個程式讀記事本看看是否可以少一個步驟好了
作者: windg (風)   2014-12-29 10:52:00
http://www.mathworks.com/help/optim/examples/using-symbolic-mathematics-with-optimization-toolbox-solvers.htm

Links booklink

Contact Us: admin [ a t ] ucptt.com