你需要确定 功率密度 很对应的频率 才能看到相应的图像 比如,如下程序可以帮到你得到些灵感,data是原始数据,a是循环变量 具体如何循环根据data的组成而定
fs=samplerate;
N = length(data);
N2 = floor(N/2);
Y = fft(data(:,1,a));
Y(N2+1:N) = 0; %滤去方程中叠加的不含frequece shift的peak
P=real(Y).^2+imag(Y).^2; %功率密度
df=fs/length(P);
f_vec=df.*[1:length(P)/2];
figure(a*10)
plot(f_vec,P(1:length(P)/2))
相位 phi = unwrap(angle(X(left_idx:right_idx)));
只供参考 希望帮到你
df=fs/length(P);