Hallar el vector X para una ecuación matricial:
A=[4 -2 -10;2 10 -12;-4 -6 16];
b=[-10 32 -16]';
x=A\b
%%Otra opcion
x1=inv(A)*b
x =
2.0000
4.0000
1.0000
x1 =
2
4
1
Ejercicio 4
Hallar los autovalores y autovectores de la matriz A:
A=[0 1 -1;-6 -11 6;-6 -11 5];
[X,D] = eig(A)
T1= A*X
T2= X*D
X =
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
T1 =
-0.7071 0.4364 0.2762
-0.0000 0.8729 1.6570
-0.7071 1.7457 2.4856
T2 =
-0.7071 0.4364 0.2762
-0.0000 0.8729 1.6570
-0.7071 1.7457 2.4856
Ejercicio 5
Para el siguiente circuito, determinar los voltajes de los nodos V1 y V2 y la potencia entregada por cada fuente:
Y=[1.5-2j -0.35+1.2j;-0.35+1.2j 0.9-1.6j];
I=[30+40j;20+15j];
V=Y\I;
S=V.*conj(I);
%El potencial V1 en voltios vale:
V1=V(1)
V1 =
3.5902 +35.0928i
%El potencial V2 en voltios vale:
V2=V(2)
V2 =
6.0155 +36.2212i
%La potencia entregada por la fuente 1 en voltiamperios (VA) es:
S1=S(1)
S1 =
1.5114e+03 + 9.0918e+02i
%La potencia entregada por la fuente 2 en voltiamperios (VA) es:
S2=S(2)
S2 =
6.6363e+02 + 6.3419e+02i
Ejercicio 7
Ajustar un polinomio de orden 2 a los siguientes datos:
x= 0:0.5:5;
y=[10 10 16 24 30 38 52 68 82 96 123];
coeficientes=polyfit(x,y,2)
yc=polyval(coeficientes,x);
plot(x,y,'o',x,yc,'r');
xlabel('x'),ylabel('y'),grid,title('Ajuste polinomico de orden 2')
legend('Datos','Ajuste polinomico')
coeficientes =
4.0233 2.0107 9.6783
Ejercicio 8
Partir la ventana Figure en cuatro particiones (2x2) y graficar las siguientes funciones para wt de 0 a 3pi en pasos de 0.05
- Graficar v = 120 seno wt e i = 100 seno(wt - pi/4 ) en función de wt en la parte superior izquierda
- Graficar p = vi en la parte superior izquierda
- Para Fm = 3.0, graficar fa = Fm seno wt, fb = Fm seno(wt – 2 pi/3) y fc = Fm seno(wt – 4 pi/3) en función de wt en la parte inferior izquierda
- Para fR = 3.0, construir un círculo de radio fR en la parte inferior derecha
wt=0:0.05:3*pi;
figure(1)
subplot(2,2,1)
v=120*sin(wt);
i=100*sin(wt-pi/4);
plot(wt,v,wt,i)
title('Tensión e Intensidad'),xlabel('\omegat, radianes')
grid on
subplot(2,2,2)
p=v.*i;
plot(wt,p)
title('Potencia'),xlabel('\omegat, radianes')
grid on
subplot(2,2,3)
Fm=3.0;
fa=Fm*sin(wt);
fb=Fm*sin(wt-2*pi/3);
fc=Fm*sin(wt-4*pi/3);
plot(wt,fa,wt,fb,wt,fc)
title('Fm trifasico'),xlabel('\omegat, radianes')
grid on
subplot(2,2,4)
fR= 3/2*Fm;
plot(-fR*cos(wt),fR*sin(wt))
title('Fm rotante')
grid on
Ejercicio 13
Tomando como base las condiciones del ejemplo de la transformada de Fourier de los apuntes (pág. 124), graficar para las siguientes señales la gráfica de la señal en el tiempo y la gráfica de la amplitud espectral en función de la frecuencia:
figure(5)
k = 5; m = 10; fo = 10; Bo = 2.5; N = 2^m; T = 2^k/fo;
ts = (0:N-1)*T/N;
df = (0:N/2-1)/T;
SampledSignal = Bo*sin(2*pi*fo*ts)+Bo/2*sin(2*pi*fo*2*ts);
An = abs(fft(SampledSignal, N))/N;
subplot(4,2,1)
plot(ts, SampledSignal)
grid on
subplot(4,2,2)
plot(df, 2*An(1:N/2))
grid on
SampledSignal1 = exp(-2*ts).*sin(2*pi*fo*ts);
An1 = abs(fft(SampledSignal1, N))/N;
subplot(4,2,3)
plot(ts, SampledSignal1)
grid on
subplot(4,2,4)
plot(df, 2*An1(1:N/2))
grid on
SampledSignal2 = sin(2*pi*fo*ts + 5*sin(2*pi*fo/10*ts));
An2 = abs(fft(SampledSignal2, N))/N;
subplot(4,2,5)
plot(ts, SampledSignal2)
grid on
subplot(4,2,6)
plot(df, 2*An2(1:N/2))
grid on
SampledSignal3 = sin(2*pi*fo*ts - 5*exp(-2*ts));
An3 = abs(fft(SampledSignal3, N))/N;
subplot(4,2,7)
plot(ts, SampledSignal3)
grid on
subplot(4,2,8)
plot(df, 2*An3(1:N/2))
grid on