conpred.m

用VC++编写的遗传算法源程序

function m = conpred(b,w,p,q,lv)

[mq,nq] = size(q);
[mw,nw] = size(w);
if nw ~= lv
  if lv > nw
    s = sprintf('Original model has a maximum of %g LVs',nw);
    disp('  '), disp(s)
	s = sprintf('Calculating vectors for %g LVs only',nw);
	disp(s), disp('  ')
	lv = nw;
  else
    w = w(:,1:lv);
	q = q(:,1:lv);
	p = p(:,1:lv);
	b = b(:,1:lv);
  end
end
m = zeros(mq*lv,mw);
if mq == 1
  m = (w*inv(p'*w)*diag(b))';
else
  mp = (w*inv(p'*w)*diag(b))';
  for i = 1:lv
    mpp = mp(i,:);
    m((i-1)*mq+1:i*mq,:) = diag(q(:,i))*mpp(ones(mq,1),:);
  end
end