Masa-Resorte Matlab
clear all,clc
E1=inline('z','t','z');
E2=inline('9.8-200*h-20*z','t','h','z');
h=0;z=(2*9.8*h)^0.5;H=0.1;
disp('t h z k1 c1 k2 c2 k3 c3 k4 c4')for t=0:H:2
fprintf('\n %7.4f %7.4f %7.4f',t,h,z)
k1=E1(t,z);
c1=E2(t,h,z);
k2=E1(t+H/2,z+c1*H/2);
c2=E2(t+H/2,h+k1*H/2,z+c1*H/2);
k3=E1(t+H/2,z+c2*H/2);c3=E2(t+H/2,h+k2*H/2,z+c2*H/2);
k4=E1(t+H,z+c3*H);
c4=E2(t+H,h+k3*H,z+c3*H);
h=h+H*(k1+2*(k2+k3)+k4)/6;
z=z+H*(c1+2*(c2+c3)+c4)/6;
fprintf(' %10.4f %10.4f %10.4f %10.4f %10.4f %10.4f%10.4f %10.4f\n',k1,c1,k2,c2,k3,c3,k4,c4)
end
t h z k1 c1 k2 c2 k3 c3 k4 c4
0.0000 0.0000 0.0000 0.00009.8000 0.4900 0.0000 0.0000 4.9000 0.4900 0.0000
0.1000 0.0245 0.3267 0.3267 -1.6333 0.2450 -3.2667 0.1633 -0.8167 0.2450 -3.2667
0.20000.0476 0.1089 0.1089 -1.9056 0.0136 -1.0889 0.0544 -0.9528 0.0136 -1.0889
0.3000 0.0519 -0.0091 -0.0091 -0.4083 -0.0295 0.0907 -0.0045 -0.2042-0.0295 0.0907
0.4000 0.0502 -0.0181 -0.0181 0.1285 -0.0117 0.1815 -0.0091 0.0643 -0.0117 0.1815
0.5000 0.0490 -0.0048 -0.0048 0.0996 0.00020.0479 -0.0024 0.0498 0.0002 0.0479
0.6000 0.0488 0.0009 0.0009 0.0153 0.0017 -0.0092 0.0005 0.0077 0.0017 -0.0092
0.7000 0.0489 0.00100.0010 -0.0087 0.0005 -0.0097 0.0005 -0.0044 0.0005 -0.0097
0.8000 0.0490 0.0002 0.0002 -0.0050 -0.0001 -0.0020 0.0001 -0.0025 -0.0001 -0.00200.9000 0.0490 -0.0001 -0.0001 -0.0005 -0.0001 0.0007 -0.0000 -0.0002 -0.0001 0.0007
1.0000 0.0490 -0.0001 -0.0001 0.0005 -0.0000 0.0005 -0.0000...
Regístrate para leer el documento completo.