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にすると上記ｴﾗｰが発生します．
初歩的な質問かと存じますがご教授いただけますと幸いです．
采纳的回答
Atsushi Ueno
2021-5-21
编辑：Atsushi Ueno
2021-5-21
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)」を入力したので三番目のパターンにも該当しないと判断され、エラーが発生します。
4 个评论
Atsushi Ueno
2021-5-21
上記プログラムで検証したのですが、どんな形にしてもエラーを再現できませんでした。
説明に誤りがありました。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の転置については、次元が合致していればどちらでも動作します(下記説明にも書いてありました)。
Atsushi Ueno
2021-5-21
修正した回答内容は間違いないですが、質問にあるエラーを再現できませんでした。
（なぜtrapz関数の2番目のパターンに該当しないのか理由が不明な為）
必要あれば当該エラーの発生が再現するプログラム(BとCの定義を含む)を添付願います。
浩祐 佐々木
2021-5-25
Atsushi Ueno
2021-5-25
permuteのエラーをみて気付きましたが、i=2:10のループにする必要があります。XY一組の数値では台形が出来ないからですね。スカラ値を入力すると座標ではなく間隔として認識するので、エラーになります。
更多回答（1 个）
浩祐 佐々木
2021-5-25
0 个投票
1 个评论
Atsushi Ueno
2021-5-25
採用ありがとうございます。
- エラー: trapz (line 47) 次元の引数は、インデックス付け範囲の中で、正の整数のスカラーでなければなりません。
- エラー: permute ORDERには無効な置き換えインデックスが含まれています
1番目は「次元：正の整数スカラー」を期待する所に「座標のベクトル」が入力されたのが原因ですが、問題が再現しませんでした。（質問者様の問題が解決したので、めでたしとしましょう）
2番目は「座標のベクトル」を期待する所に「等間隔のスカラー」が入力されたのが原因でした。trapz(1,1)を入力するとエラー無しで答え0を出してしまうtrapz関数の仕様は混乱の元ですね。
类别
在 帮助中心 和 File Exchange 中查找有关 数値積分と微分 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
