slvblk.m

演示matlab曲线拟和与插直的基本方法

% SLVBLK   僽儘僢僋懳妏慄宍僔僗僥儉偺夝朄
%
% SLVBLK(BLOKMAT,B) 偼丄傎傏僽儘僢僋懳妏側宆偺僗僾儔僀儞偱(偨偲偊偽丄
% SPCOL 偱弌椡偝傟傞傛偆側) BLOKMAT 偵奿擺偝傟偨峴楍 A 傪傕偮慄宍僔僗僥儉 
% A*X=B 偺夝傪(傕偟偁傟偽)弌椡偟傑偡丅
%
% 僔僗僥儉偑夁寛掕(偡側傢偪曽掱幃偺悢偑枹抦悢傛傝傕懡偄)偺応崌丄
% 嵟彫擇忔夝偑弌椡偝傟傑偡丅 
%
% SLVBLK(BLOKMAT,B,W) 偼丄廳傒晅偒 l_2 偺榓傪嵟彫壔偡傞儀僋僩儖 X 傪
% 弌椡偟傑偡丅
%
%      sum_j W(j)*( (A*X-B)(j) )^2 .
%
% 偙傟偼丄僔僗僥儉偑夁寛掕偺応崌丄桳岠偱偡丅
% W 偵懳偡傞僨僼僅儖僩偼丄楍 [1,1,1,...]偱偡丅
%
% 椺戣:
% 偮偓偺僗僥乕僩儊儞僩偱偼丄偁傞僲僀僘偺擖偭偨僨乕僞傪惗惉偟傑偡丅
% 偦傟偐傜丄2偮偺摍娫妘偵暲傫偩愡揰傪傕偮僉儏乕價僢僋僗僾儔僀儞偵傛傞
% 僨乕僞偵懳偟偰丄SLVBLK 傪梡偄偰暋崌戜宍懃偺廳傒偵傛偭偰廳傒晅偗傜傟偨
% 嵟彫2忔嬤帡傪寛掕偟丄寢壥傪僾儘僢僩偟傑偡丅
%   
%      rand('seed',27);
%      x = [0,sort(rand(1,31)),1]*(2*pi);
%      y = sin(x)+(rand(1,33)-.5)/10;
%      k = 4; knots = augknt(linspace(x(1),x(end),3),k);
%      dx = diff(x); w = ([dx 0] + [0 dx])/2;
%      sp = spmak(knots,slvblk(spcol(knots,k,x,'slvblk','noderiv'),y.',w).');
%      fnplt(sp,2); hold on, plot(x,y,'ok'), hold off
%
% 嶲峫 SPCOL, SPAPS, SPAPI, SPAP2.


%   Copyright 1987-2004 C. de Boor and The MathWorks, Inc.