for loo内のtrapzのエラーについて

下記のようにfor loop内でtrapzを使用しようとしています．

A = zeros(1,1039);

for i=1:1039

Cum_current(1,i)=trapz(B(1:i,1)',C(1:i,1)');

end

B,Cは1300行*1列のベクトルです

意図としては各試行iでの積分結果をCum_currentのi列目に保存して記録しておきたいというものです．

これを実行すると下記のｴﾗｰが表示されます

エラー: trapz (line 47)

次元の引数は、インデックス付け範囲の中で、正の整数のスカラーでなければなりません。

実際にiを1039と定義してtrapzを実行すると計算は回るのですが

for loopにすると上記ｴﾗｰが発生します．

初歩的な質問かと存じますがご教授いただけますと幸いです．

Answer 1

[Atsushi Ueno]

台形則による数値積分 - MATLAB trapz - MathWorks 日本

trapz関数の引数の組み合わせは下記三種類のパターンがあります。引数の意味も併記しておきます。

Q = trapz(Y) Y — 数値データ：ベクトル | 行列 | 多次元配列

Q = trapz(X,Y) X — 点の間隔：1 (既定値) | 等間隔を表すスカラー | 座標のベクトル

Q = trapz(___,dim) dim — 演算の対象の次元：正の整数スカラー

trapz(B(1:i,1)',C(1:i,1)');

上記の書き方でtrapz関数を呼び出すと、「引数が2つあって、最初の引数がXの条件に該当しないので、二番目のパターンに該当しない」と認識されます。

エラー: trapz (line 47)

次元の引数は、インデックス付け範囲の中で、正の整数のスカラーでなければなりません。

エラー文の「次元の引数」とは、trapz関数の引数dimの事です。dimと認識すべき場所にスカラでない「C(1:i,1)」を入力したので三番目のパターンにも該当しないと判断され、エラーが発生します。

配列インデックス付け - MATLAB & Simulink - MathWorks 日本

[Atsushi Ueno]

上記プログラムで検証したのですが、どんな形にしてもエラーを再現できませんでした。

説明に誤りがありました。B(1:i,1)'はベクトルになります。

誤）上記の書き方でtrapz関数を呼び出すと「引数が2つあって、最初の引数がスカラーでもベクトルでもない」ので、三番目のパターンに該当すると認識されます。「最初の引数に入力したBは1300行*1列のベクトルだ」と仰るかもしれませんが、2次元のインデックス方法「B(1:i,1)」でアクセスしているので、ベクトルではなくマトリクス(多次元配列)と解釈されます。

正）上記の書き方でtrapz関数を呼び出すと、「引数が2つあって、最初の引数がXの条件に該当しないので、二番目のパターンに該当せず、かつ三番目のパターンにも該当しないのでエラー」と認識されます。

憶測ですがお手持ちのデータB(1:i,1)がtrapz関数の引数Xの条件に該当しないのだと思います。また、XとYの転置については、次元が合致していればどちらでも動作します(下記説明にも書いてありました)。

台形則による数値積分 - MATLAB trapz - MathWorks 日本

[Atsushi Ueno]

修正した回答内容は間違いないですが、質問にあるエラーを再現できませんでした。

（なぜtrapz関数の2番目のパターンに該当しないのか理由が不明な為）

必要あれば当該エラーの発生が再現するプログラム(BとCの定義を含む)を添付願います。

[浩祐 佐々木]

返信が大変遅くなり申し訳ございません．

大変勉強になります．

すいません，簡単な例での再現ができず手こずっておりました．

ご指摘の点も踏まえまず下記のようなプログラムを書けばうまくまわりました

A=[1 2 3 4 5 6 7 8 9 10]

B=[1 4 9 16 25 36 49 64 81 100]

C=zeros(1,7)

for i=1:7

D=trapz(A(1:i),B(1:i))

end

一方，同様のA,Bをcsv形式で別のファイルで準備し（その際A,Bは縦列データなのでcsvファイルは10*2のマトリクスだと思ってください）これを読み込む形で下記のようなプログラムを作ってみました

importdata('Book1.csv') %Book1.csvは保存したcsvファイル名だとお考えください

A=ans(1:11,1)

B=ans(1:11,2)

C=zeros(10,1)

for i=1:10

C(i)=trapz(A(1:i),B(1:i))

end

このプログラムでは下記エラーでした

エラー：permute

ORDERには無効な置き換えインデックスが含まれています

前回とは異なるエラーが出て途方にくれております．．

前回投稿した１０３９行のデータは中身や行数は違うこそすれ，やっていることは上記一番初めにトライした簡単なものと同じだと自分では認識しているのですが，何がおかしいのでしょうか．．．

＊

元の「次元は正の整数で〜」のエラーでテキストファイルを同じくimportdataで読み込み

データlabelをtextファイルとして，数値データを行列として読み込んだ構造ファイルであるという違いはあります

アドバイスいただけると幸いです

[Atsushi Ueno]

permuteのエラーをみて気付きましたが、i=2:10のループにする必要があります。XY一組の数値では台形が出来ないからですね。スカラ値を入力すると座標ではなく間隔として認識するので、エラーになります。

Answer 2

[浩祐 佐々木]

ご回答ありがとうございます．

はい，ご指摘いただきました通りi=1，が問題でした．

X,Yがスカラ担っているためのエラーでした．

上記修正したら当初目的のものも無事まわりました．

ご親切にありがとうございました．

[Atsushi Ueno]

採用ありがとうございます。

1. エラー: trapz (line 47) 次元の引数は、インデックス付け範囲の中で、正の整数のスカラーでなければなりません。
2. エラー: permute ORDERには無効な置き換えインデックスが含まれています

1番目は「次元：正の整数スカラー」を期待する所に「座標のベクトル」が入力されたのが原因ですが、問題が再現しませんでした。（質問者様の問題が解決したので、めでたしとしましょう）

2番目は「座標のベクトル」を期待する所に「等間隔のスカラー」が入力されたのが原因でした。trapz(1,1)を入力するとエラー無しで答え0を出してしまうtrapz関数の仕様は混乱の元ですね。