详细平稳解

1.详细平衡

定义:一个在高斯白噪声激励下的动力学系统在状态空间中如果用如下运动方程描述:
d d t X j \frac{d}{dt}\mathbf{X}_{j} dtdXj= f j ( X ) f_{j}(\mathbf{X}) fj(X)+ ∑ l = 1 m g j l ( X ) W l ( t ) \sum_{l=1}^{m}g_{jl}(\mathbf{X})W_{l}(t) l=1mgjl(X)Wl(t),j=1,2,…n
它一般不属于平稳势类。为满足FPK方程,可将一阶导数矩分成如下两个部分:
a j ( x ) = a j R ( x ) + a j I ( x ) a_{j}(\mathbf{x})=a_{j}^{R}(\mathbf{x})+a_{j}^{I}(\mathbf{x}) aj(x)=ajR(x)+ajI(x)
其中, a j R ( x ) a_{j}^{R}(\mathbf{x}) ajR(x)为可逆分量, a j I ( x ) a_{j}^{I}(\mathbf{x}) ajI(x)为不可逆分量。不可逆部分和阻尼力相应,可逆部分与惯性力和恢复力相应 (对于拉格朗日提法和哈密顿提法也是如此)

补充:维纳过程和高斯白噪声之间的关系

考虑方程 d X ( t ) d t = W ( t ) , X ( 0 ) = 0 \frac{dX(t)}{dt}=W(t),X(0)=0 dtdX(t)=W(t),X(0)=0其中, W ( t ) W(t) W(t)是谱密度为 K K K的高斯白噪声,即 E [ W ( t ) ] = 0 , E [ W ( t ) W ( t + τ ) ] = 2 π K δ ( τ ) . E[W(t)]=0,E[W(t)W(t+\tau)]=2\pi K\delta(\tau). E[W(t)]=0E[W(t)W(t+τ)]=2π(τ).
按维纳过程的定义, X ( t ) X(t) X(t)是维纳过程,则 d B ( t ) d t = W ( t ) . ( ∗ ) \frac{dB(t)}{dt}=W(t).(*) dtdB(t)=W(t).()而维纳过程的强度 σ 2 \sigma^{2} σ2与高斯白噪声的谱密度之间的关系为 σ 2 = 2 π K . \sigma^{2}=2\pi K. σ2=2πK.
这里的(*)式只是一种形式上的关系,因为维纳过程 B ( t ) B(t) B(t) L 2 L_{2} L2意义上不可微。
作为最简单的马尔可夫扩散过程,维纳过程 B ( t ) B(t) B(t)可以通过随机微分方程用于构造其他马尔可夫过程,一个标量马尔可夫扩散过程可由下式产生: d X ( t ) = m ( X , t ) d t + σ ( X , t ) d B ( t ) dX(t)=m(X,t)dt+\sigma(X,t)dB(t) dX(t)=m(X,t)dt+σ(X,t)dB(t),其中 B ( t ) B(t) B(t)是单位维纳过程,即
E [ B ( t 1 ) B ( t 2 ) ] = m i n ( t 1 , t 2 ) , E [ d B ( t 1 ) d B ( t 2 ) ] = { 0 , t 1 ≠ t 2 d t , t 1 = t 2 = t E[B(t_{1})B(t_{2})]=min(t_{1},t_{2}),E[dB(t_{1})dB(t_{2})]=\left\{\begin{matrix} 0, &t_{1}\neq t_{2} \\ dt, &t_{1}=t_{2}=t \end{matrix}\right. E[B(t1)B(t2)]=min(t1,t2),E[dB(t1)dB(t2)]={0,dt,t1=t2t1=t2=t

2.以下举例说明:

2.1外激单自由度系统

考虑系统 X ¨ + h ( Λ ) X ˙ + u ( X ) = W ( t ) \mathbf{\ddot{X} }+h(\Lambda)\mathbf{\dot{X}}+u(\mathbf{X})=W(t) X¨+h(Λ)X˙+u(X)=W(t),式中 u ( X ) u(\mathbf{X}) u(X)是恢复力, W ( t ) W(t) W(t)是谱密度为 K K K的高斯白噪声, Λ \Lambda Λ是系统的总能量,也就是 Λ = 1 2 X ˙ 2 + ∫ 0 X u ( z ) d z . \Lambda=\frac{1}{2}\mathbf{\dot{X}}^{2}+ \int_{0}^{X}u(z)dz. Λ=21X˙2+0Xu(z)dz.
自然地,我们可以知道阻尼力为 X ˙ \mathbf{\dot{X}} X˙对应部分,惯性力和恢复力为 u ( X ) u(\mathbf{X}) u(X)

首先,我们得到系统对应的伊藤方程(令 X 1 = X \mathbf{X_{1}}=\mathbf{X} X1=X X 2 = X ˙ \mathbf{X_{2}}=\mathbf{\dot{X}} X2=X˙):
d X 1 = X 2 d t d X 2 = − [ h ( Λ ) X 2 + u ( X 1 ) ] d t + 2 π K d B ( t ) . d\mathbf{X_{1}}=\mathbf{X_{2}}dt\\ d\mathbf{X_{2}}=-[h(\Lambda)\mathbf{X_{2}}+u(\mathbf{X_{1}})]dt+\sqrt{2\pi K}dB(t). dX1=X2dtdX2=[h(Λ)X2+u(X1)]dt+2πK dB(t).
对于这里 d B ( t ) = W ( t ) d t dB(t)=W(t)dt dB(t)=W(t)dt,由于把普通随机微分方程转化成了伊藤随机微分方程,故按对于的规则,得到上式; σ ( X , t ) d B ( t ) = 2 π K d B ( t ) \sigma(X,t)dB(t)=\sqrt{2\pi K}dB(t) σ(X,t)dB(t)=2πK dB(t).

由此得到一、二阶导数矩:
a 1 = x 2 , a 2 = − h ( λ ) x 2 − u ( x 1 ) , b 11 = b 12 = b 21 = 0 , b 22 = 2 π K a_{1}=x_{2},a_{2}=-h(\lambda)x_{2}-u(x_{1}),b_{11}=b_{12}=b_{21}=0,b_{22}=2\pi K a1=x2,a2=h(λ)x2u(x1),b11=b12=b21=0,b22=2πK
其中可逆部分:
a 1 R = x 2 , a 2 R = − u ( x 1 ) a_{1}^{R}=x_{2},a_{2}^{R}=-u(x_{1}) a1R=x2,a2R=u(x1)
不可逆部分:
a 1 I = 0 , a 2 I = − h ( λ ) x 2 a_{1}^{I}=0,a_{2}^{I}=-h(\lambda)x_{2} a1I=0,a2I=h(λ)x2
将上式带入 a j I = 1 2 ∑ k = 1 n [ ∂ ∂ x k b j k ( x ) − b j k ( x ) ∂ ϕ ∂ x k ] , ∑ j = 1 n ∂ ∂ x j a j R ( x ) = ∑ j = 1 n a j R ( x ) ∂ ϕ ∂ x j . a_{j}^{I}=\frac{1}{2}\sum_{k=1}^{n}[\frac{\partial }{\partial x_{k}}b_{jk}(x)-b_{jk}(x)\frac{\partial \phi}{\partial x_{k}}],\\ \sum_{j=1}^{n}\frac{\partial }{\partial x_{j}}a_{j}^{R}(x)=\sum_{j=1}^{n}a_{j}^{R}(x)\frac{\partial \phi}{\partial x_{j}}. ajI=21k=1n[xkbjk(x)bjk(x)xkϕ],j=1nxjajR(x)=j=1najR(x)xjϕ.
(若系统属于详细平衡类,则存在一个 ϕ ( x ) \phi(x) ϕ(x)满足所有这些方程。)
得到, π K ∂ ϕ ∂ x 2 = − h ( λ ) x 2 , ( 1 ) x 2 ∂ ϕ ∂ x 1 = u ( x 1 ) ∂ ϕ ∂ x 2 . ( 2 ) \pi K \frac{\partial\phi}{\partial x_{2}}=-h(\lambda)x_{2},(1)\\ x_{2}\frac{\partial\phi}{\partial x_{1}}=u(x_{1})\frac{\partial \phi}{\partial x_{2}}.(2) πKx2ϕ=h(λ)x2,(1)x2x1ϕ=u(x1)x2ϕ.(2)
解(1)式:
请添加图片描述
带入(2)式,发现 g ( x 1 ) g(x_{1}) g(x1)必为常数,又对于一个系统属于平稳势类则其平稳概率密度表示为: p ( x ) = C e x p [ − ϕ ( x ) ] p(x)=Cexp[-\phi(x)] p(x)=Cexp[ϕ(x)]其中,C为归一化常数,于是
p ( x 1 , x 2 ) = C e x p [ − 1 π K ∫ 0 λ h ( z ) d z ] , λ = 1 2 x 2 2 + ∫ 0 x 1 u ( z ) d z . p(x_{1},x_{2})=Cexp[-\frac{1}{\pi K}\int_{0}^{\lambda} h(z)dz],\lambda=\frac{1}{2}x_{2}^{2}+\int_{0}^{x_{1}}u(z)dz. p(x1,x2)=Cexp[πK10λh(z)dz],λ=21x22+0x1u(z)dz.

故对于受外激作用的单自由度系统属于详细平衡类.
特别地,对于线性阻尼力情形: h ( Λ ) X ˙ = α X ˙ h(\Lambda)\dot{X}=\alpha \dot{X} h(Λ)X˙=αX˙, p ( x 1 , x 2 ) = C e x p { − α π K [ ∫ 0 x 1 u ( z ) d z + 1 2 x 2 2 ] } . p(x_{1},x_{2})=Cexp\{-\frac{\alpha}{\pi K}[\int_{0}^{x_{1}}u(z)dz+\frac{1}{2}x_{2}^{2}]\}. p(x1,x2)=Cexp{πKα[0x1u(z)dz+21x22]}.

2.2同受外激和参激的单自由度系统

补充:外激(施加在系统上的外力或外部扰动)参激(系统内部的反馈信号或参数输入)
X ¨ + ( α + β X 2 ) X ˙ + ω 0 2 X = X W 1 ( t ) + W 2 ( t ) \ddot{X}+(\alpha+\beta X^{2})\dot{X}+\omega_{0}^{2}X=XW_{1}(t)+W_{2}(t) X¨+(α+βX2)X˙+ω02X=XW1(t)+W2(t)
其中 W 1 ( t ) W_{1}(t) W1(t) W 2 ( t ) W_{2}(t) W2(t)是谱密度分别为 K 11 K_{11} K11 K 22 K_{22} K22的独立高斯白噪声,以 X 1 X_{1} X1 X X X X 2 X_{2} X2 X ˙ \dot{X} X˙,则对于的伊藤方程:
d X 1 = X 2 d t d X 2 = − [ ( α + β X 1 2 ) X 2 + ω 0 2 X 1 ] d t + 2 π ( K 11 X 1 2 + K 22 ) d B ( t ) dX_{1}=X_{2}dt\\ dX_{2}=-[(\alpha+\beta X_{1}^{2})X_{2}+\omega_{0}^{2}X_{1}]dt+\sqrt{2 \pi (K_{11}X_{1}^{2}+K_{22})}dB(t) dX1=X2dtdX2=[(α+βX12)X2+ω02X1]dt+2π(K11X12+K22) dB(t)

请添加图片描述
请添加图片描述

可以看到这里系统属于详细平衡解需要满足 α β = K 22 K 11 \frac{\alpha}{\beta}=\frac{K_{22}}{K_{11}} βα=K11K22

ps:可能会补充更新,仅供自己学习使用,对于其他形式的系统可作类似推导,都是书中原例。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/602830.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Open CASCADE学习|入门Hello world

目录 1、新建项目 2、写代码 3、配置 3.1配置头文件 3.2配置静态库文件 3.3配置动态库文件 4、编译运行 1、新建项目 新建一个Win32控制台应用程序,取名为HelloWorld,如下图所示: 2、写代码 测试所用的代码如下: // Use T…

Redis:原理速成+项目实战——Redis实战6(封装缓存工具(高级写法)缓存总结)

👨‍🎓作者简介:一位大四、研0学生,正在努力准备大四暑假的实习 🌌上期文章:Redis:原理速成项目实战——Redis实战5(互斥锁、逻辑过期解决缓存击穿问题) 📚订…

C# Entity Framework 中不同的数据的加载方式

延迟加载 延迟加载是指在访问导航属性时,Entity Framework 会自动查询数据库并加载相关数据。这种方式在我们需要访问导航属性时比较方便,因为我们无需手动加载相关数据,而且只会在需要时才会进行查询,从而减少了不必要的开销。但…

Python基础-07(for循环、range()函数)

文章目录 前言一、for循环1.for循环结构2.参数 end(使其输出时变为横向) 二、range()函数1.range(常数)2.range(起始值,结束值)3.range(起始值,结束值,步长)4.例子 总结 前言 此章介绍循环结构中最常用的循环&#xf…

Redisson与SQL乐观锁:实现接口幂等性的终极指南与实战演示

Redisson与SQL乐观锁:实现接口幂等性的终极指南与实战演示 Redisson与SQL乐观锁:实现接口幂等性的终极指南与实战演示 接口幂等性.md

497 蓝桥杯 成绩分析 简单

497 蓝桥杯 成绩分析 简单 //C风格解法1&#xff0c;*max_element&#xff08;&#xff09;与*min_element&#xff08;&#xff09;求最值 //时间复杂度O(n)&#xff0c;通过率100% #include <bits/stdc.h> using namespace std;using ll long long; const int N 1e4 …

java基于VUE3+SSM框架的在线宠物商城+vue论文

摘 要 信息数据从传统到当代&#xff0c;是一直在变革当中&#xff0c;突如其来的互联网让传统的信息管理看到了革命性的曙光&#xff0c;因为传统信息管理从时效性&#xff0c;还是安全性&#xff0c;还是可操作性等各个方面来讲&#xff0c;遇到了互联网时代才发现能补上自古…

【langchain】在单个文档知识源的上下文中使用langchain对GPT4All运行查询

In the previous post, Running GPT4All On a Mac Using Python langchain in a Jupyter Notebook, 我发布了一个简单的演练&#xff0c;让GPT4All使用langchain在2015年年中的16GB Macbook Pro上本地运行。在这篇文章中&#xff0c;我将提供一个简单的食谱&#xff0c;展示我们…

【Docker基础三】Docker安装Redis

下载镜像 根据自己需要下载指定版本镜像&#xff0c;所有版本看这&#xff1a;Index of /releases/ (redis.io) 或 https://hub.docker.com/_/redis # 下载指定版本redis镜像 docker pull redis:7.2.0 # 查看镜像是否下载成功 docker images 创建挂载目录 # 宿主机上创建挂…

element-ui table height 属性导致界面卡死

问题: 项目上&#xff0c;有个点击按钮弹出抽屉的交互, 此时界面卡死 原因分析: 一些场景下(父组件使用动态单位/弹窗、抽屉中使用), element-ui 的 table 会循环计算高度值, 导致界面卡死 github 上的一些 issues 和解决方案: Issues ElemeFE/element GitHub 官方讲是升…

bootstrap教程

bootstrap教程 大家好&#xff0c;我是免费搭建查券返利机器人赚佣金就用微赚淘客系统3.0的小编&#xff0c;也是冬天不穿秋裤&#xff0c;天冷也要风度的程序猿&#xff01;今天&#xff0c;让我们一同进入前端开发的世界&#xff0c;探索一款备受欢迎的前端框架——Bootstra…

修改安卓apk设置为安卓主屏幕(launcher)

修改安卓apk 将apk可以设置安卓主屏幕 原理&#xff1a; 将打包好的apk文件进行拆包增加配置文件在重新编译回apk包 需要得相关文件下载 解包 apktool :https://pan.baidu.com/s/1oyCIYak_MHDJCvDbHj_qEA?pwd5j2xdex2jar&#xff1a;https://pan.baidu.com/s/1Nc-0vppVd0G…

2024年学习计划

2024-2-29号完成 机器视觉基础知识学习&#xff0c;并可以处理视觉工作中的需求。 2024-3月份学习SCARA机械手应用开发SCARA机器人-埃斯顿自动化 - ESTUN 2024-4月份继续学习python 好了&#xff0c;今年可以完成这三个目标就满足了 好好学习&#xff0c;天天向上。每天进步…

【数据库系统概论】数据库并发控制机制——并发控制的主要技术之封锁(Locking)

系统文章目录 数据库的四个基本概念&#xff1a;数据、数据库、数据库管理系统和数据库系统 数据库系统的三级模式和二级映射 数据库系统外部的体系结构 数据模型 关系数据库中的关系操作 SQL是什么&#xff1f;它有什么特点&#xff1f; 数据定义之基本表的定义/创建、修改和…

Vue CLI组件通信

目录 一、组件通信简介1.什么是组件通信&#xff1f;2.组件之间如何通信3.组件关系分类4.通信解决方案5.父子通信流程6.父向子通信代码示例7.子向父通信代码示例8.总结 二、props1.Props 定义2.Props 作用3.特点4.代码演示 三、props校验1.思考2.作用3.语法4.代码演示 四、prop…

嵌入式培训机构四个月实训课程笔记(完整版)-Linux系统编程第四天-Linux管道练习题(物联技术666)

更多配套资料CSDN地址:点赞+关注,功德无量。更多配套资料,欢迎私信。 物联技术666_嵌入式C语言开发,嵌入式硬件,嵌入式培训笔记-CSDN博客物联技术666擅长嵌入式C语言开发,嵌入式硬件,嵌入式培训笔记,等方面的知识,物联技术666关注机器学习,arm开发,物联网,嵌入式硬件,单片机…

爬虫网易易盾滑块案例:某乎

声明&#xff1a; 该文章为学习使用&#xff0c;严禁用于商业用途和非法用途&#xff0c;违者后果自负&#xff0c;由此产生的一切后果均与作者无关 一、滑块初步分析 js运行 atob(‘aHR0cHM6Ly93d3cuemhpaHUuY29tL3NpZ25pbg’) 拿到网址&#xff0c;浏览器打开网站&#xff0…

Java十种经典排序算法详解与应用

数组的排序 前言 排序概念 排序是将一组数据&#xff0c;依据指定的顺序进行排列的过程。 排序是算法中的一部分&#xff0c;也叫排序算法。算法处理数据&#xff0c;而数据的处理最好是要找到他们的规律&#xff0c;这个规律中有很大一部分就是要进行排序&#xff0c;所以需…

python总结-装饰器

装饰器 装饰器解决日志问题&#xff08;分三个版本&#xff09;多个装饰器带参数的装饰器wraps装饰器内置装饰器property装饰器staticmethod装饰器classmethod装饰器 类装饰器 概念 装饰器来自 Decorator 的直译。什么叫装饰&#xff0c;就是装点、提供一些额 外的功能。在 pyt…

一个可能的网址服务器证书自动续期自动化实现脚本方案

需求背景&#xff1a; 目标: 您希望为您的网站启用HTTPS&#xff0c;以保护通信安全&#xff0c;并希望这个过程是免费的。 证书类型: 您需要获取SSL/TLS证书&#xff0c;并且希望证书能够自动续期&#xff0c;以确保网站不会因证书过期而停机。 服务器兼容性: 您希望生成的证…