%--------------------------------------------------
%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