分离变数法

news2024/12/23 22:48:18
  • 今天是2022年11月28号
    • 我的方程学的不太好,一些讲宇宙的,讲技术的,方程实在是看不懂
    • 很多方程的解与参数不可分割
  • 期末来了
    • 有的人回去了
    • 有的人要看光学了
    • 我呢,已经废物了,节日快乐,大家伙节日快乐啊!!
  • 分离变数法说简单,很简单,说难,肯定难
  • 话说,我看过一些视频,将方程与动画制作的。挺好的,挺好的

知识梳理

前置知识

正交函数系与傅里叶级数

\begin{Bmatrix} \varphi _n(x) \end{Bmatrix}(n=1,2,...)为定义在区间[a,b]上的函数系:

如果该函数系满足:\int_{a}^{b}\varphi_n(x)\varphi^*_m(x)dx=0,m\neq n

则称函数系\begin{Bmatrix} \varphi_n(x) \end{Bmatrix}[a,b]上的正交函数系

如果函数系还满足:\int_{a}^{b}\varphi_n(x)\varphi^*_n(x)dx=\int_{a}^{b}|\varphi_n(x)|^2dx=1

则称函数系\begin{Bmatrix} \varphi_n(x) \end{Bmatrix}[a,b]上的正交归一函数系

 正交归一关系可记为:

\int_{a}^{b}\varphi_n(x)\varphi_m^*(x)dx = \delta_{mn}=\left\{\begin{matrix} 0 &m\neq n \\ 1 & m=n \end{matrix}\right.

函数空间

  • 对建立于区间[a, b]上的函数空间,定义内积为

[\varphi_m,\varphi_n]=\int_{a}^{b}\varphi_n(x)\varphi^*_m(x)dx

  • 范数

||\varphi_n||=\sqrt{[\varphi_n,\varphi_n]}=\sqrt{\int_{a}^{b}\varphi_n(x)\varphi^*_n(x)dx}=\sqrt{\int_{a}^{b}|\varphi_n(x)|^2dx}

  • 正交(如上)
  • 归一(如上)
  • Schwartz 不等式

|\int_{a}^{b}\varphi_n(x)\varphi^*_m(x)dx|^2 \leq ||\varphi_m||^2||\varphi_n||^2

  • 正交归一函数系

\begin{Bmatrix} \frac{\varphi_n(x)}{||\varphi_n(x)||} \end{Bmatrix}(n=1,2,3,...)

正交函数系展开

[a,b]上的函数f(x),由正交函数系\begin{Bmatrix} \varphi_n(x) \end{Bmatrix}展开为:

f(x)=\sum_n C_n\varphi_n(x)

展开系数为:C_n=\frac{1}{||\varphi_n(x)||^2}\int_{a}^{b}f(x)\varphi_n^*(x)dx

傅里叶级数

  • 狄利克雷定理:略

f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_ncos\,n\omega x+b_nsin\,n\omega x)

 \left\{\begin{matrix} a_n=\frac{2}{T}\int_{-T/2}^{T/2}f(x)cos\,n\omega x\,dx\\ b_n=\frac{2}{T}\int_{-T/2}^{T/2}f(x)sin\,n\omega x\,dx \end{matrix}\right.

有限区间上函数的傅里叶展开

  • 定义在有限区间[0,l]上的函数f(x),可先将函数延拓为周期函数,再展为傅里叶级数

本征函数

  • 本征函数又称为固有函数,为特定的齐次边界条件所固有
  • 给出泛定方程和边界条件,相应的本征函数即已确定
  • 分离变量的基础, 是函数以正交的本征函数系为基函数
边界条件u|_{x=0}=u|_{x=l}=0u_x|_{x=0}=u_x|_{x=l}=0
本征值\lambda = (\frac{n\pi}{l})^2,n=1,2,...\lambda = (\frac{n\pi}{l})^2,n=0,1,2,...
本征函数\begin{Bmatrix} sin\,\frac{n\pi}{l}\,x \end{Bmatrix}\begin{Bmatrix} 1,cos\,\frac{\pi}{l}x,cos\frac{2\pi}{l}x,... \end{Bmatrix}
边界条件u|_{x=0}=u_x|_{x=l}=0u_x|_{x=0}=u|_{x=l}=0
本征值\lambda = [(n+\frac{1}{2})\frac{\pi}{l}]^2,n=0,1,2,...\lambda = [(n+\frac{1}{2})\frac{\pi}{l}]^2,n=0,1,2...
本征函数\begin{Bmatrix} sin\,\frac{(n+\frac{1}{2})\pi}{l}\,x \end{Bmatrix}\begin{Bmatrix} cos\,\frac{(n+\frac{1}{2})\pi}{l}\,x \end{Bmatrix}



齐次方程的分离变数法

  • 把偏微分方程分解成几个常微分方程,其中某些常微分方程带有附加条件,从而构成本征值问题

分离变数法的基本步骤

  1. 求满足齐次方程和齐次边界条件的变量分离形式的常微分方程和附加条件
  2. 求本征值和本征函数
  3. 得出分离变数形式的本征解
  4. 根据叠加原理求出一般解
  5. 利用初始条件求叠加系数,代入得定解问题的解

分离变数法的适用范围

  • 具有齐次线性泛定方程和齐次边界条件的定解问题
    1. 两端均为第一类齐次边界条件
    2. 两端均为第二类齐次边界条件
    3. 一端为第一类齐次边界条件,另一端为第二类齐次边界条件

非齐次边界条件问题

  • 应当让尽可能多的边界条件齐次化

非齐次振动方程和输运方程

傅里叶级数法

  • 直接求解非齐次方程的定解问题
    • 适用于:
      • 非齐次泛定方程可以分离,且有可分离的齐次边界条件
    • 试用解的技术形式由边界条件确定
    • 解决齐次边界条件的振动和输运问题

  • 齐次边界条件,分离变数法得到的解具有傅里叶级数表示形式

u(x,t)=\sum_nT_n(t)X_n(x)

  • 其中关于X的部分为傅里叶级数的基本函数族,由边界条件决定,其系数为t的函数
  • 将此试探解代入非齐次泛定方程,尝试分离出关于t的方程,结合 初值条件,求出关于t的解,最后带入试探解可得方程的解.
  • 操作步骤
    1. 根据初始条件把所求的解展开为傅里叶级数
    2. 分离出关于T的常微分方程并得出初始条件
    3. 求出关于T的常微分方程的通解初始条件
    4. 叠加求出定解问题的解

冲量定理法

  • 把非齐次方程的定解问题转化为齐次方程的定解问题后求解
  • 适用于:齐次的边界条件下的非齐次振动和输运问题

一般混合问题解的总结

例:

\left\{\begin{matrix} u_{tt}-a^2u_{xx}=f_1(x,t) &\,\,\,(0<x<l,t>0) \\ u|_{x=0}=u_1(t),u|_{x=l}=u_2(t) & \,\,\, (t \geq 0)\\ u|_{t=0}=\varphi_1(x),u_t|_{t=0}=\psi _1(x) &\,\,\,(0\leq x\leq l) \end{matrix}\right.

  1. 分离变量法中,需要齐次化的边界条件才能给出本征解
  2. 非齐次的泛定方程,可由齐次化原理,化为齐次问题求解

边界条件的齐次化

  • u(x,t)=v(x,t)+\omega(x,t)
  • v(x,t)为待定函数,并且满足非齐次边界

v|_{x=0}=u_1(t),v|_{x=l}=u_2(t)

实战

  • 周三更

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

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

相关文章

Spirng 痛苦源码学习(二)——手写spring大致总框架(一)

文章目录前言一、总体步骤如下1、spring 文件夹2、myProject 文件夹二、主要coding1、配置文件2、容器3、一些spring中的重要的注解4、项目中的使用5.重要的bean定义信息6、postProcessor重要&#xff0c;前置通知和后置主要项目的截图前言 本文主要基于spring的注解的方式完成…

Spring相关源码解读

框架1.ApplicationContext refresh的流程2.spring bean 的生命周期3.spring事务失效的几种场景以及原因4.springMVC执行流程5.一些注解&#xff08;1&#xff09;Configuration&#xff08;2&#xff09;Import&#xff08;3&#xff09;SpringBootApplication6.spring中有哪些…

BP神经网络详解,Python实现求解异或问题

BP神经网络 符号及其含义 nln_lnl​表示第lll层神经元的个数&#xff1b;f(⋅)f()f(⋅)表示神经元的激活函数&#xff1b;W(l)∈Rni∗ni−1W^{(l)}\in\mathbb R^{n_i*n_{i-1}}W(l)∈Rni​∗ni−1​表示第l−1l-1l−1层到第lll层的权重矩阵&#xff1b;wij(l)w_{ij}^{(l)}wij(l…

基于tensorflow的ResNet50V2网络识别动物

前言 之前很多人在&#xff0c;如何进行XXX的识别&#xff0c;对应的神经网络如何搭建。对应神经网络怎么搭建&#xff0c;我也是照本宣科&#xff0c;只能说看得懂而已&#xff0c;没有对这块进行深入的研究&#xff0c;但是现在tensorflow&#xff0c;paddle这些工具&#x…

长期稳定的项目—steam搬砖

大家好&#xff0c;我是阿阳 steam搬砖项目一直稳稳定定的进行着&#xff0c;有些朋友基本都观察了近2年 所以很多人问我公众号的项目是不能做了吗&#xff1f;怎么最近做新的去了&#xff1f;很明显这是几乎不可能的事情&#xff0c;steam做2年了&#xff0c;本公众号都能翻到…

这几个数据分析项目,让我看到了什么才叫专业!!

大家好&#xff0c;我是小一 新的一周又来了&#xff0c;从今天开始&#xff0c;会出一个新的系列《数分实验室》 实验室会介绍一些有内核、有科技的数据分析实战项目。 项目数据集、源代码都是公开的&#xff0c;非常适合想练手但是又没数据、没参考案例的同学 今天先热热…

ES的基础概念

1、ES是什么 Elasticsearch 是一个分布式可扩展的实时搜索和分析引擎,一个建立在全文搜索引擎 Apache Lucene(TM) 基础上的搜索引擎.当然 Elasticsearch 并不仅仅是 Lucene 那么简单&#xff0c;它不仅包括了全文搜索功能&#xff0c;还可以进行以下工作:分布式实时文件存储&am…

6-1分支限界法

6-1分支限界法 1.分支限界法与回溯法的不同 &#xff08;1&#xff09;求解目标: 回溯法的求解目标是找出解空间树中满足约束条件的所有解&#xff08;或一个最优解&#xff09;&#xff0c; 而分支限界法的求解目标则是找出满足约束条件的一个解&#xff08;或最优解&#x…

组织机器学习代码

组织机器学习代码 从note本转移到 Python 脚本时组织代码。 Intuition 有组织的代码就是有可读的、可重现的、健壮的代码。您的团队、经理&#xff0c;最重要的是&#xff0c;您未来的自己&#xff0c;将感谢您为组织工作付出的最初努力。在本课中&#xff0c;将讨论如何将代码…

pytest测试框架入门1

pytest单元测试框架 单元测试是指在软件开发当中&#xff0c;针对软件的最小单位&#xff08;函数&#xff0c;方法&#xff09;进行正确性的检查测试 单元测试框架主要做什么 测试发现&#xff1a;从多个文件里面找到我们的测试用例测试执行&#xff1a;按照一定的顺序和规则…

初学者指南: 使用NumPy数组进行图像处理

这里写自定义目录标题初学者指南: 使用NumPy数组进行图像处理1、加载图像2、裁剪图像3、分离颜色4、转换5、灰度转换6、图像分割结语初学者指南: 使用NumPy数组进行图像处理 由于图像也可以被视为由数组组成&#xff0c;因此我们也可以使用NumPy执行不同的图像处理任务。在本文…

【Lilishop商城】No2-6.确定软件架构搭建五(本篇包括定时任务xxl-job)

仅涉及后端&#xff0c;全部目录看顶部专栏&#xff0c;代码、文档、接口路径在&#xff1a; 【Lilishop商城】记录一下B2B2C商城系统学习笔记~_清晨敲代码的博客-CSDN博客 全篇只介绍重点架构逻辑&#xff0c;具体编写看源代码就行&#xff0c;读起来也不复杂~ 谨慎&#xf…

如何配置一台适合oc渲染器的电脑?

众所周知&#xff0c;Octane 是最流行的渲染引擎之一。此外&#xff0c;Octane 是一个 GPU 渲染引擎&#xff0c;它使用一种计算最终生成的图片的方法&#xff0c;试图达到照片般的真实感。Octane 是一种利用 GPU 技术的无偏渲染引擎&#xff0c;非常接近物理精度。一台好的 PC…

计算机组成原理习题课第三章-2(唐朔飞)

计算机组成原理习题课第三章-2&#xff08;唐朔飞&#xff09; ✨欢迎关注&#x1f5b1;点赞&#x1f380;收藏⭐留言✒ &#x1f52e;本文由京与旧铺原创&#xff0c;csdn首发&#xff01; &#x1f618;系列专栏&#xff1a;java学习 &#x1f4bb;首发时间&#xff1a;&…

天宇优配|平台助企“抱团出海” “小而美”中觅“先机”

天津华图轿车物流有限公司一批二手新能源车从连云港装船发往阿联酋迪拜。&#xff08;采访方针供图&#xff09; 最近&#xff0c;一笔100.8万美元的出口信誉稳妥保单融资借款&#xff0c;被划到了天津华图轿车物流有限公司的账户上。正值客户“下单”高峰期&#xff0c;这笔及…

Three.js实例详解___旋转的精灵女孩(附完整代码和资源)(一)

Three.js实例详解___旋转的精灵女孩(附完整代码和资源)&#xff08;一&#xff09; 本文目录&#xff1a; 一、【旋转的精灵女孩】案例运行效果 二、Three.js简介 三、Three.js代码正常运行显示条件 &#xff08;1&#xff09;不载入任何纹理贴图的网页 &#xff08;2&…

双十二蓝牙耳机啥牌子好?2022年度热销蓝牙耳机排名

这期双十二数码好物分享&#xff0c;工作室打算来跟大家说说蓝牙耳机这个话题&#xff0c;它已经成为出行必带的装备&#xff0c;上班族、学生党、游戏党都离不开蓝牙耳机。今年我们测评过数十款型号了&#xff0c;本期我们盘点了今年热销的蓝牙耳机排名&#xff0c;让大家直观…

【学习笔记】《Python深度学习》第五章:深度学习用于计算机视觉

文章目录1 卷积神经网络简介1.1 卷积运算1.2 最大池化运算2 在小型数据集上从头开始训练一个卷积神经网络2.1 下载数据2.2 构建网络2.3 数据预处理2.4 数据增强3 使用预训练的卷积神经网络3.1 特征提取3.2 微调模型3.3 小结4 卷积神经网络的可视化4.1 可视化中间激活4.2 可视化…

新手想开一个传奇该如何操作?开一个传奇必须掌握哪些知识要点

对于这个问题&#xff0c;近期问的人比较多&#xff0c;相比这也是热爱传奇这个游戏的朋友会问到的一个问题&#xff0c;因为喜欢玩这个游戏&#xff0c;也想要自己去开一个 经营一个 不管是电脑端也好 还是手机端也好&#xff0c;但是对于一些新手确实不知道该如何开始操作 从…

H3C opsf/rip/ftp/telent/nat/acl综合

实验拓扑 拓扑下载 https://sharewh2.xuexi365.com/share/84b85b32-acb7-4f62-a389-6188680a19f3?t3 图 1-1 注&#xff1a;如无特别说明&#xff0c;描述中的 R1 或 SW1 对应拓扑中设备名称末尾数字为 1 的设备&#xff0c;R2 或 SW2 对应拓扑中设备名称末尾数字为 2 的设备…