二自由度模型振动模型频域分析matlab程序

%%数据处理
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