一维河道水动力学建模.zip

MATLAB编程，对一维常规河道进行水动力学建模，对圣维南方程组进行Preissmann离散，给定边界条件后，可计算各时刻水位流量值。

%--------------------------------------------------
%2019.11.12最后修改
%模拟工况：一维不规则河道水动力学建模
%边界条件：上游水位为-5m，下游水位为-15m，初始水位-10m，初始流量为0。
%程序思路：先将断面数据处理成数据表，每隔0.5m确定大断面的河宽、过水面积、湿周和水力半径
%         然后对圣维南方程组进行Preissmann离散，最后求解五对角矩阵即得各断面水位和流量。
%--------------------------------------------------

%--------------------------------------------------
%定义变量及初始化

clcclear;

%导入断面原始测量数据，共115个断面，每个断面给出了若干个起点距和高程对
load（‘data.mat‘）; 

%设置大断面水位间距dhT为计算时间段（s），dt为时间步长（s）ac:连续方程离散后的系数，
%Zu、Zd分别为上下游水位，q0为河道初始流量，n0为糙率
dh=0.5;T=24*3600;dt=300;a=1;c=1;Zu=-5;Z0=-10;Zd=-15;q0=0;g=9.81;n0=0.025;

%以下为矩阵分配存储空间
%data（11）为断面总数，n为河段总数，N为断面数*2即ZiQi变量总数，t为时间分段数，x为ZiQi解向量
n=data（11）-1;N=2*data（11）;t=ceil（T/dt）;x=zeros（1N）;

%B：五对角矩阵方程的常数项，IJKMO:五对角矩阵A的五列
B=zeros（1N）;I=zeros（1N）;J=zeros（1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      210492  2018-12-04 23:00  data.mat
     文件        9467  2019-11-13 11:53  General_river.m
     文件      728799  2018-12-03 10:39  澄通河段地形资料.txt
     文件         172  2019-11-13 11:57  说明.txt