matlab编制有限元法计算管道自振频率