我们对使用 Leap-frog 方法获得一维波动方程的解感兴趣。 并且边界条件是周期性的。 然而,初始条件是T(x,0)=sin(10*pi*x); 0<= x<= 0.1 =0; 0.1<= x<= 1 u = 0.25
2023-05-22 21:07:31 2KB matlab
1
Matlab 求解偏微分的代码PyCheb 这是一个使用谱方法求解 ODE 的 Python 包 背景 微分方程用于描述状态和过程的现象。 这些问题的解解释了它们的模式,因此人们渴望寻求这些方程的解来描述状态和预测未来。 常微分方程 (ODE)是一种微分方程,其中包含一个(作为方程的变量)自变量(函数的)及其导数的函数。 求解 ODE 相对容易,但对科学家和工程师很有用。 这就是为什么我们对它感兴趣并制作这样一个 Python 包来解决它。 光谱方法 谱方法是应用数学中用于数值求解微分方程的一类技术。 这个想法是将微分方程的解写为某个“基函数”的总和(例如,作为正弦和的傅立叶级数),然后选择总和中的系数以满足微分任何给定精度的方程。 谱方法可用于求解常微分方程 (ODE)、偏微分方程 (PDE) 和涉及微分方程的特征值问题。 与传统的 ODE 求解方法相比,在目标函数足够平滑的情况下,谱方法自然具有收敛速度超快的优势。 有关光谱方法的更多详细信息,请查看 。 它列出了用于理解谱方法和 MATLAB 项目Chebfun 的参考书目,我们将在后面专门讨论。 相关作品 2002年,由牛津大学
2023-05-15 18:53:57 61KB 系统开源
1
学习偏微分方程的经典教材 evans编的 英文版
2023-05-09 09:52:16 4.67MB 偏微分方程 evans
1
求解非线性扩散方程,该方程可以线性化,如 Richtmyer & Morton [1] 中的一般非线性扩散方程所示。 该方法是将 pde 线性化并应用 Crank-Nicolson 隐式有限差分格式以数值方式求解方程。 Matlab 运行命令-------------------------- 类型: IsoFreeSurfaceSolver 解决 pde: ------------------- \frac{\partial h}{\partial t}=\frac{1}{12}\frac{\partial^2 h^4}{\partial x^2} pde 可应用于在重力作用下在水平基底上扩散的等温粘性流体流动 - Huppert [2]。 请注意,PDE 已无量纲化。 初始条件: t=0: h = (1 - x^2)_{+} + 10^-6(有预湿膜) 考虑到 x = 0
2023-05-05 10:22:48 394KB matlab
1
5.05.Multigrid1D 一维泊松方程的V周期多重网格方法
2023-04-22 13:29:30 3KB Fortran
1
作者: [苏]菲利波夫 出版社: 上海科学技术出版社 译者: 孙广成 / 张德厚 出版年: 1981-1 页数: 147 定价: 0.46 装帧: 平装 统一书号: 13119-945目录 · · · · · · 前言 目录 §1.等斜线、曲线族微分方程的建立 §2.可分离变量的方程 §3.几何与物理问题 §4.齐次方程 §5.一阶线性方程 §6.全微分方程、积分因子 §7.解的存在性与唯一性 §8.导数未解出的方程 §9.各类一阶方程 §10.可降阶的方程 §11.常系数线性方程 §12.变系数线性方程 §13.边值问题 §14.常系数线性方程组 §15.稳定性 §16.奇点 §17.相平面 §18.解对于初始条件和参数的依赖性、微分方程的近似解 §19.非线性方程组 §20.一阶偏微分方程 答案 指数函数与对数函数表
2023-04-21 07:36:50 10.61MB 孙广成 张德厚 常微分方程 习题集
1
MATLAB牛顿法求解非线性方程组 部分源码 function Newton() x0=[0.1;0.5]; x1=x0-inv(myJacobi(x0))*myfun(x0); while norm(x1-x0)>1e-3 x0=x1; x1=x0-inv(myJacobi(x0))*myfun(x0); end x1 end
2023-04-17 19:34:30 664B MATLAB 牛顿法 求解 非线性方程组
1
详细介绍了数值模拟方法、如有限差分、有限元、谱方法、谱元法。书中附带代码下载地址和视频课程地址
2023-04-17 10:11:37 17.07MB 有限元 有限差分 谱元法 微分方程
1
基于第二种压力分裂,通过 Liou 和 Steffen AUSM 通量向量分裂 (FVS) 技术求解一维欧拉方程。 添加的源项通过将参数 alpha 设置为 0(平面 1D)、1(圆柱轴对称)或 2(球对称)来解释圆柱和球对称流。 边界条件是透射类型,初始数据适用于黎曼问题 (RP)。 为了求解添加的源,在 Ut 上进行 PDE/ODE 拆分,如下所示: PDE: Ut + Fx = 0 对于 U(x,t0) 产生 Up(预测) 颂: dU/dt = S(U) for U(t0) = Up 产生 U(x,t+dt) 数字的详细信息可以在 E. Toro 的书中找到,也可以在 Liou 和 Steffen 的原始论文中找到。 适用于大多数基准测试,除了 Toro 书第 8 章中描述的测试 3,它失败了。 要确定时间步长,请运行粗网格情况,计算完整解的特征值 L1 = ua 和 L3
2023-04-14 08:40:24 2KB matlab
1
第一类边界条件下扩散方程稳态近似分析,王成善,穆小静,反应器内层状液、液界面两侧内扩散,团块、气泡、粉粒或液滴内的扩散都是有限长度区间上的传质问题;获得有限长度区间上的扩散方
2023-04-12 23:03:51 405KB 首发论文
1