spline toolbox

function splinetool(x,y) %SPLINETOOL Experiment with some spline approximation methods (GUI). % % SPLINETOOL prompts you for some data sites and data values to fit and % lets you do this in various ways. % One of the choices is to provide your own data, in which case both % the sites and values can be expressions that evaluate to vectors required % to be of the same length. % You can also specify the values by providing the name of a function whose % values at the sites are to be used as data values. % You can also provide the name of a function (like 'titanium') whose % two-argument output provides numerical values for sites and values. % In any case, there must be at least two distinct sites. % The other choices provide specific sample sites and values to illustrate % various aspects of the GUI. % % SPLINETOOL(X,Y) uses the input data sites X and data values Y, and these % must be numerical vectors of the same length > 1. The data sites need not % be distinct nor ordered, but there must be at least two distinct sites. % % Example: % x = linspace(1,pi,101); y = cos(x)+(rand(size(x))-.5)/10; % splinetool(x,y) % might start an interesting experimentation with noisy data. % % See also CSAPI, CSAPS, SPAPS, SPAP2, SPAPI. % Copyright 1987-2010 The MathWorks, Inc. % $Revision: 1.1.6.3.2.1 $ $Date: 2010/06/28 14:16:57 $ testing = 0; v = ver('curvefit'); [relYear,~,~] = datevec(v.Date); abouttext = { sprintf( '%s %s', v.Name, v.Version ) sprintf( 'Copyright 1987-%d The MathWorks, Inc.', relYear ) }; if nargin<1 % import data action = 'start'; else if ischar(x) action = x; hand = get(gcbf,'Userdata'); if isstruct(hand)&&isfield(hand,'currentline') hc = hand.currentline; % make available, for use with changes end else % check and sort the given data [x,y] = chckxy(x,y); xname = 'x'; yname = 'y'; isf = 0; action = 'start'; end end switch action case 'add_item' % add an item switch get(hand.data_etc,'Value') case 1 [xx,yy] = ask_for_add(hand); x = get(hand.dataline,'Xdata'); if xx<x(1)||xx>x(end) warndlg(['The Spline Tool will let you only', ... ' add interior data sites.'], ... 'Attempted change of interval of approximation','modal') elseif get(hand.method,'Value')==4 && any(x==xx) warndlg(['Since SPAPI would interpret repeated data sites',... ' as a request for osculatory interpolation, the Spline Tool',... ' won''t let you create repeated data sites now.'], ... 'Repeated data sites', 'modal'); else set(hand.undo,'Enable','on') insert_and_mark(x,get(hand.dataline,'Ydata'),xx,yy,hand) undoud = get(hand.undo,'Userdata'); undoud.lineud = get(hand.nameline,'Userdata'); set(hand.undo,'Userdata',undoud); get_approx(hand) end if isequal(get(hand.undo,'Enable'),'off') V = get(hand.params(1,3),'Value'); set(hand.params(2,3),'String',x(V)) if isequal(get(hand.params(2,5),'Visible'),'on') y = get(hand.dataline,'Ydata'); set(hand.params(2,5),'String',y(V)) end end case 2 switch get(hand.method,'Value') case 3 xx = ask_for_add(hand); lineud = get(hand.nameline,'Userdata'); [lineud.knots,i] = sort([lineud.knots,xx]); placed = find(i==length(lineud.knots)); if placed==1||placed==length(lineud.knots) warndlg('The Spline Tool will only let you add interior knots.',... 'The end knots are forever','modal') lineud.knots(placed)=[]; set(hand.nameline,'Userdata',lineud) set(hand.params(2,4),'String', ... lineud.knots(get(hand.params(1,4),'Value'))) else set(hand.params(1,4),'String',lineud.knots(:),'Value', placed) lineud = update(hand,'knots',lineud); lineud.kc(end+1) = ... addchange(hc,['knots = sort([knots,',num2str(xx,15),']);',... ' %% you added a knot']); set(hand.nameline,'Userdata',lineud) % update knots display set(hand.buttonm,'Visible','on') set(hand.params(4,4), 'String',... (lineud.knots(end)-lineud.knots(1))/(10*length(lineud.knots))) markb(hand), set_bdisplay(lineud.knots,hand) get_approx(hand) end end %switch get(hand.method,'Value') case 3 end %switch get(hand.data_etc,'Value') case 'align' % make sure the left and right data list are aligned % at present, this serves no function, but sits here in case clicking % anywhere on a list box activates a callback tag = get(gcbo,'Tag'); if tag(5)=='l', j=3; else j=5; end set(hand.params(1,8-j),'Value',get(hand.params(1,j),'Value')) case 'axesclick' % action when the axes is clicked % protect against clicks on Print Figure window: if isempty(hand), return, end % only react here to LEFT clicks tmp = get(hand.figure,'SelectionType'); if tmp(1)=='n' % get the location of the latest click in data coordinates: clickp = get(hand.Axes(1),'CurrentPoint'); switch get(hand.data_etc,'Value') case {1,3} % highlight the nearest data point highlightn(clickp,hand) case 2 switch get(hand.method,'Value') case {1,3,4} % find the nearest break breaks = str2num(get(hand.params(1,4),'String')); %#ok<ST2NM> breaks is a vector [~,V]=min(abs(breaks(:)-repmat(clickp(1),length(breaks),1))); if get(hand.params(1,4),'Value')~=V set(hand.params(1,4),'Value',V), markb(hand) end case 2 highlightn(clickp,hand) end %switch get(hand.method,'Value') end %switch get(hand.data_etc,'Value') end %if tmp(1)=='n' case 'change_name' % change the name of the current spline names = get_names(get(hand.list_names,'String')); V = get(hand.list_names,'Value'); oldname = names(V,:); oldname = deblank(oldname(6:end)); name = ''; while isempty(name) answer = inputdlg(... 'Give a (valid MATLAB variable) name for the current spline:',... 'Please provide a new name for the current spline ...', 1,... {oldname}); if isempty(answer)||strcmp(answer{:},oldname), return, end name = full_name(answ