load 'bneiranji';
load 'bneiranji1';
load 'bneiranji2';
load 'bneiranji3';
load 'bneiranji4';
load 'bneiranji5';
figure(1);
subplot(6,1,1);
plot(t*(0.6/h),xyz(2,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
xlim([7200 7920]);
title('額定轉速下自由端扭角位移時域圖');
subplot(6,1,2);
plot(t*(0.6/h),xyz1(2,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
xlim([7200 7920]);
title('額定轉速e=0.1自由端扭角位移時域圖');
subplot(6,1,3);
plot(t*(0.6/h),xyz2(2,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
xlim([7200 7920]);
title('額定轉速e=0.2自由端扭角位移時域圖');
subplot(6,1,4);
plot(t*(0.6/h),xyz3(2,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
xlim([7200 7920]);
title('額定轉速e=0.3自由端扭角位移時域圖');
subplot(6,1,1);
plot(t*(0.6/h),xyz4(2,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
xlim([7200 7920]);
title('額定轉速e=0.4自由端扭角位移時域圖');
subplot(6,1,1);
plot(t*(0.6/h),xyz5(2,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
xlim([7200 7920]);
title('額定轉速e=0.45自由端扭角位移時域圖');

plot(t*(0.6/h),xyz(9,:));
ylabel('扭角位移(rad)');
xlabel('曲軸轉角(度)');
plot(t*(0.6/h),xyz1(2,:));
xlim([7200 7920]);
title('額定轉速下第九軸段扭角位移時域圖');