| 查看: 3153 | 回复: 0 | ||
[求助]
matlab用fmincon函数求极大似然估计
|
|
%l1_s分别是不同实验条件下的数据，下面的函数为极大既然函数的负值，先用fmincon求解，运行出错了。 function [l1_s,l2_s,l3_s,l4_s,l5_s,l6_s] = readandevaluate(filename) global l1_s l2_s l3_s l4_s l5_s l6_s; filename = 'li_s.xlsx'; l1_s = xlsread(filename,'A:A'); l2_s = xlsread(filename,'B:B'); l3_s = xlsread(filename,'C:C'); l4_s = xlsread(filename,'D  ');l5_s = xlsread(filename,'E:E'); l6_s = xlsread(filename,'F:F'); end function f = fmin4mle(x) global l1_s l2_s l3_s l4_s l5_s l6_s; f = -(x(5)*l1_s.^(x(5)-1)*exp(-x(5)*(x(1)+x(2)/80+x(3)/20+x(4)*log(20)))*exp(-l1_s.^x(5)*exp(-x(5)*(x(1)+x(2)/80+x(3)/80+x(4)*log(20)))))... *(x(5)*l2_s.^(x(5)-1)*exp(-x(5)*(x(1)+x(2)/80+x(3)/60+x(4)*log(20)))*exp(-l2_s.^x(5)*exp(-x(5)*(x(1)+x(2)/80+x(3)/80+x(4)*log(60)))))... *(x(5)*l3_s.^(x(5)-1)*exp(-x(5)*(x(1)+x(2)/80+x(3)/40+x(4)*log(20)))*exp(-l3_s.^x(5)*exp(-x(5)*(x(1)+x(2)/55+x(3)/80+x(4)*log(40)))))... *(x(5)*l4_s.^(x(5)-1)*exp(-x(5)*(x(1)+x(2)/95+x(3)/20+x(4)*log(20)))*exp(-l4_s.^x(5)*exp(-x(5)*(x(1)+x(2)/55+x(3)/95+x(4)*log(20)))))... *(x(5)*l5_s.^(x(5)-1)*exp(-x(5)*(x(1)+x(2)/95+x(3)/60+x(4)*log(20)))*exp(-l5_s.^x(5)*exp(-x(5)*(x(1)+x(2)/70+x(3)/95+x(4)*log(60)))))... *(x(5)*l6_s.^(x(5)-1)*exp(-x(5)*(x(1)+x(2)/95+x(3)/40+x(4)*log(20)))*exp(-l6_s.^x(5)*exp(-x(5)*(x(1)+x(2)/70+x(3)/95+x(4)*log(40))))); end A = [0 0 0 0 -1]; b = [0]; x0 = [10,10,10,10,10]; [x,fval] = fmincon(@readandevaluate,x0,A,b) 提示错误如下： 错误使用 fmincon (line 637) User supplied objective function must return a scalar value. |
回复此楼» 猜你喜欢
281求调剂（0805）
已经有12人回复
-
一志愿西南交大，求调剂
已经有4人回复
-
0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287
已经有7人回复
-
297求调剂
已经有7人回复
-
材料专硕306英一数二
已经有13人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有10人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有6人回复
-
收复试调剂生
已经有4人回复
-
材料专硕274一志愿陕西师范大学求调剂
已经有6人回复
-
331求调剂（0703有机化学
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
优化问题，对于非线性约束寻优问题除了可以用fmincon函数外还可以用什么？
已经有4人回复
- Matlab中的fmincon功能求助
已经有16人回复
- matlab中fmincon的fun是矩阵运算怎么办？
已经有5人回复
- 用matlab求bessel函数的根
已经有6人回复
- matlab最小二乘去噪求助
已经有9人回复
- Matlab最小二乘参数优化
已经有3人回复
- matlab fmincon优化函数 怎样知道每次得到的迭代点
已经有3人回复
- 小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！
已经有4人回复
- 【求助】用matlab对多项式函数求导
已经有3人回复
- 用matlab对多项式函数求导
已经有1人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
