cir模型matlab代码