[メモ]本日のScilab

1/(1.0+0.01*(x-50.0)^2.0);が上に凸のグラフになってしまう罠発生

---

test=(1.0+0.01*(50-50.0)^2.0);
test =1

test=(1.0+0.01*(x-50.0)^2.0);
test(51) =1

test=1/(1.0+0.01*(x-50.0)^2.0);
test(51) =0.0000669

---

ダレカタスケテ、、、

どうやら正しくは「test=1 ./ (1.0+0.01 .* (x-50.0) .^ 2.0);」らしい、演算子を「./」にしないと行列の割り算になってしまい、数字がおかしくなるらしい。


あとあれ、比較演算子を使う（例えば「muA>=0.2」）と答えは真理値で「....FFFFTT.....」みたいな感じになるけど、グラフとして表示するためにはbool2s関数で数値に変換しないといけないらしい、、、

ファジィ理論のアルファカット(α-cut)の例をのせてみる

---

clear;

//alpha_cut

alpha=0.2

x=[0:1:100];//xの範囲

muA=1.0 ./ (1.0+0.01 .* (x-50) .^ 2.0);
alpha_cut = bool2s(muA>=alpha);//真理値を数値へ

//表示
//軸の範囲 rect=[x最小値,y最小値,x最大値,y最大値];
rect=[0,0,100,1];
//刻みの数 tics=[x大きな刻み,x小さな刻み,y大きな刻み,y小さな刻み];
tics=[10,11,1,11];

//title=["表題","x軸名","y軸名"];
title=["Membership func (muA)","x","muA"];
//フレーム作成
plotframe(rect,tics,[%f,%f],title);
//グラフ1を表示
plot(x,muA,":r");
//グラフ２を表示
plot(x,alpha_cut,"-b");
//凡例 legend("グラフ１","グラフ２",右上);
legend("Menbership","alpha cut",1);

---