Contents

Ejercicio 2

A=[4 -2 -10;2 10 -12;-4 -6 16];
C=[-10 32 -16]';
x=A\C

%%Comprobación:
comprobacion=A*x

x =

    2.0000
    4.0000
    1.0000


comprobacion =

   -10
    32
   -16

Ejercicio 4

%V,D] = eig(A) devuleve una matriz diagonal D de autovalores y una matriz V
%cuyas columnas son los correspondientes autovectores
A=[0 1 -1;-6 -11 6;-6 -11 5];
[V,D]=eig(A)

%tal que A*V = V*D

%Comprobación:

matriz1=A*V
matriz2=V*D

V =

    0.7071   -0.2182   -0.0921
    0.0000   -0.4364   -0.5523
    0.7071   -0.8729   -0.8285


D =

   -1.0000         0         0
         0   -2.0000         0
         0         0   -3.0000


matriz1 =

   -0.7071    0.4364    0.2762
   -0.0000    0.8729    1.6570
   -0.7071    1.7457    2.4856


matriz2 =

   -0.7071    0.4364    0.2762
   -0.0000    0.8729    1.6570
   -0.7071    1.7457    2.4856

Ejercicio 5

Y=[1.5-2j -.35+1.2j; -.35+1.2j 0.9-1.6j];
I=[30+40j;20+15j];

%Tensiones de los nodos V1 y V2
V=Y\I

%Potencia de la fuente
S=V.*conj(I)

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

% Vector con datos de 0 a 5 con paso de 0.5
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);

%Obtenemos el polinomio: 4.0233x^2 +2.0107x + 9.6783=0

plot(x,y,'or',x,yc);
xlabel('x'),ylabel('y'),grid,title('Ajuste polinomico')
legend('Datos ','Ajuste polinomico')

p =

    4.0233    2.0107    9.6783

Ejercicio 8

wt=0:0.05:3*pi;
v=120*sin(wt);
i=100*sin(wt-pi/4);
p=v.*i;
Fm=3.0;
fa=Fm*sin(wt);
fb=Fm*sin(wt-2*pi/3);
fc=Fm*sin(wt-4*pi/3);
fR= 3*Fm/2;

scrsz = get(groot,'ScreenSize');
figure('Position',[1 scrsz(4)/2 scrsz(3)/2 scrsz(4)/2])

subplot(2,2,1)
plot(wt,v,wt,i)
title('Tensión e Intensidad')
ylabel('v, A'),xlabel('\omegat, rad')

subplot(2,2,2)
plot(wt,p)
title('Potencia')
ylabel('Potencia, W'),xlabel('\omegat, rad')

subplot(2,2,3)
plot(wt,fa,wt,fb,wt,fc)
title('CA sinusoidal')
xlabel('\omegat, rad')

subplot(2,2,4)
plot(fR*cos(wt),fR*sin(wt))
title('Circulo de radio f_R'),axis equal

Ejercicio 13

figure('Position',[1 1 scrsz(3)/2 scrsz(4)/2])
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('Position',[1 1 scrsz(3)/2 scrsz(4)/2])
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('Position',[1 1 scrsz(3)/2 scrsz(4)/2])
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('Position',[1 1 scrsz(3)/2 scrsz(4)/2])
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

A = imread('WindTunnel.jpg', 'jpg');

figure();image(A), axis image off

scrsz = get(groot,'ScreenSize');
figure('Position',[1 scrsz(4)/2 scrsz(3)/2 scrsz(4)/2])
R = A(200, :, 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(200)]);
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(200)]);
xlabel('Ancho de la imagen (px)')

Published with MATLAB® R2017a