cell内に格納された時系列データの平均値を算出するにはどうすればいいですか？

Question

Kohei Yoshino 2024-4-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2110316-cell

评论： Kohei Yoshino 2024-4-23

采纳的回答： Kojiro Saito

在 MATLAB Online 中打开

以下のcellデータの時系列の平均を算出したいと考えています。

↓

それぞれのcellに格納された変数のうち(:,6)の列を32個抽出して行平均を出したいと考え以下のコードを作成しましたが、meandataが32列目のデータしか格納されません。いい方法はありませんでしょうか？

data = cell(1, length(A.Pelvic)) % Aに格納されているPelvicという変数を参照
for i = 1:length(A.Pelvic)
    data{i} = A.Pelvic{i}(:,6);
    meandata = arrayfun(@mean, data{i}); % cellfunだと変数が'double'なので実行できないというエラーが出るためarrayfunを使用
end
for n = 1:length(A.Pelvic);
    plot(A.Pelvic{n}(:,6), 'b')
    hold on
    plot(meandata, 'r'); % dataをあらかじめ作成し、そこにmeandataを格納するつもりでしたが、meandataが全体の平均ではなくA.Pelvicの最後の列のみが反映されており平均できていない
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Kojiro Saito 2024-4-23

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2110316-cell#answer_1446051

在 MATLAB Online 中打开

meandataが32列目のデータしか格納されないのは、for ループの meandata = arrayfun(@mean, data{i}); で同じ変数名で上書きされているので、最後のループのi=32だけが格納されているためです。

forループを使わないでcellfunで一度で格納できます。

meandata = cellfun(@(x) mean(x(:,6)), A.Pelvic); % 1x32 double

また、プロットする際もそれぞれの6列目に対する平均値を示そうとされていると思いますので、meandata(n)でやってみたらいかがでしょうか。plot よりも定数値の線を引くyline も使えます。

for n = 1:length(A.Pelvic)
    figure;
    plot(A.Pelvic{n}(:,6), 'b')
    hold on
    %plot(meandata, 'r');
    yline(meandata(n), 'r');
    hold off
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Kohei Yoshino 2024-4-23

ご回答ありがとうございます。ご提案いただいた構文ですと、1×32のデータしか得られず、やりたいこととしては、1×100（もしくは100×1）のデータにしたいです。自分がお示しした構文にて、A.Pelvic{1:n}(:,6)のグラフを以下のようにグラフ化（青線）し、その平均（赤線）を図示したいのですが、明らかに平均値ではない値が得られてしまいます。