Guias de analisis numerico

Solo disponible en BuenasTareas
  • Páginas : 3 (585 palabras )
  • Descarga(s) : 0
  • Publicado : 6 de junio de 2011
Leer documento completo
Vista previa del texto
Guía de ejercicios correspondiente a la practica 5 Métodos Numéricos y Programación Digital. Ciclo 01/2011 La siguiente guía expone los métodos de cálculo de raíces e interpolación polinomial; losarchivos de MATLAB (m-file) se proporcionan para que el estudiante los verifique en el ) laboratorio e incluso en su casa, como un complemento al trabajo numérico desarrollado en la clase. A partir deacá, MATLAB se usará para programar la solución computacional , a los distintos métodos numéricos que se vean en clase. dos Previo a esta práctica ya se han realizado cuatro laboratorios explicandoconceptos tales como: Lazos for, if, while, etc. Comandos para elaboración de gráficas y resolución de aplicaciones sencillas; además, para cada programa se especifican las entradas aciones necesariaspara poder ejecutarlo y las salidas que se esperan obtener como respuestas respuestas, como una forma de facilitar su comprensión. RAICES DE UNA ECUACIÓN. En esta sección el objetivo es encontrar lasraíces de ecuaciones de una variable. Si se ecuaciones supone que f(x) es una función continua. Cualquier numero r tal que f(r) =0 se llama raíz de la ecuación f(x) =0; también se dice que r es un cerode la función f(x). Los métodos ; que se disponen para este caso son los de la bisección y el punto fijo. Método de Bisección:
function [c,err,yc]=bisect(f,a,b,max1,delta) %Entradas %f es la función aevaluar %a y b son el extremo izquierdo y el extremo derecho %delta es la tolerancia %Salidas %c es el cero %yc=f(c) %err es el error estimado de la aproximación ya=feval(f,a); yb=feval(f,b); ifya*yb>0 end for k=1:max1 c=(a+b)/2; yc=feval(f,c); if yc==0 a=c; b=c; elseif yb*yc>0 b=c; yb=yc; else a=c; ya=yc;

end end if (b-a)> x=[4 4.2 4.5 4.7 5.1 5.5 5.9 6.3 6.8 7.1]; >> f=[102.56 113.18130.11 142.05 167.53 195.14 224.87 256.73 299.50 326.72]; >> plot(x,f,'o');grid on;hold on; >> y=@(x) 72.0845*x-194.1382; >> y1=feval(y,x1); >> plot(x1,y1,'g');hold off;

En donde se puede apreciar...
tracking img