texto: metodos numericos para ecuaciones diferenciales ordinarias
texto: metodos numericos para ecuaciones diferenciales ordinarias
texto: metodos numericos para ecuaciones diferenciales ordinarias
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
donde<br />
⎛<br />
Y ′ ˙y1<br />
⎜ ⎟<br />
⎜ ⎟<br />
⎜ ˙y2 ⎟<br />
⎜ ⎟<br />
⎜ ˙z1 ⎟<br />
= ⎜ ⎟<br />
⎜ ⎟<br />
⎜ ˙z2 ⎟<br />
⎜ ⎟<br />
⎜ ˙u1<br />
⎟<br />
⎝ ⎠<br />
˙u2<br />
⎞<br />
⎛<br />
⎜<br />
M = ⎜<br />
⎝<br />
0 1 0 0 0 0<br />
− k12<br />
M1 −(B12+B1)<br />
M1<br />
k12<br />
M1<br />
B12<br />
M1<br />
0 0<br />
0 0 0 1 0 0<br />
k12<br />
M2<br />
B12<br />
M2<br />
− (k12+k23)<br />
M2<br />
− (B23+B12+B2)<br />
M2<br />
k23<br />
M2<br />
0 0 0 0 0 1<br />
0 0<br />
⎛<br />
⎜<br />
Y = ⎜<br />
⎝<br />
y1<br />
y2<br />
z1<br />
z2<br />
u1<br />
u2<br />
⎞<br />
⎟<br />
⎠<br />
k23<br />
M3<br />
⎛ ⎞<br />
0<br />
⎜ ⎟<br />
⎜û(t)<br />
⎟<br />
⎜ M1 ⎟<br />
⎜ ⎟<br />
⎜<br />
F = ⎜<br />
0 ⎟<br />
⎜ 0<br />
⎟<br />
⎜ ⎟<br />
⎜<br />
⎝ 0 ⎟<br />
⎠<br />
0<br />
Asumimos los siguientes valores <strong>para</strong> los parámetros<br />
M2 = M3 = 2M1 = 2600kg<br />
N − s<br />
k23 = k12 = 100, 000<br />
m<br />
N − s<br />
B23 = B12 = 500<br />
m<br />
N − s<br />
m<br />
û(t) = 1300N<br />
B2 = B3 = 2B1 = 10000<br />
3.5.5. Programa del Método de Runge-Kutta<br />
Programa monorriel.m<br />
%Halla los desplazamientos y velocidades<br />
%de cada masa del monorriel de dos carros<br />
clear;<br />
clc;<br />
format long<br />
M2=2600;M3=2600;M1=1300;<br />
K23=100000;K12=100000;<br />
38<br />
B23<br />
M3<br />
B23<br />
M2<br />
− k23<br />
M3 −(B23+B3)<br />
M3<br />
⎞<br />
⎟<br />
⎠