add 传播模型matlab代码