| 查看: 1475 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】matlab运行结果如下，不知何故，请高人指路【已完结】 已有3人参与
|
|||
|
命令行： q0=[0； 0 ；0.5； 0.5 ；0； 0； 2]; aeq=[1 1 1 1 1 1 0 1 2 3 4 5 6 -0.89298]; beq=[1 0]; a=[];b=[]; vlb=[0 0 0 0 0 0 2];vub=[1 1 1 1 1 1 10]; [q,fun]=fmincon('objectfun',a,b,aeq,beq,vlb,vub,'ncon') objectfun文件： function fun=objectfun(q) n=6;s=zeros(1,n+1);g=zeros(1,n+1);h=zeros(1,41);q=[]; x=linspace(0,2,41); for k=1:41 for j=1:n for i=0:j-1 s(j)=s(j)+(q(7)*x(k))^i*exp(-q(7)*x(k))/factorial(i); end g(j)=q(j)*s(j); end h(k)=abs(exp(-x(k)^2)-sum(g)); end fun=sum(h); ncon文件： function [gc,ceq]=ncon(q) gc=[]; ceq=[2*q(1)+6*q(2)+12*q(3)+20*q(4)+30*q(5)+42*q(6)-0.902705*q(7)^2]; 运行结果如下，不知何故，请高人指路： ??? Error using ==> fmincon FMINCON cannot continue because user supplied objective function failed with the following error: Index exceeds matrix dimensions. [ Last edited by 余泽成 on 2010-4-26 at 09:00 ] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】MATLAB编程
已经有2人回复
- 【求助】matlab求助
已经有7人回复
- 【求助】Matlab
已经有4人回复
- 【求助】matlab
已经有2人回复
- 【求助】matlab资源
已经有8人回复
- 【求助】MATLAB错误
已经有6人回复
- 【求助】matlab程序
已经有1人回复
- 【求助】matlab编程求助
已经有2人回复
- 【求助】matlab
已经有0人回复
- 【求助】matlab 求助
已经有2人回复
- 【求助】matlab 求助
已经有1人回复
- 【求助】matlab求助，
已经有5人回复
- 【求助】matlab视频教程
已经有8人回复
- 【求助】Matlab程序求助
已经有2人回复
6楼2010-04-25 21:34:40
zzy870720z
荣誉版主 (文坛精英)
- 程序强帖: 1
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2010-04-22 00:00:42
3楼2010-04-22 08:04:41
4楼2010-04-24 09:18:10
30