Matlab二次规划编程-经管之家官网！

function [re1,re2]=minvar(x,y)
% calculate the optimal coefficient by minimize the variance between x and
% c*y. That is var(x-c*y).
% Input parameters:
% x: the raw data
% y: the benchmark for optimization
% Output parameters:
% re1: the R^2 (first column)
% the parameters c (rest of every row)
% re2: the ESS
% x should be a array or a matrix
% y should be a matrix,the form is either(1,b1,b2...),or(b1,b2...)
% x and y should have the same dimension
[m,n]=size(y);
[a,b]=size(x);
lb = zeros(1,n); % Set lower bounds
if y(:,1) == ones(m,1)% if b is (1,b1,b2...), the corresponding lower bound is empty
lb(1,1)=-inf;
end
ub = [ ]; % No upper bounds
Aeq= ones(1,n); % Set equality constrain
if y(:,1) == ones(m,1)% if b is (1,b1,b2...), the Aeq is (0,1,1...)
Aeq(1,1)=0;
end
beq=1;
options = optimset('LargeScale','off');
mm1 = zeros(b,n+1);
mm2 = zeros(a,b);
for i =1:b
k = x(:,i);
[v,fval] =lsqlin(y,k,[ ],[ ],Aeq,beq,lb,ub);
mm1(i,2:n+1) = v'; % calculate the optimal coefficient c
tssk=sum((k-mean(k)).^2);
essk=sum((k-y*v).^2);
mm1(i,1)=1-essk/tssk; % calculate the R^2 value
mm2(:,i)=(k-y*v).^2; % calculate the ESS
end
re1 = mm1;
re2 = mm2;
为什么运行后出现： Error: Function definitions are not permitted at the prompt or in scripts.