data=dlmread（‘I:\Protokoll\Schienenverkehr\實驗一?車輛平穩性實驗\chuixiang352.txt‘‘‘）;%讀取實驗數據
z=data（:1）;??%?導入數據
N=length（z）;??%采樣點數
figure（1）
t=1/1000:1/1000:N/1000;?%建立采樣時間數組t
plot（tz）;??%畫出以時間t為橫坐標，垂向加速度為縱坐標的圖像
xlabel（‘時間?單位：s‘）;
ylabel（‘垂向加速度?單位：m/s^2‘）;
title（‘垂向加速度-時間圖像‘）;
grid?on;
axis（[0l/1000-2.52.5]）;
axis?fill
fs=1000;???%采樣頻率fs為1000Hz
fz=fft（zN）;????%傅里葉變換
a=abs（fz）;????%各頻率點模值
a=a*2/N;????%該點的模值除以N/2就是對應該頻率下的信號的幅度
n=0:N-1;????????
f=n*fs/N;???????%各點頻率第一個點為頻率為0的分量，故N值減1
figure（2）???????%畫出垂向加速度與頻率的關系圖像
plot（fa）;
xlabel（‘頻率?單位：Hz‘）;
yla