基于第二种压力分裂,通过 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
(数值分析课程设计)Matlab求解常微分方程初值问题 欧拉方法 梯形方法 龙格-库塔方法
2023-04-12 14:05:05 398KB Matlab 常微分方程 欧拉方法 梯形方法
1
欧拉公式求长期率的matlab代码FlightGear-Space-Shuttle模拟 基于MATLAB / Simulink和FlightGear的航天飞机基于真实推力数据进行发射模拟。 这是对NASA航天飞机升空80公里的模拟。 它是用MATLAB和FlightGear编写的 为了查看模拟,您将需要: MATLAB(最好是R2014或更高版本)。 FlightGear(首选3.0版),可以在slx文件的“ FlightGear预配置6DOF动画”中进行更改。 Simulink Microsoft Excel 要运行模拟,请导航到安装文件的文件夹。 打开MATLAB运行文件“ runfirst.m”。 这会将轨迹从“ outputs.xlsx”文件扫描到Simulink中。原始的仿真代码在文件“ shuttle_model.m”中。 您可以在此处更改方程式和代码。 这些方程式是从《工程学生的轨道力学》(Howard Curtis)一书的第11章获得的。 它们可以通过通用的Euler方法(步长为0.2s)来解决。
2023-03-23 10:15:59 191KB 系统开源
1
四元数欧拉角算法库,角速度陀螺仪地磁计融合算法,包含姿态转换源码。输入加速度信号由算法库转换生成四元数和欧拉角。
2023-03-15 15:51:58 1.1MB 四元数 C语言 加速度算法 欧拉角
1
用WPF编写的欧拉角转换方向余弦矩阵的软件工具,不用手工计算
2023-03-10 17:08:29 211KB 欧拉角 方向余弦矩阵 C# WPF
1
实时欧拉视频放大 这是一个关于欧拉视频放大的程序该程序提供实时运动或彩色放大和频率分析还为视频提供运动或彩色放大该算法主要基于本文 平台和图书馆: 适用于 Linux 的 openCV2.4.2 适用于 Linux 的 Qt5.1.0
2023-02-24 11:30:59 233KB Makefile
1
matlab的欧拉方法代码 圈定种群周期的模型(在Matlab中) 该存储库包含Matlab代码,用于(大多数)描述旅居种群动态的ODE模型进行数值积分: predation文件夹中的捕食者-猎物模型。 这些包括Gilg等人的模型。 (2003)针对格陵兰岛特雷尔岛(Traill Island)和更简单变体的旅鼠种群。 图尔钦和巴茨利(Turchin&Batzli,2001)的论文集中在阿拉斯加Point Barrow的解释波动之后, vegetation文件夹中的消费者-植物模型。 host-parasite模型。 注意尽管建模人员可以方便地在单独的模型中实现这些机制,但是当然不能保证它们在现实中不会全部相互作用。 它们之间可能会相互作用,也可能与舞台结构相互作用,这在这里的大多数模型中都被忽略了。 实用笔记 大部分代码最初是由FrédéricBarraquand(CNRS;然后在特罗姆瑟大学)于2011年至2014年编写的,是根据约翰·安德烈·亨登(特罗姆瑟大学)以前的特雷尔岛模型编写的代码。 我们最初考虑将这些用于尚未实现的出版物(尚未发布),因此公开发布。 我们对特雷尔岛模型的
2023-02-19 13:06:42 7.02MB 系统开源
1
用改进欧拉方法求 dy/dx = 2/3xy^(-2),x∈[0,1], y(0) = 1 的数值解(取h = 0.1),并将计算结果与准确解y = ∛(1+x^2 )进行比较:
2023-02-17 01:03:38 3.46MB 欧拉公式
1
SpinCalc 是一个整合的 matlab 函数,它将在包含的 4 种类型之间转换任何旋转数据。 还将在 2 种不同的欧拉角集类型之间进行转换。 可以输入多个方向。 对于N个方向: DCM ===> 3x3xN多维数组EA### ===> Nx3 矩阵欧拉向量 ===> Nx4 矩阵四元数 ===> Nx4 矩阵 输入包括错误检查标志,当欧拉角接近奇点或当适当的值偏离统一时,该标志将发出警告。 因不正确的 DCM 等而发出的致命错误。 *****注意用户***** 我有很多关于转换为欧拉角集的问题。 将数据转换为欧拉角时,您必须确保要平移的方向不接近奇点。 奇异欧拉集是无法按特定旋转顺序唯一转换为 3 个变量的方向。 单数集如下: 类型 1 转数:123 - 132 - 213 - 231 - 321 - 312 如果第二个旋转角度为 -90 或 90 度,则为单数。 类型 2
2023-02-17 00:53:42 6KB matlab
1
欧拉公式求长期率的matlab代码hydro_examples 各种流体力学技术的简单一维示例 这是一些简单的python代码(加上一些Fortran代码)的集合,这些代码演示了流体力学代码中使用的一些基本技术。 所有代码都是独立的-没有相互依赖关系。 这些代码与讲义一起位于: 并使用pyro2代码: advection/ advection.py :具有多种限制器的一维二阶线性对流求解器。 fdadvect_implicit.py :使用周期边界条件的一维一阶一阶隐式有限差分线性对流求解器。 fdadvect.py :使用迎风微分的一维一阶显式有限差分线性对流求解器。 fv_mol.py :一维线法二阶精确对流求解器。 Fortran/ : advect.f90 :二阶线性对流的Fortran实现。 此版本执行分段常数,分段线性和分段抛物线(PPM)重建。 basic-numerics orbit-converge.py (和orbit.py ):演示了各种ODE积分方法对太阳绕地球旋转问题的收敛性。 burgers/ burgers.py :不粘的Burgers方程的一维二阶求解器
2023-02-13 20:45:06 346KB 系统开源
1