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/1000);
fs1=ns*fs/Ns;
ft1=nt*ft/Nt;
distance=fs1*1000/(1/(2*pi*R/1000));%若分析的是径向电磁力空间FFT,转换为空间阶次
xiebo=ft1/(P*zhuansu/60);
fn=zeros(Ns,Nt);%创建一个行数为空间采样数,列数为时间采样数
mag=zeros(Ns,Nt);%创建一个行数为空间采样数,列数为时间采样数的幅度矩阵
1