Matlab矩阵处理小结（2）-读数据求均值和中位数 绘制errorbar图

files = dir('*.log');

for i=1:length(files)

File =files(i).name;

a=importdata(File);

b=reshape(a,6,12)';%calculate tramsportation to make std calculation easy

kth(:,:,i)=b;

end

for i=1:24

mymean=mean(kth(:,:,i));

mystd=std(kth(:,:,i));

mymedian=median(kth(:,:,i));

% myresult(:,:,1)=cat(2,myresult(:,:,1),mymean)

% myresult(:,:,2)=cat(2,myresult(:,:,2),mystd)

% a=cat(2,a,mymean);

% b=cat(2,b,mystd);

y(i,:)=mymean;

e(i,:)=mystd;

m(i,:)=mymedian;

end

y

e

m

fid=fopen('v1mean.txt','wt');

fprintf(fid,'%g\n',y);

fclose(fid);

fid1=fopen('v1mid.txt','wt');

fprintf(fid1,'%g\n',m);

fclose(fid1);

x=[1:1:24];

for i=1:6

subplot(2,3,i)

% figure('hello')

% text('Message transformation of Node v1')

% scatter(x,m(:,i),10,'Red');

plot(x,m(:,i),'RED')

hold on

errorbar(x,y(:,i),e(:,i));

% errorbar(x,y(:,i),1)

%title('type"i"')

xlabel('time');

ylabel('message');

switch i

case 1

title('type1')

case 2

title('type2')

case 3

title('type3')

case 4

title('type4')

case 5

title('type5')

case 6

title('type6')

end

end