MATLABの積分について

質問

MATLABである離散値の集合で信号の波形をあらわしているのですが、その波形を積分したいのですが方法がわかりません。
減衰している波形なので関数の式はわからないため、どうすればいいか迷っています。
お願いいたします。

回答No.2

rabbit_cat

関数系は分かってないんですか。
となると、

素直に台形公式で計算するか（trapz関数）
http://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/trapz.shtml

シンプソン則で計算するか
なぜかmatlabにはシンプソンで計算する組み込み関数はないみたいなので自前でやることになりますが。
http://sach1o.blog80.fc2.com/blog-entry-20.html
こことか参考

あるいは、quadlを使うことにして、中間値は、interp1を使って補完するか
http://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/interp1.shtml

xxという名前のベクトルにxの値、yyという名前のベクトルに関数の値が入っているとして、
quadl(inline('interp1(xx,yy,x)','x','xx',yy'), x0, x1, [], [], xx, yy)
みたいな感じで。

回答No.1

rabbit_cat

quadl 関数で数値積分すればいいです。
http://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/quadl.shtml

補足

ありがとうございます。
これだと関数funの形が式で表せないと使えないように思ってしまったのですが、そんなことはないのでしょうか？よろしくお願いします。