用MATLAB求解求解定谔代码我的薛定谔方程求解器代码 出版商 作者：Tsogbayar Tsednee（博士）联系人： 介绍： Matlab 代码 H_atom_DC_Stark_resonance.m 使用伪光谱方法结合复数吸收电位计算原子氢的 DC Stark 共振参数。 Matlab 代码 H2plus_eig_values_for_sigma_states.m 使用伪光谱方法计算 H2+ 离子的 sigma 状态的特征能。 要求： 任何版本的 Matlab 软件。 实现细节和运行 H_atom_DC_Stark_resonance.m 和 H2plus_eig_values_for_sigma_states.m 代码都使用代码 legDC2.m，该代码使用 Legendre-Gauss-Lobatto 节点和相应权重的搭配计算勒让德微分矩阵。 您可以下载它们并直接运行代码。 版权/许可 这些代码作为免费软件发布。 基本上，您有权使用它们并修改它们中的任何一个。 在以下位置找到 GNU 通用公共许可证：

Matlab代码sqrt BESOneo2 快速有效的BESO拓扑优化功能，用于受体积限制的2D最小符合性拓扑优化。 关于 二维双向进化拓扑优化（BESO），用于计算受给定体积约束的最佳最小顺应性结构。 BESO基于的代码，通过减少索引和有效的矩阵构造实现了显着的加速。 入门 可以使用[x, obj] = BESOneo2(400,200,0.3,0.02,3);来调用一个基本的悬臂示例[x, obj] = BESOneo2(400,200,0.3,0.02,3); BESOneo2.m中包含有关框架加固问题，L型支架和MBB梁的BESOneo2.m并在主要代码后进行注释。 这些可以替换C部分，也可以在此处定义自定义结构。 输入项 范围 描述 nx ny 在x和y方向上设计域大小 volfrac 体积约束，其中0 < volfrac <= 1 er 进化率 rmin 灵敏度滤波器半径解决网格相关性问题 特征 可以通过在代码部分C的pasS和pasV集中添加元素来实现非设计区域的指定。这将从优化过程中排除指定的元素，并将它们永久设置为实心或空元素。 示例：不可设计的区域 在C部分pasV

CFD2D是用于Linux的开源软件，用于求解任意二维域内的无量纲不可压缩Navier-Stokes方程（NSE），该二维域刻在具有Dirichlet和“不做任何事情”边界条件的单位正方形中。 空间离散化基于有限元方法（FEM），使用近似均匀的三角形网格。 提供了两种FE空间选择，分别是所谓的MINI元素和Taylor-Hood元素。 第一个元素由具有速度气泡的三次气泡函数的连续分段线性组成，第二个元素由连续分段线性组成。 在这两种情况下，都通过分段线性来近似压力场。 固定和时间相关的制度都支持。 提供了基本的绘图工具。 GMRES和CG迭代算法用于求解线性系统。 软件“ Triangle”用于网格生成。

二维稳态求解器，基于 SIMPLE 算法。 使用 10x10 交错网格。 采用盖子驱动的腔边界条件，但右手边的速度不为零。 如果您发现任何错误，请联系。 您可以根据自己的喜好添加任何可视化图。 轮廓非常粗糙（因为它是一个 10x10 的网格）。 网格可以做得更细，以获得更好的可视化。

这些函数求解周期性 LQ 状态反馈设计的离散时间周期 Riccati 方程 (DPRE)。 这些函数计算离散时间周期 Riccati 方程的唯一稳定解 X{k} 并返回状态反馈中的增益矩阵 K{k} u{k} = -K{k}x{k}，其中k = 1:P。 m文件“dpre”通过循环QZ或牛顿反向迭代法解决离散时间周期最优控制问题。 这些不是可用的最快方法，但效果很好。 mex 文件“dprex”通过周期性 QR（使用来自 matlab 内部 slicot 库的函数）或复杂的周期性 QC 方法（使用从 pqzschur 库中转换为 c 代码的 fortran 来解决离散时间周期性最优控制问题）。 mex文件的实现要快得多，但是需要编译mex文件，这可以通过运行make_dprex.m来完成。

Dirichlet问题的边界积分方程法 使用边界积分方程方法求解平面域上拉普拉斯方程的狄利克雷问题，其中边界是具有 C^2 参数化的平滑简单闭合曲线。

QMPython Python 中的量子力学和薛定谔方程求解器。 这些计算的起点是笔记本“Solving_the_Schrodinger_Equation_Numerically.ipynb”。 在本笔记本中，您将找到有关如何为无限方形井设置哈密顿量的说明。 很容易将其更改为任何形状势，限制在某个域 [xlow,xhigh] 上。 可以在 Harmonic Oscillator.ipynb 中找到以这种方式求解薛定谔方程的这种方法的示例用法，其中显示了 HO 的谐波状态如何随时间演变。 这个问题使用梯形算子的代数是完全可以解决的，并且与这里的数值解一致。 您可以使用这款笔记本玩一些游戏，例如在井的中心添加一个尖峰。 在 TimePropagation_of_WF.ipynb 中可以找到一个不同的解决方案，直接求解瞬态薛定谔方程。 这仍然使用哈密顿量，但现在直接解决时间相关性，步入时

为了协调Richards方程数值模型在精度、稳定性和计算效率之间的矛盾，对时间步长进行优化控制，发展了一种变步长的计算模型。采用van Genu chten模型，将描述土壤水分运动的Richards方程（为偏微分方程）在空间上半离散后得到常微分方程组，借助于该求解器求解。针对以含水率θ为变量的Richards方程进行了试验，常微分方程组求解采用CVODE求解器，并对土壤水分特征曲线和水量平衡进行了检验。结果表明，该模型具有较高的计算精度、求解效率和稳定性。

matlab发布代码LocalPeaksSearchMethod 程序简介：已在Matlab版本R2014a和R2018a上进行了测试 程序： 局部峰搜索方法（LPS方法）-色散方程求解方法。 适用于以下型号： 自由或流体加载的单层或双层弹性或粘弹性板； 自由或流体加载的单层或双层弹性或粘弹性圆柱壳。 开发人员：龚家园， 2009年在中国哈尔滨哈尔滨工程大学（HEU）最初编程； 2011年-2016年，中国青岛市中国科学院声学研究所修订 于湖北汽车技术大学（HUAT）修改，中国十堰，2018 版权所有（c）2009-2018，龚嘉元，隶属于HEU，IACAS和HUAT。 参考： 嘉苑宫水中的粘弹性阻尼材料涂覆的弹性板和圆柱壳中泄漏兰姆波的传播研究。 哈尔滨工程大学硕士学位论文，2010。 笔记： 圆柱壳的规范存在一些错误。 任何优化程序的人，请与开发者分享，非常感谢。 电子邮件： 示范方案 Fluid 1 [Fl, Va]: row1,c1 Material 2 [So]: rowvm, Evm0, ytavm, sigmavm Material 1 [So]: rowem, Eem0

matlab的欧拉方法代码可压缩的二维Euler方程求解器 该存储库包含MATLAB代码，用于使用磁通分解方法求解2D可压缩的Euler方程。 它目前使用Steger-Warming方案（1981）。 资源