domingo, 2 de febrero de 2014

MATLAB Matriz exponencial via laplace

Este script es para encontrar la matriz exponencial de una matriz cuadrada via la transformación de lapace utilizando MATLAB.

% ****************************************************************************
% * Este script encuentra la matriz exponencial y la solución de la matriz A *
% *            Hecho por Jonathan Israel Hernández Hernández                 *
% *                   Nospherathuxxx.blogspot.mx                       *
% ****************************************************************************
clc;clear all;close all; 
disp('Matriz exponencial y la solución de la matriz A')
disp(' ')
q=input('Ingrese la dimención de la matriz A: ');
A=zeros(q);
for i = 1:q       
    for j = 1:q   
        A(i,j) = input(strcat('Ingrese valor numerico de A (',num2str(i),',',num2str(j),') : '));
    end 
    xo(i,1) = input(strcat('Ingrese valor numerico de Xo (',num2str(i),',1) : '));
end
syms s
s1=inv(s.*(eye(q))-A);
disp('La matriz (SI-A)^-1 es: ')
pretty(simplify(s1))
s2=ilaplace(s1);
disp('La matriz e(At) es: ')
pretty(s2)
s3=s2*xo;
for k = 1:q
   strcat('El valor de X(', num2str(k),') = ')
   pretty(s3(k,1))
end
clear all;close all;

No hay comentarios:

Publicar un comentario