Trabajo Matlab PDF

Preview

Summary

Porgramación en matlab desde cero, funciones basicas como lo son el plot....

Description

1. Escriba un programa de funciones bien comentado que grafica las funciones sin x, sin 2x, sin 3x, sin 4x, sin 5x y sen 6x en el intervalo [0, 2π] en un gráfico. Use un tamaño de paso suficientemente pequeño para suavizar todos los gráficos. Entregue el programa y el gráfico. Código: clc clear all x=[0:pi/100:2*pi];%definimos a x %nombramos cada grafica de manera diferente g1=sin(x);%definimos cada grafica g2=sin(2*x); g3=sin(3*x); g4=sin(4*x); g5=sin(5*x); g6=sin(6*x); %utilizamos subplot para poder analizar una cantidad n de graficas en una %sola figura subplot(2,3,1),plot(x,g1)%pemite graficar g1(x) xlabel('x'); %esta opcion me permite nombrar el eje x ylabel('sin(x)');%esta opcion me permite nombrar el eje y subplot(2,3,2),plot(x,g2) xlabel('x'); ylabel('sin(2x)'); subplot(2,3,3),plot(x,g3) xlabel('x'); ylabel('sin(3x)'); subplot(2,3,4),plot(x,g4) xlabel('x'); ylabel('sin(4x)'); subplot(2,3,5),plot(x,g5) xlabel('x'); ylabel('sin(5x)'); subplot(2,3,6),plot(x,g6) xlabel('x'); ylabel('sin(6x)');

Gráficas:

2. Use MATLAB para generar una gráfica de esta función de t = 0 a 0.8, dado que q0 = 10, R =60 60, L = 9 y C = 0.00005. Código: clc clear all hold off %CARGA EN UN CONDENSADOR r=60 ;%ohms q=10 ; t=[0:0.001:0.8];%definimos t L=9;% Henrios C=1/20000;% farads %separamos la funcion q(t) en varias funciones exponencial=exp((-r*t)/(2*L)); %exp(x)simula la funcion de euler y el argumento es el exponente al cual %queriamos elevar euler raiz=sqrt((1/(L*C))-((r^2)/(2*L)));%sqrt es la funcion de raiz cuadratica carga=q*exponencial.*cos(raiz*t); plot(t,carga)%Gráfica t vs carga title(' Proceso de carga de un condensador');%ponemos titulo a la figura xlabel('Tiempo (t)');%Ponemos titulo en el eje X e Y ylabel('Carga q(t)') Gráficas:

3. Use MATLAB para almacenar F y x como vectores y luego calcule los vectores de las constantes de resorte y las energías potenciales. Use la función max para determinar la energía potencial máxima. Código: clc clear all K=[];%defino un espacio vacio para el vector de las constantes Distancia=[0.013 0.02 0.009 0.010 0.012] F=[14 18 8 9 13] %definimos la distancia y fuerza como un vector X=Distancia.^2; for a=1:5 %utilizamos un for para que realice la division para encontrar K K(a)=F(a)/Distancia(a); end U=0.5*[K.*X] MaxE=max(U); texto1=['El vector de la energia potencial es U=[ ' num2str(U) ' ]']; %la funcion num2str convierte a un numero en una cadena de caracteres texto2=['El vector de las constantes de los resortes es K=[ ' num2str(K) ']']; texto3=['La energia potencial maxima obtenida fue ' num2str(MaxE)]; disp(texto1)%la funcion disp muestra como cadena de caracteres un texto disp(texto2)

disp(texto3)

4. Genere valores de x e y para valores de t de 0 a 100 con ∆t = 1/16 Construya gráficas de: a. x e y versus t. b. y versus x. Use subplot mostrar ambas gráficas en una sola. Incluya títulos y etiquetas de eje. Código: clc clear all t=[0:1/16:100]%Definimos a t X= sin(t).*(exp(cos(t))-2*cos(4*t)-(sin(t/12)).^5); Y= cos(t).*(exp(cos(t))-2*cos(4*t)-(sin(t/12)).^5); plot(t,X)%graficamos a X(t) %la funcion subplot permite graficar en una sola figura varias funciones %como lo vemos en esta ultimo graficamos a x e y al mismo tiempo vs t subplot(2,2,1) title('Subplot 1: x=[sin(t).*(exp(cos(t))-2*cos(4*t)-(sin(t/12)).^5)]') xlabel('Tiempo (t)');%nombramos el eje x e y ylabel('Función X'); subplot(2,2,2) plot(t,Y)%graficamos a Y(t) title('Subplot 2: Y=[cos(t).*(exp(cos(t))-2*cos(4*t)-(sin(t/12)).^5)]') xlabel('Tiempo (t)'); ylabel('Función Y'); subplot(2,2,[3,4]); plot(t,X,'b',t,Y,'g'); title('Subplot 3 and 4: Y e X vs t')%tittle permite nombrar cada grafica xlabel('Tiempo (t)'); ylabel('Función X e Y'); Gráficas:...