summaryrefslogtreecommitdiff
path: root/fuzzylite/src/defuzzifier/Bisector.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'fuzzylite/src/defuzzifier/Bisector.cpp')
-rw-r--r--fuzzylite/src/defuzzifier/Bisector.cpp59
1 files changed, 21 insertions, 38 deletions
diff --git a/fuzzylite/src/defuzzifier/Bisector.cpp b/fuzzylite/src/defuzzifier/Bisector.cpp
index ee4d2fc..ed61b6f 100644
--- a/fuzzylite/src/defuzzifier/Bisector.cpp
+++ b/fuzzylite/src/defuzzifier/Bisector.cpp
@@ -1,76 +1,60 @@
/*
- Author: Juan Rada-Vilela, Ph.D.
- Copyright (C) 2010-2014 FuzzyLite Limited
- All rights reserved
+ fuzzylite (R), a fuzzy logic control library in C++.
+ Copyright (C) 2010-2017 FuzzyLite Limited. All rights reserved.
+ Author: Juan Rada-Vilela, Ph.D. <jcrada@fuzzylite.com>
This file is part of fuzzylite.
fuzzylite is free software: you can redistribute it and/or modify it under
- the terms of the GNU Lesser General Public License as published by the Free
- Software Foundation, either version 3 of the License, or (at your option)
- any later version.
+ the terms of the FuzzyLite License included with the software.
- fuzzylite is distributed in the hope that it will be useful, but WITHOUT
- ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
- FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License
- for more details.
-
- You should have received a copy of the GNU Lesser General Public License
- along with fuzzylite. If not, see <http://www.gnu.org/licenses/>.
-
- fuzzyliteâ„¢ is a trademark of FuzzyLite Limited.
+ You should have received a copy of the FuzzyLite License along with
+ fuzzylite. If not, see <http://www.fuzzylite.com/license/>.
+ fuzzylite is a registered trademark of FuzzyLite Limited.
*/
#include "fl/defuzzifier/Bisector.h"
-#include "fl/term/Accumulated.h"
#include "fl/term/Term.h"
namespace fl {
Bisector::Bisector(int resolution)
- : IntegralDefuzzifier(resolution) {
- }
+ : IntegralDefuzzifier(resolution) { }
- Bisector::~Bisector() {
-
- }
+ Bisector::~Bisector() { }
std::string Bisector::className() const {
return "Bisector";
}
+ Complexity Bisector::complexity(const Term* term) const {
+ return Complexity().comparison(1).arithmetic(1 + 2 + 5) +
+ term->complexity().comparison(1).arithmetic(1 + 5).multiply(getResolution());
+ }
+
scalar Bisector::defuzzify(const Term* term, scalar minimum, scalar maximum) const {
- if (not fl::Op::isFinite(minimum + maximum)) {
- return fl::nan;
- }
- if (maximum - minimum > _resolution) {
- FL_DBG("[accuracy warning] the resolution <" << _resolution << "> "
- "is smaller than the range <" << minimum << ", " << maximum << ">. In order to "
- "improve the accuracy, the resolution should be at least equal to the range.");
- }
- scalar dx = (maximum - minimum) / _resolution;
+ if (not Op::isFinite(minimum + maximum)) return fl::nan;
- int counter = _resolution;
+ const scalar dx = (maximum - minimum) / getResolution();
+ int counter = getResolution();
int left = 0, right = 0;
scalar leftArea = 0, rightArea = 0;
scalar xLeft = minimum, xRight = maximum;
while (counter-- > 0) {
- if (fl::Op::isLE(leftArea, rightArea)) {
+ if (Op::isLE(leftArea, rightArea)) {
xLeft = minimum + (left + 0.5) * dx;
leftArea += term->membership(xLeft);
- left++;
+ ++left;
} else {
xRight = maximum - (right + 0.5) * dx;
rightArea += term->membership(xRight);
- right++;
+ ++right;
}
}
-
//Inverse weighted average to compensate
- scalar bisector = (leftArea * xRight + rightArea * xLeft) / (leftArea + rightArea);
- return bisector;
+ return (leftArea * xRight + rightArea * xLeft) / (leftArea + rightArea);
}
Bisector* Bisector::clone() const {
@@ -81,5 +65,4 @@ namespace fl {
return new Bisector;
}
-
}