| 查看: 1456 | 回复: 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 ] |
回复此楼» 猜你喜欢
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】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人回复
3楼2010-04-22 08:04:41
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
4楼2010-04-24 09:18:10
5楼2010-04-24 18:08:16
20