自由始端和终端的动态规划,求指标函数最小值的逆序算法递归 % 计算程序。x是状态变量,一列代表一个阶段状态;M-函数 % DecisFun(k,x)由阶段k的状态变量x求出相应的允许决策变量 % M-函数ObjFun(k,x,u)是阶段指标函数,M-函数TransFun(k,x,u) % 是状态转移函数,其中x是阶段k的某状态变量,u是相应的决策变量; % 输出p_opt由4列构成,p_opt=[序号组 最优策略组 最优轨线组 % 指标函数值组];fval是一个列向量,各元素分别表示p_opt各 % 最优策略组对应始端状态x的最优函数值; %
上传时间: 2016-12-20
上传用户:wangzhen1990
MATLAB 数学手册复合版内容详尽 包括矩阵及其基本运算,特征值与二次型数值计算与数据分析 插值、拟合与查表 常微分方程数值解偏微分方程的数值解 符号运算积分变换Taylor级数 概率统计 随机数的产生 随机变量的概率密度计算 随机变量的累积概率值(分布函数值)正整数的频率表 4.6.2 经验累积分布函数图形 最小二乘拟合直线 绘制正态分布概率图形 绘制威布尔(Weibull)概率图形 样本数据的盒图 给当前图形加一条参考线 在当前图形中加入一条多项式曲线 样本的概率图 附加有正态密度曲线的直方图在指定的界线之间画正态密度曲 假设检验 方差分析 foptions函数 非线性规划问题 隶属函数 模糊推理结构FIS 绘图与图形处理 含有基本实例代码
上传时间: 2017-05-14
上传用户:磊子226
有限元求解柏松方程。本文采用FORTRAN语言编制程序。程序中大部分变量采用有名公共区存储方式存储,这样可以减少内存占用量。 IFG:生成有限元网格信息,即元素节点局部编码与总体编码对照表,节点实际坐标,边界节点编码与边界点上的已知值 GKD:生成总刚一维存储对角元的地址,计算总刚一维存储长度 FIXP:设置已知节点函数值 GK(NI,NJ,ADJ,AIJ):单元刚度矩阵计算 GF(NI,N,M,LE,YI,FE):单元列阵的计算 AK(I,J,AIJ):总刚度矩阵元素迭加 QEB:总刚度矩阵和总列阵合成 BDE:边界条件处理 SOLGS:Gauss-Seidel迭代法求解方程组 UDIFF(NI,NFLAG,UDIF,LE,ADJ):标准元素内形状函数导数计算 DIFF:节点上 , 加权平均
上传时间: 2017-09-12
上传用户:erkuizhang
有限元求解柏松方程。本文采用FORTRAN语言编制程序。程序中大部分变量采用有名公共区存储方式存储,这样可以减少内存占用量。 IFG:生成有限元网格信息,即元素节点局部编码与总体编码对照表,节点实际坐标,边界节点编码与边界点上的已知值 GKD:生成总刚一维存储对角元的地址,计算总刚一维存储长度 FIXP:设置已知节点函数值 GK(NI,NJ,ADJ,AIJ):单元刚度矩阵计算 GF(NI,N,M,LE,YI,FE):单元列阵的计算 AK(I,J,AIJ):总刚度矩阵元素迭加 QEB:总刚度矩阵和总列阵合成 BDE:边界条件处理 SOLGS:Gauss-Seidel迭代法求解方程组 UDIFF(NI,NFLAG,UDIF,LE,ADJ):标准元素内形状函数导数计算 DIFF:节点上 , 加权平均
上传时间: 2017-09-12
上传用户:问题问题
1. 给定函数 ,. (1) 1.在定义区间上等间隔地取5个节点,并计算节点处的函数值(取6位有效数字). 借助Matlab 求该函数的Lagrange 插值基函数以及差值多项式的表达式. 2.利用图示的方法将插值多项式曲线与曲线(1)进行比较. 3.求该函数在区间端点处的一阶和二阶导数. 4.利用Matlab 函数csape和1.所得数据求该函数的三次样条插值函数,其中,边界条件分别为固定边界条件和自然边界条件。以图示的方式将所求样条函数曲线和被插值函数曲线进行比较。 5. 写出所求三次样条函数在各个小区间上的表达式(系数取2位有效数字)。 6.利用Matlab函数fnval计算所求三次样条函数在点处的值,画图并与被插值曲线作比较。
上传时间: 2015-03-29
上传用户:kmyilang
遗传算法为群体优化算法,也就是从多个初始解开始进行优化,每个解称为一个染色体,各染色体之间通过竞争、合作、单独变异,不断进化。 优化时先要将实际问题转换到遗传空间,就是把实际问题的解用染色体表示,称为编码,反过程为解码,因为优化后要进行评价,所以要返回问题空间,故要进行解码。SGA采用二进制编码,染色体就是二进制位串,每一位可称为一个基因;解码时应注意将染色体解码到问题可行域内。 遗传算法模拟“适者生存,优胜劣汰”的进化机制,染色体适应生存环境的能力用适应度函数衡量。对于优化问题,适应度函数由目标函数变换而来。一般遗传算法求解最大值问题,如果是最小值问题,则通过取倒数或者加负号处理。SGA要求适应度函数>0,对于<0的问题,要通过加一个足够大的正数来解决。这样,适应度函数值大的染色体生存能力强。 遗传算法有三个进化算子:选择(复制)、交叉和变异。 SGA中,选择采用轮盘赌方法,也就是将染色体分布在一个圆盘上,每个染色体占据一定的扇形区域,扇形区域的面积大小和染色体的适应度大小成正比。如果轮盘中心装一个可以转动的指针的话,旋转指针,指针停下来时会指向某一个区域,则该区域对应的染色体被选中。显然适应度高的染色体由于所占的扇形区域大,因此被选中的几率高,可能被选中多次,而适应度低的可能一次也选不中,从而被淘汰。算法实现时采用随机数方法,先将每个染色体的适应度除以所有染色体适应度的和,再累加,使他们根据适应度的大小分布于0-1之间,适应度大的占的区域大,然后随机生成一个0-1之间的随机数,随机数落到哪个区域,对应的染色体就被选中。重复操作,选出群体规模规定数目的染色体。这个操作就是“优胜劣汰,适者生存”,但没有产生新个体。 交叉模拟有性繁殖,由两个染色体共同作用产生后代,SGA采用单点交叉。由于SGA为二进制编码,所以染色体为二进制位串,随机生成一个小于位串长度的随机整数,交换两个染色体该点后的那部分位串。参与交叉的染色体是轮盘赌选出来的个体,并且还要根据选择概率来确定是否进行交叉(生成0-1之间随机数,看随机数是否小于规定的交叉概率),否则直接进入变异操作。这个操作是产生新个体的主要方法,不过基因都来自父辈个体。 变异采用位点变异,对于二进制位串,0变为1,1变为0就是变异。采用概率确定变异位,对每一位生成一个0-1之间的随机数,看是否小于规定的变异概率,小于的变异,否则保持原状。这个操作能够使个体不同于父辈而具有自己独立的特征基因,主要用于跳出局部极值。 遗传算法认为生物由低级到高级进化,后代比前一代强,但实际操作中可能有退化现象,所以采用最佳个体保留法,也就是曾经出现的最好个体,一定要保证生存下来,使后代至少不差于前一代。大致有两种类型,一种是把出现的最优个体单独保存,最后输出,不影响原来的进化过程;一种是将最优个体保存入子群,也进行选择、交叉、变异,这样能充分利用模式,但也可能导致过早收敛。 由于是基本遗传算法,所以优化能力一般,解决简单问题尚可,高维、复杂问题就需要进行改进了。 下面为代码。函数最大值为3905.9262,此时两个参数均为-2.0480,有时会出现局部极值,此时一个参数为-2.0480,一个为2.0480。算法中变异概率pm=0.05,交叉概率pc=0.8。如果不采用最优模式保留,结果会更丰富些,也就是算法最后不一定收敛于极值点,当然局部收敛现象也会有所减少,但最终寻得的解不一定是本次执行中曾找到过的最好解。
标签: 遗传算法
上传时间: 2015-06-04
上传用户:芃溱溱123
为了系统深入地研究MoS2的电子能带结构和光电性质,基于密度泛函理论的第一性原理平面波赝势方 法,计算和分析了材料MoS2的电子结构及其光学性质,给出了MoS2 的能带结构、光吸收谱、电子态密度、能 量损失谱、反射谱、介电函数谱等光学性质。计算结果表明:体材料MoS2的电子跃迁形式是非垂直跃迁,具有 间接带隙的半导体材料,带隙宽为1.126 eV;价带和导带的形成是由Mo和S 的价电子起作用产生的。通过分析 其光学性质,发现MoS2的介电函数的实部和虚部的峰值都出现在低能区,当光子能量的升高,介电函数值会缓 慢降低;材料MoS2对可见紫外区域的光子具有很强的吸收,最大吸收系数为3.17×105cm-1;MoS2在能量为18.33eV 位置出现了共振现象,其它区域内能量的损失值都趋于为0,说明电子之间共振非常微弱。这些光学性质奠定了 该材料在制作微电子和光电子器件方面的作用,尤其是在紫外探测器应用方面有着潜在的应用前景,为未来对 MoS2材料的进一步研究提供理论参考。
上传时间: 2020-11-08
上传用户:
Mathlab发行的图形计算器应用于安卓设备的高品质显示屏上,对用户来说,计算更加清晰易懂且一目了然。这个程序有两大优势:首先,它不仅是一个精细的科学计算器,而且更重要的是,它在您输入过程中显示计算步骤,可以让学生观看和学习如何得出最终答案。第二,它的图形显示能力超乎寻常!不仅计算器显示图精美,而且会自动并显示生成x和y的值。本软件适用于Android平台专业版的功能* 3D图形* 全屏* 9的工作区域* 保存常量和函数库* 不要求因特网* 没有广告科学计算器* 算术表达式 +, - ,*,/,÷* 平方根,立方和多次方根 (保持‘√’密钥)* 指数,对数 (ln,log)* 三角函数sin π/2,cos 30°,...* 双曲线函数:正弦,余弦,正切,...(按“e”键切换)* 反函数(按直接功能键)* 复数,所有功能都支持复数* 导数 sin x' = cos x,... (按 x^n 键)* 科学记数法(在菜单中启用)* 百分比模式* 保存/载入历史图形计算器* 多种功能绘图* 隐函数的第二度(椭圆 2x^2+3y^2=1,等等)* 极性图 (r=cos2θ)* 参数函数,输入新线 (x=cos t,y=sin t)* 功能根和交叉点的图表,请点选的传说开启和关闭(左上角),使用菜单显示為一个列表* 图交叉口 (x^2=x+1)* 跟踪函数值和斜坡* 滚动和缩放图表* 捏放大* 横向全屏图* 函数表* 保存為图像图形* 表保存為 CSV分数计算器* 简单和复杂的分数 1/2+1/3=5/6* 混合数字时,使用空格输入值 3 1/2代数计算器* 线性方程 x+1=2 -> x=1* 二次方程 x^2-1= 0 -> x=-1,1* 较高多项式近似根* 系统线性方程组,每行写一个方程式,x1+x2=1,x1-x2=2* 多项式长除法* 多项式展开,多项式展开,因式分解矩阵计算器* 矩阵和向量运算* 点击点积 (按住*),多种功能的图形* 行列式,逆,规范,移调,跟踪库自定义* 用户定义的常量和函数* 保存/加载表达式
上传时间: 2021-12-12
上传用户:XuVshu
数值泛函与小波理论,学习小波分析的一本很好的书籍,是某高校研究生的教材,是PDF版本
上传时间: 2016-04-15
上传用户:songrui
电磁场数值计算与电磁铁设计
上传时间: 2013-08-05
上传用户:eeworm