版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

如瑜?运�

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 600.9
帖子: 27
在线: 12.4小时
虫号: 2019575
注册: 2012-09-21
性别: GG
专业: 制造系统与自动化

[求助] matlab优化设计目标函数复杂，求帮助已有2人参与

我优化程序的本意是想在两条曲线取j个点，在y相等的情况下，两条曲线对应的x的差的总和值最小，但是好像目标函数太复杂，计算出的结果就迭代一次，根本就不是想要的结果，求高手帮忙呀，看看程序的问题再哪里，我用的是fmincon的优化，下面是我的程序
主函数

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

下面这个是目标函数

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;

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab最优化问题求解中，复杂的目标函数和约束条件无法直接表示的问题已经有3人回复
matlab求函数积分出错已经有12人回复
一个复杂函数的matlab绘图已经有6人回复
隐函数积分问题的matlab编程求助已经有3人回复
求matlab高手帮忙画复杂函数的函数图像已经有6人回复
Matlab 已知函数y值求对应自变量x值已经有11人回复
求助matlab的发fgoalattain函数已经有4人回复
求助如何用Matlab求解复杂的不规则积分区间的二重定积分的数值解？已经有3人回复
matlab高手帮忙做个自定义复杂函数拟合，含有复数，已经有25人回复
求助MATLAB隐函数等已经有7人回复
帮忙给看看这个matlab优化函数问题已经有8人回复
求matlab画一个函数的图已经有6人回复
利用粒子群算法做多目标优化问题已经有7人回复
matlab求解多元函数微分学问题已经有4人回复
求助matlab求解三角函数方程已经有9人回复
matlab代码定义函数问题已经有3人回复
问个关于matlab里find函数的问题！已经有4人回复
matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助已经有3人回复
MATLAB求解一个比较复杂的方程已经有14人回复
【求助】关于有MATLAB求复杂的目标函数的最优解的问题？万分感谢！已经有6人回复
【求助】关于matlab中的函数代换已经有10人回复

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

如瑜?运�

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 600.9
帖子: 27
在线: 12.4小时
虫号: 2019575
注册: 2012-09-21
性别: GG
专业: 制造系统与自动化

怎么没人理呀，求帮助

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

把你的问题用文字图片简明扼要描述清楚，已知参数、数据、目标函数、约束条件、待求参数等，直接放Matlab代码不知道有几个人愿意根据代码去还原原问题？

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

2楼: Originally posted by 如瑜?运� at 2014-12-08 15:14:59

怎么没人理呀，求帮助

用遗传算法函数ga试试。
或者换用1stopt软件试试。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

如瑜?运�

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 600.9
帖子: 27
在线: 12.4小时
虫号: 2019575
注册: 2012-09-21
性别: GG
专业: 制造系统与自动化

引用回帖:

3楼: Originally posted by dingd at 2014-12-08 16:45:23
把你的问题用文字图片简明扼要描述清楚，已知参数、数据、目标函数、约束条件、待求参数等，直接放Matlab代码不知道有几个人愿意根据代码去还原原问题？

您好，优化程序的本意是曲线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个点，使用以下的循环嵌套，判断在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的值作为目标函数，但是好像目标函数太复杂，计算出的结果就迭代一次，这个您看懂了吗，求高手帮忙

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

如瑜?运�

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 600.9
帖子: 27
在线: 12.4小时
虫号: 2019575
注册: 2012-09-21
性别: GG
专业: 制造系统与自动化

引用回帖:

4楼: Originally posted by 月只蓝 at 2014-12-08 17:05:49
用遗传算法函数ga试试。
或者换用1stopt软件试试。...

我试过用ga算法了，但是结果依旧不理想，

，另外那个软件我没用过呀

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主如瑜?运� 的主题更新

返回列表