版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

qiaoye32803

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 27.5
帖子: 7
在線: 2.3小時
蟲號: 2124141
注冊: 2012-11-13
專業(yè): 制造系統(tǒng)與自動化

[求助] 函數(shù)調用

這是我的目標函數(shù)，里面有隨機生成數(shù):
function opforpaper=lq_paper(x)
syms theta;
mu_ra=0.0645;
mu_rb=0.0889;
mu_rc=0.043;
sigma_ra=0.008;
sigma_rb=0.07;
sigma_rc=0.004;

for i=1:100
f_1(i)=cos(theta)*(x(1)+x(3)*normrnd(0,1)+normrnd(mu_rb,sigma_rb)*cos(unifrnd(0,2*pi))-normrnd(mu_ra,sigma_ra)*cos(unifrnd(0,2*pi)))+sqrt((x(2)+x(4)*normrnd(0,1)+normrnd(mu_rc,sigma_rc)*cos(unifrnd(0,2*pi))).^2-(x(1)+x(3)*normrnd(0,1)+normrnd(mu_rb,sigma_rb)*cos(unifrnd(0,2*pi))-normrnd(mu_ra,sigma_ra)*sin(theta)).^2);
f_2(i)=cos(theta+0.0063)*(x(1)+x(3)*normrnd(0,1)+normrnd(mu_rb,sigma_rb)*cos(unifrnd(0,2*pi))-normrnd(mu_ra,sigma_ra)*cos(unifrnd(0,2*pi)))+sqrt((x(2)+x(4)*normrnd(0,1)+normrnd(mu_rc,sigma_rc)*cos(unifrnd(0,2*pi))).^2-(x(1)+x(3)*normrnd(0,1)+normrnd(mu_rb,sigma_rb)*cos(unifrnd(0,2*pi))-normrnd(mu_ra,sigma_ra)*sin(theta+0.0063)).^2);
f3(i)=abs(f_1(i)-f_2(i));
end
opforpaper=mean(f3);
下面的是約束，
function [c,ceq]=consfun(x)
M=384;
c=[(M*sqrt(M^2+f3^2)-M)/(2*sqrt(M^2+f3^2))-0.05;f3*M/(2*sqrt(M^2+f3^2))-0.5];
ceq=[];
下面的是所用的計算函數(shù)：
clear
x0=[50,335,0.005,0.05];
lb=[45,330,0.002,0.007];
ub=[55,340,0.02,0.07];
options=optimset('largescale','off','Algorithm','active-set');
[x,fval]=fmincon('lq_paper',x0,[],[],[],[],lb,ub,'consfun',options)
現(xiàn)在有個問題：在約束中的f3必須是目標函數(shù)中生成的那個f3，因為目標函數(shù)中的f3含有隨機數(shù)，請問哪位大神能幫我下，在約束中調用目標函數(shù)中的那個f3！

回復此樓

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有254人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。已經(jīng)有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經(jīng)有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人已經(jīng)有0人回復
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8 已經(jīng)有0人回復

1樓 2012-11-19 21:15:33

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

感謝參與，應助指數(shù) +1

用全局變量試試，應該可以的吧

贊一下

回復此樓

showmethemoney

2樓2012-11-20 12:01:26