function T=phispan(x,phitype,dilation,uo,a,bord1,bord2,Family,par);% USAGE T=phispan(x,phitype,dilation,uo,a,bord1,bord2) % % This function returns a matrix Tij that contains elements% \phi_j(x_i). The size of the matrix is N*M where N is the number% of data x and M the number of parametric functions.%% phitype   is a string containing the type of parametric functions desired.%           'sin', 'sin2','unite','wavelet',sin2D    %% dilation, uo,a, bord1, bord2 are parameters that should be used for wavelet span%%% Octobre 2000 - ARswitch lower(phitype)    case 'linear'        %[L,d]=size(x)        %x0=ones(size(x))        %x1=[x(:,1), zeros(L,1)]        %x2=[zeros(L,1),x(:,2)]        %x3=x1.*x2        %T=[x0 x1 x2 x3]                        [L,d]=size(x);        x0=ones(L,1);        x1=x(:,1);        x2=x(:,2);        x3=x1.*x2;        T=[x0 x1 x2 x3];           case 'quadra'        %[L,d]=size(x)        %x0=ones(size(x))        %x1=x(:,1)        %x2=x(:,2)        %x3=x1.*x2        case 'sin'    phi0=ones(size(x));    phi1=sin(x);    phi2=cos(x);    T=[phi0 phi1 phi2];case  'sin2'    phi0=ones(size(x));    phi3=sin(x);    phi4=cos(x);    phi1=sin(2*x);    phi2=cos(2*x);    T=[phi0 phi1 phi2 phi3 phi4];    case 'unite'    phi0=ones(size(x));    T=phi0;case 'sinc_sin'    phi0=ones(size(x));    phi3=sin(x);    phi4=cos(x);    phi5=sinc(x-1);    phi6=sinc(x-2);	    phi7=sinc(x-3);    phi8=sinc(x-4);    phi9=sinc(x-5);    phi10=sinc(x-6);    phi11=sinc(x-7);    phi12=sinc(x-8);    phi13=sinc(x-9);    T=[phi0 phi3 phi4 phi5 phi6 phi7 phi8 phi9 phi10 phi11 phi12 phi13];case 'sinc13_sin'    phi0=ones(size(x));    phi3=sin(x);    phi4=cos(x);    phi5=sin(1*pi*(x-1))./(1*pi*(x-1));    phi6=sin(1*pi*(x-2))./(1*pi*(x-2));	    phi7=sin(1*pi*(x-3))./(1*pi*(x-3));    phi8=sin(1*pi*(x-4))./(1*pi*(x-4));    phi9=sin(1*pi*(x-5))./(1*pi*(x-5));    phi10=sin(1*pi*(x-6))./(1*pi*(x-6));    phi11=sin(1*pi*(x-7))./(1*pi*(x-7));    phi12=sin(1*pi*(x-8))./(1*pi*(x-8));    phi13=sin(1*pi*(x-9))./(1*pi*(x-9));    phi14=sin(3*pi*(x-1))./(3*pi*(x-1));    phi15=sin(3*pi*(x-2))./(3*pi*(x-2));    phi16=sin(3*pi*(x-3))./(3*pi*(x-3));    phi17=sin(3*pi*(x-4))./(3*pi*(x-4));    phi18=sin(3*pi*(x-5))./(3*pi*(x-5));    phi19=sin(3*pi*(x-6))./(3*pi*(x-6));    phi20=sin(3*pi*(x-7))./(3*pi*(x-7));    phi21=sin(3*pi*(x-8))./(3*pi*(x-8));    phi22=sin(3*pi*(x-9))./(3*pi*(x-9));    T=[phi0 phi3 phi4 phi5 phi6 phi7 phi8 phi9 phi10 phi11 phi12 phi13 phi14 phi15 phi16 phi17 phi18 phi19 phi20 phi21 phi22];case 'sinc13'    phi4=sin(2*x);    phi5=sin(1*pi*(x-1))./(1*pi*(x-1));    phi6=sin(1*pi*(x-2))./(1*pi*(x-2));	    phi7=sin(1*pi*(x-3))./(1*pi*(x-3));    phi8=sin(1*pi*(x-4))./(1*pi*(x-4));    phi9=sin(1*pi*(x-5))./(1*pi*(x-5));    phi10=sin(1*pi*(x-6))./(1*pi*(x-6));    phi11=sin(1*pi*(x-7))./(1*pi*(x-7));    phi12=sin(1*pi*(x-8))./(1*pi*(x-8));    phi13=sin(1*pi*(x-9))./(1*pi*(x-9));    phi14=sin(3*pi*(x-1))./(3*pi*(x-1));    phi15=sin(3*pi*(x-2))./(3*pi*(x-2));    phi16=sin(3*pi*(x-3))./(3*pi*(x-3));    phi17=sin(3*pi*(x-4))./(3*pi*(x-4));    phi18=sin(3*pi*(x-5))./(3*pi*(x-5));    phi19=sin(3*pi*(x-6))./(3*pi*(x-6));    phi20=sin(3*pi*(x-7))./(3*pi*(x-7));    phi21=sin(3*pi*(x-8))./(3*pi*(x-8));    phi22=sin(3*pi*(x-9))./(3*pi*(x-9));    T=[phi5 phi6 phi7 phi8 phi9 phi10 phi11 phi12 phi13 phi14 phi15 phi16 phi17 phi18 phi19 phi20 phi21 phi22];    case 'wavelet';    T=waveletspan(x,dilation,uo,a,bord1,bord2);    case 'orthowavelet';    i=1;    nx=length(x);       n=1024;    vector=linspace(0,1,n);        for p=1:length(dilation)        for k=1:2^dilation(p)            aux=makewavelet(dilation(p),k,Family,par,'Father',n);            for m=1:nx;                T(m,i)=valinterp(aux,vector,x(m));            end;                        i=i+1;        end;    end;    case 'sin2d';        phi0=ones(length(x),1);    T=[phi0 sin(x(:,1)).*sin(x(:,2))];    end;