質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
for

for文は、様々なプログラミング言語で使われている制御構造です。for文に定義している条件から外れるまで、for文内の命令文を繰り返し実行します。

MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

Q&A

解決済

2回答

595閲覧

matlabでfor文を使って計算した結果を一つの行列にまとめる方法をご教示ください

tohakujp

総合スコア10

for

for文は、様々なプログラミング言語で使われている制御構造です。for文に定義している条件から外れるまで、for文内の命令文を繰り返し実行します。

MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

0グッド

0クリップ

投稿2022/06/09 10:24

編集2022/06/10 21:50

質問内容の詳細

matlabで加速度の時間変化データを処理したいと考えております。matlabのヘルプページを読み、for文を使い、一定区間毎にpsdを計算するところまではできました。
1回の計算の度出た結果を列に追加していきたいのですが、方法が分かりません。
ご教示のほど、よろしくお願い申し上げます。
数日間悩んでおりますが、解決方法が見つかりません。そこでこちらのページにたどり着きました。初めて質問するもので、作法に従っていないところ多々あるかと思います。何卒ご容赦ください。

実現したいこと

  • for文で計算した結果を1つずつ新規の列に追加して1つの行列として保存したい

該当のソースコード

#wsize = 4096(ウィンドウサイズ) #fs = 5000(サンプリング周波数) #psdE(処理対象のデータ。0.0002[s]毎に計測された加速度[G]の1次元のデータ。実際は237568×1 double) #size1:処理対象のデータの総数(psdEのデータ数。実際は237568) for k = 0:wsize:(size1/wsize-1)*wsize [pxx,f] = pwelch(psdE(k+1:k+wsize),wsize,4000,fs,fs,"onesided","psd","mean"); end

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

jbpb0

2022/06/09 10:35

「size1」とは何でしょうか? 「psdE」のデータ数でしょうか? 「psdE」は1次元のデータでしょうか?
tohakujp

2022/06/09 11:31

ご指摘ありがとうございます。 情報が抜けており、大変失礼いたしました。 size1:処理対象のデータの総数(psdEのデータ数) psdE:0.0002[s]毎に計測された加速度[G]の1次元のデータ 以上となります。 よろしくお願いいたします。
guest

回答2

0

Signal Processing Toolboxを使える環境が今手元に無いので動作確認できませんが、質問に掲載されてるコードがエラー出ずに実行できるなら、たぶんこんな感じでいけるのではないかと思います

matlab

1[pxx,f] = pwelch(psdE(1:wsize),wsize,4000,fs,fs,"onesided","psd","mean"); 2pxx2 = zeros(length(pxx), length(0:wsize:(size1/wsize-1)*wsize)); 3ind = 0; 4for k = 0:wsize:(size1/wsize-1)*wsize 5 ind = ind + 1; 6 [pxx,f] = pwelch(psdE(k+1:k+wsize),wsize,4000,fs,fs,"onesided","psd","mean"); 7 pxx2(:, ind) = pxx; 8end

投稿2022/06/10 12:41

jbpb0

総合スコア7651

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

jbpb0

2022/06/10 12:44

「自己解決」に気付かずに回答してしまいました 失礼しました
guest

0

自己解決

kが1ずつ増加する変数になっていなかったのが原因のようでした。
そこで、全データ数をウィンドウサイズで割ると58であるため(237568/wsize = 58)、これを用いて、

for j = 1:1:58
[pxx,f] = pwelch(psdE(4096j-4095:4096j),4096,4000,5000,fs,"onesided","psd","mean");
psdT(:,j) = pxx;
end

とすることで解決できました(解決できたと思っています)。

投稿2022/06/10 10:54

tohakujp

総合スコア10

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

jbpb0

2022/06/10 12:50

> 全データ数をウィンドウサイズで割ると58であるため(237568/wsize = 58)、これを用いて、 みたいな数値をコード中に入れると、データの長さが変わったり、wsizeを変えたりした時に変更するのを忘れて、ハマりますよ
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問