利用辛方法及虚拟激励法计算轮轨力的功率谱密度函数(轨道为无砟轨道)