• Author / Uploaded
• Paloma Almodovar

Citation preview

1. Realizar un programa en Matlab que calcule el centro de masas de un sistema de partículas en 3D Interés: Sumatorios y medias SOLUCION:  % Script para calcular y representar el centro de masas de un % sistema de partículas n = input('Número de partículas: '); mt = 0; mx = 0; my = 0; mz = 0; for i=1:n fprintf(1,'Partícula %d\n',i); x (i)= input ('Coordenada x: '); y (i)= input ('Coordenada y: '); z (i)= input ('Coordenada z: '); m (i)= input ('Masa: '); mt = mt + m(i); mx = mx + x(i)*m(i); my = my + y(i)*m(i); mz = mz + z(i)*m(i); end xcm = mx / mt; ycm = my / mt; zcm = mz / mt; % Dibuja la posición de las partículas y su centro de masas figure(1) clf hold on for i = 1:n plot3 (x(i),y(i),z(i),'bo'); text (x(i)+.1,y(i)+.1,z(i)+.1,num2str(m(i))); end plot3 (xcm,ycm,zcm,'rx'); text (xcm+.1,ycm+.1,zcm+.1,num2str(mt)); grid on hold off

 

2. Realizar un programa en Matlab que calcule la solución de una ecuación de forma iterativa utilizando el método de Newton-Raphson Interés: Uso de funciones y cálculos iterativos SOLUCION: % Script para aplicar Newton-Raphson a mifuncion % Lectura de datos % Tengo que leer la cota de error 'cota', el número de iteraciones máximo 'n'

% y el valor inicial 'xi' cota = input ('Cota de error objetivo: '); n = input ('Nº máximo de iteraciones: '); xi = input ('Valor inicial de la x: '); xanterior = xi; x = xanterior; eanterior = inf; e = eanterior; i = 1; % contador de iteraciones % Voy a dibujar las tangentes mientras calculo figure(1) clf hold on fplot('mifuncion(x)',[-2 2],'b'); grid % Debo realizar un procedimiento iterativo mientras el error sea mayor que la cota % y el número de iteración sea menor que n while (e>cota) & (i