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

銀蟲 (初入文壇)

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

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

我優(yōu)化程序的本意是想在兩條曲線取j個(gè)點(diǎn)，在y相等的情況下，兩條曲線對(duì)應(yīng)的x的差的總和值最小，但是好像目標(biāo)函數(shù)太復(fù)雜，計(jì)算出的結(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

下面這個(gè)是目標(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ù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(fù)

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

matlab最優(yōu)化問題求解中，復(fù)雜的目標(biāo)函數(shù)和約束條件無法直接表示的問題已經(jīng)有3人回復(fù)
matlab求函數(shù)積分出錯(cuò) 已經(jīng)有12人回復(fù)
一個(gè)復(fù)雜函數(shù)的matlab繪圖已經(jīng)有6人回復(fù)
隱函數(shù)積分問題的matlab編程求助已經(jīng)有3人回復(fù)
求matlab高手幫忙畫復(fù)雜函數(shù)的函數(shù)圖像已經(jīng)有6人回復(fù)
Matlab 已知函數(shù)y值求對(duì)應(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高手幫忙做個(gè)自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，已經(jīng)有25人回復(fù)
求助MATLAB隱函數(shù)等已經(jīng)有7人回復(fù)
幫忙給看看這個(gè)matlab優(yōu)化函數(shù) 問題已經(jīng)有8人回復(fù)
求matlab畫一個(gè)函數(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ù)
問個(gè)關(guān)于matlab里find函數(shù)的問題！已經(jīng)有4人回復(fù)
matlab 求解隱式常微分方程時(shí)用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復(fù)
MATLAB求解一個(gè)比較復(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的回帖

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

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 600.9
帖子: 27
在線: 12.4小時(shí)
蟲號(hào): 2019575
注冊: 2012-09-21
性別: GG
專業(yè): 制造系統(tǒng)與自動(dò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小時(shí)
蟲號(hào): 291104
注冊: 2006-10-28

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

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

贊一下

回復(fù)此樓

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

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

月只藍(lán)

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

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

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

???????:

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

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

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

贊一下

回復(fù)此樓

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

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

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

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

銀蟲 (初入文壇)

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

引用回帖:

3樓: Originally posted by dingd at 2014-12-08 16:45:23
把你的問題用文字圖片簡明扼要描述清楚，已知參數(shù)、數(shù)據(jù)、目標(biāo)函數(shù)、約束條件、待求參數(shù)等，直接放Matlab代碼不知道有幾個(gè)人愿意根據(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個(gè)點(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ù)雜，計(jì)算出的結(jié)果就迭代一次，這個(gè)您看懂了嗎，求高手幫忙

贊一下

回復(fù)此樓

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

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

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

銀蟲 (初入文壇)

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

引用回帖:

4樓: Originally posted by 月只藍(lán) at 2014-12-08 17:05:49
用遺傳算法函數(shù)ga試試。
或者換用1stopt軟件試試。...

我試過用ga算法了，但是結(jié)果依舊不理想，

，另外那個(gè)軟件我沒用過呀

贊一下

回復(fù)此樓

6樓2014-12-09 14:38:54

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

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

[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +11	手機(jī)用戶 2026-03-17	12/600	2026-03-23 23:51 by 熱情沙漠
[考研] 求調(diào)劑 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 熱情沙漠
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng