clear
clc

D=load('Z_Beschleunigung.mat');
D = D.Z_vibration_Zeit;
Z_a_vibration=D(:,2)./10;
Z_vibration_Zeit(:,1)=D(:,1);
Z_vibration_Zeit(:,2)=Z_a_vibration;%%Z_vibration_Zeit 真实数据
Ts=1e-4; %sample time 
Fs=1/Ts; %sample frequenz(Hz)
x=D(:,2);
t1=D(:,1);
Ndata=length(x); %sample punkte
t=[0:1/Fs:(Ndata-1)/Fs];
N=2.^10;
N=1000;
n=0:N-1; %plus
f=n*Fs/N;
subplot(2,2,1) %plus
plot(t1,Z_a_vibration,'g','Linewidth',1.0);
title('signal_original');
xlabel('Zeit');
ylabel('Beschleunigung');
Y_z = fft(Z_a_vibration,N); %fft translation for z-achse fft变化
Ayy_z =(abs(Y_z));
subplot(2,2,2); %plus
plot(f,Ayy_z,'g');
title('FFT-Modulus');
xlabel('sample-punkt');
ylabel('FFT-Modulus');
Ayy_z=Ayy_z/(N/2);
Ayy_z(1)=Ayy_z(1)/2;
subplot(2,2,3);
plot(f(1:N/2),Ayy_z(1:N/2),'g'); 
title('Amplitude-frequenz');
xlabel('Frequenz');
ylabel('Betrag');
%% phase
Pyy=[1:N/2];
for i=1:N/2
    Pyy(i)=phase(Y_z(i));   % Phase
    Pyy(i)=Pyy(i)*180/pi; % as Winkel
end
%figure;
subplot(2,2,4);
plot(f(1:N/2),Pyy(1:N/2));
title('Phase-frequency');

grid on;