代码是快速傅里叶算法,将时域信号转化成频域信号,但是出来的图像上没有曲线
ppgwatch=evalin('base','ppgwatch');
y=ppgwatch(:,2);
Fs=10;
T=1/Fs;
N=height(y);
t=(0:1:N-1)*T;
t=t';
y = table2array(ppgwatch(:, 2));
Y = fft(y);
Y=Y(1:N/2+1);
A=abs(Y);
f=(0:1:N/2)*Fs/N;
f=f';
A_adj=zeros(N/2+1,1);
A_adj(1)=A(1)/N;
A_adj(end)=A(end)/N;
A_adj(2:end-1)=2*A(end-1)/N;
subplot(2,1,1);
plot(f,A_adj);
xlabel('hz');
ylabdel=('幅值');