find関数とfor​関数をもちいた全デー​タの取得について

21 views (last 30 days)
ryu
ryu on 28 Feb 2022
Commented: ryu on 1 Mar 2022
find関数とfor関数を用いて以下のような計算を行いたいです(重心動揺外周面積の計算)
ちなみにデータは重心動揺計のデータでx軸とy軸のデータがあります
theta = mod(atan2d(X,Y)+360,360) ;
X(x軸方向のデータ),Y(y軸方向のデータ)に対して原点からの角度thetaを求めました。
for i = 0 :3: 360 ; %3°ずつずらす
a = find(theta > i & theta < i+3) ; %i<theta<i+3の範囲にあるデータ数を抽出
end
その後上記のようにfind関数を使用して0から3°ずつ範囲を移動(合計120区間=360°/3°)させてその範囲内のデータをaに取得していきたいです。
しかしfindで抽出されるデータ数がそれぞれの区間によって異なるためaに代入することができません。
(ex.0-3°の範囲ではデータ数10個,3-6°の範囲ではデータ数15個→aに代入できない)
このようにデータ数が違うものをforを使用して取得する方法はございますでしょうか。
またもっといい方法がございましたら教えていただいたいです。

Accepted Answer

Hernia Baby
Hernia Baby on 28 Feb 2022
cell配列 に格納してください
for ii = 1:10
a{ii} = randi(10,ii,ii);
end
a
a = 1×10 cell array
{[8]} {2×2 double} {3×3 double} {4×4 double} {5×5 double} {6×6 double} {7×7 double} {8×8 double} {9×9 double} {10×10 double}
cellfunを用いて一気にそれぞれの要素サイズを見ていきます。
cellfun(@size,a,'UniformOutput',false)
ans = 1×10 cell array
{[1 1]} {[2 2]} {[3 3]} {[4 4]} {[5 5]} {[6 6]} {[7 7]} {[8 8]} {[9 9]} {[10 10]}
サイズが異なるものが格納できました
  2 Comments
ryu
ryu on 1 Mar 2022
非常にわかりやすい説明ありがとうございます。
無事にプログラムを完成させることができました。

Sign in to comment.

More Answers (1)

Atsushi Ueno
Atsushi Ueno on 1 Mar 2022
X = rand(10000,1)*2-1; % X: -1.0~1.0の乱数10000個
Y = rand(10000,1)*2-1; % Y: -1.0~1.0の乱数10000個
theta = mod(atan2d(X,Y)+360,360);
[N,edges] = histcounts(theta, 0:3:360)
N = 1×120
58 71 58 54 57 89 56 85 70 97 90 103 82 103 139 136 137 109 106 91 86 75 65 75 85 71 68 64 45 53
edges = 1×121
0 3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48 51 54 57 60 63 66 69 72 75 78 81 84 87
histogram(theta, 0:3:360);
  1 Comment
ryu
ryu on 1 Mar 2022
このような関数があるとは知りませんでした。
ありがとうございます。使用したいと思います。

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!