采用“四阶龙格—库塔法求解"求解公式为:
式中的h为时间步长。编写程序时需要将逐时的太阳辐射、气温、相对湿度、风速及空调负荷数据用txt文件的形式储存,程序运行时会自动调用这些数据。采用上式计算K2、K3时,各气象参数及空调负荷均采用前后两个时刻的算术平均值。
计算参数准备
水温模拟计算需要干球温度、相对湿度、辐射强度、风速、全年负荷等
相关数据可以到气象局网站查询,或者从DeST、鸿业EP等全年负荷模拟计算软件中导出相关气象参数。
matlab模型的建立
start.m	//启动函数 main	//计算函数文件夹 ff.m fun.m psie.m psig.m psis.m QL.m runge.m testrk.m Tetens.m public	//计算数据文件夹 avgt.txt //逐时的平均气温 fengsu.txt //全年逐时风速 fuhe.txt //全年逐时负荷 fushe.txt	//全年逐时辐射 qiwen.txt	//全年逐时干球气温 shidu.txt	全年逐时相对湿度
四阶龙格—库塔法的matlab实现
runge函数为四阶龙格—库塔法的计算函数,需要输入x0,x1,y0,h等四项参数。