大家帮我分析分析
>> function = 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=;
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中各小时的拟合数据————————————
=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 = chebshev2(sp3,n)
|
Error: Function definitions are not permitted in this context.
求高手帮我分析分析问题出在哪，怎么解决啊 给您跪下了
页:
[1]