| 查看: 3163 | 回复: 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. |
回复此楼» 猜你喜欢
293调剂
已经有3人回复
-
机械工程264学硕求调剂
已经有3人回复
-
264求调剂
已经有9人回复
-
求调剂
已经有6人回复
-
22408 266求调剂
已经有7人回复
-
求调剂
已经有14人回复
-
273求调剂
已经有43人回复
-
307求调剂
已经有13人回复
-
327求调剂
已经有7人回复
-
338求调剂
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
优化问题，对于非线性约束寻优问题除了可以用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
