%功能:画火箭弹性振动的时域相应
%作者:何骁
%时间:2018/09/19
%%静态的模态图
clc; clear;
x0=[0 0.335 0.835 1.335 1.835 2.335 2.791 3.291 3.791 4.291 4.662 5.162 5.662 6.162 6.662 7.162 7.662 8.162 8.662 9.153 9.653 10.153 10.353 10.853 11.353 11.853 12.353 12.853 13.353 13.853 14.353 14.693 15.193 15.693 16.193 16.693 17.193 17.693 18.193 18.693 19.193 19.693 20.193 20.693 21.193 21.693 21.803 21.811 22.411 22.911 23.411 23.911 24.411 24.911 25.411 25.711 26.211 26.711 27.211 27.711 28.211 28.711 28.911 29.411 30.011 30.511 31.011 31.511 32.011 32.511 33.011 33.511 34.011 34.209 34.709 35.359 35.859 36.359 36.859 37.359 37.859 38.359 38.859 39.359 39.859 40.359 40.859 41.359 41.859 42.342 42.842 43.342 43.842 44.342 44.842 45.342 45.842 46.342 46.842 47.048 47.848 48.348 48.848 49.348 49.848 50.348 50.848 51.348 51.848 52.348 52.848 53.348 53.62 54.12 54.398 54.998 55.498 55.998 56.498 56.998 57.498 57.998 58.498 58.998 59.498 59.998 60.498 60.998 61.498 61.998 62.498 62.998 63.498 63.998 64.498 64.998 65.498 65.998 66.498 66.998 67.498 67.998 68.498 68.998 69.498 69.998 70.498 70.998 71.498 71.998 72.498 72.998 73.481 73.931 74.431 74.931 75.431 75.492 75.942 76.442 76.942 77.442 77.942 78.442 78.942 79.442 79.942 80.442 80.942 81.442 81.942 82.442 82.942 83.442 83.942 84.442 84.942 85.134 85.734 86.234 86.734 87.234 87.534 88.034 88.534 89.034 89.534 90.034 90.085 47.842 48.342 48.842 49.342 49.842 50.342 50.842 51.342 51.842 52.342 52.678 53.178 53.619 54.119 54.398 54.998 55.498 55.998 56.498 56.998 57.498 57.998 58.498 58.998 59.498 59.998 60.498 60.998 61.498 61.998 62.498 62.998 63.498 63.998 64.498 64.998 65.498 65.998 66.498 66.998 67.498 67.998 68.498 68.998 69.498 69.998 70.498 70.998 71.498 71.998 72.498 72.998 73.481 73.931 74.431 74.931 75.431 75.492 75.942 76.442 76.942 77.442 77.942 78.442 78.942 79.442 79.942 80.442 80.942 81.442 81.942 82.442 82.942 83.442 83.942 84.442 84.942 85.134 85.734 86.234 86.734 87.234 87.534 88.034 88.534 89.034 89.534 90.034 90.085 47.842 48.342 48.842 49.342 49.842 50.342 50.842 51.342 51.842 52.342 52.678 53.178 53.619 54.119 54.398 54.998 55.498 55.998 56.498 56.998 57.498 57.998 58.498 58.998 59.498 59.998 60.498 60.998 61.498 61.998 62.498 62.998 63.498 63.998 64.498 64.998 65.498 65.998 66.498 66.998 67.498 67.998 68.498 68.998 69.498 69.998 70.498 70.998 71.498 71.998 72.498 72.998 73.481 73.931 74.431 74.931 75.431 75.492 75.942 76.442 76.942 77.442 77.942 78.442 78.942 79.442 79.942 80.442 80.942 81.442 81.942 82.442 82.942 83.442 83.942 84.442 84.942 85.134 85.734 86.234 86.734 87.234 87.534 88.034 88.534 89.034 89.534 90.034 90.085 47.842 48.342 48.842 49.342 49.842 50.342 50.842 51.342 51.842 52.342 52.678 53.178 53.619 54.119 54.398 54.998 55.498 55.998 56.498 56.998 57.498 57.998 58.498 58.998 59.498 59.998 60.498 60.998 61.498 61.998 62.498 62.998 63.498 63.998 64.498 64.998 65.498 65.998 66.498 66.998 67.498 67.998 68.498 68.998 69.498 69.998 70.498 70.998 71.498 71.998 72.498 72.998 73.481 73.931 74.431 74.931 75.431 75.492 75.942 76.442 76.942 77.442 77.942 78.442 78.942 79.442 79.942 80.442 80.942 81.442 81.942 82.442 82.942 83.442 83.942 84.442 84.942 85.134 85.734 86.234 86.734 87.234 87.534 88.034 88.534 89.034 89.534 90.034 90.085 47.842 48.342 48.842 49.342 49.842 50.342 50.842 51.342 51.842 52.342 52.678 53.178 53.619 54.119 54.398 54.998 55.498 55.998 56.498 56.998 57.498 57.998 58.498 58.998 59.498 59.998 60.498 60.998 61.498 61.998 62.498 62.998 63.498 63.998 64.498 64.998 65.498 65.998 66.498 66.998 67.498 67.998 68.498 68.998 69.498 69.998 70.498 70.998 71.498 71.998 72.498 72.998 73.481 73.931 74.431 74.931 75.431 75.492 75.942 76.442 76.942 77.442 77.942 78.442 78.942 79.442 79.942 80.442 80.942 81.442 81.942 82.442 82.942 83.442 83.942 84.442 84.942 85.134 85.734 86.234 86.734 87.234 87.534 88.034 88.534 89.034 89.534 90.034 90.085]';
y0=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';
z0=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 4.1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1 -4.1]';
load('u_7.txt')
omega=2.337304E+00;
lmd=1e3;%放大倍数
%%
%T2/T1决定播放的快慢
% T1=0.1;T2=0.01;
T1=1;t1=0.028;
tic
figure
for t=0:T1:10
T2=T1-t1;%计算暂停时间
t0=cputime;
x=x0+lmd*u_7(1:545,2)*sin(omega*t);
y=y0+lmd*u_7(1:545,3)*sin(omega*t);
z=z0+lmd*u_7(1:545,4)*sin(omega*t);
plot3(x0,y0,z0,x,y,z,'r.','MarkerSize',5);
grid on
axis([-10 100 -10 10 -10 10]);%限定显示的范围
xlabel('x轴');%x轴坐标
ylabel('y轴');%y轴坐标
zlabel('z轴');%z轴坐标
title('振型图');%标题
view(15,15)
% drawnow
t1=cputime-t0;
pause(T2);
end
toc