如瑜?運?/a>

銀蟲 (初入文壇)

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

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

我優(yōu)化程序的本意是想在兩條曲線取j個點，在y相等的情況下，兩條曲線對應的x的差的總和值最小，但是好像目標函數(shù)太復雜，計算出的結果就迭代一次，根本就不是想要的結果，求高手幫忙呀，看看程序的問題再哪里，我用的是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

下面這個是目標函數(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;

回復此樓

» 猜你喜歡

0703化學調(diào)劑，求導師收已經(jīng)有3人回復
生物學學碩求調(diào)劑已經(jīng)有6人回復
305分求調(diào)劑（食品工程）已經(jīng)有7人回復
材料與化工考研調(diào)劑已經(jīng)有5人回復
299求調(diào)劑已經(jīng)有6人回復
材料292調(diào)劑已經(jīng)有5人回復
085404電子信息284分求調(diào)劑已經(jīng)有4人回復
收08調(diào)劑生已經(jīng)有3人回復
081700 調(diào)劑 267分已經(jīng)有9人回復
293求調(diào)劑已經(jīng)有5人回復

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

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

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

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

如瑜?運?/a>

銀蟲 (初入文壇)

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

怎么沒人理呀，求幫助

贊一下

回復此樓

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

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

dingd

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

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

【答案】應助回帖

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

贊一下

回復此樓

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

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

月只藍

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

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

【答案】應助回帖

???????:

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

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

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

贊一下

回復此樓

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

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

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

如瑜?運?/a>

銀蟲 (初入文壇)

應助: 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ù)、目標函數(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個點，使用以下的循環(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的值作為目標函數(shù)，但是好像目標函數(shù)太復雜，計算出的結果就迭代一次，這個您看懂了嗎，求高手幫忙

贊一下

回復此樓

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

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

如瑜?運?/a>

銀蟲 (初入文壇)

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

引用回帖:

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

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

，另外那個軟件我沒用過呀

贊一下

回復此樓

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

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

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

[考研] 305分求調(diào)劑（食品工程） +5	Sxy112 2026-03-21	7/350	2026-03-24 12:27 by 544594351
[考研] 086003食品工程求調(diào)劑 +4	淼淼111 2026-03-24	4/200	2026-03-24 11:53 by 544594351
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 材料專業(yè)求調(diào)劑 +11	hanamiko 2026-03-18	11/550	2026-03-23 23:12 by peike
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[論文投稿] 急發(fā)核心期刊論文 +3	賢達問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 316求調(diào)劑 +7	梁茜雯 2026-03-19	7/350	2026-03-23 16:21 by lingjue
[考研] 招08考數(shù)學 +6	laoshidan 2026-03-20	14/700	2026-03-23 14:37 by 15614359529
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 一志愿北京化工大學070300 學碩336求調(diào)劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 一志愿蘇州大學材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 一志愿南京航空航天大學大學，080500材料科學與工程學碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 材料學碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 一志愿福大288有機化學，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨