function [w,wbias] = solveout2(g,t)
%隱層到輸出層的權(quán)值
% g為np×m的隱含層的輸出矩陣
% t為np×no的期望輸出矩陣
% w為m×no的權(quán)值矩陣
% wbias為1×no的偏值向量
if nargin < 2,error('Not enough arguments.'),end

if nargout <= 1
  w= inv(g'*g)*g'*t;
else
  [np,m] = size(g);
  y = [g'; ones(1,np)]';
  x=inv(y'*y)*y'*t;
  w=x(1:m,:);
  wbias=x(m+1,:);
end