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

登录注册

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

24小时热门版块排行榜

返回列表

bingbing5380

铁虫 (正式写手)

应助: 1 (幼儿园)
金币: 72.7
散金: 450
红花: 3
帖子: 331
在线: 144.1小时
虫号: 640275
注册: 2008-10-29
专业: 宏观经济管理与战略

[求助] 亲，毕业论文急用，调半天不知道错在哪？！求点醒啊~~（别被 F 吓到，只是有点长）

程序目的：设置一个扰动参数e，如果不满足退出条件，就把e变小（即：e乘以一个小于1的常数c）。知道满足退出条件位置。其中，x0为初始点，x1为新的迭代点，两个点之间的关系就是利用fmincon函数。

function x=FISHER(x0)
e=0.5;
c=0.2;
eps=1.0e-6; %可接受误差。
while 1
fx=@(x)[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9),x(10),x(11),x(12),x(13),x(14),x(15),x(16),x(17),x(18),x(19),x(20),x(21),x(22),x(23),x(24),x(25),x(26),x(27),x(28),x(29),x(30),x(31),x(32),x(33),x(34),x(35),x(36),x(37),x(38),x(39),x(40),x(41),x(42),x(43),x(44),x(45),x(46),x(47),x(48),x(49),x(50),x(51),x(52),x(53)]*[4/5;4/5;4/5;2/5;2/5;2/5;2/5;2/5;6/5;6/5;-2;-2;-2;2/5;2/5;2/5;2/5;6/5;6/5;6/5;-2;0;-2;-1;-1;-1;-1;-1;-1;zeros(24,1)];

options=optimset('LargeScale','off');
x1=fmincon(fx,x0,[],[],[],[],zeros(1,53),[],@mycon_fish,options); %x0与x1之间的迭代关系。

F=@(x)[3.8*x(1)+1.36*x(2)+3.2-0.6*x(34)+0.4*x(47)-0.4*x(37);3.8*x(2)+1.36*x(1)+3.2-0.6*x(35)+0.4*x(48)-0.4*x(38);3.08*x(3)+6-0.6*x(36)+0.4*x(49)-0.4*x(39);x(4)+3.5+x(34)+1.92*x(4)+1.92*x(6)+1.28*x(5)+1.28*x(7)-0.8*x(40)+0.2*x(50)-0.2*x(43);x(5)+3.5+x(34)+1.92*x(5)+1.92*x(7)+1.28*x(4)+1.28*x(6)-0.8*x(41)+0.2*x(51)-0.2*x(44);2*x(6)+2+x(35)+1.92*x(6)+1.92*x(4)+1.28*x(5)+1.28*x(7)-0.8*x(40)+0.2*x(50)-0.2*x(43);2*x(7)+2+x(35)+1.92*x(7)+1.92*x(5)+1.28*x(4)+1.28*x(6)-0.8*x(41)+0.2*x(51)-0.2*x(44);x(8)+3+x(36)+3.2*x(8)+4-0.8*x(42)+0.2*x(52)-0.2*x(45);2*x(9)+2+0.8*x(9)+0.8*x(10)-0.4*x(42)+0.6*x(52)-0.6*x(45);x(10)+2.5+0.8*x(9)+0.8*x(10)-0.4*x(42)+0.6*x(52)-0.6*x(45);2*x(11)+3+2*x(11)+2*x(12)+2*x(13)+0.4-0.8*x(46)-x(47)+x(37)-x(53);4*x(12)+1+2*x(11)+2*x(12)+2*x(13)+0.4-0.8*x(46)-x(48)+x(38)-x(53);x(13)+5+2*x(11)+2*x(12)+2*x(13)+0.4-0.8*x(46)-x(49)+x(39)-x(53);4*x(14)+2.5+2.56*x(14)+0.64*x(15)-0.8*x(36)+0.2*x(49)-0.2*x(39);3*x(15)+2.5+2.56*x(15)+0.64*x(14)-0.8*x(36)+0.2*x(49)-0.2*x(39);3*x(16)+3+3.2*x(16)+1.2-0.8*x(34)+0.2*x(47)-0.2*x(37);4*x(17)+2+3.2*x(17)+1.2-0.8*x(35)+0.2*x(48)-0.2*x(38);1.6*x(18)+1.6-0.4*x(34)+3*x(18)+2.5+x(46)+0.6*x(47)-0.6*x(37);1.6*x(19)+1.6-0.4*x(35)+5*x(19)+2+x(46)+0.6*x(48)-0.6*x(38);1.28*x(20)+1.2-0.4*x(36)+6*x(20)+1+x(46)+0.6*x(49)-0.6*x(39);2*x(21)+3+2*x(21)+2*x(22)+2*x(23)+0.8-0.6*x(46)-x(50)+x(43)-x(53);x(22)+2.5+2*x(21)+2*x(22)+2*x(23)+0.8-0.6*x(46)-x(51)+x(44)-x(53);3*x(23)+2+2*x(21)+2*x(22)+2*x(23)+0.8-0.6*x(46)-x(52)+x(45)-x(53);x(24)+3+x(40)+x(24)+5-x(30);3*x(25)+2.5+x(41)+x(25)+5-x(30);5*x(26)+1+x(42)+x(26)+5-x(31);x(27)+3+x(40)+x(27)+5-x(32);3*x(28)+2.5+x(41)+x(28)+5-x(32);5*x(29)+1+x(42)+x(29)+5-x(33);x(24)+x(25)+2*x(30)+1.5*x(32)-1000;x(27)+x(28)+2*x(32)+1.5*x(30)-1000;x(26)+2.5*x(31)+2*x(33)-1000;x(29)+2.5*x(33)+2*x(31)-1000;0.6*x(1)+0.8*x(16)+0.4*x(18)-x(4)-x(5);0.6*x(2)+0.8*x(17)+0.4*x(19)-x(6)-x(7);0.6*x(3)+0.8*x(14)+0.8*x(15)+0.4*x(29)-x(8);0.4*x(1)+0.2*x(16)+0.6*x(18)-x(11);0.4*x(2)+0.2*x(17)+0.6*x(19)-x(12);0.4*x(3)+0.2*x(14)+0.2*x(15)+0.6*x(20)-x(13);0.8*x(4)+0.8*x(6)-x(24)-x(27);0.8*x(5)+0.8*x(7)-x(25)-x(28);0.4*x(9)+0.4*x(10)+0.8*x(8)-x(26)-x(29);0.2*x(4)+0.2*x(6)-x(21);0.2*x(5)+0.2*x(7)-x(22);0.6*x(9)+0.6*x(10)+0.2*x(8)-x(23);0.8*x(11)+0.8*x(12)+0.8*x(13)+0.6*x(21)+0.6*x(22)+0.6*x(23)-x(18)-x(19)-x(20)];

Fx=F(x1);%把x1带入F求得函数值F(x1).
nv=zeros(46,1);
for i=1:46
      nv(i)=x1(i)+Fx(i)-sqrt(x1(i)^2+Fx(i)^2);
end
tol=norm(nv);
   if tol<=eps
      x=x1;
      break;
   else
      x0=x1;
      e=e*c;%如果不满足退出，则把x0更新为x1,且把e变小。
   end
end
x=x0;

另附上：
function [c,ceq]=mycon_fish(x,e)
c=(0.6*x(18)+0.2*x(16)+0.4*x(1)-x(11))*x(47)+(0.6*x(19)+0.2*x(17)+0.4*x(2)-x(12))*x(48)+(0.6*x(20)+0.2*x(14)+0.2*x(15)+0.4*x(3)-x(13))*x(49)+(0.2*x(40)+0.2*x(6)-x(21))*x(50)+(0.2*x(49)+0.2*x(7)-x(22))*x(51)+(0.2*x(8)+0.6*x(9)+0.6*x(10)-x(23))*x(52)+(x(11)+x(12)+x(13)+x(21)+x(22)+x(23))*x(53)-10000;
x=[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9),x(10),x(11),x(12),x(13),x(14),x(15),x(16),x(17),x(18),x(19),x(20),x(21),x(22),x(23),x(24),x(25),x(26),x(27),x(28),x(29),x(30),x(31),x(32),x(33),x(34),x(35),x(36),x(37),x(38),x(39),x(40),x(41),x(42),x(43),x(44),x(45),x(46),x(47),x(48),x(49),x(50),x(51),x(52),x(53)];
F=[3.8*x(1)+1.36*x(2)+3.2-0.6*x(34)+0.4*x(47)-0.4*x(37);3.8*x(2)+1.36*x(1)+3.2-0.6*x(35)+0.4*x(48)-0.4*x(38);3.08*x(3)+6-0.6*x(36)+0.4*x(49)-0.4*x(39);x(4)+3.5+x(34)+1.92*x(4)+1.92*x(6)+1.28*x(5)+1.28*x(7)-0.8*x(40)+0.2*x(50)-0.2*x(43);x(5)+3.5+x(34)+1.92*x(5)+1.92*x(7)+1.28*x(4)+1.28*x(6)-0.8*x(41)+0.2*x(51)-0.2*x(44);2*x(6)+2+x(35)+1.92*x(6)+1.92*x(4)+1.28*x(5)+1.28*x(7)-0.8*x(40)+0.2*x(50)-0.2*x(43);2*x(7)+2+x(35)+1.92*x(7)+1.92*x(5)+1.28*x(4)+1.28*x(6)-0.8*x(41)+0.2*x(51)-0.2*x(44);x(8)+3+x(36)+3.2*x(8)+4-0.8*x(42)+0.2*x(52)-0.2*x(45);2*x(9)+2+0.8*x(9)+0.8*x(10)-0.4*x(42)+0.6*x(52)-0.6*x(45);x(10)+2.5+0.8*x(9)+0.8*x(10)-0.4*x(42)+0.6*x(52)-0.6*x(45);2*x(11)+3+2*x(11)+2*x(12)+2*x(13)+0.4-0.8*x(46)-x(47)+x(37)-x(53);4*x(12)+1+2*x(11)+2*x(12)+2*x(13)+0.4-0.8*x(46)-x(48)+x(38)-x(53);x(13)+5+2*x(11)+2*x(12)+2*x(13)+0.4-0.8*x(46)-x(49)+x(39)-x(53);4*x(14)+2.5+2.56*x(14)+0.64*x(15)-0.8*x(36)+0.2*x(49)-0.2*x(39);3*x(15)+2.5+2.56*x(15)+0.64*x(14)-0.8*x(36)+0.2*x(49)-0.2*x(39);3*x(16)+3+3.2*x(16)+1.2-0.8*x(34)+0.2*x(47)-0.2*x(37);4*x(17)+2+3.2*x(17)+1.2-0.8*x(35)+0.2*x(48)-0.2*x(38);1.6*x(18)+1.6-0.4*x(34)+3*x(18)+2.5+x(46)+0.6*x(47)-0.6*x(37);1.6*x(19)+1.6-0.4*x(35)+5*x(19)+2+x(46)+0.6*x(48)-0.6*x(38);1.28*x(20)+1.2-0.4*x(36)+6*x(20)+1+x(46)+0.6*x(49)-0.6*x(39);2*x(21)+3+2*x(21)+2*x(22)+2*x(23)+0.8-0.6*x(46)-x(50)+x(43)-x(53);x(22)+2.5+2*x(21)+2*x(22)+2*x(23)+0.8-0.6*x(46)-x(51)+x(44)-x(53);3*x(23)+2+2*x(21)+2*x(22)+2*x(23)+0.8-0.6*x(46)-x(52)+x(45)-x(53);x(24)+3+x(40)+x(24)+5-x(30);3*x(25)+2.5+x(41)+x(25)+5-x(30);5*x(26)+1+x(42)+x(26)+5-x(31);x(27)+3+x(40)+x(27)+5-x(32);3*x(28)+2.5+x(41)+x(28)+5-x(32);5*x(29)+1+x(42)+x(29)+5-x(33);x(24)+x(25)+2*x(30)+1.5*x(32)-1000;x(27)+x(28)+2*x(32)+1.5*x(30)-1000;x(26)+2.5*x(31)+2*x(33)-1000;x(29)+2.5*x(33)+2*x(31)-1000;0.6*x(1)+0.8*x(16)+0.4*x(18)-x(4)-x(5);0.6*x(2)+0.8*x(17)+0.4*x(19)-x(6)-x(7);0.6*x(3)+0.8*x(14)+0.8*x(15)+0.4*x(29)-x(8);0.4*x(1)+0.2*x(16)+0.6*x(18)-x(11);0.4*x(2)+0.2*x(17)+0.6*x(19)-x(12);0.4*x(3)+0.2*x(14)+0.2*x(15)+0.6*x(20)-x(13);0.8*x(4)+0.8*x(6)-x(24)-x(27);0.8*x(5)+0.8*x(7)-x(25)-x(28);0.4*x(9)+0.4*x(10)+0.8*x(8)-x(26)-x(29);0.2*x(4)+0.2*x(6)-x(21);0.2*x(5)+0.2*x(7)-x(22);0.6*x(9)+0.6*x(10)+0.2*x(8)-x(23);0.8*x(11)+0.8*x(12)+0.8*x(13)+0.6*x(21)+0.6*x(22)+0.6*x(23)-x(18)-x(19)-x(20)];
ceq=zeros(46,1);
for i=1:46
      ceq(i)=x(i)+F(i)-sqrt(x(i)^2+F(i)^2+e);
end

回复此楼

» 猜你喜欢

求调剂已经有3人回复
265求调剂已经有4人回复
085700资源与环境308求调剂已经有6人回复
一志愿吉林大学材料学硕321求调剂已经有12人回复
286分人工智能专业请求调剂愿意跨考！已经有3人回复
329求调剂已经有5人回复
申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复已经有4人回复
材料学硕318求调剂已经有5人回复
一志愿中国海洋大学，生物学，301分，求调剂已经有6人回复
081700化工学硕调剂已经有3人回复

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

请问大家哪个分析类杂志接受综述容易点啊，明年要毕业了，急用！已经有4人回复
求助想发有机氯农药分析类的文章，毕业急用，不知道核心期刊既快且接受容易已经有12人回复

1楼 2014-06-24 15:49:36

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

baobiao007

木虫 (职业作家)

中国特色

应助: 201 (大学生)
金币: 6482.7
散金: 557
红花: 40
帖子: 3050
在线: 1009.9小时
虫号: 505962
注册: 2008-02-18
专业: 应用地球物理学

.....................................想不吓到？不可能

赞一下

回复此楼

我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦

2楼2014-06-24 16:55:43

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

qingyi12075277

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 2.6
散金: 1359
帖子: 188
在线: 56.7小时
虫号: 672736
注册: 2008-12-11
专业: 微生物生态学

被吓到了

回复此楼

3楼2014-06-25 00:01:45

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

fatcat916

至尊木虫 (著名写手)

金肥猫

应助: 1 (幼儿园)
金币: 15360.2
帖子: 1168
在线: 375小时
虫号: 1852958
注册: 2012-06-09
性别: GG
专业: 计算机软件

我Kao，还有这么写代码的？截屏供观摩。

赞一下

回复此楼

4楼2014-06-26 01:29:36

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

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

这是完全不理解Matlab的语法; 应该从头学习矩阵和向量的用法

赞一下

回复此楼

ส็็็็็็็็็็็็็็็็็็็็

5楼2014-06-26 07:32:02

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

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

对大多数熟悉Matlab的人来说, 解决你原始问题, 比调试你的代码要容易

赞一下

回复此楼

ส็็็็็็็็็็็็็็็็็็็็

6楼2014-06-26 07:33:11

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

bingbing5380

铁虫 (正式写手)

引用回帖:

6楼: Originally posted by cooooldog at 2014-06-26 07:33:11
对大多数熟悉Matlab的人来说, 解决你原始问题, 比调试你的代码要容易

7楼2014-07-01 07:48:25

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

bingbing5380

铁虫 (正式写手)

应助: 1 (幼儿园)
金币: 72.7
散金: 450
红花: 3
帖子: 331
在线: 144.1小时
虫号: 640275
注册: 2008-10-29
专业: 宏观经济管理与战略

引用回帖:

5楼: Originally posted by cooooldog at 2014-06-26 07:32:02

这是完全不理解Matlab的语法; 应该从头学习矩阵和向量的用法

初级新手啊，别打击俺啊

，问题自己已解决。

赞一下

回复此楼

8楼2014-07-01 07:49:13

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

lyf1850

祝福祝福

9楼2014-07-01 08:29:54

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

相关版块跳转我要订阅楼主 bingbing5380 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 329求调剂 +5	想上学吖吖 2026-03-19	5/250	2026-03-19 23:56 by 23Postgrad
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +3	sbdksD 2026-03-19	3/150	2026-03-19 23:21 by fmesaito
[考研] 294求调剂材料与化工专硕 +14	陌の森林 2026-03-18	14/700	2026-03-19 22:38 by 学员8dgXkO
[考研] 梁成伟老师课题组欢迎你的加入 +9	一鸭鸭哟 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4	脱颖而出 2026-03-16	12/600	2026-03-19 16:17 by 脱颖而出
[考研] 一志愿北京化工大学0703化学318分，有科研经历，求调剂 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 尽舜尧1
[考研] 328求调剂，英语六级551，有科研经历 +4	生物工程调剂 2026-03-16	12/600	2026-03-19 11:10 by 生物工程调剂
[考研] 一志愿华中科技大学，080502，354分求调剂 +4	守候夕阳CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 【同济软件】软件（085405）考研求调剂 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏击518
[考研] 277调剂 +5	自由煎饼果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[硕博家园] 湖北工业大学生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 材料与化工专硕调剂 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太监王公公
[考研] 11408 一志愿西电，277分求调剂 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 机械专硕325，寻找调剂院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0854控制工程 359求调剂可跨专业 +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 321求调剂 +5	大米饭！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 297求调剂 +4	学海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 热情沙漠
[考研] 招收0805（材料）调剂 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 290求调剂 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike