版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

makechange74

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 322.2
帖子: 70
在線: 84.2小時(shí)
蟲號: 3646387
注冊: 2015-01-15
專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)

[求助] matlab用fmincon函數(shù)求極大似然估計(jì)

%l1_s分別是不同實(shí)驗(yàn)條件下的數(shù)據(jù)，下面的函數(shù)為極大既然函數(shù)的負(fù)值，先用fmincon求解，運(yùn)行出錯(cuò)了。
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)
提示錯(cuò)誤如下：
錯(cuò)誤使用 fmincon (line 637)
User supplied objective function must return a scalar value.

回復(fù)此樓

» 猜你喜歡

0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有14人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有12人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有13人回復(fù)
288求調(diào)劑已經(jīng)有14人回復(fù)
307求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有5人回復(fù)
288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有5人回復(fù)

1樓 2015-03-10 10:38:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

優(yōu)化問題，對于非線性約束尋優(yōu)問題除了可以用fmincon函數(shù)外還可以用什么？已經(jīng)有4人回復(fù)
Matlab中的fmincon功能求助已經(jīng)有16人回復(fù)
matlab中fmincon的fun是矩陣運(yùn)算怎么辦？已經(jīng)有5人回復(fù)
用matlab求bessel函數(shù)的根已經(jīng)有6人回復(fù)
matlab最小二乘去噪求助已經(jīng)有9人回復(fù)
Matlab最小二乘參數(shù)優(yōu)化已經(jīng)有3人回復(fù)
matlab fmincon優(yōu)化函數(shù) 怎樣知道每次得到的迭代點(diǎn) 已經(jīng)有3人回復(fù)
小弟新學(xué)matlab，想借用 fmincon求個(gè)函數(shù)的最小值，求助錯(cuò)誤該法，感謝！已經(jīng)有4人回復(fù)
【求助】用matlab對多項(xiàng)式函數(shù)求導(dǎo) 已經(jīng)有3人回復(fù)
用matlab對多項(xiàng)式函數(shù)求導(dǎo) 已經(jīng)有1人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 makechange74 的主題更新