A novel matrix model of district heating network was proposed.
Thermal conditions can be solved analytically and efficiently by matrix operation.
The number of uncertain parameters was reduced by an aggregation operation.
Model parameters were calibrated by an optimized method using genetic algorithms.
After calibration, predicted temperatures fitted well with measured data.