Calculo Numerico - Algoritmia, Apuntes de Cálculo para Ingenierios

¡Descarga Calculo Numerico - Algoritmia y más Apuntes en PDF de Cálculo para Ingenierios solo en Docsity!

1

Dra. Ing. Selva Soledad Rivera

Prof. Adjunta

ALGORITMIA

Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

for i=a:b:c

sentencia

sentencia n

end

for i=1:

disp('El valor ingresado es ');

disp(i)

end

 El valor ingresado es  1  El valor ingresado es  2  El valor ingresado es  3 i disp(….) disp(i) 1 'El valor ingresado es

2 'El valor ingresado es

3 'El valor ingresado es

4 Fin del lazo 4 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

Ej: Calcular la sumatoria de los números desde el 1 al 10. %%%SUMATORIA%%%%% suma=0; for i=1: suma=suma+i; end disp('La sumatoria es igual a = '); disp(suma) 5

Leer componentes de un vector for j=1: disp('Ingrese la componente') vec(j)= input(' '); end (^) j disp(…) vec(j)

1 Ingrese la componente vec(1)

2 Ingrese la componente vec(2)

3 Ingrese la componente vec(3)

4 Ingrese la componente vec(4)

5 Ingrese la componente vec(5)

7 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

Escribir componentes de un vector

for j=1:

vec(j)= input(‘Ingrese una componente: ’);

end

vec = [1 2 3 4 5] vector fila vec = [1; 2; 3 ;4 ;5] vector columna 8 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

for j=1:

vec(j)= j;

end

a=0; for i=1: vec(i)=a; a=a+0.5; end vec=0:0.1:20;

EJERCICIO Una estación meteorológica recibe datos de temperatura en cada hora del día. Desarrollar un algoritmo que reciba dichos datos, y al final del día calcule el promedio e imprima el resultado. 10

 % Estación meteorológica

 disp(‘Ingrese una temperatura por cada hora del día')

 N=24; % Asigna en N las 24 horas de un día

 for j=1:N % Varia j desde 1 hasta N

 disp('Ingrese la componente') % escribe

 vec(j)= input(' '); % lee desde teclado la temperatura

 end % fin del bloque variar

 sum=0; % Hace la suma de las componentes

 for j=1:N % Varia j desde 1 hasta N

 sum = sum + vec(j); % suma

 end % fin del bloque variar

 vm = sum/N; % Calcula el promedio

 disp ('El promedio es')

 disp (vm)

j vec(j) sum 1 vec(1) 0+vec(1) 2 vec(2) vec(1)+vec(2) 3 vec(3) vec(1)+vec(2)+vec(3) … …. 11 …….. Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

EJ: Sumar 1+2+3 usando while clc;clear all; i=0; suma=0; while i<= i=i+ suma=suma+i end suma 13

Ej: Sume valores ingresados por teclado mientras la suma sea menor que 200. %%%USO de WHILE flag=1; suma=0; while (flag==1) aux=input('Ingrese un valor < 200 : '); suma=suma+aux; if suma>= flag=0; end end disp(suma) 14

 tol=0.0001; % Valor debajo del cual se considera CERO  ym=1;  iter=0; % contador de iteraciones  while (abs(ym)> tol) % Bloque MIENTRAS  iter = iter + 1; % acumula el número de iteraciones  % Bloque ingreso de datos  disp('Ingrese x1');x1=input(' '); disp('Ingrese y1') ; y1=input(' ');  disp('Ingrese x2'); x2=input(' '); disp('Ingrese y2'); y2=input(' ');  % Bloque de cálculo del Punto Medio  xm=(x1+x2)/2; ym=(y1+y2)/2; disp('El valor de ym es'), ym  end % FIN del Bloque MIENTRAS  % Entrega de resultados  disp ('El valor de ym es prácticamente nulo')  disp ('El número de iteraciones fue'), iter  disp ('El primer punto es'), x1, y  disp ('El segundo punto es'),x2, y 16 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

 tol=0.0001; % Valor debajo del cual se considera CERO  flag=true % control de while  iter=0; % contador de iteraciones  while (flag==true) % Bloque MIENTRAS  iter = iter + 1; % acumula el número de iteraciones  % Bloque ingreso de datos  disp('Ingrese x1');x1=input(' '); disp('Ingrese y1') ; y1=input(' ');  disp('Ingrese x2'); x2=input(' '); disp('Ingrese y2'); y2=input(' ');  % Bloque de cálculo del Punto Medio  xm=(x1+x2)/2; ym=(y1+y2)/2; disp('El valor de ym es'), ym  if (abs(ym)<tol)  flag=false;  end  end % FIN del Bloque MIENTRAS  % Entrega de resultados  disp ('El valor de ym es prácticamente nulo')  disp ('El número de iteraciones fue'), iter  disp ('El primer punto es'), x1, y  disp ('El segundo punto es'),x2, y 17 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

Multiplicar dos vectores  a = [a1 a2 a3] b = [b1 b2 b3]  c = a x b T = a1 x b T 1 + a2 x b T 2 + a3 x b T 3  i=1 i=2 i= c=0; for i=1:  Con matlab: c= a(i) x b(i)+c; end 19 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación

Multiplicar dos vectores

 % Multiplica dos vectores de 5 componentes

 clc;clear;

 disp('Ingrese el vector a');

 for i=1:

 a(i)=input('');

 end

 disp('Ingrese el vector b');

 for i=1:

 b(i)=input('');

 end

 c=0;

 for i=1:

 c=a(i)*b(i)+c;

 end

 disp(c)

20 Dra. Ing. Selva S. Rivera - Prof. Adjunta - Cálculo Numérico y Computación