close all;clear all;clc;
data = [];
N=200;%采样点个数
fs = 50; %采样频率
t = 0:1/fs:(N-1)/fs;  %生成信号时间
data = [data csvread('D:\matlab_notebook\PD表\data\data_4s\196.csv')];%导入数据
% plot(data)
x=data(1:200,1);%提取第一列数据 1：N  (1:200,1)
% x=data(1:200,4);%提取第一列数据 1：N  (1:200,1)

%FFT
figure(1);subplot(211);plot(t,x);% 画出原始时域信号
% title('The combined acceleration');%score:2
ylabel('Acc'); %Acceleration
xlabel('Time/s');

signalFFT=abs(fft(x,N));%真实的幅值
Y=2*signalFFT/N;
f=(0:N/2)*(fs/N);
subplot(212);
plot(f,Y(1:N/2+1));
xlim([0 15]);
ylabel('Amp'); %Amplitude
xlabel('Frequency');
title('Spectrum');%Spectrum输入信号的频谱

%功率谱

%调用periodogram函数
[Pxx,f] = periodogram(x,rectwin(length(x)), N, fs);
subplot(212);
freq = 0:fs/length(x):fs/2;%频率刻度
train_y = 10*log10(Pxx)
for i= 1:size(train_y,2)
    if train_y(1,i)< -40
       train_y(1,i)==-40
    end
end
figure(2);plot(freq,train_y);
grid on;xlim([0 fs/2]);
title('Power spectral density');%用periodogram函数的周期图
xlabel('Frequency');
ylabel('PSD');%功率谱密度/db