連番名のcsvファイ​ルからforループ等​を使って標準偏差をも​とめたい

18 次查看(过去 30 天)
H.O
H.O 2020-8-12
评论: H.O 2020-8-16
プログラミング初心者です.
①ファイルは:N-1.csv, N-2.csv, N-3.csv, , N-4.csv,, N-5.csv, N-6.csv・・・・N-9.csv
  の連番で,すべて29X29の行列です.
②1セットを3枚(N)とし,5セット(cyc)の平均(Dave)と偏差(平均との差),標準偏差を求めようとしています.
③統計処理で対応するファイルは,3枚とびなのでそれぞれ
 N-1.csv, N-4.csv, N-7.csv,  ・・・1セット目.これに平均・標準偏差の処理をして,29X29の行列Na-1.csv として保存
N-2.csv, N-5.csv, N-8.csv, ・・・2セット目.これに平均・標準偏差の処理をして,29X29の行列Na-2.csv として保存
N-3.csv, N-6.csv, N-9.csv, ・・・3セット目.これに平均・標準偏差の処理をして,29X29の行列Na-3.csv として保存
 ・・・5セットまで実行します.
④”平均”や標準偏差は行列 要素(mean、std)でなく,同じサイズの複数の行列間 との演算です.
⑤標準偏差は,エクセル関数STDEV.Pに該当するものとします.
平均処理ができたプログラムは以下です.※実際には扱うファイル数が多いため変数を使用して作成しています.
同様に偏差,標準偏差の処理を行いたいのでご教示頂けませんでしょうか.
matlab2013で機能が乏しい場合は2019でも大丈夫です.
宜しくお願い致します.
data = zeros(29,29);
n= N*cyc
for a=1:N % N=3 ...3枚飛びを指定
Dsum=0;
for number = a:N:n
num = num2str(number);
data= xlsread(['絶対パス\N-' num '.csv']);
%1セット分の平均を求める
Dsum = Dsum + data ;
Dave=Dsum/cyc ; %平均
No=num2str(a);
xlswrite(['絶対パス\Na-' No '.csv'],Dave); %CSVファイルとして連番保存
contourf(Dave,200,'LineStyle','none'); %カラーの分布画像にする
end
end
  2 个评论
Kenta
Kenta 2020-8-14
偏差や標準偏差を求める関数があります。ネットで、「matlab 標準偏差」などと調べてみてはいかがでしょうか。
H.O
H.O 2020-8-14
閲覧頂きありがとうございます。検討継続致します。

请先登录,再进行评论。

采纳的回答

Kenta
Kenta 2020-8-15
clear;clc
cyc=5;
N=3;
n=N*cyc;
for a=1:N % N=3 ...3枚飛びを指定
Dsum = [];
for number = a:N:n
num = num2str(number);
% data= xlsread(['絶対パス\N-' num '.csv']);
data=rand(29); % データがないので乱数で29×29の行列を作成
%1セット分の平均を求める
Dsum = cat(3,Dsum,data);
end
% Dave=Dsum/cyc ; %平均
Dave=mean(Dsum,3);
Dstd=std(Dsum,[],3);
% No=num2str(a);
% xlswrite(['絶対パス\Na-' No '.csv'],Dave); %CSVファイルとして連番保存
% contourf(Dave,200,'LineStyle','none'); %カラーの分布画像にする
end
こんにちは、少しコードを修正しました。5セット分繰り返すのならはじめのfor a=1:Nは、for a=1:cycな気がしますがいかがでしょうか?ひとまずそのままにしています。ただ、
  • N-1.csv, N-4.csv, N-7.csv,  ・・・1セット目.これに平均・標準偏差の処理をして,29X29の行列Na-1.csv として保存
  • N-2.csv, N-5.csv, N-8.csv, ・・・2セット目.これに平均・標準偏差の処理をして,29X29の行列Na-2.csv として保存
  • N-3.csv, N-6.csv, N-9.csv, ・・・3セット目.これに平均・標準偏差の処理をして,29X29の行列Na-3.csv として保存
とありますが、これを5セット繰り返すと例えば、N-7は、2度使われますが大丈夫なのでしょうか?
Daveも、奥行方向(Z方向)にデータを重ねて、最後にチャンネル方向に平均(または標準偏差)の計算をしました。あまり、質問の条件を汲み取れてないかもしれませんので、適宜変更を加えてください。
  1 个评论
H.O
H.O 2020-8-16
ありがとうございます。 意図するものができました。質問して良かったです。
>N-7は、2度使われますが大丈夫なのでしょうか? 失礼しました、例が悪かったです。 実際にはcycよりNの方が数が大きいため同じファイルが2回以上使われることはありません。
とても助かりました。ありがとうございました。

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Matrix Indexing 的更多信息

产品


版本

R2013a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!