krylov ayuda!!!!
Publicado por Juan Manuel (3 intervenciones) el 13/10/2015 21:20:31
Estoy haciendo el código del método de krylov en matlab, pero no he podido avanzar de este punto, básicamente me falta que muestre la ecuación característica con el carácter s como sym, pero no me sale, les anexo el código, si me pueden ayudar lo agradecería de verdad...
si lo copian y pegan en Matlab les mostrará paso por paso el proceso, pero se supone que debe llegar a una ecuación, de la forma
s^3+s^2+s+1, solo es un ejemplo, los coeficientes varian...
les dejo la matriz con la que lo estuve intentando ejecutar... [1 -1 0;-2 0 1;-1 3 -1]
si lo copian y pegan en Matlab les mostrará paso por paso el proceso, pero se supone que debe llegar a una ecuación, de la forma
s^3+s^2+s+1, solo es un ejemplo, los coeficientes varian...
les dejo la matriz con la que lo estuve intentando ejecutar... [1 -1 0;-2 0 1;-1 3 -1]
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
%Metodo de Krylov
close all
clear all
clc
%matriz cuadrda
N=input('Introduce el tamaño de matriz:');
A=input('Introduce la matriz cuadrada:');
% Aqui inicializamos las matrices a rellenar con valores calculdos
C=zeros(N,N)
y=zeros(N,1)
y(1)=1
C(:,N)=y
%Aqui va el ciclo
for i=(N-1):-1:1
C(:,i)=A*C(:,(i+1))
end
a=-(A^N)*y
%Calcular los coeficientes b1, b2,... bn
b=inv(C)*a
s=sym('s')
I=s*eye(size(A))
b2=I*b
%Aqui falta anotar el polinomio en forma symbolic
p=[1 b(1) b(2) b(3)];
%Calculamos las raices del polinomio
r = roots(p)
Valora esta pregunta
0
