home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 5.ddi / SPLINES.DI$ / AVEKNT.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  876 b   |  31 lines
  1. function tstar=aveknt(t,k)
  2. % AVEKNT    Provide knot averages.
  3. %
  4. %        tstar=aveknt(t,k)
  5. %
  6. %  returns the averages of successive  k-1  knots, i.e., the points
  7. %
  8. %           tstar(i) = (t_{i+1} + ... + t_{i+k-1})/(k-1)
  9. %
  10. %  recommended as good interpolation point choices when interpolating from
  11. %
  12. %           S_{k,t} .
  13.  
  14. % C. de Boor / latest change: May 25, 1989
  15. % C. de Boor / latest change: March 22, 1991 (program around a MATLAB
  16. %                                                            discontinuity)
  17. % Copyright (c) 1990-92 by Carl de Boor and The MathWorks, Inc.
  18.  
  19. t=t(:);lt=length(t);
  20. if (k<2), tstar=[], error('second argument must be at least 2');
  21. elseif (k==2), tstar=t(2:lt-1)';
  22. else,
  23.    n=lt-k;
  24.    t*ones(1,k-1);
  25.    junk=ans(:);
  26.    try=zeros(n+k+1,k-1);
  27.    try(:)=[junk;zeros(k-1,1)];
  28.    sum(try')/(k-1);
  29.    tstar=ans(1+[1:n]);
  30. end
  31.