Contents
clear all
clc
Ejercicio 2
A=[4 -2 -10;2 10 -12;-4 -6 16]
B=[-10;32;-16]
X=A\B
A =
4 -2 -10
2 10 -12
-4 -6 16
B =
-10
32
-16
X =
2.0000
4.0000
1.0000
Ejercicio 4
A=[0 1 -1;-6 -11 6;-6 -11 5];
[auto_vectores auto_valores]=eig(A)
auto_vectores =
0.7071 -0.2182 -0.0921
0.0000 -0.4364 -0.5523
0.7071 -0.8729 -0.8285
auto_valores =
-1.0000 0 0
0 -2.0000 0
0 0 -3.0000
Ejercicio 5
Y=[1.5-2j -.35+1.2j;-.35+1.2j 0.9-1.6j];
I=[30+40j;20+15j]
V=Y\I
S=V.*conj(I)
I =
30.0000 +40.0000i
20.0000 +15.0000i
V =
3.5902 +35.0928i
6.0155 +36.2212i
S =
1.0e+03 *
1.5114 + 0.9092i
0.6636 + 0.6342i
Ejercicio 6
Hanoi(5,'A','B','C')
Mueva el disco 1 de A a C
Mueva el disco 2 de A a B
Mueva el disco 1 de C a B
Mueva el disco 3 de A a C
Mueva el disco 1 de B a A
Mueva el disco 2 de B a C
Mueva el disco 1 de A a C
Mueva el disco 4 de A a B
Mueva el disco 1 de C a B
Mueva el disco 2 de C a A
Mueva el disco 1 de B a A
Mueva el disco 3 de C a B
Mueva el disco 1 de A a C
Mueva el disco 2 de A a B
Mueva el disco 1 de C a B
Mueva el disco 5 de A a C
Mueva el disco 1 de B a A
Mueva el disco 2 de B a C
Mueva el disco 1 de A a C
Mueva el disco 3 de B a A
Mueva el disco 1 de C a B
Mueva el disco 2 de C a A
Mueva el disco 1 de B a A
Mueva el disco 4 de B a C
Mueva el disco 1 de A a C
Mueva el disco 2 de A a B
Mueva el disco 1 de C a B
Mueva el disco 3 de A a C
Mueva el disco 1 de B a A
Mueva el disco 2 de B a C
Mueva el disco 1 de A a C
Ejercicio 7
x=0:0.5:5;
y=[10 10 16 24 30 38 52 68 82 96 123];
p=polyfit(x,y,2)
yc=polyval(p,x);
figure(1);
plot(x,y,'*',x,yc);
xlabel('x'),ylabel('y'),grid,title('Ajuste polinomico')
legend('Datos ','Ajuste polinomico',2)
p =
4.0233 2.0107 9.6783
Ejercicio 8
omt=0:0.05:3*pi;
v=120*sin(omt);
i=100*sin(omt-(pi/4));
figure(2);
subplot(2,2,1)
plot(omt,v,omt,i)
title('Grafica Intensidad y Tensión'),xlabel('\omt(radianes)')
p=v.*i;
subplot(2,2,2)
plot(omt,p)
title('Potencia'),xlabel('\omt(radianes)'),ylabel('watios')
Fm=3.0;
fa=Fm*sin(omt);
fb=Fm*sin(omt-2*pi/3);
fc=Fm*sin(omt-4*pi/3);
subplot(2,2,3)
plot(omt,fa,omt,fb,omt,fc)
title('Trifasico'),xlabel('\omt(radianes)')
subplot(2,2,4)
fr=3.0;
plot(-fr*cos(omt),fr*sin(omt))
title('Radio fr')
Ejercicio 13
figure(3);
k = 5; m = 10; fo = 10; Bo = 2.5; N = 2^m; T = 2^k/fo;
t = (0:N-1)*T/N; df = (0:N/2-1)/T;
g1 = Bo*sin(2*pi*fo*t)+Bo/2*sin(2*pi*2*fo*t);
An = abs(fft(g1, N))/N;
subplot(2,1,1)
plot(t, g1)
title('g_1(t)=B_o·sin(2·\pi·f_o·t) + B_o/2·sin(2\pi·2f_o·t)');
subplot(2,1,2);
plot(df, 2*An(1:N/2))
title('Amplitud espectral de g_1(t) en función de f');
figure(4);
g2 = exp(-2*t).*sin(2*pi*fo*t);
An1 = abs(fft(g2, N))/N;
subplot(2,1,1)
plot(t, g2)
title('g_2(t)=e^{-2t}·sin(2\pi f_ot)');
subplot(2,1,2)
plot(df, 2*An1(1:N/2))
title('Amplitud espectral de g_2(t) en función de f');
figure(5);
g3 = sin(2*pi*fo*t + 5*sin(2*pi*fo/10*t));
An2 = abs(fft(g3, N))/N;
subplot(2,1,1)
plot(t, g3)
title('g_3(t)=sin(2\pi f_ot + 5 sin(2\pi 2 f_o/10t))');
subplot(2,1,2)
plot(df, 2*An2(1:N/2))
title('Amplitud espectral de g_3(t) en función de f');
figure(6);
g4 = sin(2*pi*fo*t - 5*exp(-2*t));
An3 = abs(fft(g4, N))/N;
subplot(2,1,1)
plot(t, g4)
title('g_4(t)=sin(2\pi f_o t - 5 e^{-2t})');
subplot(2,1,2)
plot(df, 2*An3(1:N/2))
title('Amplitud espectral de g_4(t) en función de f');
Ejercicio 14
v = imread('WindTunnel.jpg');
figure(7), image(v), image(v), axis image off
fila=200;
figure(8);
R = v(fila, :, 1);
s1=subplot(2,1,1);s1.YGrid='on';
stem(R,'Marker','none','LineWidth',2,'Color',[1 0 0]);
xlim([0 275]);ylim([0 255]);
title(['Niveles de color rojo en la fila ',num2str(fila)]);
xlabel('Ancho de la imagen (px)')
s2=subplot(2,1,2);
hist(R,0:5:255);xlim([0 275]);
s2.XTick=0:25:275;
title(['Histograma del color rojo en la fila ',num2str(fila)]);
xlabel('Ancho de la imagen (px)')