function [f, y] = zfft(x, fi, fa, fs)
% x为採集的数据
% fi为分析的起始频率
% fa为分析的截止频率
% fs为採集数据的採样频率
% f为输出的频率序列
% y为输出的幅值序列(实数）
R=210.2%
P=50;
zhuansu=220;
Nt=51;
Ns=51;
ns=0:Ns-1;
nt=0:Nt-1;
% ns=0:1.3207:132;
% nt=0:109000:5450000;
fs=1/1.3207;
ft=1/(109000/1
