function [jx,SNR,result]=myesprit(x,M,k,dt) %实现《基于快速TLS_ESPRIT的间谐波检测算法》张滨生1.1 1.2 1.3 % input x:原始数据% input M:时间窗宽度% input dt:采样间隔% output jx:拟合数据% output Ad：振幅% output Qd：相位% output ad：衰减因子% output fd：振荡频率% output result=[Ad,Qd,ad,fd]