大家帮我分析分析

飘零de叶@_@

>> function [ddr,hour1,hour2,hour3] = chebshev2(sp3,n)

x0=sp3(:,4)*60+sp3(:,5);y0=sp3(:,7:9);        %x0为以分钟为单位的时间序列，y0为坐标
x=min(x0):0.5:max(x0);
nn=length(x0);
m=length(x);
c=x0(n)-x0(1);

r0=2*(x0-x0(1))/c-1;         %变量变为（-1，1）区间
T=zeros(n,n);
for i=1:nn
        T(i,1)=1;
        T(i,2)=r0(i);
        for j=3:n
                T(i,j)=2*r0(i)*T(i,j-1)-T(i,j-2);
        end
end

C=inv(T'*T)*T'*y0;       %最小二乘求切比雪夫多项式系数
%*————————求未知历元拟合坐标——————————
X=zeros(1,n);y=[];

for i=1:m
    k=2*(x(i)-x0(1))/c-1;
        X(1)=1;
        X(2)=k;
        for j=3:n
                X(j)=2*k*X(j-1)-X(j-2);
    end       
    y=[y;X*C];
end
%————————求切比雪夫拟合精度（内符合精度‘中误差’）————————
zz=T*C;cz=zz-y0;
ddrx=sqrt(cz(:,1)'*cz(:,1)/(nn-n-1));
ddry=sqrt(cz(:,2)'*cz(:,2)/(nn-n-1));
ddrz=sqrt(cz(:,3)'*cz(:,3)/(nn-n-1));   %   xyz值拟合中误差
ddr=sqrt(ddrx^2+ddry^2+ddrz^2);    %坐标拟合精度
ddr
%————————取出y中各小时的拟合数据————————————
[mm1,nn1]=size(y);
code=121:120:mm1;hour1=[];hour2=[];hour3=[];   %定义三个小时的矩阵
hour1=y(1:code(1),;
hour2=y(code(1):code(2),;
hour3=y(code(2):code(3),;


format long g


function [ddr,hour1,hour2,hour3] = chebshev2(sp3,n)
|
Error: Function definitions are not permitted in this context.
求高手帮我分析分析问题出在哪，怎么解决啊   给您跪下了