Solución numérica y analítica de problema matemático (página 2)
Tabla 1.
Resultados de la Solución Analítica de la Ecuación.
Tiempo | X=0.05 | x=0.10 | X=0.15 | X=0.20 | X=0.25 | X=0.30 | X=0.35 | X=0.40 | X=0.45 | X=0.50 | X=0.55 | X=0.6 | X=0.65 | X=0.70 | X=0.75 | X=0.80 | X=0.85 | X=0.90 | X=0.95 |
0.00 | 0.0023 | 0.0095 | 0.0233 | 0.0390 | 0.0638 | 0.0884 | 0.1245 | 0.1577 | 0.2052 | 0.2468 | 0.3062 | 0.3556 | 0.4277 | 0.4838 | 0.5702 | 0.6302 | 0.7357 | 0.7899 | 0.9426 |
0.05 | 0.0047 | 0.0130 | 0.0242 | 0.0435 | 0.0637 | 0.0941 | 0.1230 | 0.1648 | 0.2022 | 0.2558 | 0.3012 | 0.3670 | 0.4197 | 0.4989 | 0.5570 | 0.6530 | 0.7101 | 0.8392 | 0.3950 |
0.10 | 0.0103 | 0.0195 | 0.0333 | 0.0489 | 0.0739 | 0.0983 | 0.1345 | 0.1676 | 0.2153 | 0.2567 | 0.3165 | 0.3653 | 0.4382 | 0.4929 | 0.5738 | 0.6368 | 0.7339 | 0.3151 | -0.1034 |
0.15 | 0.0147 | 0.0305 | 0.0442 | 0.0636 | 0.0836 | 0.1142 | 0.1429 | 0.1850 | 0.2219 | 0.2760 | 0.3207 | 0.3877 | 0.4385 | 0.5210 | 0.5817 | 0.6850 | 0.2419 | -0.1862 | -0.0799 |
0.20 | 0.0203 | 0.0394 | 0.0609 | 0.0788 | 0.1040 | 0.1282 | 0.1647 | 0.1973 | 0.2457 | 0.2860 | 0.3474 | 0.3940 | 0.4705 | 0.5184 | 0.6204 | 0.1778 | -0.2415 | -0.1531 | -0.0828 |
0.25 | 0.0247 | 0.0506 | 0.0741 | 0.1012 | 0.1234 | 0.1545 | 0.1826 | 0.2254 | 0.2613 | 0.3170 | 0.3592 | 0.4301 | 0.4738 | 0.5739 | 0.1234 | -0.2953 | -0.2172 | -0.1540 | -0.0732 |
0.30 | 0.0303 | 0.0594 | 0.0910 | 0.1187 | 0.1517 | 0.1778 | 0.2152 | 0.2466 | 0.2967 | 0.3345 | 0.3997 | 0.4392 | 0.5335 | 0.0788 | -0.3687 | -0.2716 | -0.2148 | -0.1373 | -0.0712 |
0.35 | 0.0347 | 0.0707 | 0.1039 | 0.1415 | 0.1731 | 0.2125 | 0.2419 | 0.2864 | 0.3199 | 0.3795 | 0.4144 | 0.5032 | 0.0442 | -0.4091 | -0.3161 | -0.2653 | -0.1917 | -0.1320 | -0.0641 |
0.40 | 0.0404 | 0.0792 | 0.1212 | 0.1583 | 0.2022 | 0.2371 | 0.2837 | 0.3151 | 0.3693 | 0.3998 | 0.4829 | 0.0195 | -0.4392 | -0.3508 | -0.3056 | -0.2363 | -0.1825 | -0.1185 | -0.0607 |
0.45 | 0.0446 | 0.0909 | 0.1336 | 0.1819 | 0.2224 | 0.2735 | 0.3104 | 0.3664 | 0.3949 | 0.4726 | 0.0047 | -0.459 | -0.3755 | -0.3360 | -0.2709 | -0.2229 | -0.1630 | -0.1112 | -0.0544 |
0.50 | 0.0505 | 0.0990 | 0.1516 | 0.1977 | 0.2532 | 0.2956 | 0.3562 | 0.3902 | 0.4699 | 0 | -0.4699 | -0.390 | -0.3562 | -0.2956 | -0.2532 | -0.1977 | -0.1516 | -0.0990 | -0.0505 |
1- Gráfica de la Función con la Solución Analítica.
x=0:0.05:1;
t=0:0.05:0.5;
[T X]=meshgrid (t,x);
suma=0;
N=5000;
k=size(T);
Z=zeros(k(1),k(2));
for n=1:N
Z= Z + (2*cos(n*pi)/(n*pi)*(2/(n*pi)^2 – 1) – 4/(n*pi)^3)*cos(n*pi*T).*sin(n*pi*X);
mesh(T,X,Z)
hold on
end
xlabel('Eje t');
ylabel('Eje x');
zlabel('Eje z');
print -dpdf grafica_ninamaria.pdf
Donde, la i y la j representan enteros que interceptan rectas
horizontales y verticales llamadas puntos de malla.
Se procede a resolver las ecuaciones (ec.20) y (ec.21) con las condiciones
de frontera e iniciales en una región rectangular en el plano x,
t definido por las desigualdades
donde:
La a es el extremo derecho del intervalo de la variable x, t es
el extremo derecho del intervalo de la variable t, T es el valor de
tiempo especificado, n es un entero positivo, al igual que la m;
h representa la unidad de separación entre las líneas verticales;
k es la unidad de separación que existe en las líneas
horizontales. Xi son las retículas de líneas verticales;
tj son retículas de líneas horizontales. depende
de los enteros n y m, y su valor nos da la precisión
de los resultados con las técnicas de diferencias finitas. De lo anterior
se desprende que:
Continuando así con las sustituciones de j en la ecuación
(ec.23) hasta llegar al 9, haciendo uso tanto de las condiciones de fronteras
como de las condiciones iniciales y de los resultados de la (ec.22), encontramos
en cada línea de tiempo todos los valores aproximados de la ecuación
que se está resolviendo.
Ver resultados de la solución numérica de la ecuación
(Tabla 2).
Página anterior | Volver al principio del trabajo | Página siguiente |