MATLAB TRABAJOS COMANDO

% PREGUNTA N°01

% PROGRAMA QUE DETERMINA SI UN NÚMERO ES O NO PERFECTO

n= input('Ingrese un número entero positivo =')

S= 0;

for i=1:n-1

    if rem(n,i)== 0

        S= S+i;

    end

end

if S==n

    disp('El número es perfecto')

else

    disp('El número no es perfecto')

end

% PROGRAMA GUARDADO CON EL NOMBRE DE:'NUMPER'

% PREGUNTA N°02

L= input('Ingrese la longitud del alambre,L=')

if L>0

x= 0:0.1:L;

    v= '((L-3*x)/4).^2 + (sqrt(3)*x.^2)/4';

    A= eval(v);

    B= max(A);

    C= min(A);

fprintf('Máxima área:%4.5f\n',B)

fprintf('Mínima área:%4.5f\n',C)

end

% PROGRAMA GUARDADO CON EL NOMBRE DE:'MAXMINAR'

% PREGUNTA N°02

L= input('Ingrese la longitud del alambre,L=')

if L>0

x= 0:0.1:L;

    v= '((L-3*x)/4).^2 + (sqrt(3)*x.^2)/4';

    A= eval(v);

    B= max(A);

    C= min(A);

fprintf('Suma máxima de las áreas:%2.2f\n',B)

fprintf('Suma mínima de las áreas:%2.2f\n',C)

end

% PROGRAMA GUARDADO CON EL NOMBRE DE:'MAXMINAR'

%PREGUNTA N°03

% PROGRAMA QUE CALCULA LA DISTANCIA ENTRE UN PUNTO Y UNA RECTA

function [d] = DistPal( xo,yo,A,B,C )

%d: Distancia entre el punto y la recta

%xo:Punto en el eje de abcisas 'x'

%yo:Punto en el eje de ordenadas 'y'

%A: Coeficiente de x

%B: Coeficiente de y

%C: Término independiente

d=abs(A*xo+B*yo+C)/sqrt(A^2+B^2);

end

% PROGRAMA GUARDADO CON EL NOMBRE DE:'DistPal'

% PREGUNTA N°04

% PROGRAMA QUE DETERMINA LA VELOCIDAD Y ÁNGULO DE LANZAMIENTO

% DE DOS MISILES

% Para el misil 1

% LEYENDA

% tetha1: Ángulo 1 (°)

% tetha2: Ángulo 2 (°)

% vo1: Velocidad inicial 1(m/s)

% vo2: velocidad inicial 2(m/s)

% vo1x: vo1 en el eje de abscisas  

% vo2x: vo2 en el eje de abscisas

% vo1y: vo1 en el eje de ordenadas

% vo2y: vo2 en el eje de ordenadas

% a1: tetha1 (rad)

% a2: tetha2 (rad)

% tt1: Tiempo total del primer misil

% hmax1: Altura máxima 1

% thmax1: Tiempo hasta la altura máxima 1

% dAB: Distancia entre el punto A y B (m)

% dAC: Distancia entre el punto A y C (m)

% dCD: Distancia entre el punto C y D (m)

g=9.8

theta1= input('Ingrese el valor del ángulo 1 en °:')

vo1= input('Ingrese la velocidad inicial 1 (m/s):')

dAC= input('Ingrese la distancia entre el punto A y C (m):')

a1= theta1*pi/180;

vo1x= vo1*cos(a1);

vo1y= vo1*sin(a1);

tt1= 2*vo1y/g;

fprintf('Tiempo total hasta que el misil 1 impacta con el suelo=%2.2f\n', tt1)

thmax1= vo1y/g;

fprintf('Tiempo hasta la altura máxima del misil 1=%2.2f\n',thmax1)

dAB= vo1x*tt1;

fprintf('Distancia entre el punto A y B del misil 1=%2.2f\n', dAB)

hmax1=(vo1y^2)/(2*g);

t1=0:0.1:tt1;

x=vo1x*t1;

y=vo1y*t1-(0.5*g.*t1.^2);

plot(x,y,'m')

xlabel('EJE X')

ylabel('EJE Y')

text(8000,4250,'TRAYECTORIA MISIL 1')

grid on

hold on 

% Para el misil 2

% RECORDAR

% vo2x: vo2 en el eje de abscisas

% vo2y: vo2 en el eje de ordenadas

% tetha2: Ángulo 2 (°)

% a2: tetha2 (rad)

% LEYENDA

% tc1: Tiempo de choque entre el misil 1 y 2

% ti1: Tiempo desde la altura máxima hasta el el tiempo de choque del

%        misil 1

% vf1y: Velocidad final 1 en el de ordenadas en el punto de choque

% dAi: Distancia desde el punto A hasta el punto de impacto

% dCi: Distancia desde el punto C hasta el punto de impacto

% D: Punto de impacto del misil 2 con el suelo

% tt2: Tiempo total del proyectil 2

% Por dato: Altura de choque= 0.8*hmax1

% Además se sabe que la velocidd en la altura máxima es 0 m/s

vf1y=sqrt(2*g*0.8*hmax1);

ti1=vf1y/g;

tc1=thmax1+ti1;

fprintf('Tiempo de choque entre el misil 1 y 2=%2.2f\n',tc1)

dAi=vo1x*tc1;

fprintf('Distancia desde el punto A hasta el punto de impacto=%2.2f\n',dAi)

dCi=dAC-dAi;

fprintf('Distancia desde el punto C hasta el punto de impacto=%2.2f\n',dCi)

vo2x=dCi/(tc1+5);

fprintf('Velocidad inicial 2 en el eje de abscisas=%2.2f\n',vo2x)

vo2y= (0.8*hmax1+(0.5*g*tc1^2))/(tc1+5);

...