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

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

應(yīng)助: 0 (幼兒園)
金幣: 600.9
帖子: 27
在線: 12.4小時(shí)
蟲(chóng)號(hào): 2019575
注冊(cè): 2012-09-21
性別: GG
專(zhuān)業(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é)果，求高手幫忙呀，看看程序的問(wèn)題再哪里，我用的是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)

%畫(huà)圖

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ù)此樓

» 猜你喜歡

一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專(zhuān)業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)
327求調(diào)劑已經(jīng)有3人回復(fù)

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

matlab最優(yōu)化問(wèn)題求解中，復(fù)雜的目標(biāo)函數(shù)和約束條件無(wú)法直接表示的問(wèn)題已經(jīng)有3人回復(fù)
matlab求函數(shù)積分出錯(cuò) 已經(jīng)有12人回復(fù)
一個(gè)復(fù)雜函數(shù)的matlab繪圖已經(jīng)有6人回復(fù)
隱函數(shù)積分問(wèn)題的matlab編程求助已經(jīng)有3人回復(fù)
求matlab高手幫忙畫(huà)復(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ù) 問(wèn)題已經(jīng)有8人回復(fù)
求matlab畫(huà)一個(gè)函數(shù)的圖已經(jīng)有6人回復(fù)
利用粒子群算法做多目標(biāo)優(yōu)化問(wèn)題已經(jīng)有7人回復(fù)
matlab求解多元函數(shù)微分學(xué)問(wèn)題已經(jīng)有4人回復(fù)
求助matlab求解三角函數(shù)方程已經(jīng)有9人回復(fù)
matlab代碼定義函數(shù)問(wèn)題已經(jīng)有3人回復(fù)
問(wèn)個(gè)關(guān)于matlab里find函數(shù)的問(wèn)題！已經(jīng)有4人回復(fù)
matlab 求解隱式常微分方程時(shí)用ode15i 函數(shù)出現(xiàn)的問(wèn)題求助已經(jīng)有3人回復(fù)
MATLAB求解一個(gè)比較復(fù)雜的方程已經(jīng)有14人回復(fù)
【求助】關(guān)于有MATLAB求復(fù)雜的目標(biāo)函數(shù)的最優(yōu)解的問(wèn)題？萬(wàn)分感謝！已經(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的回帖

dingd

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

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

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

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

贊一下

回復(fù)此樓

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

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

查看全部 6 個(gè)回答

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

銀蟲(chóng) (初入文壇)

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

怎么沒(méi)人理呀，求幫助

贊一下

回復(fù)此樓

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

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

月只藍(lán)

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

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

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

???????:

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

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

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

贊一下

回復(fù)此樓

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

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

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

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

銀蟲(chóng) (初入文壇)

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

引用回帖:

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

您好，優(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的差值加起來(lái)最小
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的回帖

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

[考研] 085600 286分材料求調(diào)劑 +10	麻辣魷魚(yú) 2026-03-27	11/550	2026-03-29 21:52 by nanaliuyun
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 281求調(diào)劑 +4	亞克西good 2026-03-26	6/300	2026-03-29 21:06 by 啊李999
[考研] 283求調(diào)劑（080500） +7	A child 2026-03-27	7/350	2026-03-29 20:46 by 唐沐兒
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +8	舉個(gè)栗子oi 2026-03-24	9/450	2026-03-29 13:08 by peike
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 一志愿廈門(mén)大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 311求調(diào)劑 +9	lin0039 2026-03-26	9/450	2026-03-28 13:05 by 唐沐兒
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚(yú)的貓
[考研] 324求調(diào)劑 +8	hanamiko 2026-03-26	10/500	2026-03-27 08:06 by hypershenger
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無(wú)際的草原
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 材料專(zhuān)碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt