| 查看: 2386 | 回复: 37 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xizi19880211新虫 (初入文坛)
|
[求助]
1stopt求最小值 已有1人参与
|
||
| 如果是非线性约束求最值，并且约束条件里面的量是变化，比如是随时间变化的，可以用1stopt来实现么？MATlab的GA工具箱可以实现么？求实例代码。 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有185人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
xizi19880211
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 83.5
- 散金: 5
- 帖子: 42
- 在线: 17.6小时
- 虫号: 2752848
- 注册: 2013-10-25
- 专业: 电路与系统
|
是这个意思不？ Constant Pw=[25.33 25.33 25.67 25.67 33.33 34.67 38.33 33 28 21 14.67 14.33 4.43 5.32 10.33 9.88 8.53 4.43 3.33 3.33 5.03 8.01 5 9.03 10.99 10.53 12.52 20 23.98 30.99 30 28.33 27.53 21.56 27.01 30.32 46.98 46.07 55.65 51.55 59.88 53.66 52.03 47.55 40 53.09 54.99 53.22 50.02]; Pl=[15 15 15 15 15 15 28.99 28.99 28.99 28.99 28.99 28.99 16.55 16.55 16.55 16.55 9.59 9.59 9.59 9.59 9.59 9.59 9.59 9.59 9.59 23.33 23.33 23.33 23.33 23.33 23.33 23.33 23.33 45.88 45.88 45.88 45.88 45.88 45.88 45.88 45.88 45.88 49.55 49.55 49.55 49.55 30.89 30.89 30.89]; Parameter x(1)[50,],x(2)[100,]; MinFunction y=0.288*x(1)+0.000735*x(2); StartProgram; var i: Integer; Eb(0)=0.0012*x(1); Ec(0)=0.0000030375*x(2); Pb(1)=0.0012*x(1); Pc(1)=0.0000030375*x(2); for i:=1 to 49 do begin if Pw(i)-Pl(i)>0; Eb(i)=Eb(k-1)+Pb(i)*i*0.7; Ec(i)=Ec(i-1)+Pc(i)*i*0.98; else Eb(i)=Eb(i-1)-Pb(i)*i*0.8; Ec(i)=Ec(i-1)-Pc(i)*i*0.98; end Pb(i)=Eb(i)/i; Pc(i)=Ec(i)/i; end for i=1:49 if dp=Pw(i)*0.95-Pl(i)>0; Elsp(i)=0; else Elsp(i)=[Pl(i)-Pw(i)]*0.95*i; end LPSP++=Elsp(k)/Pl(k); end; Eb(k)*1000/1.2<=x(1)<=Eb(k)*1000/0.72; Ec(k)*100000/0.3037<=x(2)<=Ec(k)*1000000/0.267; x(1)<=0.7*[Pb(k)+Pc(k)]*1000/1.2; EndProgram; |
19楼2014-01-01 15:44:27
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2013-12-30 09:12:53
xizi19880211
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 83.5
- 散金: 5
- 帖子: 42
- 在线: 17.6小时
- 虫号: 2752848
- 注册: 2013-10-25
- 专业: 电路与系统
|
也就是说把约束条件那一块用循环来实现？ 代码在下面，帮我看看如何循环，谢谢 Pw=[25.3300000000000 25.3300000000000 25.6700000000000 25.6700000000000 33.3300000000000 34.6700000000000 38.3300000000000 33 28 21 14.6700000000000 14.3300000000000 4.43000000000000 5.32000000000000 10.3300000000000 9.88000000000000 8.53000000000000 4.43000000000000 3.33000000000000 3.33000000000000 5.03000000000000 8.01000000000000 5 9.03000000000000 10.9900000000000 10.5300000000000 12.52 20 23.98 30.99 30 28.33 27.53 21.56 27.01 30.32 46.98 46.07 55.65 51.55 59.88 53.66 52.03 47.55 40 53.09 54.99 53.22 50.02]; Pl=[15 15 15 15 15 15 28.9900000000000 28.9900000000000 28.9900000000000 28.9900000000000 28.9900000000000 28.9900000000000 16.5500000000000 16.5500000000000 16.5500000000000 16.5500000000000 9.59000000000000 9.59000000000000 9.59000000000000 9.59000000000000 9.59000000000000 9.59000000000000 9.59000000000000 9.59000000000000 9.59000000000000 23.3300000000000 23.3300000000000 23.3300000000000 23.3300000000000 23.3300000000000 23.3300000000000 23.3300000000000 23.3300000000000 45.8800000000000 45.8800000000000 45.8800000000000 45.8800000000000 45.8800000000000 45.8800000000000 45.8800000000000 45.8800000000000 45.8800000000000 49.5500000000000 49.5500000000000 49.5500000000000 49.5500000000000 30.8900000000000 30.8900000000000 30.8900000000000]; Eb(0)=1.2e-3*x(1); Ec(0)=3.0375e-6*x(2); for k=1:49 if Pw(k)-Pl(k)>0; Eb(k)=Eb(k-1)+Pb(k)*k*0.7; Ec(k)=Ec(k-1)+Pc(k)*k*0.98; else Eb(k)=Eb(k-1)-Pb(k)*k*0.8; Ec(k)=Ec(k-1)-Pc(k)*k*0.98; end end for k=1:49 if dp=Pw(k)*0.95-Pl(k)>0; Elsp(k)=0; else Elsp(k)=[Pl(k)-Pw(k)]*0.95*k; end LPSP++=Elsp(k)/Pl(k); end; MinFunction 0.288*x(1)+7.35e-4*x(2); Eb(k)*1000/1.2<=x(1)<=Eb(k)*1000/0.72; Ec(k)*100000/0.3037<=x(2)<=Ec(k)*1000000/0.267;x(1)<=0.7*[Pb(k)+Pc(k)]*1000/1.2; |
3楼2013-12-30 10:33:28
xizi19880211
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 83.5
- 散金: 5
- 帖子: 42
- 在线: 17.6小时
- 虫号: 2752848
- 注册: 2013-10-25
- 专业: 电路与系统
？ 5楼2013-12-30 20:11:44
