Capitulo 5 de Matlab

Diego Fernando reyes Triana 1110138

Electiva
x=0:10; %Declaramos otra con la exponencial de x
y1=exp(x);
plot(x,y1);
title('exponencial');
xlabel('x');
ylabel('y1=exp(x)');
grid
%%
x=0:10;

y2=sin(x);

plot(x,y2);%une los puntos de la grafica desde -pi hasta pi

title('sinusoidal wave');

xlabel('x');

ylabel('y2=sin(x)');

grid

%%
x=0:10;
a=5;
b=2;
c=4;
y3=(a*x.^2)+(b.*x)+c;
plot(x,y3);
title('quadratic wave');
xlabel('x');
ylabel('y3=ax^2+bx+c');
grid
%%
x=0:10;

y4=sqrt(x);

plot(x,y4);

title('sqrt(x) wave');

xlabel('x');
ylabel('y4=sqrt(x)');

grid

%%
y=[12,14,12,22,8,9];

plot(y)

title('plot of the data');

xlabel('values set by MATLAB');

ylabel('given data y');

%%
x=-pi:0.01*pi:pi;

y1=sin(x);

plot(x,y1);

title('Figures for y1,y2 and y3')

xlabel('x in radians');

ylabel('sin(x)/ sin(2x)/sin(3x)');

hold on

y2=sin(2*x);

plot(x,y2);

hold on

y3=sin(3*x);

plot(x,y3);

hold off

grid

%%
x=-pi:0.01*pi:pi;

y1=sin(x);

hold on

plot(x,y1,'--r');

title('Figures for y1,y2 and y3');% usamos el title para poner titulo

xlabel('x en raddianes');
ylabel('sin(x)/ sin(2x)/sin(3x)');

y2=sin(2*x);

hold on

plot(x,y2,'-b');%Declaramos otra con la exponencial de x

y3=sin(3*x);

hold on

plot(x,y3,':g');

hold off
%%
%La distancia que recorre un proyectil cuando se dispara a un ángulo u es
función
%del tiempo y se puede dividir en distancias horizontal y vertical de acuerdo
con las
%fórmulas

%Grafica distancia horizontal contra tiempo.

g=9.81;
vo=input("ingrese velocidad incial ");
angulo=input("ingrese el angulo");
Vhorizontal=((vo)*cos(angulo));
H=input("ingrese altura");
teta=teta*pi/180;
t=((vo*sin(teta))+((vo*sin(teta))^2+2*H*g)^0.5)/g;

xMax=vo*cos(teta)*teta;
yMax=H((vo*sin(teta)^2))/(2*g);

Times=[0:2:20];
x=vo*cos(teta).*Times;
y=h+(vo*sin(teta).*Times)

%%

g=9.8;
t=(0:20);
v0=100;
A1=pi/2;
A2=pi/4;
A=pi/4;% las 7 primras lineas elejimos nuestras variables
Horizon1=(t*v0*cos(A1)) %
Vert1=(t*v0*sin(A1)-1/2*g*t.^2)%cogemos
subplot(2,2,1) ; % usamos suplot para unir las graficsa en un mismo cuadro
plot(Horizon1,Vert1)% usamos plot para graficar nuestras anteriores variables
con respecto al tiempo
title('vertical vs horizontal ');% usamos el title para poner titulo
xlabel('holizontal'); % usamos el xlabel para llamar el eje x
ylabel('vertical'); % usamos el ylabel para llamar y
legend('vertical')% usamos leged para marcar las lineas

Horizon=(t*v0*cos(A)) %cogemos la formula fiica para

Vert=(t*v0*sin(A)-1/2*g*t.^2)%cogemos
subplot(2,2,2); % usamos suplot para unir las graficsa en un mismo cuadro
plot(Horizon,Vert)% usamos plot para graficar nuestras anteriores variables con
respecto al tiempo
title('vertical vs horizontal ');% usamos el title para poner titulo
xlabel('holizontal'); % usamos el xlabel para llamar el eje x
ylabel('vertical'); % usamos el ylabel para llamar y
legend('vertical')% usamos leged para marcar las lineas

Horizon2=(t*v0*cos(A2))
Vert2=(t*v0*sin(A2)-1/2*g*t.^2)%cogemos
subplot(2,2,3) % usamos suplot para unir las graficsa en un mismo cuadro
plot(Horizon2,Vert2)% usamos plot para graficar nuestras anteriores variables
con respecto al tiempo
title('vertical vs horizontal ');% usamos el title para poner titulo
xlabel('holizontal'); % usamos el xlabel para llamar el eje x
ylabel('vertical'); % usamos el ylabel para llamar y
legend('vertical')% usamos leged para marcar
grid

%%

%5.14
p0=1000;
r=0.08;
t=(1:30); % se definen los valors de las variables
P=p0*(exp(r*t))
ingreso=table(P) %creamos una tabla con la formula P
% punto A:
polarplot(t,P) % graficamos convirtiendo a polar
%punto B:
subplot(3,3,1) % une las graficas en un mismo apartado
semilogy(t,P) % pone una parte del plano en unidad log
grid
%punto C:
subplot(3,3,3)
semilogx(t,P)

grid

%%
%5.15
% punto A:
tiempo=(0:45);
trasnsistor=30*(2.^tiempo/1.5); % definimos las variables y hacemos la formula
tabladetiempo=table(tiempo,trasnsistor); % creamos una tabla con la palabra
table
% Punto B:
subplot(2,2,1) % con la palabra suplot unimos las graficas en un mismo apartado
semilogy(tiempo,trasnsistor) %semilog hace que el eje x o y sea en formato log
subplot(2,2,2)
semilogx(tiempo,trasnsistor)
subplot(2,2,3)
loglog(tiempo,trasnsistor) %loglog crea una tabla con formato log

%% %5.17(1)
% Definir el vector de calificaciones
G=[68,83,61,70,75,82,57,5,76,85,62,71,96,78,76,68,72,75,83,93];

% Ordenar el vector de calificaciones

G_ordenado = sort(G);

% Crear la gráfica de barras

bar(G_ordenado);

% Agregar etiquetas al eje x y al eje y

xlabel('Calificaciones');
ylabel('Frecuencia');

%% %5.17(2)
% Definir el vector de calificaciones
G = [68,83,61,70,75,82,57,5,76,85,62,71,96,78,76,68,72,75,83,93];

% Crear histograma
histogram(G);

% Agregar etiquetas al eje x y al eje y

xlabel('Calificaciones');
ylabel('Frecuencia');

%%
%5.18
X=[2,4,8,4,2];%variable
figure
pie(X)
comentario={'A','B','C','D','E'};
legend(comentario,'Location','southoutside','Orientation','horizontal')
pie3(X)
comentario={'A','B','C','D','E'};
legend(comentario,'Location','southoutside','Orientation','horizontal')

%%

%5.22
g=9.8;
t=(0:20);
v0=100;
a=0;
A=pi/4;% las 4 primras lineas elejimos nuestras variables
Vert=(t*v0*sin(A)-1/2*g*t.^2)%cogemos la formula fiica para una vertical
ver=(v0*sin(a)-g*t)
yyaxis right
plot(t,ver)
hold on
plot(t,Vert)

%%
%5.19(1)
% Crear datos
cantidad = [2345 4363 3212 4565 8776 7679 6532 2376 2238 4509 5643 1137 2343
5766 4534 4719 3422 2200 3454 7865 6543 4508 2312 4566];
% Crear gráfico de barras para la primera serie de datos
bar(cantidad)

% Etiquetar los ejes y la leyenda

xlabel('meses')
ylabel('cantidad')

%%

%5.20(1)
% Definir la media y la desviación estándar
media = 70;
desviacion_estandar = 3.5;

% Generar los datos aleatorios

datos = media + desviacion_estandar * randn(1, 1000);

% Crear histograma
histogram(datos);

% Agregar etiquetas al eje x y al eje y

xlabel('Datos');
ylabel('Frecuencia');

%%
%5.26
x=linspace(-5,5,21);
y=x;
[X,Y]= meshgrid(x,y);
z=(sin(sqrt(X.^2+Y.^2)));%la funcion
mesh(z)
hold on %grafica en 3D para ver la funcion con cada limite
contour3(z)
hold off
surf(X,Y,z)
colormap winter
%% %5.24(2)
p = pi; % Definimos p como pi
x = linspace(0, 20*p, 201); % Vector x con 201 elementos
y = x .* sin(x);
polar(x, y)
title('Gráfica polar de x y y')%geneneramos la grafica
%%
%5.24(3)
p = pi; % Definimos p como pi
x = linspace(0, 20*p, 201); % Vector x con 201 elementos
y = x .* sin(x);
z = x .* cos(x);
plot3(x, y,z)%creamos esta funcion ya que es para
xlabel('x')
ylabel('y')
zlabel('z')
title('Gráfica de x y y z')%cramos para que sea vean en la graficada cada una de
ellas

capitulo 4

a =[15,3,22;3,8,5;14,3,82];
b=[1,;5,;6];%estas son variables
c=[12 18 5 2];

d= a(:,3);%cogemos la fila de la columna

f = [b ; d];%cogemos estos vectores
e=[b,d];
g=[a;c(1,1),c(1,2),c(1,3)];%las agremos apara crera un vector bidemensional
h=[a(1,3),c(1,2),b(2,1)]

%%

time=[0:1:19]
sensor1=[70.6432,73.2823,64.1909,67.6970,68.6878,63.9342,63.4028,74.6561,70.0562
,66.7743,74.0286,71.1581,65.0512,76.6979,71.4475,77.3946,75.6901,66.5793,63.5403
,69.6354];
sensor2=[68.3470,65.7819,72.4888,77.4425,67.2676,65.7662,69.7683,73.3151,65.7290
,63.9934,69.4007,69.6735,72.4265,67.0225,69.2517,67.8262,69.6033,77.6758,66.9676
,63.2632];
sensor3=[72.3469,65.4822,70.1794,66.8623,72.6770,2.7644,68.9815,59.7284,70.6628,
77.9647,75.0921,62.0980,69.6067,66.5917,64.8772,63.8282,71.440,67.8535,70.2790,6
8.1606];
sensor4=[67.6751,71.8548,73.6414,80.5608,63.2135,64.8869,75.1892,68.0510,63.0937
,71.5777,77.7662,73.5395,79.7869,72.5227,79.3226,68.3009,64.3011,68.9444,75.9512
,64.4190];%todas la varaibles de los sensores
sensor5=[73.1764,66.9929,72.7559,64.5008,70.4300,59.9772,67.5346,72.3102,68.3950
,76.1828,16.8436,58.3739,63.8418,75.2782,69.4339,71.8961,74.7210,59.3979,66.7766
,66.4785];

valor_maxtimes=max(sensor1)% esto son los numero maximo la sensor

valor_mintimes=min(sensor1)%estos son los minimos de los sensor1

%%
%
times=[0:2:24]
Thermocouple_1=[84.3,86.4,85.2,87.1,83.5,84.8 85.0 85.3 85.3 85.2 82.3 84.7
38.6];
Thermocouple_2=[90.0 89.5 88.6 88.9 88.9 90.4 89.3 89.5 88.9 89.1 89.5 89.4
89.8];
Thermocouple_3=[86.7 87.6 88.3 85.3 80.3 82.4 83.4 85.4 86.3 85.3 89.0 87.3
87.2];
% combine el vector times que creo en la parte a, con los datos de
%thermocouple para crear unamatriz que corresponda a la tabla del problema.
A=[times,Thermocouple_1,Thermocouple_2,Thermocouple_3]

Valor_maxtimes=max(times)
Valor_mintimes=min(times)
%%

%%
disp ('Problemas de dos variables')
disp('ejecrcio 4.5')
a= (1/2)*(2)*2;
b= (1/2)*(4)*3;
C=(1/2)*(6)*4;
d=(1/2)*(8)*(5);%esta en hallar ma intervalos de los numero que nos pierdon
crear un matriz bidemensional con varables hallada

e=(1/2)*(10)*(6);
A=[a,b;C,d;d,e]
%%
disp('ejercicio 4.6')
disp('calculos para el agua')
z=2000/((1000)*(9.81))
x=4000/((1000)*(9.81))
c=6000/((1000)*(9.81));%fromulas para hallar el agua de cada uno
v=8000/((1000)*(9.81));
b=10000/((1000)*(9.81));
disp('calculos para mercurio')
n=1000/((13560)*(9.81));
m=2000/((13560)*(9.81));
l=4000/((13560)*(9.81));
k=6000/((13560)*(9.81));
j=8000/((13560)*(9.81));
h=10000/((13560)*(9.81));

F=[0,z;x,c;v,b]%hallar con resultados un matriz

d=[0,m;l,k;j,h]
%%
disp('La ley del gas ideal Py = RT describe el comportamiento de muchos gases.')
q=(8.31)*100/100E3
w=(8.31)*(200)/100E3
e=(8.31)*(300)/100E3
r=(8.31)*(400)/100E3

%%
%Problema 4.8
%matrices especiales
%cree una matriz de ceros del mismo tamano de las matrices a,b,c
%delproblema4.1
a=[15 3 22;3 8 5;14 3 82];
b=[1;5;6];
c=[12 18 5 2];
Tamano_a=size(a)
Tamano_b=size(b)%tamaño de
Tamano_c=size(c)
%crear una matriz cero del mismo tamano que las matrices a,b,c zeros(3,3)
zeros(3,2)
zeros(1,4)
%%
%Problema 4.9
%cree una matriz magica de 6x6
A=magic(6)
suma_fila=sum(A')%la suma invertida de la matriz la fila
suma_columna=sum(A)%columna la puso en uno
suma_diagonal_1=sum(diag(A))%con hallar la diogonal de suman de la matriz
B=fliplr(A)%la inverte la
suma_diagonal_2=diag(B)
sum(suma_diagonal_2)%la suman digonal la digoanles