www.gusucode.com > mbcmodels 工具箱 matlab 源码程序 > mbcmodels/@localpspline/symmetric.m

    function [x,y]= symmetric(ps,x,y)
% localpspline/SYMMETRIC

%  Copyright 2000-2008 The MathWorks, Inc. and Ford Global Technologies, Inc.



%[x,y]= NewSym(ps,x,y);
%return

m = (x > ps.knot);
nm= sum(m);
if nm<ps.order(1)
   % rhs
   x= [x ; ps.knot-(x(~m)-ps.knot)];
   if nargin==3
      y= [y ; y(~m)];
   end
elseif length(x)-nm<ps.order(2)
   % lhs
   x= [x ; ps.knot-(x(m)-ps.knot)];
   if nargin==3
      y= [y ; y(m)];
   end
end