代码出错
clc
clear
close all
Untitled =xlsread('E:\试验数据\风场实测\单天时程数据\1112\1112超声波.xlsx');
%fprintf("load data successfully\n");
X=floor(size(Untitled,1)/6000)-1;
result=cell(X,1);
Gust=cell(X,1);
I_result=cell(X,1);
L_result=cell(X,1);
Corrcoef_result=cell(X,1);
A_result=cell(X,1);
B_result=cell(X,1);
m_result=cell(X,1);
Sse_result=cell(X,1);
Rsquare_result=cell(X,1);
Rmse_result=cell(X,1);
Dfe_result=cell(X,1);
Adjrsquare_result=cell(X,1);
%--架构：对所有样本进行循环计算；进行折叠五分钟的样本wind
for i=1:X
a=6000*(i-1)+1;
b=6000*i;
wind=Untitled(a:b,:);
end
%*对单个10分钟样本即6000个数据进行分析
%--剔除僵值（考虑到下一步命令，将僵值修改为50，便于下一步处理）
%--定义：以连续二十个相同的数为僵值.
%-第二列数据僵值处理
for ii=1:5980
if wind.VarName2(ii)==wind.VarName2(ii+1)&wind.VarName2(ii)==wind.VarName2(ii+2)&wind.VarName2(ii)==wind.VarName2(ii+3)&wind.VarName2(ii)==wind.VarName2(ii+4)&wind.VarName2(ii)==wind.VarName2(ii+5)...
&wind.VarName2(ii)==wind.VarName2(ii+6)&wind.VarName2(ii)==wind.VarName2(ii+8)&wind.VarName2(ii)==wind.VarName2(ii+9)&wind.VarName2(ii)==wind.VarName2(ii+10)&&wind.VarName2(ii)==wind.VarName2(ii+11)...
&wind.VarName2(ii)==wind.VarName2(ii+12)&wind.VarName2(ii)==wind.VarName2(ii+13)&wind.VarName2(ii)==wind.VarName2(ii+14)&wind.VarName2(ii)==wind.VarName2(ii+15)&wind.VarName2(ii)==wind.VarName2(ii+16)...
&wind.VarName2(ii)==wind.VarName2(ii+17)&wind.VarName2(ii)==wind.VarName2(ii+18)&wind.VarName2(ii)==wind.VarName2(ii+19)&abs(wind.VarName2(ii+20)-wind.VarName2(ii))>0.5;
wind.VarName2(ii)=50;wind.VarName2(ii+1)=50;wind.VarName2(ii+2)=50;wind.VarName2(ii+3)=50;wind.VarName2(ii+4)=50;wind.VarName2(ii+5)=50;
wind.VarName2(ii+6)=50;wind.VarName2(ii+7)=50;wind.VarName2(ii+8)=50;wind.VarName2(ii+9)=50;
wind.VarName2(ii+10)=50;wind.VarName2(ii+11)=50;wind.VarName2(ii+12)=50;wind.VarName2(ii+13)=50;wind.VarName2(ii+14)=50;wind.VarName2(ii+15)=50;
wind.VarName2(ii+16)=50;wind.VarName2(ii+17)=50;wind.VarName2(ii+18)=50;wind.VarName2(ii+19)=50;
end
end