close all

clear all

fidin=fopen('20131118B.txt');                 %打开文件one minute

fidout=fopen('data.txt','w');                   %创建文件，用于写进数据

while ~feof(fidin)                              %判断是否为文件末尾

tline=fgetl(fidin);                             %从文件中读取行

if double(tline(1))>=48&&double(tline(1))<=57   %判断首字符是否是数值

fprintf(fidout,'%s\n\n',tline);                  %如果是数字行，把此行数据写入文件data中

continue                                         %如果是非数字继续下一次循环

end

end

fclose(fidout);

M=importdata('data.txt');                   %将生成的data.txt文件导入到工作空间。

% M=importdata('one minute.txt');

%%%%%%%%%%%%%%%%%%数据处理%%%%%%%%%%%%%%%%%%%

t=M.textdata(:,:);

p=M.data(:,4);        %取倒数第4列数据，即第6列数据。

DD=[t num2cell(M.data)];

xlswrite('C:\Users\Administrator\Desktop\2345.xlsx',DD,'Sheet1');

%%%%%%%%%%%%%%%%%%%%求标准差和斜率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

b=p;

M=30;                        %M为周期数，自己可以改

M1=60;

for i=1:length(b)-M+1

ysum(i)=sum(b(i:M+i-1));

xsum(i)=sum(i:i+M-1);

xxsum(i)=0;

xysum(i)=0;

for j=i:M+i-1

xxsum(i)=xxsum(i)+j*j;

xysum(i)=xysum(i)+j*b(j);

end

K(i)=(xysum(i)-((ysum(i)/M)*xsum(i)))/(xxsum(i)-((xsum(i)/M)*xsum(i)));

end

K1=[ones(1,M-1)*K(1),K];

xlswrite('C:\Users\Administrator\Desktop\2345.xlsx',K1','Sheet1','J');

for i=1:length(b)-M1+1

ysum(i)=sum(b(i:M1+i-1));

xsum(i)=sum(i:i+M1-1);

xxsum(i)=0;

xysum(i)=0;

for j=i:M1+i-1

xxsum(i)=xxsum(i)+j*j;

xysum(i)=xysum(i)+j*b(j);

end

K(i)=(xysum(i)-((ysum(i)/M1)*xsum(i)))/(xxsum(i)-((xsum(i)/M1)*xsum(i)));

end

K2=[ones(1,M1-1)*K(1),K];

xlswrite('C:\Users\Administrator\Desktop\2345.xlsx',K2','Sheet1','K');

错误如图：

数据见附件！

2013-11-19 18:05 上传

k列多出的数据是j列最后的一部分数据，不知道为什么会这样？

2013-11-19 18:07 上传

点击文件名下载附件

12.7 KB, 下载次数: 5