求高手帮忙改动一下小程序

zhangshixiang

Mu''+Cu'+Ku'=F(t)这是方程，激励列向量取F(t)={100sin(pi*t);100sin(pi*t);100sin(pi*t)}程序是计算顶层的水平位移，我现在想得到顶层水平位移的时域图，本人水平太菜了，请求高手帮忙！搞完寒假早点回家了！

1. clear
   
2. M=diag([100;100;100]);
   
3. C=[170 -85 0;-85 170 -85;0 -85 170];
   
4. K=[2800 -1400 0;-1400 2800 -1400;0 -1400 2800];
   
5. H=[-inv(M)*C/2 inv(M);C*(inv(M))*C/4-K -C*(inv(M))/2];
   
6. I=eye(size(H));
   
7. N=20;
   
8. delt_t=0.5;
   
9. tao=delt_t/20^N;
   
10. T1=H*tao;
    
11. for i=1:4
    
12. T1=T1+(H*tao)^(i+1)/prod(1:i+1);
    
13. end
    
14. clear i
    
15. for i=1:N
    
16. T_N=2*T1+T1*T1;
    
17. T1=T_N;
    
18. end
    
19. T=eye(size(T_N))+T_N;
    
20. v=[0;0;0;0;0;0];
    
21. r1=[0;0;0;100;100;100];
    
22. r2=[0;0;0;0;0;0];
    
23. A=(eye(6)*pi+H^2/pi)^-1*(r2-H*r1/pi);
    
24. B=(eye(6)*pi+H^2/pi)^-1*(-r1-H^2*r2/pi);
    
25. clear i
    
26. for i=1:1:(20/delt_t+1)
    
27. v=T*v+(eye(6)*sin(pi*delt_t*i)-T*sin(pi*delt_t*(i-1)))*A+(eye(6)*cos(pi*delt_t*i)-T*cos(pi*delt_t*(i-1)))*B+T*(inv(H))*[0;0;0;2.725;2.725;2.725]+(inv(H))*[0;0;0;2.725;2.725;2.725];
    
28. V(i,:)=v;
    
29. end
    

复制代码

anbcjys

画图用plot

liushuangq05

Mu''+Cu'+Ku'=F(t)
你这个公式是不是没写对?是Mu''+Cu'+Ku=F(t)吧?如果是的话,那么一般的,这个公式的u是u(x,y,z,t)的函数形式,你把F(t)={100sin(pi*t);100sin(pi*t);100sin(pi*t)}代入,计算出每个离散t值对应的u(x,y,z),那么你绘制u(x,y,z,t)随时间t的变化图是很容易的.
顶层水平位移可能是对z的值做了限制.