MATLABで外れ値をプロットする方法を教えてくだ

質問

MATLABで外れ値をプロットする方法を教えてください。
　具体的には下記のとおりです。

　MATLABで、中央値＋中央値から±20％離れた時系列データをプロットするにはどうすればよいでしょうか？

冗長なprogramしか作成できないため、よい関数か、方法があれば教えていただきたく思い、
質問しました。

まず、時系列データですが、数十万点数からなる信号になっています。
この信号の中央値をだし、その中央値から±20%の範囲でひっかかるデータをplotしたく思っています。

中央値算出は、関数を知っていますので処理が早いのですが中央値から±20%の範囲でひっかかるデータをプロットするのにいまは、以下のような考えで行っています。

行データ、size(ｘ（：））　ans＝：1000000　1とおきます。

まず、
１）中央値を求める。
２）median+20%、median-20%を求める。
３）
hold on;
　for i=1:1:data-length
　　if（x(i)<=median-20% & x(i)>=median+20%） plot x(i) endif
　end

というようにしています。しかし、これではデータ一点一点検討する必要があるので冗長で、処理に時間が
かかります。

そのため、for文を使わずに、
medi-20(:)=function(x(:),median-20)
medi+20(:)=function(x(:),median+20)
で計算できるような、関数functionを教えていただけると幸いです。

よろしくお願いします。

回答No.1

m0r1_2006 · ベストアンサー

配列 x にデータが入っていて，
そこから，A < x <B の範囲のデータを取り出したければ，

index = find((x>A) & (x<B));
y = x(index);

とすれば，配列 y には A<x<B を満たす x のみ入る．

例
x = rand(1,10000);
index = find((x>0.2) & (x< 0.8));
y = x(index);
whos

お礼

大変ありがとうございます。
助かりました。

また御礼が遅くなり申し訳ありませんでした。

matlab likeのオープンソースでもこのやり方が使えるようでしたので、大変助かりました。

ありがとうございます

回答No.2

Kules

両方の回答合わせると完全なものになってますよね。
まあ関数を作るのであれば
function [m_minus20 m_plus20]=myfunction(x,percent)
みたいなものを作って、その中でmedianから外れ値の計算まで全てする感じでしょうか。

参考になれば幸いです。

お礼

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

自作関数についてお教えいただきありあとうございました。

この場合は、外から呼び出すこともできるのですね。
なるほどです。

ありがとうございます。

補足

大変ありがとうございます。回答に気づくのが遅く、質問を締め切ってなく申し訳ありません。

心としてベストアンサーに選ばせていただきたいのですが、、ベストアンサーに選べるのはお一人だけなので、、、先に回答いただいた方にお礼のベストアンサーをさせていただきます。

処理としては、おっしゃられる関数を作るほうが長大なプログラムを作る際には便利だと思いました。

回答ありがとうございました