%%數據處理
clc
clear
t=0:1/（10*20）:50;
q=xlsread（‘參數‘‘B3:B10003‘）;????%讀取路面數據
y1=xlsread（‘參數‘‘C3:C10003‘）;???%讀取被動懸架車身加速度數據
y2=xlsread（‘參數‘‘D3:D10003‘）;???%讀取被動懸架懸架動行程數據
y3=xlsread（‘參數‘‘E3:E10003‘）;???%讀取被動懸架輪胎動載荷數據
y4=xlsread（‘參數‘‘F3:F10003‘）;???%讀取半主動懸架車身加速度數據
y5=xlsread（‘參數‘‘G3:G10003‘）;???%讀取半主動懸架懸架動行程數據
y6=xlsread（‘參數‘‘H3:H10003‘）;???%讀取半主動懸架輪胎動載荷數據
q1=sqrt（sum（q.*q）/（size（q1）））
y1_RMS=sqrt（sum（y1.*y1）/（size（y11）））;?%車身加速度均方根值
y2_RMS=sqrt（sum（y2.*y2）/（size（y21）））;?%懸架動行程
y3_RMS=sqrt（sum（y3.*y3）/（size（y31）））;?%輪胎動載荷
y4_RMS=sqrt（sum（y4.*y4）/（size（y41）））;?%優化后的車身加速度均方根值
y5_RMS=sqrt（sum（y5.*y5）/（size（y51）））;?%優化后的懸架動行程
y6_RMS=sqrt（sum（y6.*y6）/（size（y61）））;?%優化后的輪胎動載荷
%%?
fs=10*20;??????????????????????????????%采樣頻率
Nt=length（q）;
dt=1/（10*20）;
tend=（Nt-1）*dt