math_chap10 (1)

星期二,2013-5-21,05:50:51,第10章数学问题的非传统解法,高等应用数学问题的MATLAB求解,清华大学出版社2008,CAI课件开发：薛定宇、刘莹莹、董雯彬,星期二,2013-5-21,05:50:51,第10章数学问题的非传统解法,集合论、模糊集与模糊推理粗糙集理论与应用人工神经网络及其在数据拟合中的应用进化算法及其在最优化问题中的应用小波变换及其在数据处理中的应用分数阶微积分学问题求解及应用,星期二,2013-5-21,05:50:51,10.1集合论、模糊集与模糊推理,经典可枚举集合论问题及MATLAB求解模糊集合与隶属度函数模糊推理系统及其MATLAB求解,星期二,2013-5-21,05:50:51,10.1.1经典可枚举集合论问题及MATLAB求解,集合论是现代数学的基础集合：一些事物的全体，而其中每一个事物均称为集合中的一个元素。可枚举集合：该集合中的所有元素均可以一一列出的集合在MATLAB中，可枚举集合用向量或单元数组的形式就可以表示这样的集合,星期二,2013-5-21,05:50:51,集合表示举例,数字构成的集合，可以有重复元素上述集合的单元数组表示方法，二者等价字符串集合，可以为人名等,星期二,2013-5-21,05:50:51,MATLAB下集合运算的函数,常用的集合运算函数适用于几何变量,星期二,2013-5-21,05:50:51,例10.1,给定，对这些集合进行各种运算，并验证：进行各种集合运算：,星期二,2013-5-21,05:50:51,验证交换律：使用ismember()函数求集合的交,星期二,2013-5-21,05:50:51,例10.2,给定试求集合A与B的并集与交集，并验证结合律,星期二,2013-5-21,05:50:51,进行集合运算：验证结合律：,星期二,2013-5-21,05:50:51,求解的函数调用格式求解的函数调用格式,集合“属于”关系判定,星期二,2013-5-21,05:50:51,例10.3,给定验证是否验证集合A的自反律亦即,星期二,2013-5-21,05:50:51,例10.4,考虑哥德巴赫猜想，即，任何大于2的偶数均能分解两个质数的和该猜想是Euler时代以来，尚未严格证明的最古老的数论问题计算机现在证明从4到1019，猜想都是正确的本例不是试图证明该猜想，而是由MATLAB对小偶数验证该猜想,星期二,2013-5-21,05:50:51,用集合运算的方法验证小于2000的偶数均满足该猜想由某范围内的两质数所有的可能的和构造出一个集合，然后判定是否有限偶数均属于该集合,哥德巴赫猜想的验证,星期二,2013-5-21,05:50:51,10.1.2模糊集合与隶属度函数,经典集合：一个事物a，要么就属于集合A，要么就不属于集合A美国控制论专家LotfiAZadeh教授，1965拓展了此集合，提出了模糊集合fuzzyset的概念，说明a在一定程度下属于AZadeh教授指出，当问题的复杂性增加时，精确的描述将失去意义，而有意义的描述将失去精度。,星期二,2013-5-21,05:50:51,星期二,2013-5-21,05:50:51,例10.5,年老与年轻的模糊表示及隶属度函数，假设论域U=0,120，则,星期二,2013-5-21,05:50:51,模糊集合论域隶属度a在一定程度下属于A绘制出隶属度函数：,星期二,2013-5-21,05:50:51,钟形隶属度函数,数学描述：函数调用格式：其中，x是独立变量的值，a,b为参数y是隶属度函数的值,星期二,2013-5-21,05:50:51,例10.6,绘制出不同参数a,b,c组合下的钟形隶属度函数曲线MATLAB求解语句：,星期二,2013-5-21,05:50:51,Gauss隶属度函数,数学描述：函数调用格式：其中，x是独立变量的值y是隶属度函数的值,星期二,2013-5-21,05:50:51,例10.7,绘制不同参数c,s组合下Gauss隶属度函数MATLAB求解语句：,星期二,2013-5-21,05:50:51,Sigmoid型隶属度函数,数学描述：函数调用格式：其中，x是独立变量的值y包含隶属度函数的值,星期二,2013-5-21,05:50:51,例10.8,绘制不同参数a,c组合下Gauss隶属度函数MATLAB求解语句：,星期二,2013-5-21,05:50:51,隶属度函数的图形编辑界面,使用GUI编辑隶属度函数打开编辑界面的函数调用格式修改隶属度函数的参数添加隶属度函数,星期二,2013-5-21,05:50:51,10.1.3模糊推理系统及其MATLAB求解,模糊化模糊规则解模糊化,星期二,2013-5-21,05:50:51,模糊逻辑工具箱中提供的函数可以构建出模糊推理系统的数据结构，函数的调用格式其中，name为字符串，表示模糊推理系统的名称，通过该函数可以建立起结构体fis，其内容包括模糊的与、或运算，解模糊算法等,建立模糊推理系统,星期二,2013-5-21,05:50:51,添加系统的输入和输出变量的调用格式定义一个输入变量iname定义一个输出变量oname其中，及为输入或输出变量的取值范围每个变量的隶属函数可以用addmf()函数定义，也可以用mfedit()定义,定义输入、输出变量,星期二,2013-5-21,05:50:51,某模糊推理系统有两个输入变量，并有一个输出变量，又有，分为3个区间，隶属函数选择为钟形函数，和，分为3个区间，隶属函数选择为Gauss型函数；，隶属函数为Sigmoid型函数，则可以用下面的语句构造模糊推理系统,例10.9,星期二,2013-5-21,05:50:51,MATLAB求解语句：建立模糊推理系统定义模糊变量在GUI中用交互方式定义隶属函数,星期二,2013-5-21,05:50:51,10.1.3.1模糊化,若将某信号用三个隶属函数表示，则一般对应的物理意义是“很小”、“中等”与“较大”，若分为5段，则可以表示为“很小”、“较小”、“中等”、“较大”和“很大”，一个精确的信号可以通过这样一组隶属函数模糊化，变成模糊信号,星期二,2013-5-21,05:50:51,10.1.3.2模糊规则,模糊推理关系可以表示成规则矩阵不加到模糊推理系统fis,星期二,2013-5-21,05:50:51,10.1.3.3解模糊化,通过模糊推理可以得出模糊输出量op，再通过指定的算法精确化，亦称解模糊化模糊逻辑工具箱提供了多种解模糊化的算法求解模糊推理问题的函数调用格式其中，X为矩阵，其各列为各个输入信号的精确值，变量y是输出信号的精确值,星期二,2013-5-21,05:50:51,例10.10,某模糊推理系统有两个输入变量，并有一个输出变量。又有，分为3个区间，隶属函数选择为钟形函数，和，分为3个区间，隶属函数选择为Gauss型函数；，隶属函数为Sigmoid型函数，用此模糊推理系统绘制出输出的三维曲面,星期二,2013-5-21,05:50:51,MATLAB求解语句：,星期二,2013-5-21,05:50:51,10.2粗糙集理论与应用,粗糙集理论简介粗糙集的基本概念信息决策系统粗糙集数据处理问题的MATLAB求解粗糙集约简的MATLAB程序界面,星期二,2013-5-21,05:50:51,10.2.1粗糙集理论简介,粗糙集(roughset)是波兰数学家ZPawlak为开发自动规则生成系统及研究软计算问题于1982年提出的粗糙集理论是一种处理不精确、不确定与不完全数据的新的数学方法能有效地分析和处理不精确、不一致、不完整等各种不完备信息，并从中发现隐含的知识，揭示潜在的规律,星期二,2013-5-21,05:50:51,10.2.2粗糙集的基本概念,下近似集：根据现有知识判断肯定属于X的对象组成的最大的集合，称为正区，记为POS(X)上近似集：那些可能属于X的对象组成的最小集合,星期二,2013-5-21,05:50:51,边界集c：上近似集与下近似集的差集如果BND(X)是空集，则称X是清晰的若BND(X)非空，则称X为粗糙集,星期二,2013-5-21,05:50:51,例10.11,玩具积木集合具有“颜色R1”、“形状R2”、“体积R3”属性分别对应红、黄和绿色分别对应方、圆、三角形分别对应于“大”和“小”的物体集合举例红色的积木、绿色的积木蓝色的积木粗糙集表示,星期二,2013-5-21,05:50:51,10.2.3信息决策系统,信息决策系统可以表示为U是对象的集合，即论域A是属性集合，C条件属性集D决策属性集其中，CD=A并且CD=,星期二,2013-5-21,05:50:51,粗糙集理论中使用决策表来描述论域中对象,星期二,2013-5-21,05:50:51,例10.12,积木集合的4个相关属性：（颜色,形状,大小,价位）黄色，红色，绿色；方形，圆形，三角形；大，其余小；价低，价中，价高销售情况：好，一般，较差,星期二,2013-5-21,05:50:51,信息决策系统表,星期二,2013-5-21,05:50:51,函数调用格式Sl求出下近似集求出上近似集Su利用MATLAB的差集函数可以求出边界集Sd,近似集的MATLAB计算,星期二,2013-5-21,05:50:51,编辑函数rslower(),下近似集函数,星期二,2013-5-21,05:50:51,编辑函数rsupper(),上近似集函数,星期二,2013-5-21,05:50:51,不可分辨关系的定义：构造MATLAB函数,不可分辨关系,星期二,2013-5-21,05:50:51,例10.13,论域关系为且若，试求出集合X的上近似集和下近似集,星期二,2013-5-21,05:50:51,旋转的信息系统决策表MATLAB求解语句：,星期二,2013-5-21,05:50:51,接上页,星期二,2013-5-21,05:50:51,10.2.4粗糙集数据处理问题的MATLAB求解,利用粗糙集理论的约简粗糙集理论在信息约简中的应用举例,星期二,2013-5-21,05:50:51,10.2.4.1利用粗糙集理论的约简,条件约简：不含多余属性并保证分类正确的最小条件属性集找出从C属性中选定条件推出D的最小集合核集：所有不可约去的关系求取从C属性中选定条件推出D的核集,星期二,2013-5-21,05:50:51,例10.14,试用粗糙集的方法对进行约简，找出不必要的数码管,星期二,2013-5-21,05:50:51,数码管显示真值表,星期二,2013-5-21,05:50:51,MATLAB求解语句：,星期二,2013-5-21,05:50:51,例10.15,SARS是2003年给全球带来恐慌的疾病，其准确诊断是很困难的。这里给出从报刊提取出的一些数据，如下表所示，试利用粗糙集理论对给出的12个条件进行约简，找出辅助诊断的最主要的条件这里的例子只用于演示粗糙集的约简，数据也说不全面的，不能用于临床诊断,星期二,2013-5-21,05:50:51,SARS患者和正常人若干检测指标表,星期二,2013-5-21,05:50:51,MATLAB求解语句：最后得出的条件为Y=3,4，表示第3和4列是诊断SARS的重要因素，亦即“血液检测呈阳性”和“高烧28o”,星期二,2013-5-21,05:50:51,10.2.5粗糙集约简的MATLAB程序界面,启动粗糙集约简界面的函数调用格式用户可以由其中的Browse按钮读入信息系统决策表，给出C属性和D属性所需的列号，则可以进一步进行分析,星期二,2013-5-21,05:50:51,10.3人工神经网络及其在数据拟合中的应用,神经网络基础知识径向基网络结构与应用神经网络界面人工神经网络是在对复杂的生物神经网络研究和理解的基础上发展起来的。人工神经元互联构成的网络,星期二,2013-5-21,05:50:51,10.3.1神经网络基础知识,前馈神经网络结构与应用神经网络的训练与泛化,星期二,2013-5-21,05:50:51,10.3.1.1前馈神经网络结构与应用,单个人工神经元的数学表示形式,星期二,2013-5-21,05:50:51,Sigmoid函数对数Sigmoid函数也可以使用简单的饱和函数和阶跃函数等作为传输函数,常用的传输函数,星期二,2013-5-21,05:50:51,例10.16,试绘制出各种常用的传输函数曲线。MATLAB求解语句：结果：,星期二,2013-5-21,05:50:51,接上页：,星期二,2013-5-21,05:50:51,BP神经网络结构示意图,星期二,2013-5-21,05:50:51,两层网络示意图,两层网络,星期二,2013-5-21,05:50:51,两层网络的数学表示,隐层节点在传输函数前后的值分别为输出层传输函数前后的信号分别为训练的目标函数残差最小,星期二,2013-5-21,05:50:51,训练的样本,随机选择训练用样本通过训练更新权值,星期二,2013-5-21,05:50:51,早期版本新版本隐层节点个数向量传输函数f按需求使用fitnet()、patternfit(),构造BP网络的函数调用格式,星期二,2013-5-21,05:50:51,神经网络对象的常用属性,星期二,2013-5-21,05:50:51,例10.17,假设输入信号为2路，其信号范围分别为0,1和(-1,5)，且输出信号为单路信号，试建立所需的前馈神经网络对象MATLAB求解语句：设置网络结构设置其他参数可以直接设置,星期二,2013-5-21,05:50:51,10.3.1.2神经网络的训练与泛化,神经网络训练函数调用格式其中，变量X为nM矩阵，n为输入变量的路数，M为样本的组数，Y为mM矩阵，m为输出变量的路数目标值曲线绘制神经网络的仿真或泛化的函数调用格式,星期二,2013-5-21,05:50:51,例10.18,用下面语句生成一组数据，试用神经网络对其进行拟合，样本点生成选择2个隐层，隐层节点选择为5：,星期二,2013-5-21,05:50:51,显示出神经网络的权值选择不同的训练算法,星期二,2013-5-21,05:50:51,增加节点个数能否改善拟合效果？,星期二,2013-5-21,05:50:51,例10.19,二元函数拟合：样本数据生成选择两个隐层，每层10个节点,星期二,2013-5-21,05:50:51,单个隐层三个隐层,星期二,2013-5-21,05:50:51,10.3.2径向基网络结构与应用,径向基函数(RBF)的数学描述其中，c为聚类中心点，而b0为调节聚类效果的参数神经网络工具箱中，radbas()函数可以计算出标准径向基函数的曲线参数,星期二,2013-5-21,05:50:51,若隐层的传输函数F1(x)为径向基函数，输出层的传输函数F2(x)为线性函数，则此结构的网络称为径向基网络,星期二,2013-5-21,05:50:51,例10.20,试绘制不同参数(c,b)下的径向基函数曲线取中心点c=-2,0,2，并假设b=1,星期二,2013-5-21,05:50:51,选择中心点c=0，并假设b=0.1,1,5径向基网络的训练方式不是采用反向误差传播实现的，故不是BP网络用函数newrbe()和sim()可以实现神经网络的建立、训练和泛化全过程,星期二,2013-5-21,05:50:51,例10.21,用径向基神经网络重新拟合一元函数径向基网络的拟合,星期二,2013-5-21,05:50:51,例10.22,用径向基神经网络重新拟合二维曲面生成样本数据曲面拟合,星期二,2013-5-21,05:50:51,10.3.3神经网络界面,启动神经网络工具箱的图形用户界面的函数调用格式该界面可以建立所需的神经网络模型，并可以由已知数据对该网络进行训练、仿真,星期二,2013-5-21,05:50:51,例10.23,用神经网络工具箱的界面对例10.21中的数据进行拟合MATLAB求解语句：,星期二,2013-5-21,05:50:51,单击按钮Import来导入数据单击主界面中的NewNetwork按钮来选择神经网络的结构单击按钮Train来训练神经网络单击按钮Export将训练得出的网络模型输出到MATLAB环境中绘制出曲线拟合与泛化结果,星期二,2013-5-21,05:50:51,10.4进化算法及其在最优化问题中的应用,遗传算法的基本概念介绍及MATLAB实现遗传算法在求解最优化问题中的应用举例遗传算法在有约束最优化问题中的应用粒子群优化算法与求解其他全局优化算法求取精确的全局最优解基于遗传算法的混合整数规划求解,星期二,2013-5-21,05:50:51,10.4.1遗传算法的基本概念介绍及MATLAB实现,遗传算法是基于进化论，在计算机上模拟生命进化机制而发展起来的一门新学科，它根据适者生存、优胜劣汰等自然进化规则搜索和计算问题的解美国Michigen大学的JohnHolland于1975年提出的更可能获得全局最优解,星期二,2013-5-21,05:50:51,遗传算法工具箱,MATLAB7.0的遗传算法与直接搜索工具箱现改名全局优化工具箱用过Sheffield大学的遗传算法工具箱MATLABCentral下有大量工具箱GAOT：最早的遗传算法最优化工具箱,星期二,2013-5-21,05:50:51,遗传算法的基本思想,从一个代表最优化问题解的一组初值开始进行搜索，这组解称为一个种群，这里种群由一定数量的、通过基因编码的个体组成，其中每一个个体称为染色体，不同个体通过染色体的复制、交叉或变异又生成新的个体，依照适者生存的规则，个体也在一代一代进化，通过若干代的进化最终得出条件最优的个体。,星期二,2013-5-21,05:50:51,选择n个个体构成初始种群P0，求出各个个体的函数值，P0可以随机生成设置代数为i=1，即设置其为第一代。计算选择函数的值，所谓选择即通过概率的形式从种群中选择若干个个体的方式通过染色体个体基因的复制、交叉、变异等创造新的个体，构成新的种群Pi+1i=i+1，若终止条件不满足，则继续进化,简单遗传算法的一般步骤,星期二,2013-5-21,05:50:51,不同于从一个点开始搜索最优解的传统的最优化算法，遗传算法从一个种群开始对问题的最优解进行并行搜索，所以更利于全局最优化解的搜索。遗传算法并不依赖于导数信息或其他辅助信息来进行最优解搜索。遗传算法采用的是概率型规则而不是确定性规则，所以每次得出的结果不一定完全相同，有时甚至会有较大的差异。,遗传算法和传统优化算法比较,星期二,2013-5-21,05:50:51,10.4.2遗传算法在求解最优化问题中的应用举例,GAOT工具箱中的函数调用格式MATLAB全局优化工具箱打开遗传算法程序界面：gatool(),星期二,2013-5-21,05:50:51,例10.24,给定试求出f(x)取最大值时x的值绘制目标函数曲线：试测不同初值：,星期二,2013-5-21,05:50:51,不同初值下搜索到的“最优解”及目标函数值,星期二,2013-5-21,05:50:51,编写目标函数：调用gaopt()函数：与用最优化工具箱得出的结果进行比较：,遗传算法求解,星期二,2013-5-21,05:50:51,更大的求解空间，假设改变遗传算法的求解区间,星期二,2013-5-21,05:50:51,例10.25,给定试求其的最小值描述最优化问题的目标函数：,星期二,2013-5-21,05:50:51,假设为自变量的求解范围,星期二,2013-5-21,05:50:51,gaopt()函数稍复杂一点的调用格式其中，p为附加参数，v为精度及显示控制向量，P0为初始种群，fun1为终止函数的名称，默认值为maxGenTerm，n为最大的允许代数该函数还有更复杂的调用格式,函数其他调用格式,星期二,2013-5-21,05:50:51,例10.26,给定试设置更多的允许代数，观察寻优的结果，并和最优化搜索算法得出结果在时间、精度上的差异进行比较,星期二,2013-5-21,05:50:51,利用ga()函数，编写目标函数指定一些搜索参数,星期二,2013-5-21,05:50:51,用传统的无约束最优化方法求解得出更高的求解精度，求解时间大大减少,星期二,2013-5-21,05:50:51,10.4.3遗传算法在有约束最优化问题中的应用,GAOT工具箱不能直接用于有约束最优化问题的求解对于有约束最优化问题，应该将其转变为无约束最优化问题；GAOT工具箱求解的是最大值不等式约束则可以通过惩罚函数方法转换到目标函数中等式约束，则可以通过等式求解的方式将其中若干个自变量用其他自变量表示,星期二,2013-5-21,05:50:51,例10.27,试用遗传算法求解线性规划问题对于等式，用x1和x2表示x3,星期二,2013-5-21,05:50:51,描述目标函数：运用遗传算法求解问题：,星期二,2013-5-21,05:50:51,用线性规划方法求解：linprog()得出的结果更精确线性规划是凸问题，没有必要用遗传算法建议求解方法：用GA找出全局最优解的大致位置以其为初值调用最优化函数求精确解,星期二,2013-5-21,05:50:51,全局优化工具箱函数求解,有约束优化问题的直接求解等式约束效果似乎不佳求解的是最小值，而不像gaopt的最大值得出的是可行解,星期二,2013-5-21,05:50:51,例10-28,用ga函数直接求解线性规划问题可以用ga函数得到可行解，但不一定是真正的最优解也可以使用结构体型调用格式,星期二,2013-5-21,05:50:51,求解语句结构体求解,星期二,2013-5-21,05:50:51,10.4.4粒子群优化算法与求解,粒子群优化算法是一种进化算法，该算法是受生物界鸟群觅食的启发而提出的搜索食物假设某个区域内有一个全局最优点，和位于随机初始位置的粒子，每一个粒子有到目前为止自己的个体最优值pi,b，整个粒子群有到目前为止群体的最优值gb,星期二,2013-5-21,05:50:51,粒子按下式更新自己的速度和位置其中，g1a，g2a为0,1区间内均匀分布的随机数f(k)为惯量函数a1，a2为加速常数,星期二,2013-5-21,05:50:51,求解粒子群优化的函数调用格式其中，fun为目标函数（不支持匿名函数），n是x向量的维数，这两个量是必须提供的该函数在随书的程序包中提供支持向量化运算，目标函数的描述有所不同,星期二,2013-5-21,05:50:51,例10.29,给定用粒子群算法求其最小值编辑向量化的目标函数：求解该最优化问题,星期二,2013-5-21,05:50:51,例10.30,试用粒子群优化算法，求解线性规划问题需要用惩罚方式表示等式约束目标函数的向量描述,星期二,2013-5-21,05:50:51,编写向量化的模式的目标函数得出精确的最优解,星期二,2013-5-21,05:50:51,10.4.5其他全局搜索方法,MATLAB全局优化工具箱其他方法模拟退火方法无约束最优化问题求解模式搜索方法和fmincon调用格式一致patternsearch函数,星期二,2013-5-21,05:50:51,例10.31,线性规划问题模式搜索方法求解,星期二,2013-5-21,05:50:51,例10.32,试求出下面最优化问题的精确全局最优解绘制目标函数的三维曲面,星期二,2013-5-21,05:50:51,引入决策变量向量x1=x,x2=y直接求解模式搜索算法求解,星期二,2013-5-21,05:50:51,10.4.6求取精确的全局最优解,传统的最优化方法可能得出精确的最优解，然而，经常可能得到的是局部最优解而不是全局最优解，而进化方法通常能得出全局最优解，但解的精度可能很低在实际最优化问题求解中可以考虑将二者的优势结合起来，得到精确的全局最优解,星期二,2013-5-21,05:50:51,重新求解模式搜索方法求初值，传统方法求精确解,例10-33,星期二,2013-5-21,05:50:51,用遗传算法求解问题：将遗传算法得出的求解最优解作为搜索起点，采用传统最优化方法,星期二,2013-5-21,05:50:51,10.3.7基于遗传算法的混合整数规划求解,求解函数如果采用结构体求解，则使用IntCon局限性：目前只能求解不含有等式约束的问题非线性约束f1函数也不能含有等式约束,星期二,2013-5-21,05:50:51,例10-35,用遗传算法重新求解x1是0.25的整数倍，x2是0.1的整数倍引入，整数规划,星期二,2013-5-21,05:50:51,非线性约束基于遗传算法的求解,星期二,2013-5-21,05:50:51,10.5小波变换及其在数据处理中的应用,小波变换及基小波波形小波变换技术在信号处理中的应用小波问题的程序界面,星期二,2013-5-21,05:50:51,10.5.1小波变换及基小波波形,连续小波变换离散小波变换小波工具箱中提供的基小波函数,星期二,2013-5-21,05:50:51,10.5.1.1连续小波变换,连续小波变换的变换公式为其中，且(t)称为基小波，为基小波通过平移、比例缩放构成的小波信号,星期二,2013-5-21,05:50:51,例10.35,假设“墨西哥帽”基小波函数由下式给出试绘制出不同a,b值变换下的小波函数。,星期二,2013-5-21,05:50:51,小波反变换的数学描述为：其中，计算连续小波变换的系数的函数调用格式计算小波系数矩阵小波系数绝对值图其中，fun为基小波名称,小波反变换求解,星期二,2013-5-21,05:50:51,例10.36,试对信号f(t)=sint2进行连续小波分解，并绘制出其系数图绘制出时域信号曲线：选择mexh基小波作为模板绘制出小波系数的三维表面图,星期二,2013-5-21,05:50:51,10.5.1.2离散小波变换,离散信号的小波变换的数学描述：函数调用格式离散小波反演公式的数学描述函数调用格式,星期二,2013-5-21,05:50:51,例10.37,给定信号f(t)=sint2叠加均值0，标准差0.1的Gaussian白噪声MATLAB求解语句：小波反变换，误差,星期二,2013-5-21,05:50:51,10.5.1.3小波工具箱中提供的基小波函数,生成和测试显示基小波类型的函数调用格式计算基小波函数的函数调用格式Gauss、墨西哥帽Daubechies/SymletsBior等,星期二,2013-5-21,05:50:51,例10.38,绘制出不同阶次下的Daubechies6基小波波形MATLAB求解语句：,星期二,2013-5-21,05:50:51,例10.39,试绘制出常用基小波波形MATLAB求解语句：,星期二,2013-5-21,05:50:51,10.5.2小波变换技术在信号处理中的应用,小波分解的过程,星期二,2013-5-21,05:50:51,一维信号的小波分解的函数调用格式分解结果及存储方式：,一维小波处理,星期二,2013-5-21,05:50:51,提取系数的函数调用格式提取近似系数提取第i段细节系数信号重建的函数调用格式,小波系数的提取与小波重建,星期二,2013-5-21,05:50:51,例10.40,给定信号f(t)=sint2叠加均值0，标准差0.1的Gaussian白噪声对数据进行3次小波分解，各种基小波函数对其进行降噪处理，比较降噪效果绘制原信号波形曲线：,星期二,2013-5-21,05:50:51,对给定的数据进行3次小波分解：,星期二,2013-5-21,05:50:51,采用db6基小波采用基小波bior2.6和coil4,星期二,2013-5-21,05:50:51,例10.41,如下数字滤波问题试用小波对其进行滤波，并比较滤波效果,星期二,2013-5-21,05:50:51,MATLAB求解语句：,星期二,2013-5-21,05:50:51,10.5.3小波问题的程序界面,启动小波变换问题的图形用户界面他可以用于求解一维、二维小波变换问题,星期二,2013-5-21,05:50:51,10.6分数阶微积分学问题求解及应用,分数阶微积分的定义不同分数阶微积分定义的关系与性质分数阶微积分的计算方法分数阶线性微分方程的求解方法基于框图的非线性分数阶微分方程近似解法分数阶传递函数模型与分析,星期二,2013-5-21,05:50:51,分数阶微积分学概述,第3章的整数阶微积分问题1695年法国数学家GuillaumeFranoisAntoineLHpital问过微积分学创造者之一GottfriedWilhelmLeibniz的问题，n=1/2?分数阶微积分理论已经有300年的历史了早期主要侧重于理论研究近年来在很多领域都已经开始应用自动控制领域出现了分数阶控制理论等新的分支,星期二,2013-5-21,05:50:51,10.6.1分数阶微积分的定义,四种常用的定义、相互关联分数阶Cauchy积分公式Grnwald-Letnikov分数阶微积分定义Riemann-Liouville分数阶微积分公式Caputo分数阶微分定义,星期二,2013-5-21,05:50:51,分数阶Cauchy积分公式,数学描述：其中，C为包围f(t)单值与解析开区域的光滑曲线,星期二,2013-5-21,05:50:51,Grnwald-Letnikov分数阶微积分定义,数学描述：其中，为二项式系数,星期二,2013-5-21,05:50:51,Riemann-Liouville分数阶微积分公式,积分的数学描述：其中，0a1，且a为初始时间分数阶微分定义的数学描述：,星期二,2013-5-21,05:50:51,Caputo分数阶微分定义,微分的数学描述：其中，a=m+g，m为整数，0g1当g0时，Caputo分数阶积分的数学描述,星期二,2013-5-21,05:50:51,10.6.2不同分数阶微积分定义的关系与性质,Caputo定义和Riemann-Liouville定义0a1更一般地,星期二,2013-5-21,05:50:51,分数阶微积分的性质,解析函数f(t)的分数阶导数对t和a都是解析的。a为整数时，分数阶微分与整数阶微分的值完全一致，且分数阶微积分算子为线性的，对常数a、b分数阶微积分算子满足交换律,星期二,2013-5-21,05:50:51,函数分数阶导数的Laplace变换,函数的分数阶积分表达式的Laplace变换Riemann-Liouville微分的Laplace变换为Caputo/Riemann-Louiville定义函数积分的Laplace完全一致，微分的Laplace变换满足,星期二,2013-5-21,05:50:51,10.6.3分数阶微积分的计算方法,用Grnwald-Letnikov定义求解分数阶微分Caputo微积分定义的数值计算Oustaloup滤波算法改进的Oustaloup滤波器方法分数阶线性微分方程的求解方法零初值分数阶线性微分方程的解法非零初值Caputo微分方程的数值求解,星期二,2013-5-21,05:50:51,10.6.3.1用Grnwald-Letnikov定义求解分数阶微分,数学描述其中，,星期二,2013-5-21,05:50:51,构造一个MATLAB函数函数调用格式,分数阶微分数值计算,星期二,2013-5-21,05:50:51,例10-42,阶跃函数的导数和积分是什么？整数阶的积分与微分对不同的阶次,星期二,2013-5-21,05:50:51,例10.43,考虑一个初值非零的函数试求出其分数阶导数。MATLAB求解语句：,星期二,2013-5-21,05:50:51,对不同的g选值，可以调用下面的语句绘制出分数阶导函数的三维图,星期二,2013-5-21,05:50:51,例10.44,试用不同定义求取函数为f(t)=sin(3t+1)的0.75阶微分，并比较得出的结果。Cauchy公式结果MATLAB求解语句：,星期二,2013-5-21,05:50:51,8.6.3.2Caputo微积分定义的数值计算,Caputo分数阶积分与Grnwald-Letnikov定义完全一致Caputo/Grnwald-Letnikov微分关系MATLAB函数实现,星期二,2013-5-21,05:50:51,MATLAB代码,星期二,2013-5-21,05:50:51,例10-45,函数为f(t)=sin(3t+1)，求0.3、1.3、2.3阶的Caputo导数0.3阶两个定义之差为绘制0.3阶导数,星期二,2013-5-21,05:50:51,其余两个导数曲线绘制,星期二,2013-5-21,05:50:51,10.6.3.3Oustaloup滤波算法,如果函数f(t)的表达式不是事先已知连续滤波器传递函数模型：Oustaloup算法感兴趣区间,星期二,2013-5-21,05:50:51,构造MATLAB函数，设计连续滤波器函数调用格式,Oustaloup滤波器求解,星期二,2013-5-21,05:50:51,例10.46,函数感兴趣频率区间分数阶阶次0.5，5阶滤波器MATLAB求解语句：绘制分数阶微分曲线：,星期二,2013-5-21,05:50:51,更大的频率范围不同的滤波器阶次,星期二,2013-5-21,05:50:51,10.6.3.4改进的Oustaloup滤波器方法,数学描述：其中，通常可以选择加权参数b=10，d=9,星期二,2013-5-21,05:50:51,编写出下面的MATLAB函数函数调用格式：,星期二,2013-5-21,05:50:51,例10.47,给定函数相同频率范围0.5阶导数MATLAB求解语句：,星期二,2013-5-21,05:50:51,10.6.4分数阶微分方程求解,一类分数阶线性系统时域响应解析解方法零初值分数阶线性微分方程的解法非零初值Caputo微分方程的数值求解非线性微分方程求解非线性非零初值Caputo微分方程求解,星期二,2013-5-21,05:50:51,分数阶线性微分方程,一般形式初值为零，Riemann-Liouville微分方程初值不为零，Caputo微分方程传递函数,星期二,2013-5-21,05:50:51,10.6.4.1一类分数阶线性系统时域响应解析解方法,如果方程右侧只有u(t)，微分方程的解析解,Mittag-Leffler函数，第8章,星期二,2013-5-21,05:50:51,左侧只含有三项的方程解析解,传递函数其中,星期二,2013-5-21,05:50:51,例10-48,零初值微分方程输入u(t)=1得出直接求解此方法适用面窄，计算量大,星期二,2013-5-21,05:50:51,10.6.4.2零初值分数阶线性微分方程的解法,分数阶线性微分方程如果右侧不是u(t)，可以通过下式计算假设如果有负阶次，应该事先变换,星期二,2013-5-21,05:50:51,由Grnwald-Letnikov定义方程左侧每项均含y(t)，其余项移到等号右侧，则微分方程的闭式数值解为,星期二,2013-5-21,05:50:51,分数阶微分方程求解函数,星期二,2013-5-21,05:50:51,例10-49,求解零初值分数阶线性微分方程选择计算步长h=0.002求解，选择不同的h检验,星期二,2013-5-21,05:50:51,10.6.4.3非零初值Caputo微分方程的数值求解,微分方程中是Caputo分数阶导数如果，应该已知m个初始值引入辅助信号，将原方程变换为零初值,星期二,2013-5-21,05:50:51,辅助信号变换函数新微分方程,星期二,2013-5-21,05:50:51,星期二,2013-5-21,05:50:51,例10-50,微分方程初始条件MATLAB求解,星期二,2013-5-21,05:50:51,10.6.5基于框图的非线性分数阶微分方程近似解法,基于框图解法的必要性复杂非线性方程没有通用的数值解法如果微分方程是整个系统的一部分求解方法设计一个分数阶微分算子，可以采用Oustaloup滤波器或改进形式可以考虑初值，处理Caputo微分方程由分数阶微分算子搭建系统仿真模型,星期二,2013-5-21,05:50:51,10.6.5.1零初值非线性分数阶微分方程的求解,微分算子的封装初始化栏目语句,星期二,2013-5-21,05:50:51,例10-51,框图求解零初值线性分数阶微分方程引入变换后,星期二,2013-5-21,05:50:51,星期二,2013-5-21,05:50:51,例10-52,非线性分数阶微分方程（零初值）变换形式,星期二,2013-5-21,05:50:51,星期二,2013-5-21,05:50:51,10.6.5.2非零初值的Caputo微分方程数值解法,需要对初值补偿构造封装模块,星期二,2013-5-21,05:50:51,取定义辅助信号由poly2caputo()填写Fcn模块的补偿函数可以封装Caputo微分算子模块仿照前面例子建模,星期二,2013-5-21,05:50:51,例10.53,非线性分数阶微分方程已知初值解析解辅助信号等效于,星期二,2013-5-21,05:50:51,微分方程变换,直接得出化简新微分方程,星期二,2013-5-21,05:50:51,星期二,2013-5-21,05:50:51,10.6.6对象编程实例分数阶传递函数模型,分数阶传递函数（FOTF）对象FOTF对象的创建FOTF短信的重载函数编写面向对象的分析与设计方法,星期二,2013-5-21,05:50:51,10.6.6.1FOTF对象定义,FOTF类用于描述分数阶传递函数，其作用类似于MAT