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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

如瑜?運(yùn)?/a>

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 600.9
帖子: 27
在線: 12.4小時
蟲號: 2019575
注冊: 2012-09-21
性別: GG
專業(yè): 制造系統(tǒng)與自動化

[求助] matlab優(yōu)化設(shè)計目標(biāo)函數(shù)復(fù)雜，求幫助已有2人參與

我優(yōu)化程序的本意是想在兩條曲線取j個點(diǎn)，在y相等的情況下，兩條曲線對應(yīng)的x的差的總和值最小，但是好像目標(biāo)函數(shù)太復(fù)雜，計算出的結(jié)果就迭代一次，根本就不是想要的結(jié)果，求高手幫忙呀，看看程序的問題再哪里，我用的是fmincon的優(yōu)化，下面是我的程序
主函數(shù)

CODE:

clc;clear;

%format long;

A=[-1 0;0 -1];b=[0;0];

%A=[0 -zp/rp;0 zp/rp;1 0;-1 0];b=[-c1;c2;rp*(sin(pi/zp))/d1;-rp*(sin(pi/zp))/d2];

%lb=[1;0.6];ub=[7;3];

[x,fval,exitflag,output]=fmincon(@xiuxingOptimFun,[0.1 0.1],A,b,[],[],[0;0],[1.5;1.5],@xiuxingConFun)

%畫圖

rp=72.5;

rrp=4;

a=1.4;

zc=43;

zp=44;

data=0.14*pi/180;

kd1=a*zp/(rp);

for pha=0:0.002*pi/zc:pi/zc

s1=1+kd1^2-2*kd1.*cos(pha);

i=zp/zc;

xd5=(rp-rrp.*s1.^(-1/2)).*cos((1-i).*pha-data)-(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*cos(i.*pha+data);

yd5=(rp-rrp.*s1.^(-1/2)).*sin((1-i).*pha-data)+(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*sin(i.*pha+data);

xd7=(rp-rrp.*s1.^(-1/2)).*cos((1-i).*pha)-(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*cos(i.*pha);

yd7=(rp-rrp.*s1.^(-1/2)).*sin((1-i).*pha)+(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*sin(i.*pha);

%plot(x,y,'-','linewidth',3);

%plot(x,y,'-');

%hold on

kd2=a*zp/(rp+x(2));

s2=1+kd2^2-2*kd2.*cos(pha);

xd6=(rp+x(2)-(rrp+x(1)).*s2.^(-1/2)).*cos((1-i).*pha)-(a/(rp+x(2))).*(rp+x(2)-zp*(rrp+x(1)).*s2.^(-1/2)).*cos(i.*pha);

yd6=(rp+x(2)-(rrp+x(1)).*s2.^(-1/2)).*sin((1-i).*pha)+(a/(rp+x(2))).*(rp+x(2)-zp*(rrp+x(1)).*s2.^(-1/2)).*sin(i.*pha);

plot(xd5,yd5,'r');

hold on

plot(xd6,yd6);

hold on

plot(xd7,yd7,'g');

hold on

end

下面這個是目標(biāo)函數(shù)

CODE:

function y=xiuxingOptimFun(x)

rp=72.5;

rrp=4;

a=1.4;

zc=43;

zp=44;

P=zeros(1,10000);

xd1=P;

yd1=P;

xd2=P;

yd2=P;

yd3=P;

yd4=0;

m=1;

j=1;

data=0.14*pi/180;

kd1=a*zp/(rp);

for pha=0:0.0002*pi/zc:pi/zc

s1=1+kd1^2-2.*kd1.*cos(pha);

i=zp/zc;

xd1(j)=(rp-rrp.*s1.^(-1/2)).*cos((1-i).*pha-data)-(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*cos(i.*pha+data);

yd1(j)=(rp-rrp.*s1.^(-1/2)).*sin((1-i).*pha-data)+(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*sin(i.*pha+data);

%plot(x,y,'-','linewidth',3);

%plot(x,y,'-');

%hold on

kd2=a*zp./(rp+x(2));

s2=1+kd2^2-2.*kd2.*cos(pha);

xd2(j)=(rp+x(2)-(rrp+x(1)).*s2.^(-1/2)).*cos((1-i).*pha)-(a/(rp+x(2))).*(rp+x(2)-zp*(rrp+x(1)).*s2.^(-1/2)).*cos(i.*pha);

yd2(j)=(rp+x(2)-(rrp+x(1)).*s2.^(-1/2)).*sin((1-i).*pha)+(a/(rp+x(2))).*(rp+x(2)-zp*(rrp+x(1)).*s2.^(-1/2)).*sin(i.*pha);

j=j+1;

end

for k=1:1:j-1

    for o=1:1:j-1

       if(yd1(k)-yd2(o)<=0.01)

          yd3(m)=abs(xd1(k)-xd2(o));

          yd4=yd4+yd3(m);

          m=m+1;

       end

    end

end

%plot(x,y,'-','linewidth',3);

%plot(x,y,'-');

%hold on

y=yd4;

約束條件

CODE:

function [c,ceq]=xiuxingConFun(x)

c=[];

ceq=x(1)-x(2)-0.02;

回復(fù)此樓

» 猜你喜歡

317分一志愿南理工材料工程本科湖工大求調(diào)劑已經(jīng)有11人回復(fù)
332求調(diào)劑已經(jīng)有12人回復(fù)
材料專碩 085600求調(diào)劑已經(jīng)有5人回復(fù)
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可已經(jīng)有14人回復(fù)
283求調(diào)劑（080500）已經(jīng)有14人回復(fù)
材料工程專碩求調(diào)劑已經(jīng)有6人回復(fù)
材料專碩調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑有機(jī)化學(xué)考研356分已經(jīng)有11人回復(fù)
環(huán)境科學(xué)與工程334分求調(diào)劑已經(jīng)有8人回復(fù)
考研調(diào)劑已經(jīng)有5人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab最優(yōu)化問題求解中，復(fù)雜的目標(biāo)函數(shù)和約束條件無法直接表示的問題已經(jīng)有3人回復(fù)
matlab求函數(shù)積分出錯已經(jīng)有12人回復(fù)
一個復(fù)雜函數(shù)的matlab繪圖已經(jīng)有6人回復(fù)
隱函數(shù)積分問題的matlab編程求助已經(jīng)有3人回復(fù)
求matlab高手幫忙畫復(fù)雜函數(shù)的函數(shù)圖像已經(jīng)有6人回復(fù)
Matlab 已知函數(shù)y值求對應(yīng)自變量x值已經(jīng)有11人回復(fù)
求助matlab的發(fā)fgoalattain函數(shù) 已經(jīng)有4人回復(fù)
求助如何用Matlab求解復(fù)雜的不規(guī)則積分區(qū)間的二重定積分的數(shù)值解？已經(jīng)有3人回復(fù)
matlab高手幫忙做個自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，已經(jīng)有25人回復(fù)
求助MATLAB隱函數(shù)等已經(jīng)有7人回復(fù)
幫忙給看看這個matlab優(yōu)化函數(shù) 問題已經(jīng)有8人回復(fù)
求matlab畫一個函數(shù)的圖已經(jīng)有6人回復(fù)
利用粒子群算法做多目標(biāo)優(yōu)化問題已經(jīng)有7人回復(fù)
matlab求解多元函數(shù)微分學(xué)問題已經(jīng)有4人回復(fù)
求助matlab求解三角函數(shù)方程已經(jīng)有9人回復(fù)
matlab代碼定義函數(shù)問題已經(jīng)有3人回復(fù)
問個關(guān)于matlab里find函數(shù)的問題！已經(jīng)有4人回復(fù)
matlab 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復(fù)
MATLAB求解一個比較復(fù)雜的方程已經(jīng)有14人回復(fù)
【求助】關(guān)于有MATLAB求復(fù)雜的目標(biāo)函數(shù)的最優(yōu)解的問題？萬分感謝！已經(jīng)有6人回復(fù)
【求助】關(guān)于matlab中的函數(shù)代換已經(jīng)有10人回復(fù)

1樓 2014-12-05 09:48:07

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

???????:

2?: Originally posted by ?????? at 2014-12-08 15:14:59

????????????????

???????????ga?????
???????1stopt?????????

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2014-12-08 17:05:49

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

查看全部 6 個回答

如瑜?運(yùn)?/a>

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 600.9
帖子: 27
在線: 12.4小時
蟲號: 2019575
注冊: 2012-09-21
性別: GG
專業(yè): 制造系統(tǒng)與自動化

怎么沒人理呀，求幫助

贊一下

回復(fù)此樓

2樓2014-12-08 15:14:59

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

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

把你的問題用文字圖片簡明扼要描述清楚，已知參數(shù)、數(shù)據(jù)、目標(biāo)函數(shù)、約束條件、待求參數(shù)等，直接放Matlab代碼不知道有幾個人愿意根據(jù)代碼去還原原問題？

贊一下

回復(fù)此樓

3樓2014-12-08 16:45:23

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

如瑜?運(yùn)?/a>

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 600.9
帖子: 27
在線: 12.4小時
蟲號: 2019575
注冊: 2012-09-21
性別: GG
專業(yè): 制造系統(tǒng)與自動化

引用回帖:

3樓: Originally posted by dingd at 2014-12-08 16:45:23
把你的問題用文字圖片簡明扼要描述清楚，已知參數(shù)、數(shù)據(jù)、目標(biāo)函數(shù)、約束條件、待求參數(shù)等，直接放Matlab代碼不知道有幾個人愿意根據(jù)代碼去還原原問題？

您好，優(yōu)化程序的本意是曲線1(xd1,yd1)，與曲線2（xd2，yd2），方程為
xd1(j)=(rp-rrp.*s1.^(-1/2)).*cos((1-i).*pha-data)-(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*cos(i.*pha+data);
yd1(j)=(rp-rrp.*s1.^(-1/2)).*sin((1-i).*pha-data)+(a/rp).*(rp-zp*(rrp).*s1.^(-1/2)).*sin(i.*pha+data);
xd2(j)=(rp+x(2)-(rrp+x(1)).*s2.^(-1/2)).*cos((1-i).*pha)-(a/(rp+x(2))).*(rp+x(2)-zp*(rrp+x(1)).*s2.^(-1/2)).*cos(i.*pha);
yd2(j)=(rp+x(2)-(rrp+x(1)).*s2.^(-1/2)).*sin((1-i).*pha)+(a/(rp+x(2))).*(rp+x(2)-zp*(rrp+x(1)).*s2.^(-1/2)).*sin(i.*pha);

這兩段曲線上面分別有j個點(diǎn)，使用以下的循環(huán)嵌套，判斷在yd1=yd2相等的情況下，xd1-xd2的差值加起來最小
for k=1:1:j-1
for o=1:1:j-1
   if(yd1(k)-yd2(o)<=0.01)
      yd3(m)=abs(xd1(k)-xd2(o));
      yd4=yd4+yd3(m);
      m=m+1;
   end
end
end
最后輸出yd4的值作為目標(biāo)函數(shù)，但是好像目標(biāo)函數(shù)太復(fù)雜，計算出的結(jié)果就迭代一次，這個您看懂了嗎，求高手幫忙

贊一下

回復(fù)此樓

5樓2014-12-09 14:36:42

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

查看全部 6 個回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 297求調(diào)劑 +13	田洪有 2026-03-26	13/650	2026-03-30 11:33 by 替代品000
[論文投稿] arxiv背書人郵件求助 20+3	redcom 2026-03-25	7/350	2026-03-30 10:56 by bobvan
[考研] 材料與化工272求調(diào)劑 +21	阿斯蒂芬2004 2026-03-28	21/1050	2026-03-30 10:52 by 晴空210210
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 材料與化工328分調(diào)劑 +9	。，。，。，。i 2026-03-23	9/450	2026-03-30 07:30 by 無際的草原
[考研] 327求調(diào)劑 +4	小卡不卡. 2026-03-29	4/200	2026-03-30 06:15 by wxiongid
[考研] 086000生物與醫(yī)藥調(diào)劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 0856求調(diào)劑 +13	zhn03 2026-03-25	14/700	2026-03-29 08:13 by fmesaito
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 311（085601）求調(diào)劑 +4	liziyeyeye 2026-03-28	4/200	2026-03-28 18:50 by 535743368
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8	曼111 2026-03-24	9/450	2026-03-28 07:58 by YYYYX1234
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4	手機(jī)用戶 2026-03-26	4/200	2026-03-28 07:17 by mmm just
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

如瑜?運(yùn)?/a>

[求助] matlab優(yōu)化設(shè)計目標(biāo)函數(shù)復(fù)雜，求幫助 已有2人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

月只藍(lán)

【答案】應(yīng)助回帖

如瑜?運(yùn)?/a>

dingd

【答案】應(yīng)助回帖

如瑜?運(yùn)?/a>

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] matlab優(yōu)化設(shè)計目標(biāo)函數(shù)復(fù)雜，求幫助已有2人參與