clear
k=400;m=2;A=0.1;t=linspace(0,0.9,1000);
x=dsolve('D2x+w^2*x=0','Dx(0)=0,x(0)=A','t')    %Çó½âÕñ×ÓÕñ¶¯·½³Ì
v=diff(x,'t')
a=diff(x,'t',2)
w=sqrt(k/m);
V1=0.5*m*w^2*x^2;
T=0.5*m*v^2;
E=T+V1;
x=eval(x);
v=eval(v);
a=eval(a);
V1=eval(V1);
T=eval(T);fontweight取值
E=eval(E);
%pretty(x);
%pretty(v);
%pretty(a);
%pretty(V1);
subplot(3,2,1);plot(t,x,'-
'),xLabel('t','fontsize',12,'fontweight','bold','color','b'),yLabel('x', 'fontsize',12,'fontweight','bold','color','b');
subplot(3,2,2);plot(t,v,'-
'),xLabel('t','fontsize',12,'fontweight','bold','color','b'),yLabel('v', 'fontsize',12,'fontweight','bold','color','b');
subplot(3,2,3);plot(t,a,'-
'),xLabel('t','fontsize',12,'fontweight','bold','color','b'),yLabel('a', 'fontsize',12,'fontweight','bold','color','b');
subplot(3,2,4);plot(x,T,'-
'),xLabel('x','fontsize',12,'fontweight','bold','color','b'),yLabel('T', 'fontsize',12,'fontweight','bold','color','b');
subplot(3,2,5);plot(t,E,'-
'),xLabel('x','fontsize',12,'fontweight','bold','color','b'),yLabel('E(x )','fontsize',12,'fontweight','bold','color','b');
subplot(3,2,6);plot(x,V1,'-
'),xLabel('x','fontsize',12,'fontweight','bold','color','b'),yLabel('V(x )','fontsize',12,'fontweight','bold','color','b');
x =
(A*(1/exp(t*w*i)))/2 + (A*exp(t*w*i))/2
v =
- (A*w*(1/exp(t*w*i))*i)/2 + (A*w*exp(t*w*i)*i)/2
a =
- (A*w^2*(1/exp(t*w*i)))/2 - (A*w^2*exp(t*w*i))/2

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。