版本：R2016a

利用SCE-UA调用新安江模型进行率定时，显示参数太多，求大神指导。

clear;

clc;

spath=cd

P0=(textread([spath,'\matlabxingshanjiangyu.txt']))';

P=P0(2,:);%降水量(mm)

EI0=(textread([spath,'\matlabxingshanzhengfa.txt']))';

EI=EI0(2,:);%蒸发皿蒸发量(mm)

cszt=(textread([spath,'\matlabxingshanchushizhuangtai.txt']));

T=cszt(2:cszt(2):cszt(1)); %计算时间para(1)，计算时段para(2)

%初始状态

F=cszt(3);%流域面积(km2)

WU=cszt(4);WL=cszt(5);WD=cszt(6);%上下深层流域蓄水量(mm)

FR=cszt(7);%产流面积比重

S=cszt(8);%自由水蓄水量(mm)

QRSS0=cszt(9);%壤中流初始流量(m3/s)

QRG0=cszt(10);%地下水初始流量(m3/s)

DT=cszt(2);%计算时段

para=(textread([spath,'\matlabxingshanpara.txt']));

%基本参数

WUM=para(1);WLM=para(2);WDM=para(3);%上下深层流域蓄水容量(mm)

K=para(4);%蒸发皿系数

C=para(5);%深层蒸散发系数

B=para(6);%蓄水容量曲线指数

IMP=para(7);%不透水面积比重

SM=para(8);%自由水蓄水库容量(mm)

EX=para(9);%自由水蓄水容量曲线指数

KG=para(10);%地下水从自由水库的出流系数

KSS=para(11);%壤中流出流系数

KKG=para(12);%地下水退水系数

KKSS=para(13);%壤中流退水系数

UH = [para(14),para(15),para(16)]; %无因次单位线时段

F0=(textread([spath,'\matlabxingshanjingliu.txt']))';

f0=F0(2,:);%jingliu

x0=[WUM,WLM,WDM,K,C,B,IMP,SM,EX,KG,KSS,KKG,KKSS,UH(1),UH(2),UH(3)];

f0=xinanjiang1(P,EI,cszt,x0)

2018-12-14 10:41 上传