用于计算雷诺数（Re）和相对粗糙度系数（epsilon）给定值下管道内摩擦系数的MATLAB代码。语法：f=colebrook（Re，epsilon）示例1:Single Re，Single epsilon Re=1e5；ε=1e-4；f=colebrook（Re，ε）示例2：多个Re，单个εRe=5000:1000:100000；ε=1e-4；f=科尔布鲁克（Re，ε）；图（Re，f）例3：单个Re，多个εRe=1e5；ε=linspace（1e-4,1e-1100）；f=科尔布鲁克（Re，ε）；plot（epsilon，f）示例4:Multiple Re，Multiple epsilon Re=logspace（4,8100）；ε=linspace（1e-4,1e-1100）；[RE，EPSILON]=meshgrid（RE，EPSILON）；F=科尔布鲁克（RE，ε）；surf（RE，EPSILON，F）参考文献：[1]Colebrook，C.F.，&；怀特，C.M.（1937年）。粗糙管道中流体摩擦的实验。伦敦皇家学会会刊。A辑-数学和物理科学，161（906），367-381。[2] 科尔布鲁克，C.（1939年）。管道中的紊流，特别是指光滑管道和粗糙管道之间的过渡区规律。土木工程师学会杂志，11（4），133-156。

MATLAB类计算密度、粘度和体积模量随温度的变化，从1摄氏度到99摄氏度。

MATLAB代码绘制Moody图，显示管道中不同粗糙度系数的摩擦系数和雷诺数之间的关系。

FIGURE2PDF将当前的MATLAB图形导出为PDF文件。此函数用于为研究论文准备高质量的图形文件。

GETFIGDATA从FIG文件中提取xyz数据；[x，y，z]=GETFIGDATA（文件名）；[x，y，z]=GETFIGDATA（文件名，n，类型）；如果图形包含多个轴，则必须指定要使用的轴编号（n）。默认情况下，使用第一个轴。如果图形包含多个绘图类型（直线、散点等），则必须指定所需的类型，否则默认情况下将选择“直线”。可以使用GETFIGDATA检查绘图类型，无输出参数。每个数据系列都是返回单元格数组的一个元素：[x{i}，y{i}，z{i}]必须用于访问第i个数据系列。例如：myfigfile='simulation_output.fig'；GETFIGDATA（myfigfile）[x，y]=GETFIGDATA（myfigfile）；情节（x{1}，y{1}）作者：Ildeberto de los Santos Ruizidelossantos@ittg.edu.mxSee还有OPENFIG，SAVEFIG。

epanet_地图使用地理轴和卫星地图在新的MATLAB图形中绘制epanet网络。语法：默认情况下，epanet_映射（输入文件名、utm_区域、语言、特殊_节点），utm_区域为'14Q'，语言为'ENG'（英语），特殊_节点为空。示例：epanet_地图（'Madrid.inp'，'30T'）epanet_地图（'Guanajuto.inp'，'14Q'，'SPA'，[15,30,60]）作者：Ildeberto de los Santos Ruizidelossantos@ittg.edu.mx需要互联网连接和MATLAB R2018b或更高版本。

[H，delta]=hurwitz（p）返回多项式p的hurwitz矩阵。可选的输出参数delta包含所有主要子变量。示例：symskp=[1，K，2，5]；[H，delta]=hurwitz（p）

非线性动态模型的线性近似矩阵。在x=x0，u=u0时，从xdot=f（x，u）转换为xdot=a*x+B*u。

根据雷诺数和相对粗糙度系数，采用不同的方法（隐式和显式）计算达西-魏斯巴赫摩擦系数。要使用的计算方法由输入参数指定。此外，还提供了一个子程序，用于根据管道不同运行条件下的雷诺数和摩擦系数测量值估算粗糙度系数。http://doi.org/10.5281/zenodo.1481992

采用梯形积分法和恒定时间步长模拟线性动态系统。系统模型必须表示为：xdot=A*x+B*uy=C*x+D*u

恒定输入下线性动态系统的稳态（最终）响应。系统必须稳定。系统模型必须表示为：xdot=A*x+B*u y=C*x+D*u

非线性动力系统的平衡点。xdot=f（x，u）的平衡搜索从x=x0，u=u0开始。