Matlab矩陣處理小結(2)-讀資料求均值和中位數 繪製errorbar圖
阿新 • • 發佈:2019-02-09
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
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