悬赏50金币求解如何用matlab编写如下程序

邹华伟

2186

收藏 2016-04-13

悬赏 1 个论坛币已解决

附件: 您需要登录才可以下载或查看附件。没有帐号？我要注册

最佳答案

Weierstrass. 查看完整内容

函数定义中最后以 end 结束最后一行加 end 这个是每个MATLAB函数定义必须的。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

Weierstrass.

2016-4-13 15:43:41

函数定义中最后以 end 结束
最后一行加 end 这个是每个MATLAB函数定义必须的。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

成联正

2016-4-13 15:59:55

MATLAB Optimization Techniques或者搜索这方面的书籍，很多的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-13 16:14:30

成联正发表于 2016-4-13 15:59
MATLAB Optimization Techniques或者搜索这方面的书籍，很多的

新手找不到相应的程序请问可以麻烦您帮我写一下吗或者说一下在哪里有类似的情况

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

seaknight

2016-4-13 18:53:52

首先建立一个M文件Fun75.m
Function y = fun75(X)
Y=-350*(1-exp-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7)))-200-400+x(1)+x(2)+x(3)+x(4)+x(5)+x(6)+x(7)
储存为fun75.m文件

A=[1,0,0,0.65,0.16,0,0.16;0,-1,0,-0.65,0,-0.34,-0.16;0,0,-1,0,-0.16,-0.34,-0.16]
B=[4.47;-3.91;-4.6]
Aeq=[];Beq=[];Lb=[0;0;0;0;0;0;0]
[x,g]=fmincon(‘myfun75’,[],A,B,Aeq,Beq,Lb)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

seaknight

2016-4-13 18:55:38

差不多是这样的，主要学过matlab，再去网上找相同的问题解答，就可以照葫芦画瓢地便出代码。
记得给金币喔！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

seaknight

2016-4-13 18:56:22

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-13 19:14:19

seaknight 发表于 2016-4-13 18:56
首先建立一个M文件Fun75.m
Function y = fun75(X)
Y=-350*(1-exp-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))) ...

为什么出现了这个我错在哪里了

附件: 您需要登录才可以下载或查看附件。没有帐号？我要注册

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

gydgzd

2016-4-13 19:31:59

你好像用了中文的引号

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-13 20:06:47

gydgzd 发表于 2016-4-13 19:31
你好像用了中文的引号

恩的确是可是我改了之后提示我缺少初始值我该怎么办

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

joylk

2016-4-13 20:09:24

楼主可以参考这篇文章，有详细的举例说明。建议自己编写一遍。这样才能快速掌握。
http://blog.sina.com.cn/s/blog_62e4f3830102v11n.html

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

seaknight

2016-4-13 21:18:34

[] 就是初始变量，因为你的题目中没有设置初始变量，所以这是用[]来代替

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

seaknight

2016-4-13 21:20:49

如果不行的话，你可以加入一个X0=[1;1;1;1;1;1;1],然后就把A之前的[]，变成X0

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

shuxue_ch

2016-4-13 21:32:56

1.线性规划，首选Lingo啊。
2. MATLAB教程，照猫画虎

附件: 您需要登录才可以下载或查看附件。没有帐号？我要注册

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-13 21:43:51

seaknight 发表于 2016-4-13 21:18
[] 就是初始变量，因为你的题目中没有设置初始变量，所以这是用[]来代替

多谢回复如果加上一个限制条件x2=0 能不能再把详细的步骤在写一次

给你100金币最好麻烦你先跑出来我看看新手不太懂改起来好麻烦多谢啊

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Weierstrass.

2016-4-13 22:05:17

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Weierstrass.

2016-4-13 22:22:43

// 存储为myfun.m 切记，切记
function y= myfun(x)
y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))))-200-400+x(1)+x(2)+x(3)+x(4)+x(5)+x(6)+x(7);
end
A=[1, 0, 0, 0.65, 0.16, 0, 0.16;
0, -1, 0, -0.65,0,-0.34,-0.16;
0, 0, -1, 0, -0.16, -0.34, -0.16];
B=[4.47;-3.91;-4.6];
x0=zeros(1,7);
Aeq=[];Beq=[];Lb=[0;0;0;0;0;0;0];
[x,fval]=fmincon(@myfun,x0,A,B,Aeq,Beq,Lb);
xmin=x
fmax=-fval
结果为------------xmin = 0.5600 0.0000 4.6000 6.0154 0.0000 0.0000 0.0000-----对应于x1,x2...x7
-------------fmax = 934.8181

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-14 00:35:41

Weierstrass. 发表于 2016-4-13 22:22
// 存储为myfun.m 切记，切记
function y= myfun(x)
y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x ...

嗯你是对的，你重发一张随意图片设置50金币我买吧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-14 00:36:32

seaknight 发表于 2016-4-13 18:53
首先建立一个M文件Fun75.m
Function y = fun75(X)
Y=-350*(1-exp-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))) ...

感谢热情的网友，你重发任意一张图片设置五十金币我购买

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Weierstrass.

2016-4-14 00:39:46

附件: 您需要登录才可以下载或查看附件。没有帐号？我要注册

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-14 21:04:56

Weierstrass. 发表于 2016-4-14 00:39
// 存储为myfun.m 切记，切记
function y= myfun(x)
y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x ...

你好你的答案是正确的十分感谢我想问一下如果将目标函数改为y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))))-200-400+x(1)+c*(x(2)+x(3)+x(4)+x(5)+x(6)+x(7));c为不定常数，取值为（0,100）那么该如何编写程序在此感谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-14 21:05:43

seaknight 发表于 2016-4-13 21:20
如果不行的话，你可以加入一个X0=[1;1;1;1;1;1;1],然后就把A之前的[]，变成X0

我想问一下如果将目标函数改为y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))))-200-400+x(1)+c*(x(2)+x(3)+x(4)+x(5)+x(6)+x(7));c为不定常数，取值为（0,100）那么该如何编写程序在此感谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Weierstrass.

2016-4-14 23:23:11

在函数定义中加入循环
function.....
for c=(0:100)
.....其他不变。最后x0 初始值很重要的，确定只限制x(2)=0?

function y= myfun(x)
for c=(1:100)
y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))))-200-400+x(1)+c*(x(2)+x(3)+x(4)+x(5)+x(6)+x(7));
end

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-15 10:19:25

Weierstrass. 发表于 2016-4-14 23:23
在函数定义中加入循环
function.....
for c=(0:100)

??? Error using ==> fmincon
FMINCON cannot continue because user supplied objective function failed with the following error:
Error: File: F:\MATLAB7\myfun.m Line: 3 Column: 102
This statement is incomplete.
这是咋的了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Weierstrass.

2016-4-15 17:03:22

函数定义不完整，具体看myfun.m的第3行。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-15 17:10:29

Weierstrass. 发表于 2016-4-15 17:03
函数定义不完整，具体看myfun.m的第3行。

function y=myfun(x)
for c=(0:100)
y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))))-200-400+c*(x(1)+x(2)+x(3)+x(4)+x(5)+x(6)+x(7));
错了吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-15 17:33:21

Weierstrass. 发表于 2016-4-15 17:03
函数定义不完整，具体看myfun.m的第3行。

我的具体程序是这样的  你看看哪里出错了？
function y=myfun(x)
for c=(0:100)
y=-350*(1-exp(-(x(1)+0.65*x(4)+0.16*x(5)+0.16*x(7))))-200-400+c*(x(1)+x(2)+x(3)+x(4)+x(5)+x(6)+x(7));

A=[1, 0, 0, 0.65, 0.16, 0, 0.16;
            0, -1, 0, -0.65,0,-0.34,-0.16;
         0, 0, -1, 0, -0.16, -0.34, -0.16];
      B=[4.47;-3.91;-4.6];
Aeq=[]; Beq=[];
Lb=[0;0;0;0;0;0;0];
[x,fval]=fmincon(@myfun,  0.000000001*ones(1,7)', A,B,Aeq,Beq,Lb);
xmin=x
fmax=-fval
提示我??? Error using ==> fmincon
FMINCON cannot continue because user supplied objective function failed with the following error:
Error: File: F:\MATLAB7\myfun2.m Line: 3 Column: 102
This statement is incomplete.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

邹华伟

2016-4-15 18:37:33

Weierstrass. 发表于 2016-4-15 17:03
函数定义不完整，具体看myfun.m的第3行。

十分感谢！的确是少end。我想研究的是c取不同值时候原函数的最优化问题：就是说当c=0，c=0.5，c=1....c=100，分别求出函数的最大值然后画出c的取值与最大值的图像以及最大值时候对应的x1，x2....x7的取值图像

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Weierstrass.

2016-4-15 23:06:29

含参数的非线性规划，感觉c 如果定义为 x8可能更好一点。或者可以重新考虑模型。
附件无法上传，（仅从数学角度看发现结果不太好，请讨论

）

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群