原文:
先说一下,为什么题目是简单实现,因为我实在没有弄出好的例子。
我原来用AForge.net做的项目中的模糊系统融入了神经网络和向量机,没法抽出来当例子,就用了个最老的自动车辆的例子。
模糊系统(Fuzzy System)架构
前面零零散散说来一下有关模糊系统和模糊理论的东西,这里来个总结。
模煳系统的基本架构如图8.1所示,其中主要的功能方块包括:(1)模煳化机构、(2)模煳规则库、(3)模煳推论引擎、以及(4)去模煳化机构。
模糊化机构是有关模糊集合和隶属度函数的内容。
模糊规则前文也有提及,去模糊化只涉及了重心法一种,这里简单说一下模糊推理引擎。
模糊推理引擎是模煳系统的核心,它可以藉由近似推论或模煳推论的进行,来模拟人类的思考决策模式,以达到解决问题的目地。
比如:
前提(premise)一:x is A’
前提(premise)二:if x is A,y is B
结论:y is B’
模糊推理系统的简单实现
在AForge.Net中的对应类是InferenceSystem。
在程序实现中一个模糊推理系统(Fuzzy Inference System )由数据库(Database)和规则库(Rulebase)组成,一般操作如下:
1.获取数值输入
2.通过数据库(Database)将数值输入转为语意含义
3.验证规则库(Rulebase)中的哪些规则被输入激活
4.组合被激活的规则,得到模糊输出(Fuzzy Output)
5.去模糊化(实现IDefuzzifier接口)
下面看看例子,这个例子是关于控制车辆避免正面冲撞的系统。
输入为距离,论域[0,120],隶属度函数
输出为角度,论域[-10,50],隶属度函数
核心代码:
// 隶属度函数(距离) FuzzySet fsNear = new FuzzySet("Near", new TrapezoidalFunction(15, 50, TrapezoidalFunction.EdgeType.Right)); FuzzySet fsMedium = new FuzzySet("Medium", new TrapezoidalFunction(15, 50, 60, 100)); FuzzySet fsFar = new FuzzySet("Far", new TrapezoidalFunction(60, 100, TrapezoidalFunction.EdgeType.Left)); //绘制图像 int length1 =(int) (chart1.RangeX.Max - chart1.RangeX.Min); double[,] NearValues = new double[length1, 2]; for (int i = (int)chart1.RangeX.Min; i < chart1.RangeX.Max; i++) { NearValues[i, 0] = i; NearValues[i, 1] = fsNear.GetMembership(i); } double[,] MediumValues = new double[length1, 2]; for (int i = (int)chart1.RangeX.Min; i < chart1.RangeX.Max; i++) { MediumValues[i, 0] = i; MediumValues[i, 1] = fsMedium.GetMembership(i); } double[,] FarValues = new double[length1, 2]; for (int i = (int)chart1.RangeX.Min; i < chart1.RangeX.Max; i++) { FarValues[i, 0] = i; FarValues[i, 1] = fsFar.GetMembership(i); } chart1.UpdateDataSeries("Near", NearValues); chart1.UpdateDataSeries("Medium", MediumValues); chart1.UpdateDataSeries("Far", FarValues); // 距离(输入) LinguisticVariable lvFront = new LinguisticVariable("FrontalDistance", 0, 120); lvFront.AddLabel(fsNear); lvFront.AddLabel(fsMedium); lvFront.AddLabel(fsFar); // 隶属度函数 FuzzySet fsZero = new FuzzySet("Zero", new TrapezoidalFunction(-10, 5, 5, 10)); FuzzySet fsLP = new FuzzySet("LittlePositive", new TrapezoidalFunction(5, 10, 20, 25)); FuzzySet fsP = new FuzzySet("Positive", new TrapezoidalFunction(20, 25, 35, 40)); FuzzySet fsVP = new FuzzySet("VeryPositive", new TrapezoidalFunction(35, 40, TrapezoidalFunction.EdgeType.Left)); //绘制图像 int length2 = (int)(chart2.RangeX.Max - chart2.RangeX.Min); double[,] ZeroValues = new double[length2, 2]; for (int i = (int)chart2.RangeX.Min; i < chart2.RangeX.Max; i++) { ZeroValues[i + 10, 0] = i; ZeroValues[i + 10, 1] = fsZero.GetMembership(i); } double[,] LittlePositiveValues = new double[length2, 2]; for (int i = (int)chart2.RangeX.Min; i < chart2.RangeX.Max; i++) { LittlePositiveValues[i + 10, 0] = i; LittlePositiveValues[i + 10, 1] = fsLP.GetMembership(i); } double[,] PositiveValues = new double[length2, 2]; for (int i = (int)chart2.RangeX.Min; i < chart2.RangeX.Max; i++) { PositiveValues[i + 10, 0] = i; PositiveValues[i + 10, 1] = fsP.GetMembership(i); } double[,] VeryPositiveValues = new double[length2, 2]; for (int i = (int)chart2.RangeX.Min; i < chart2.RangeX.Max; i++) { VeryPositiveValues[i + 10, 0] = i; VeryPositiveValues[i + 10, 1] = fsVP.GetMembership(i); } chart2.UpdateDataSeries("Zero", ZeroValues); chart2.UpdateDataSeries("LittlePositive", LittlePositiveValues); chart2.UpdateDataSeries("Positive", PositiveValues); chart2.UpdateDataSeries("VeryPositive", VeryPositiveValues); //角度 LinguisticVariable lvAngle = new LinguisticVariable("Angle", -10, 50); lvAngle.AddLabel(fsZero); lvAngle.AddLabel(fsLP); lvAngle.AddLabel(fsP); lvAngle.AddLabel(fsVP); //设置数据库 Database fuzzyDB = new Database(); fuzzyDB.AddVariable(lvFront); fuzzyDB.AddVariable(lvAngle); // 创建模糊推理系统 InferenceSystem IS = new InferenceSystem(fuzzyDB, new CentroidDefuzzifier(1000)); // 直行规则 IS.NewRule("Rule 1", "IF FrontalDistance IS Far THEN Angle IS Zero"); // 左转规则 IS.NewRule("Rule 2", "IF FrontalDistance IS Near THEN Angle IS Positive"); //开始推理 // 设定输入 IS.SetInput("FrontalDistance",float.Parse(inputBox.Text)); //打印输出 try { float newAngle = IS.Evaluate("Angle"); outputBox.Text = newAngle.ToString(); } catch (Exception ex) { MessageBox.Show("error ===> "+ex.Message); }
效果:
如果要输出模糊结论可以使用:
FuzzyOutput fuzzyOutput = IS.ExecuteInference("Angle"); foreach (FuzzyOutput.OutputConstraint oc in fuzzyOutput.OutputList) { Console.WriteLine(oc.Label + " - " + oc.FiringStrength.ToString()); }
模糊推理系统的扩展和完善
借助AForge.Net我们可以很快构建一个系统,但是AForge.Net也不是十全十美的。
首先它的合成运算子没有实现完全,有些常用的,比如最大边界算子就没有。
其次其去模糊化的实现只有一种,可以考虑实现自己补充一下算法:
1.最大平均法 (modified mean of maxima defuzzifier)
2.修正型最大平均法 (modified mean of maxima defuzzifier)
3.中心平均法 (modified center average defuzzifier)
4.修正型重心法 (modified center average defuzzifier)
5.权重式平均法(weighted average method)结合了启动强度,使用更为广泛
不过AForge.Net的接口和设计比较好,扩展也是很方便的。
其实模糊系统的话,matlab也是可以做到的,但是一来matlab不开源,二来扩展确实不方便,所以没有使用,如果确实没有替代品,可以考虑混编。
模糊系统的使用本身也有一些不方便,主要是规则库的建立。
第一种也是最直接的方式就是经由询问人类专家而得。但是人类专家往往无法完整地提供所有必需的语意式模煳规则,以致于规则库的不完全。
而模糊系统的效果主要受规则和隶属度函数影响,所以有时候规则库的不全对系统的效果是毁灭性伤害。
第二种取得语意式模煳规则的方式,则是经由训练法则,从数值型资料(numerical data)中取得模煳规则,此种作法往往牵涉如何分割输入及输出变数空间。常用的就是神经网络。
AForge.Net也有神经网络相关的东西,所以使用起来很方便的。
最后附上整个项目的下载地址: