algoritmo en matlab

“Año de la Integración Nacional y el Reconocimiento de Nuestra Diversidad”
UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERIA INDUSTRIAL Y DE SISTEMAS
ESCUELA PROFESIONAL DE INGENIERIA DE SISTEMAS

PROGRAMACIÓN EN MATLAB 7.8.0
“MÉTODOS NUMÉRICOS”

ALUMNOS: CÓDIGO:
BENITES ESPINOZA, Marco Antonio 092812H
TATAJE ROJAS, Axhel Joan 092829H

ASIGNATURA:
ANÁLISISNUMÉRICO
PROFESOR:
Sergio Huaranca Tanta

FEBRERO, 2012


TEMA Nº 1
MÉTODO DE LA BISECCIÓN
ALGORITMO EN MATLAB 7.8.0 (2009)
disp(' METODO DE BISECCION ');
disp(' ------------------- ');

f=input('INGRESE FUNCION: ','s');
ezplot(f);
gridon;
A=input('INGRESE LIMITE INFERIOR DEL INTERVALO: ');
B=input('INGRESE LIMITE SUPERIOR DEL INTERVALO: ');
error=input('INGRESE ERROR: ');f=inline(f);
i=1;
E(1)=1;

if f(A)*f(B)0
signo_a='+';
else
signo_a='-';
end
ifsign(f(x(1)))>0
signo_x='+';
else
signo_x='-';
end
ifsign(f(b(1)))>0
signo_b='+';
else
signo_b='-';
end

fprintf('\n i a f(a) x f(x) b f(b) Error aprox \n');
fprintf('%2d \t %11.7f \t %s \t %11.7f \t %s \t %11.7f \t %s\n',i,a(1),signo_a,x(1),signo_x,b(1),signo_b);while E(i)>error

ifsign(f(a(i)))>0
signo_a='+';
else
signo_a='-';
end
ifsign(f(x(i)))>0
signo_x='+';
else
signo_x='-';
end
ifsign(f(b(i)))>0
signo_b='+';
else
signo_b='-';
end

if f(a(i))*f(x(i))>Metodo_Biseccion
METODO DE BISECCION
-------------------
INGRESE FUNCION: sin(log(x)+sqrt(x))
INGRESE LIMITE INFERIOR DEL INTERVALO: 3
INGRESE LIMITE SUPERIOR DEL INTERVALO: 4INGRESE ERROR: 0.00001

i a f(a) x f(x) b f(b) Error aprox
1 3.0000000 + 3.5000000 + 4.0000000 -
2 3.5000000 + 3.7500000 + 4.0000000 - 0.0666667
3 3.5000000 + 3.6250000 - 3.7500000 - 0.0344828
4 3.5000000 + 3.5625000 - 3.6250000 - 0.0175439
5 3.5000000 + 3.5312500 - 3.5625000 - 0.00884966 3.5312500 + 3.5468750 + 3.5625000 -0.0044053
7 3.5312500 + 3.5390625 - 3.5468750 -0.0022075
8 3.5312500 + 3.5351563 - 3.5390625 -0.0011050
9 3.5312500 + 3.5332031 - 3.5351563 -0.0005528
10 3.5312500 + 3.5322266 - 3.5332031 - 0.0002765
11 3.5322266 + 3.5327148 + 3.5332031 - 0.0001382
12 3.5322266 + 3.5324707 - 3.5327148-0.0000691
13 3.5324707 + 3.5325928 + 3.5327148 -0.0000346
14 3.5325928 + 3.5326538 + 3.5327148 -0.0000173

La raíz existente en el intervalo es aprox: 3.532654

Ejemplo 2:Aproxima la primera raíz de la función , mediante el método de la bisección, para un error menor o igual a .
Operando en MATLAB

METODO DE BISECCION
-------------------
INGRESE FUNCION:-2*x^5-6*x^3+2*x+4
INGRESE LIMITE INFERIOR DEL INTERVALO: 0
INGRESE LIMITE SUPERIOR DEL INTERVALO: 1
INGRESE ERROR: 0.0001

i a f(a) x f(x) b f(b) Error aprox
1 0.0000000 + 0.5000000 + 1.0000000 -
2 0.5000000 + 0.7500000 + 1.0000000 - 0.3333333
3 0.7500000 + 0.8750000 + 1.0000000 - 0.1428571
4 0.8750000 + 0.9375000 + 1.0000000- 0.0666667
5 0.8750000 + 0.9062500 - 0.9375000 - 0.0344828
6 0.9062500 + 0.9218750 + 0.9375000 - 0.0169492
7 0.9062500 + 0.9140625 - 0.9218750 - 0.0085470
8 0.9062500 + 0.9101563 - 0.9140625 - 0.0042918
9 0.9101563 + 0.9121094 + 0.9140625 - 0.0021413
10 0.9121094 + 0.9130859 + 0.9140625 - 0.0010695
11 0.9121094+ 0.9125977 - 0.9130859 - 0.0005350
12 0.9121094 + 0.9123535 - 0.9125977 - 0.0002676
13 0.9123535 + 0.9124756 + 0.9125977 - 0.0001338

La raiz existente en el intervalo es aprox: 0.912476

TEMA Nº 2
MÉTODO DE LA FALSA POSICIÓN
ALGORITMO EN MATLAB 7.8.0 (2009)
disp('METODO DE FALSA POSICIÓN');
disp('========================');

f=input('Ingrese la...