FFT.rar_fft

clc; set(0,'defaultfigurecolor','w') %横向力 [SH]=xlsread('100直线参数.xlsx','横向力'); t01=SH(:,1); h01=SH(:,2);h1 = max(abs(h01)) h02=SH(:,3);h2 = max(abs(h02)) h03=SH(:,4);h3 = max(abs(h03)) h04=SH(:,5);h4 = max(abs(h04)) figure(1) plot(t01,h01,t01,h02,t01,h03,t01,h04) grid on xlabel('时间/t') ylabel('横向力/N') title('横向力') legend('一位轮对','二位轮对','三位轮对','四位轮对') %标注曲线名称 hold on %垂向力 [SC]=xlsread('100直线参数.xlsx','垂向力'); t02=SC(:,1); c01=SC(:,2);c1 = max(abs(c01)) c02=SC(:,3);c2 = max(abs(c02)) c03=SC(:,4);c3 = max(abs(c03)) c04=SC(:,5);c4 = max(abs(c04)) figure(2) plot(t02,c01,t02,c02,t02,c03,t02,c04) grid on xlabel('时间/t') ylabel('垂向力/N') title('垂向力') legend('一位轮对','二位轮对','三位轮对','四位轮对') hold on %轮重减载率 [SL]=xlsread('100直线参数.xlsx','轮重减载率'); t03=SL(:,1); l01=SL(:,2);l1 = max(abs(l01)) l02=SL(:,3);l2 = max(abs(l02)) l03=SL(:,4);l3 = max(abs(l03)) l04=SL(:,5);l4 = max(abs(l04)) figure(3) plot(t03,l01,t03,l02,t03,l03,t03,l04) grid on xlabel('时间/t') ylabel('轮重减载率') title('轮重减载率') legend('一位轮对','二位轮对','三位轮对','四位轮对') hold on %脱轨系数 [ST]=xlsread('100直线参数.xlsx','脱轨系数'); t04=ST(:,1); tt01=ST(:,2);tt1 = max(abs(tt01)) tt02=ST(:,3);tt2 = max(abs(tt02)) tt03=ST(:,4);tt3 = max(abs(tt03)) tt04=ST(:,5);tt4 = max(abs(tt04)) figure(4) plot(t04,tt01,t04,tt02,t04,tt03,t04,tt04) grid on xlabel('时间/t') ylabel('脱轨系数') title('脱轨系数') legend('一位轮对','二位轮对','三位轮对','四位轮对') hold on %倾覆系数 [SQ]=xlsread('100直线参数.xlsx','倾覆系数'); q01=SQ(:,1); v1 = max(abs(q01)) q02=SQ(:,2); v2 = max(abs(q02)) q03=SQ(:,3); v3 = max(abs(q03)) q04=SQ(:,4); v4 = max(abs(q04)) [SHJ]=xlsread('100直线参数.xlsx','横向振动加速度'); t01=SHJ(:,1); hj01=SHJ(:,2);hj1 = max(abs(hj01)) hj02=SHJ(:,3);hj2 = max(abs(hj02)) vhj=(hj1+hj2)/2 [SCQ]=xlsread('100直线参数.xlsx','垂向振动加速度'); t01=SCQ(:,1); cq01=SCQ(:,2);cq1 = max(abs(cq01)) cq02=SCQ(:,3);cq2 = max(abs(cq02)) cqv=(cq1+cq2)/2