这个摩根平台通过调用摩根。米功能：

R=摩根（network_id，scenario_id，model_id，solver_id，reductor_id，varargin）

并且有四个强制参数：

• 网络id(一串)网络标识符
• 场景_id(一串)场景标识符
• 模型id(一串)模型标识符
• 解算器id(一串)解算器标识符

以及可选参数和其他可变长度参数列表：

• 还原器ID(细胞)还原器标识符数组（也可以为空）
• 变精蛋白变量参数列表，每个变量包含一个字符串（见下文）

所有可接受的附加(一串)参数描述如下：

• dt=X-用X覆盖配置中的时间步长（以秒为单位）
• tf=X-用X覆盖配置中的可调效率因子（正实数）
• y=X-使用10^X强制误差图的最小y刻度（默认值：-16）
• ord=X-用X（自然数）覆盖配置中的最大简化顺序
• pid=X-将自定义字符串标识符添加到绘图文件（默认值：“”）
• 笔记-不要测试降阶模型
• 契约-在一个图形中显示所有绘图

这个摩根。米函数返回结构对成员取决于参数。如果只计算降阶模型：

• .减速器(细胞)带有减速机名称的字符串数组
• .脱机(细胞)减速机离线时间数组

如果计算和测试降阶模型：

• .name（名称）(一串)实验的输出名称（保存的图和分数使用的名称）
• .减速器(细胞)包含还原器名称的字符串数组
• .订单(矢量)测试的减少订单
• .l0错误,.l1错误,.l2错误,.l8错误(细胞)每降阶平均误差数组
• .l0分,.l1芯,.l2分,.l8分(细胞)每降阶平均数数组MORscore公司秒
• .脱机(细胞)减速机离线时间数组
• .在线(细胞)减速机的平均相对在线时间数组
• .breven（布列文）(细胞)平均相对离线/在线盈亏平衡数数组

如果只运行模拟，对是矩阵，并包含离散输出轨迹，其维数为outputs-time-steps。

网络由（有向）图描述，通过边列表给出，它还指定了其边类型及其物理尺寸和属性。

• 网络必须至少有一个供应节点！
• 所有边界节点（供应或需求）必须正好通过一条边连接！
  • 可以插入短管来强制执行此操作。
• 供应节点的边缘必须远离它！
  • 因此，不能直接连接两个供应节点。
• 需求节点的边缘必须指向它！
  • 因此，没有两个需求节点可以直接连接。

列出了所有可用的网络数据集，其中包含网络编号

• 内部连接节点(编号0),
• 供应边界节点(新南威尔士州)、和
• 需求边界节点(nD公司).

• 钻石-钻石网络(n0=8，nS=1，nD=1，nC=0)
• 叉1-分叉式管道(n0=12，nS=1，nD=2，nC=0)
• 分叉2-分叉式管道(n0=12，nS=2，nD=1，nC=0)
• 比较测试-压缩机测试(n0=1，nS=1，nD=1，nC=1)
• 并列测试-平行管道测试(n0=2，nS=1，nD=1，nC=0)
• 手册DB16-三角形网络(n0=0，nS=1，nD=2，nC=0)

• MORGEN公司-小型网络(n0=27，nS=2，nD=4，nC=1)
• 阿泽J07-小型网络(n0=5，nS=1，nD=2，nC=1)
• GruHKetal13集团-小型网络(n0=11，nS=1，nD=8，nC=0)
• 九四-小型网络(n0=8，nS=1，nD=14，nC=0)
• GruJ金属14-中型网络(n0=45，nS=4，nD=2，nC=0)
• 液化天然气11-中型网络(n0=6，nS=3，nD=3，nC=2)
• 液化天然气24-中型网络(n0=14，nS=3，nD=5，nC=3)
• 液化天然气40-中型网络(n0=40，nS=3，nD=29，nC=6)
• 天然气Lib135-中型网络(n0=135，nS=3，nD=45，nC=29)
• 佩尔17a-中型网络(n0=41，nS=1，nD=15，nC=5)

• 管道-管道(n0=0，nS=1，nD=1，nC=0)
• Cha09号-管道(n0=0，nS=1，nD=1，nC=0)
• 杆S18-树(n0=6，nS=1，nD=4，nC=0)
• 盖67-树(n0=8，nS=1，nD=8，nC=0)
• H67a标段-管道(n0=0，nS=1，nD=1，nC=0)
• H67b标段-管道(n0=0，nS=1，nD=1，nC=0)
• H67c标段-树(n0=6，nS=2，nD=2，nC=2)
• H67d标段-树(n0=4，nS=2，nD=2，nC=1)

• AzePA19公司-葡萄牙(n0=0，nS=1，nD=1，nC=0)
• 第19页-西班牙(n0=6，nS=1，nD=5，nC=0)
• DeWS00公司-比利时(n0=20，nS=6，nD=9，nC=0)
• EkhD详细信息19-爱尔兰(n0=26，nS=3，nD=10，nC=0)
• 天然气Lib134-希腊(n0=134，nS=3，nD=45，nC=1)
• 气体Lib582-德国(n0=582，nS=31，nD=129，nC=5)
• 天然气Lib4197-德国(n0=4197，nS=11，nD=1009，nC=12)
• 科学网格编号-挪威(n0=44，nS=11，nD=9，nC=0)
• JinW公司-中国(n0=45，nS=5，nD=3，nC=38)

GasLib网络数据集源自：

M.Schmidt、D.Aßmann、R.Burraca、J.Humpola、I.Joorman、N.Kanelakis、，T.Koch、D.Oucherif、M.E.Pfetsch、L.Schewe、R.Schwarz、M.Sirvent：GasLib-天然气网络实例库；数据2（4）：2017年40日。

并根据CC-BY 3.0版，请参见：https://gaslib.zib.de

SciGrid网络数据集源自：

J.Dasenbrock、J.Diettrich、A.Pluta、W.Medjroubi：科学GRID_gas NO_Raw；Zenodo:10.5281/Zenodo.39852682020年。

并根据CC-BY 4.0版，请参见：https://www.gas.scigrid.de

网络编码在CSV公司文件扩展名为的文件.net文件。第一行是带有列、其含义和单位的描述。

#类型、标识输入、标识输出、管道长度[m]、管道直径[m]，高度差[m]和管道粗糙度[m]

第一行下方的每一行都包含一个带有列的边定义：

• 边缘类型(P（P）：管道，S公司：短管，C类：压缩机，V（V）：阀门）
• 开始节点标识符（正整数）
• 结束节点标识符（正整数）
• 管道长度[米]（正实数）
• 管道直径[米]（正实数）
• 管道高度差[米]（正实数）
• 管道粗糙度[米]（正实数）

因此，天然气网络的有向图表示为边缘列表，而除了边界节点外，边缘方向与流动方向不对应。注意，当前只有正整数可以用作开始和结束标识符。

解析的网络.net文件文件作为网络带构件的结构：

• 网络(结构)
  • .长度(矢量)管道长度
  • .倾斜(矢量)管道倾斜
  • .直径(矢量)管道直径
  • .粗糙度(矢量)管道粗糙度
  • .nomLen（标称长度）(矢量)每根管道长度
  • .A0(矩阵)供应节点减少的关联矩阵
  • 交流(矩阵)仅压缩机出口节点的关联矩阵
  • .B格式(矩阵)供应节点关联矩阵
  • .Bd（.Bd）(矩阵)需求节点关联矩阵
  • .Fc（预测值）(矢量)压缩机负荷矢量
  • .边缘(标量)边的数量
  • .n供应(标量)供应节点数量
  • .需求(标量)需求节点数
  • .内部(标量)内部节点数
  • .压缩机(标量)压缩机数量

场景数据集描述了天然气管网的边界值和外部不均匀性。供给和需求函数的瞬态行为仅通过标记变化以压缩形式表示为阶跃函数。每个网络都有一个培训场景(培训.ini)，对于降阶模型组件具有恒定的边界值。

场景编码在INI公司文件，带分机.ini文件。每行包含一个键值对，用于以下键：

• T0类-平均环境温度[C类]
• 卢比-平均比气体常数[焦耳/（千克*K）]
• t小时-时间范围[秒]
• 与-阀门设置[1]（管道分隔列表{0,1}）{UNDER CONSTRUCTION，当前被视为短管}
• 内容提供商-压缩机（输出）压力[酒吧]（管道分隔列表）
• 向上的-气源压力变化[酒吧]（分号分隔序列的管道分隔列表）
• 单位（uq）-需求流变化[千克/秒]（分号分隔系列的管道分隔列表）
• 美国犹他州-中更改的时间标记向上的和单位（uq）[秒]

解析的网络.ini文件文件如下：

• 脚本(结构)
  • .T0型(标量)全球平均温度
  • .卢比(标量)全球平均比气体常数
  • .tH（.tH）(标量)时间范围
  • .我们(矢量)稳态输入
  • .ut型(手柄)签名为u_t=ut（t）的函数句柄
  • .cp（.cp）(矢量)压缩机出口压力

模型对空间离散输入输出系统进行编码，形式如下：

E（p）x'（t）=A x（t）+B u（t）+F c_p+F（x（ty（t）=C x（t）

它由一个隐式非线性常微分方程组成，和（单向耦合代数）输出方程。

离散=模型（网络、配置）

• 网络(结构)解析的网络结构
• 配置(结构)配置结构

• 离散的(结构)（半）离散模型结构
  • .nP（当前位置）(标量)压力状态数
  • .nQ文件(标量)质量流量状态数
  • .nPorts（端口）(标量)端口数
  • .E型(手柄)质量矩阵函数句柄E_rtz=E（rtz）
  • .A型(矩阵)系统矩阵
  • .B类(矩阵)输入矩阵（模型边界节点）
  • .F型(矩阵)源矩阵（模拟压缩机动作）
  • .C类(矩阵)输出矩阵（边界节点处的传感器）
  • .f型(手柄)非线性向量场x=f（x，x，us，u，rtz）
  • .J（日本）(手柄)雅可比（Jacobian）x=J（x，x，u，rtz）
  • .双重(布尔)这只是一个成员（任何价值），如果它是一个双重模式！

• 模式（_M）-使用中点离散化的ODE模型
• 结束（_E）-使用端点离散化的ODE模型（端口哈密顿量）

• 争论X轴是在解算器（包装器）中计算的稳态。
• 争论x个非线性（f）和雅可比J型指与稳态的差值。这意味着xs+x生成实际状态。
• 只有组件E类,（f）和J型参数化。尤其，A类和B类不要依赖于参数。
• 压缩机只能在排放压力控制模式下运行。
• 争论无线电频率是产品卢比*T0*z0在解算器中形成。

解算器是模拟模型和场景轨迹的时间步进器。先决稳态初始值是根据场景的边界值计算得出的。

解决方案=求解器（离散、场景、配置）

• 离散的(结构)离散模型结构
• 脚本(结构)场景结构
• 配置(结构)配置结构

• 解决方案(结构)
  • t吨(矢量)时间步长矢量
  • u个(矩阵)离散输入-时间-步长轨迹
  • 年(矩阵)离散输出-时间-步长轨迹
  • steady_z0(标量)全球平均压缩系数
  • 稳定错误(标量)稳定状态错误
  • steady_iter1(标量)代数稳态迭代
  • steady_iter2(标量)微分稳态迭代
  • 运行时(标量)瞬态解算器运行时

• imex1型-一阶隐式显式解算器

• imex2-二阶隐式显式Runge-Kutta解算器

• 中国国家银行2-二阶曲柄-尼科尔森-亚当斯-巴什福斯解算器

• rk4号机组-四阶“经典”显式Runge-Kutta解算器（不稳定，仅用于测试）

• rk2型-二阶显式Runge-Kutta解算器（增加了双曲线稳定性）

• rk4hyp-四阶显式Runge-Kutta解算器（增加了双曲线稳定性）

• 通用的-二阶隐式自适应代码23Rosenbrock解算器

• 线性_导出-线性化和导出状态空间模型（包装imex1型解算器）

这个线性_导出“解算器”不是实际的解算器，但导出具有固定参数化的线性化：

E x’（t）=A x（t）+B u（t）+F，y（t）=C x（t），

选定网络场景的（E、A、B、C、F）状态空间模型，具有荷载向量F类共同描述压缩机、稳态和稳态效应：

F:=F*c_p+A*xs+B*us。

这些稀疏系统矩阵存储在.垫文件并命名网络_id-scenario_id-IySxOy.mat，其中年是边界端口数（输入和输出），以及x个是离散化的状态空间维度。

• 在内部，线性_导出调用imex1型solver返回解决方案。

• 对于线性_导出“解算器”模型_重力配置应为设置为没有人.

• 要获得（E、A、B、C）系统，负载向量F类可以连接到输入矩阵B类即：B'：=[B，F]，增加输入数量。

约简器计算降阶离散模型，旨在近似输入-输出（边界-感兴趣的数量）行为。

[proj，name]=还原器（求解器，离散，场景，配置）

• 解算器(手柄)解算器过程句柄
• 离散的(结构)离散模型结构
• 脚本(结构)场景结构
• 配置(结构)配置结构

• 项目(细胞)投影仪阵列{LP，RP；LQ，RQ}（双正交/斜向）或{LP；LQ}（LP；LQ）（正交）
• 名称(一串)减速器详细名称

这些结构式减压器近似于压力和质量流量分量单独（“结构化”缩写为“结构”）：

• 吊舱_r 结构。本征正交分解
• 编辑（_R）/eds_ro_l 结构。经验主导子空间
• 编辑_wx/eds_wx_l 结构。基于经验交叉矩阵的优势子空间
• 编辑_编辑/eds_wz_l 结构。经验非对称交叉广义优势子空间
• mpod_ro公司/mpod_ro_l 结构。改进的本征正交分解
• mpod_wx/mpod_wx_l 结构。改进的本征正交分解
• mpod_wz/mpod_wz_l 结构。改进的本征正交分解
• bpod_ro/bpod_ro_l 结构。经验平衡本征正交分解
• 电子束反辐射/ebt_ro_l 结构。经验平衡截断
• 电子束宽x/ebt_wx_l 结构。基于经验交叉矩阵的平衡截断
• ebt_wz（电子束宽度）/ebt_wz_l公司 结构。基于经验非对称交叉矩阵的平衡截断
• 袋鼠 结构。面向目标的本征正交分解
• ebg_ro（R）/ebg_ro_l公司 结构。经验平衡收益
• ebg_wx公司/ebg_wx_l 结构。基于经验的跨部门平衡收益
• ebg_wz公司/ebg_wz_l 结构。基于经验非对称交叉矩阵的平衡收益
• dmd_反应堆 结构。动态模式分解Galerkin

所有利用可观测性信息的减速器都有两种变体。默认情况下，使用非线性变量（无后缀）。这个_我后缀表示减速机的“线性”变体，假设有双系统可用。虽然这两种方法都可以适用，模式（_M）和结束（_E）、型号、，理论建议只对端口哈密顿量使用线性变量结束（_E）.

测试定义了一个实验，该实验被实现为一个脚本，其文件名为由前缀和被测试网络的名称组成。两种类型的实验是当前实施：

• 模拟实验的前缀为模拟_并执行摩根。米针对固定网络和脚本。

• 模型降阶实验的前缀是莫尔_并执行停尸房。米模型、求解器和约简器组合的功能使用时的网络和场景培训.ini计算场景降阶模型。

请注意，测试只能从摩根基本目录，运行后设置脚本，或手动添加测验文件夹到路径：

addpath（“测试”）；

通过运行设置脚本，其中列出了测验文件夹。

降阶模型保存在z_roms（_R）文件夹（或指定的文件夹由摩根（_R）配置条目）。

通过存储投影仪并对关联：文件名中的网络、模型和还原器，如下所示：

网络-模型-电抗器-pid.rom

具有pid控制器可以通过可选参数配置自定义标识符。

要加载简化模型，请提供已保存简化模型的文件名而不是减速机标识符。

绘图和MOR得分计算公式摩根。米存储在z_插槽文件夹（或由摩根广场配置条目）。

这个MORscore公司是基准指数测量模型简化误差图上方的面积，并绘制该图。该分数共同评估最小尺寸和最大精度的模型缩减目标。不稳定降阶模型计算为相对误差为1.

这个摩根平台假定配置INI公司-文件在名为morgen.ini，如果未找到硬编码默认值使用。

• 摩根广场(字符串)存储绘图的文件夹，默认：z_地块

• 摩根（_R）(字符串)用于存储简化订单模型的文件夹，默认：z_roms（_R）

• 网络dt(正浮点数)请求的时间步长（秒），默认值：60

• 网络_vmax(正浮点数)气体的最大速度（米/秒），默认值：20

• 网络cfl(正浮点数)空间离散化的目标CFL常量，默认值：0.5

• 模型_调整(正浮点数)可调效率系数缩放摩擦项，默认值：1

• 模型雷诺数(正浮点数)估计雷诺数，默认值：1000000

• 模型摩擦(字符串)摩擦系数模型，从中选择霍弗,尼古拉兹,阿尔舒尔,希弗里森,下午1025,伊格特，默认值：希弗里森

• 模型压缩性(字符串)压缩因子模型，选择：理想的,数字视频录像机,aga88型,教皇，默认值：aga88型

• 模型_压缩(字符串)压缩性参考：稳定的,正常的，默认值：稳定的

• 模型_重力(字符串)重力计算：没有人,静止的,动态，默认值：静止的

• steady_maxiter_lin(正整数)优化稳态估计的最小范数迭代次数，默认值：20

• steady_maxiter_non(正整数)细化稳态估计的时间步长迭代次数，默认值：1000

• steady_maxerror（稳定_最大误差）(正浮点)优化稳态的最大误差，默认值：1e-6

• steady_Tc（稳定_Tc）(浮子)临界温度（摄氏度），默认值：-82.595

• steady_pc(浮子)临界压力（单位：巴），默认值：45.988

• steady_pn(浮子)正常压力（单位：巴），默认值：1.01325

• 解算器松弛(浮点输入（0,1])IMEX解算器松弛，默认：1

• 解算器_rk2类型(正整数)二阶双曲龙格-库塔级数5,6,7,8,9,10,11,12，默认值：11

• 求解器rk4类型(字符串)四阶双曲Runge-Kutta型MeaR99a公司,测量99b卢比,TseS05型，默认值：MeaR99a公司

• T0_分钟(浮子)最低环境温度（摄氏度），默认值：0

• T0最大值(浮子)最高环境温度（摄氏度），默认值：20

• Rs_最小值(浮子)最小比气体常数[J/（kg*K）]，默认值：500

• Rs_最大值(浮子)最大比气体常数，单位为[J/（kg*K）]，默认值：600

• 莫尔激励(字符串)通用培训输入类型，从以下选项中选择：冲动,步,随机二进制,白噪声，默认值：步

• 最大mor(正整数)最大降序，默认值：200

• mor_参数(布尔值)使用参数化模型降阶，从中选择真的,假，默认值：真的

• mor网格(正整数)稀疏参数网格细化级别，默认值：1

• 评估_参数(浮子)参数规范：1,2,Inf公司，默认值：2

• 评估跳过(正整数)评估每个第n个降阶模型，默认值：三

• 评估最大值(正整数)要评估的最大简化顺序，默认值：200（使用Inf公司最大可能）

• 评估_参数化(布尔值)参数化降阶模型评估：真的,假，默认值：真的

• 评估_测试(正整数)测试参数数量，默认值：5

• 评估收益(布尔值)使用增益校正：真的,假，默认值：真的

在内部，配置存储在结构的结构中，如下所示：

• 配置(结构)
  • .网络(结构)成员：.日期,.vmax（最大值）,.cfl（立方英尺）
  • .型号(结构)成员：.调谐,.雷诺,.摩擦,.重力
  • .稳定(结构)成员：.日期,.maxiter_lin（最大值）,.maxiter_non（最大值）,.maxerror（最大误差）,.Tc（变速箱）,.pc版,.pn码,.压缩性
  • .solver解决方案(结构)成员：.日期,.放松,.rk2类型,.rk4类型
  • .mor(结构)成员：.rom最大值,.参数化,.solver解决方案,.激励,.T0_分钟,.T0最大值,.Rs_最小值,.Rs_最大值,.pgrid
  • .评估(结构)成员：.参数化,.p测试,.T0_分钟,.T0最大值,.Rs_最小值,.Rs_最大值,.skip（跳过）,最大值,.pnorm文件,.增益

所有输入温度，即：

• 摩根伊尼配置
• XXXXXX.ini版脚本

所有输出温度均为摄氏度.内部所有温度都在开尔文.

• xml2net.xsl转换天然气 .net文件XML网络定义到MORGEN公司-相容的.net文件通过XSLTproc的CSV网络定义：

xsltproc-o GasLib-X.csv xml2net.xsl GasLib-X.xml

• csv2net。米转换SciGRID_气体 .csv文件CSV网络定义到MORGEN公司-相容的.net文件CSV网络定义：

csv2net（'X_Y_PipeSegments.csv'，'myX_Y'）

• json2net。米转换数学能量 .json文件网络定义到MORGEN公司-相容的.net文件CSV网络定义：

json2net（'X.json'，'myX'）

• vf2千克。米将体积流量转换为kg/s（默认气体密度为0.7）

vf2kgs（值、体积单位、时间单位、密度）

• psi2巴。米将psi转换为bar

b=磅/平方英寸巴（p）

• 兰斯奇。米生成给定网络的随机场景，隐式定义训练场景

randscen（网络，场景名称）

• cmp_摩擦。米比较摩擦系数

cmp_摩擦（Re、D、k）

• cmp压缩性。米比较压缩系数

cmp_压缩性（p，T，pc，Tc）