谱减法原理及程序实现