matlab悬链线方程的求解,Matlab建模教程-变分法简介.doc
Matlab建模教程-变分法简介
§1 变分法简介
作为数学的一个分支,变分法的诞生,是现实世界许多现象不断探索的结果,人们可以追寻到这样一个轨迹:
约翰·伯努利(Johann Bernoulli,1667-1748) 这就是著名的“最速降线”问题(The Brachistochrone Problem)Guillaume Francois Antonie de l'Hospital 1661-1704)、雅可比·伯努利(Jacob Bernoulli 1654-1705)Gottfried Wilhelm Leibniz,1646-1716)Isaac Newton1642—1727)Euler Lonhard,1707~1783)Lagrange, Joseph Louis,1736-1813约翰·伯努利(The Hanging Chain Problem)向界求答案固定的端,在重力中它自然垂下,的曲方程是什在大自然中,除了的外,我們可以察到吊上方的索,水珠以及根之所架的,些都是(catenary)。
伽利略(Galileo, 1564~1643)惠更斯(Huygens, 1629~1695)在1646年(17岁),物理的,得知伽利略的猜不,但他求不出答案。约翰·伯努利
解此方程并适当选取参数,得
(1)
即为悬链线。
悬链线问题本身和变分法并没有关系,然而这和最速降线问题一样都是贝努利兄弟间的相互争强好胜、不断争吵的导火索,虽然雅可比·贝努利在解决悬链线问题时略占下风,但他随后所证明的“悬挂于两个固定点之间的同一条项链,在所有可能的形状中,以悬链线的重心最低,具有最小势能”,算是扳回了一局,俩兄弟扯平了!之所以提到悬链线问题,有两方面考虑,其一,这是有关数学史上著名的贝努利家族内的一个趣闻,而这是一个在变分法乃至整个数学物理领域有着巨大贡献的家族,其二,有关悬链线的得几个结论,可以用变分法来证明!
现实中很多现象可以表达为泛函极小问题,我们称之为变分问题。求解方法通常有两种:古典变分法和最优控制论。我们这儿要介绍的基本属于古典变分法的范畴。
变分法的基本概念
泛函的概念
设为一函数集合,若对于每一个函数有一个实数与之对应,则称是定义在上的泛函,记作。称为的容许函数集。
例如,在上光滑曲线y(x)的长度可定义为
(2)
考虑几个具体曲线,取,
若,则
若y(x)为悬链线,则
对应中不同的函数y(x),有不同曲线长度值J,即J依赖于y(x),是定义在函数集合上的一个泛函,此时我们可以写成
我们称如下形式的泛函为最简泛函
(3)
被积函数包含自变量,未知函数(t)及导数(t)。上述曲线长度泛函即为一最简泛函。
泛函极值问题
考虑上述曲线长度泛函,我们可以提出下面问题:
在所有连接定点的平面曲线中,试求长度最小的曲线。
即,求,使
取最小值。此即为泛函极值问题的一个例子。以极小值为例,一般的泛函极值问题可表述为,
称泛函在取得极小值,如果对于任意一个与接近的,都有。所谓接近,可以用距离来度量,而距离可以定义为
泛函的极大值可以类似地定义。其中称为泛函的极值函数或极值曲线。
泛函的变分
如同函数的微分是增量的线性主部一样,泛函的变分是泛函增量的线性主部。作为泛函的自变量,函数在的增量记为
也称函数的变分。由它引起的泛函的增量记作
如果可以表为
其中为的线性项,而是的高阶项,则称为泛函在的变分,记作
。用变动的代替,就有。
泛函变分的一个重要形式是它可以表为对参数的导数:
(4)
这是因为当变分存在时,增量
根据和的性质有
所以
泛函极值的相关结论
泛函极值的变分表示
利用变分的表达式(4),可以得到有关泛函极值的重要结论。
泛函极值的变分表示:若在达到极值(极大或极小),则
(5)
证明:对任意给定的,是变量的函数,该函数在处达到极值。根据函数极值的必要条件知
再由(4)式,便可得到(5)式。
变分法的基本引理:,,,有
,
则。
证明略。
泛函极值的必要条件
考虑最简泛函(3),其中F具有二阶连续偏导数,容许函数类S取为满足端点条件为固定端点(6)的二阶可微函数。
, (6)
泛函极值的必要条件:设泛函(3)在x(t)∈S取得极值,则x(t)满足欧拉方程
matlab悬链线方程的求解,Matlab建模教程-变分法简介.doc相关推荐
- matlab二元方程区间求解,matlab求解二元方程组
陈星似 魔法师 matlab求解二元方程组 悬赏分:0 提问时间:2010-11-30 23:29回答数:1浏览量:241问题指向:全国 t1=(q1+q2+q3+q4-q5-q6-q7)/g1/c1 ...
- matlab狐狸与野兔数学模型,数学建模狐狸野兔问题.doc
狐狸野兔问题 摘要:封闭自然环境中的狐狸和野兔存在捕食与被捕食关系,本题旨在通过对自然状态下两物种数量变化规律的分析,推测加入人类活动(即人工捕获)时两物种数量的变化,进而得出人类活动对自然物种的影响 ...
- matlab 画三维花瓶,精美花瓶建模教程
1.首先,草图单位为mm,进入前视图绘制如图草图,花瓶的基本形状轮廓 2.然后对草图进行旋转 3.旋转出曲面后,在顶部边线新建一个基准面 4.继续在前视图绘制草图,如图绘制一弧线 5.然后进行旋转 6 ...
- matlab曲面方程,如何用Matlab拟合数据求解曲面方程 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...
我不知道你的方程,我粗略得到一个R2 = 0.99以上的方程,代码及结果如下 高次多项式拟合代码 : load shuju.mat B=fit([Dotx,Dotx],Dotz,'poly43') 结 ...
- 基于matlab弹道方程代码,基于Matlab/Simulink的弹丸外弹道6自由度运动仿真
第 23卷 第 4期 Vol. 23 No. 4 重 庆 工 学 院 学 报 (自然科学 ) Journal of Chongqing Institute of Technology(Natural ...
- matlab解方程的程序,matlab算法程序解方程.ppt
<matlab算法程序解方程.ppt>由会员分享,可在线阅读,更多相关<matlab算法程序解方程.ppt(10页珍藏版)>请在人人文库网上搜索. 1.一 般 的 代 数 方 ...
- 永磁直流电机 matlab仿真,永磁直流电动机的Simulink建模仿真优秀教学.doc
个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 永磁直流电动机地Simulink建模仿真教学-机电论文 永磁直流电动机地Simuli ...
- matlab复数方程的根,matlab解超越方程的复数根
本帖最后由 zhaocb 于 2015-7-18 16:22 编辑 各位大神,我需要解一个超越方程,'0=(r^2+(50*m*(2*pi*k*C)^2))*tanh(r*0.04)+i*r*(50- ...
- matlab怎么选择高频变压器,高频变压器的建模与仿真定.doc
高频变压器的建模与仿真定 毕业设计(论文) 课题名称 高频变压器的建模与仿真 学生姓名 学 号 系 . 专 业 电气工程系.0级电气工程及其自动化 指导教师 职 称 讲 师 2013 年 5月 15日 ...
最新文章
- 2020年度“中国神经科学重大进展”获奖名单【附成果介绍】
- Android控件美化Shape
- python版本选择-【小白学python】之一:版本选择
- vivox21升级鸿蒙,vivo X21i相机规格再升级,加持AI成又一拍照神器
- Python2与Python3区别
- 计算机专业指南考试题6,计算机专业指南考试题.doc
- 系统动力学软件vensim学习之一阶负反馈
- Windows 下编译 OpenSSL
- 使用SHP数据批量裁剪TIF图像并转为JPG格式
- 老毛桃装ubuntu
- 支付宝运营10元花呗红包活动被薅羊毛了
- 程序员需知的11个在线教程网站,建议收藏!
- oracle—ebs_采购功能点操作手册,oracle—EBS_采購功能点操作手册.doc
- linux上数据库账号被锁定,ORA-28000账户被锁和解锁
- 安全生产双重预防体系建设数字化解决方案
- 忆鲁迅《故乡》中的一句话
- Fighting!Hack for wuhan 黑客马拉松等你来战
- 51nod 1298
- 安徽涉外经济职业学院有计算机网络,校级精品课程——会计基础
- java 判断文件类型是否是音频_用java流方式判断文件类型