matlab编写的半功率法求阻尼因子+++++++++++++++++++++++++++++++