登录论坛

查看完整版本 : [MATLAB混合编程] 大侠帮看看这段matlab程序错在哪里


zhaoguoxi
2010-04-22, 09:04
一直出现提示
??? Error using ==> fmincon
FMINCON cannot continue because user supplied objective function
failed with the following error:

Index exceeds matrix dimensions.

我的程序如下:
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')

function fun=obfun(q)
clear;
n=6;s=zeros(n+1,1);g=zeros(n+1,1);h=zeros(41,1);
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);

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];