Hamilton力学的辛算法简介

已于 2023-09-03 11:16:36 修改
阅读量575 收藏
点赞数
分类专栏： # Hamilton 系统的辛算法 文章标签： 数学建模 抽象代数 python numpy mathematica
于 2023-09-02 23:15:45 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Chandler_river/article/details/132642858
版权
  • Hamilton力学的辛算法简介
  • 冯康
  • 我的熟人和我
  • 我的朋友和我
  • 只能说 嗨 蛮好的

外微分形式与辛几何

外微分形式

  • 1-形式
  • 2-形式
  • 闭2-形式（辛构造）

Euclid Space

  • 符合如下内积定义的线性空间V称为Euclid空间
    • 对称性 (a, b) = (b, a)
    • 线性
      • (a, kb) = k(a, b)
      • (a+c, b) = (a, b) + (c, b)
    • 非简并性
      • (a, a)>= 0 当且仅当 a=0 时 (a, a) = 0

Simplectic Space

defination 

  • 具有如下内积定义的线性空间W称为 辛空间
    • 反对称性 (a, b) = -(b, a)
    • 双线性
      • (a_1+a_2, b) = (a_1, b) + (a_2, b)
      • (a, b_1+b_2) = (a, b_1) + (a, b_2)
    • 非简并性
      • 若向量a 对于W中的任意向量b均匀(a, b)=0,则a=0

Simplectic Space 

  • 度量
    • 做功
    • 面积 体积
    • 流量
  • 辛内积
    • 单位辛矩阵 J_{2n}\equiv \begin{bmatrix} 0 &1_n \\ -1_n & 0 \end{bmatrix}
    • <a, b>\equiv (a,J_{2n}b)\equiv \sum_{i=1}^n (a_ib_{n+i}-a_{n+i}b_i)

单位辛矩阵的性质

\begin{matrix} J^2=-1\\ J^T=J^{-1}=-J\\ \end{matrix}

decision theorem of simplectic matrix 

Euclid space and Simplectic space 

  • 内积
  • 单位阵与单位辛矩阵
  • 正交与辛正交
  • 正交归一基于共轭正交归一基
  • 对称变换与Hamilton变换
    • 实对称矩阵的本征值为实数
    • 实对称矩阵不同本征值的本征向量正交
    • 实对称矩阵所有的本征向量构成一组正交归一基

Simplectic structure of Simplectic space 

q\equiv [q_1,q_2,...,q_f]^T\,\,\,p \equiv [p_1,p_2,...,p_f]^T

\dot{q_i}=\frac{\partial H}{\partial p_i}\,\,\, \dot{p_i}=\frac{\partial H}{\partial q_i}

z\equiv \begin{bmatrix} p\\ q \end{bmatrix}\,\,\, \dot{z}\equiv \begin{bmatrix} \dot{p}\\ \dot{q} \end{bmatrix}

J^{-1}\frac{\partial H}{\partial z}=\dot{z}

Simplectic structure of canonical transformation

basic concept

  1. J^{-1}\frac{\partial H}{\partial z}=\dot{z}
  2. y=Mz\,\,\, \dot{y}=M\dot{z}
  3. J^{-1}\frac{\partial H}{\partial z}=\dot{z}  J^{-1}\frac{\partial H}{\partial y}=\dot{y}
  4. simplectic transformation MJM^T=J

Properties of Simplectic transformation 

\left\{\begin{matrix} MJM^T=J\\ (det \,M)^2=1 \end{matrix}\right.

infinitesimal Simplectic matrix

define

  • 若B^TJ_{2n}+J_{2n}B=0，则该阵为无穷小辛阵

Dcision Theorem

  1. 若C为对称阵，当且仅当B = J_{2n}C 时，B为无穷小辛阵

Advantage of Simplectic algorithm 

  • 保证相空间体积不变
  • 保证力学守恒量守恒

River Chandler
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
Hamilton力学的辛算法PPT教案.pptx
10-20
Hamilton力学的辛算法PPT教案.pptx
Hamilton力学的辛算法.ppt
05-25
Hamilton力学的辛算法.ppt
SGMD(辛几何分解)+峭度值+能量熵+近似熵+模糊熵+排列熵+多尺度排列熵+样本熵​
2401_83250053的博客
03-02 1016
它计算了时间序列中相似模式出现的概率，但与近似熵不同的是，样本熵在计算过程中去除了自身匹配的情况。通过对时间序列进行多尺度分析（如小波变换或滑动窗口方法），可以计算每个尺度上的排列熵，并综合这些熵值以评估整个时间序列的复杂性。对于正态分布，峭度值为3。SGMD的核心思想是在信号的辛空间中寻找一个合适的辛变换，使得信号在该变换下呈现出最简单的结构。模糊熵是近似熵的一种扩展，允许在比较时间序列中的模式时引入一定的模糊性。各分解分量峭度值、能量熵、近似熵、模糊熵、排列熵、多尺度排列熵、样本熵，程序实用性高，
弹性力学极坐标辛体系Hamilton函数的守恒律
12-30
弹性力学极坐标辛体系Hamilton函数的守恒律，朱炳麒，卓家寿，由弹性力学的Hellinger-Reissner变分原理，将极坐标问题导向辛体系，对极坐标辛体系中Hamilton函数的守恒性进行了研究，在径向和环向两种�
2017辛算法课程.pdf
11-27
中科院数学所 唐贻发教授的哈密顿系统辛几何算法讲义 HamiltonX dZ dt = J−1∇H(Z), Z ∈R2n, J = O In −In O, (♠) "A{Gτ : Z →e Z "∂e Z ∂Z#T J"∂e Z ∂Z#= J, (♥) Euler¥:“(x, 1984) e Z = Z + τJ−1∇H e Z + Z 2 ! (♣) ·_"{" ♦ ºÆ [1] V.I. Arnold, Mathematical Methods of Classical Mechanics, Springer-Verlag, 1978. [2] x§ 5x'86(II), I§fi§1995. [3] x!§ 5MX"A{6§œ˘E §2004. ='§Springer, 2013. [4] E. Hairer, Ch. Lubich, G. Wanner, Geometric Numerical Integration, Springer, 2002.
一个求简单图中所有Hamilton回路的算法* (2005年)
06-05
从Hamilton回路的定义和图的邻接矩阵的定义入手，建立了图中的初级通路的关联关系。利用长度为k的初级通路及其关联关系逐步求长度为k+1的初级通路及其关联关系的方法。求得图的所有Hamilton回路。通过理论分析，说明该算法比已有的求图的所有的Hamilton回路的算法降低了算法的复杂度，为求解Hamilton回路问题提供了新思路。
基于FPGA的Hamiton方程--辛几何算法实现（全网唯一）
qq_37912811的博客
11-25 898
本实验做的是简单谐振子运动方程组的运算，会给出matlab代码以及相应的FPGA仿真截图。
自适应辛普森算法的代码实现
qq_41734963的博客
11-06 264
#include <iostream> using namespace std; typedef long long ll; const int inf = 0x3f3f3f3f; const double eps = 1e-6; double a, b, c, d, l, r; // 要进行积分的原函数 double f(double x) { return ; } //计算辛普森积分 double sim(double l, double r) { return
辛普森算法
weixin_30699443的博客
10-20 312
double func(double a, double b, double x) { double r = a * exp(- x*x) + b * sqrt(x); return r*r;//return的积分常数，例如这里f（x*x）dx|[0,h] } double integrate(double a, double b, double h) { unsig...
计算机科学中最重要的32个算法
热门推荐
掉下个小石头
08-06 1万+
英文原址：http://www.risc.jku.at/people/ckoutsch/stuff/e_algorithms.html奥地利符号计算研究所（Research Institute for Symbolic Computation，简称RISC）的Christoph Koutschan博士在自己的页面上发布了一篇文章，提到他做了一个调查，参与者大多数是计算机科学家，他请这些科学家投票选出
辛算法的发展历史与现状 (2004年)
05-27
Hamilton系统是用来描述无耗散的物理过程与物理现象的一种力学系统.辛几何算法是保结构算法中的一种,国内外学者在这一领域的研究,取得了丰硕的成果 .文中介绍针对 Hamilton系统的辛几何算法发展的简要历史、研究现状...
哈密顿插值c语言,哈密顿系统有限元方法的长时间性质研讨——冯康猜想的证明.pdf...
weixin_29560063的博客
05-19 192
哈密顿系统有限元方法的长时间性质研讨——冯康猜想的证明摘 要哈密顿系统是最重要的动力系统之一，它有两个重要的性质：辛结构和能量守恒．此外，在一定条件下还具有周期性．优秀的算法应当尽可能的保持系统原有的性质．本文将保能量的连续有限元法用于求解哈密顿系统，并研究它的长时间性质．主要结果及创新点如下．本文综合考察了辛算法和有限元方法的长时间性质．通过分析大量的数值结果提出了关于轨道，能量...
计算物理专题：薛定谔方程的辛算法
Chandler_river的博客
06-24 139
计算物理专题：薛定谔方程的辛算法
前置知识-辛 Runge-Kutta 方法
图灵的猫的博客
02-28 279
简单介绍辛几何和辛代数的基础上, 介绍相空间、Hamilton 系统及其两类辛格式: 线性 Hamilton 系统的中心 Euler 格式和一般的辛 Runge-Kutta 方法.
(转）网上摘抄：计算数学研究方向及网上资料
09-26 1744
计算数学目的为物理学和工程学作计算。主要研究方向包括： 数值泛函分析；连续计算复杂性理论；数值偏微与有限元；非线性数值代数及复动力系统； 非线性方程组的数值解法；数值逼近论；计算机模拟与信息处理等；工程问题数学建模与计算等等。 目前发展最好的方向已经与应用数学的CAGD方向合二为一。现在最热的方向应该是微分方程的数值求解、数值代数和流形学习，数值计算名校：西安交通大学、北京...
谈谈计算数学-
guo_shaokun的专栏
08-14 3976
谈谈计算数学        从计算数学的字面来看，应该与计算机有密切的联系，也强调了实践对于计算数学的重要性。也许Parlett教授的一段话能最好地说明这个问题：How could someone as brilliant as von Neumann think hard about a subject as mundane as triangular factoriz-ation o
欧氏空间内积定义_MP5：内积、外积、面积、Hermite内积、辛内积
weixin_32233909的博客
12-31 2502
我们发现，内积和外积都是和相对夹角相关，而和一对向量的整体刚体变换无关。本讲用一种特别的角度，从勾股定理出发，把两个向量长度构成的矩形面积分解称内积和外积两个部分。两个向量的夹角，在复数上可以表达为一个向量和另一个向量共轭的乘积，于是我们可以使用Hermite内积的结构。进一步，我们发现Hermite内积将内积和外积分解到实部和虚部，从而统一了二者。这样的结构在复Hilbert空间和量子力学中有重...
2024年北京高校数学建模校际联赛竞赛赛题浅析
最新发布
qq_33690821的博客
04-29 1566
这个问题涉及对气候变化下冰川动态的建模，主要目的是分析冰川末端位置的变化、表面物质平衡，以及预测冰储量的变化。线路质量评估：利用空间分析技术按公里评价线路质量，找出线路质量最好与最差的区间，这可能涉及到数据聚类或时间序列分析。归因分析：使用多元回归或机器学习方法来分析各种因素（如表面物质平衡、海水温度等）对冰川末端位置变化的贡献。案例分析：应用模型到实际订单数据，分析现有包装选择的效率，并提出改进建议。模拟和优化：通过模拟不同的包装策略，计算空间利用率和成本，以找到最佳方案。C题：订单商品包装问题。
广义Hamilton鲁棒控制器如何优化
08-13
### 回答1： 广义Hamilton鲁棒控制器可以通过使用多样化和自适应控制算法来优化性能，例如时域滤波器、自适应调节器、模糊控制、模型预测控制（MPC）和优化控制（OC）。此外，它还可以使用状态空间模型对系统进行模型预测，以优化系统的性能。 ### 回答2： 广义Hamilton鲁棒控制器是一种基于广义Hamilton函数的控制策略，它可以在系统具有不确定性和扰动的情况下实现控制目标。为了优化这种控制器，我们可以采取以下几种方法。 第一，通过选择合适的权重矩阵来优化广义Hamilton函数。权重矩阵可以用来平衡系统的性能和鲁棒性。通过调整权重矩阵的值，可以使广义Hamilton函数在不同的工况下具有更好的控制性能和鲁棒性。 第二，引入自适应控制技术来增强广义Hamilton鲁棒控制器的自适应性。自适应控制技术可以根据系统的实时状态和参数信息来调整控制器的参数，以便应对系统的变化和不确定性。通过引入自适应控制技术，可以使广义Hamilton鲁棒控制器更加适应不确定的系统环境，提高控制性能。 第三，利用模型预测控制技术来进一步提高广义Hamilton鲁棒控制器的性能。模型预测控制技术可以建立系统的数学模型，并基于该模型进行预测和优化控制。通过结合广义Hamilton函数和模型预测控制技术，可以实现更精确、更快速的控制响应，提高系统的控制精度和鲁棒性。 综上所述，通过选择合适的权重矩阵、引入自适应控制技术和结合模型预测控制技术，可以对广义Hamilton鲁棒控制器进行优化。这些方法将有助于提高控制性能、增强系统的鲁棒性和自适应性，使广义Hamilton鲁棒控制器在实际应用中更加有效和可靠。 ### 回答3： 广义Hamilton鲁棒控制器是一种在控制系统中使用的强鲁棒控制方法，旨在解决参数不确定、外部扰动和建模误差等问题。为了进一步优化广义Hamilton鲁棒控制器，可以采取以下措施。 首先，可以采用模糊PID控制器来改进广义Hamilton鲁棒控制器。模糊PID控制器结合了模糊逻辑和PID控制器的优点，能够更好地处理参数不确定和系统非线性。通过模糊化控制输入和输出，可以根据实际系统状态调整模糊规则库，从而优化控制性能。 其次，引入自适应控制方法优化广义Hamilton鲁棒控制器。自适应控制可以自动调整控制参数，以适应系统变化和不确定性。可以根据系统的特性和性能要求设计适应性算法，不断在线更新控制参数，实现自适应鲁棒控制。 此外，可以考虑采用鲁棒最优控制器优化广义Hamilton鲁棒控制器。鲁棒最优控制器综合了鲁棒控制和最优控制的优点，能够在系统不确定性下实现最优性能。通过最小化控制过程中的鲁棒度和性能指标，可以得到最优的控制输入，提高系统的鲁棒性和控制性能。 最后，引入先进的优化算法来进一步优化广义Hamilton鲁棒控制器。例如，粒子群算法、遗传算法和模拟退火算法等都可以用于寻找控制参数的最优解。通过这些优化算法，可以快速且准确地找到最优的控制参数，提高广义Hamilton鲁棒控制器的性能。 综上所述，通过采用模糊PID控制器、自适应控制方法、鲁棒最优控制器和先进的优化算法，可以对广义Hamilton鲁棒控制器进行优化，实现更好的控制性能和鲁棒性。这些方法可以根据具体系统的要求和特性进行选择和组合使用，以满足实际应用中的控制需求。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

River Chandler

谢谢,我会更努力学习工作的!!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值