楼主: diane1005
11174 1

[问答] 求助关于matlab中fminunc的用法 [推广有奖]

  • 2关注
  • 0粉丝

本科生

91%

还不是VIP/贵宾

-

威望
0
论坛币
415 个
通用积分
0.9000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
513 点
帖子
32
精华
0
在线时间
213 小时
注册时间
2011-12-11
最后登录
2024-3-6

楼主
diane1005 发表于 2014-5-7 16:22:51 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
我想用matlab中的fminunc命令做最大似然估计,请问我这样写的错误在哪里?
mylog_lik.m里的code是:

function val=mylog_lik(B,omega)

t = 200;

T = 280;

cov1 = [5.2486 1.9729 0.7363;1.9729 3.0687 0.9202;0.7363 0.9202 5.3797];

cov2 = [37.7688 8.0768 8.5434;8.0768 6.4612 6.9920;8.5434 6.9920 38.9579];

val=(t-1)*(log(det(B*B'))+trace(inv(cov1*(B*B'))))/2+(T-t+1)*(log(det(B*omega*B'))+trace(cov2*(B*omega*B')^(-1)))/2;


maxlike.m里是

B0=[0.1 0.1 0.1;0.1 0.1 0.1;0.1 0.1 0.1];

omega0 = [0.1 0 0 ; 0 0.1 0; 0 0 0.1];

x0 = [B0,omega0];


[x,fval,exitflag,output,grad,hessian]=fminunc('mylog_lik',x0);


初始值应该根据什么设置呢?初始值应该是向量吗?



二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:fminunc MATLAB matla atlab FMI matlab

luetkepohl(2007) .pdf
下载链接: https://bbs.pinggu.org/a-1544879.html

330.58 KB

我的likelihood function是根据这篇文章编的

沙发
Xaero 发表于 2014-5-8 16:33:41
把第一行全行改为:

function val=mylog_lik(X)
B=X(1);omega=X(2);
十年一觉扬州梦。
智不足以Academy,才尚不够Industry,[情无力于Life]。

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-21 13:01