多自由度,中心差分法,采用matlab软件编写的程序