Contents

Ejercicio 1

%TABLA FORMATEADA DE LOS NUMEROS 10,20,40,60 Y 80

x=[10;20;40;60;80];
y=[x,log(x)];

fprintf('\n Numero Natural \t log\n'), fprintf('\t%4i\t\t%8.5f\n',y')

 Numero Natural 	 log
	  10		 2.30259
	  20		 2.99573
	  40		 3.68888
	  60		 4.09434
	  80		 4.38203

Ejercicio 2

%HALLAR X PARA EL VECTOR MATRICIAL A:

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

x =

    2.0000
    4.0000
    1.0000

Ejercicio 3

%A PARTIR DEL EJERCICIO ANTERIOR OBTENER LA FACTORIZACION LU Y APLICAR

A=[4 -2 -10;2 10 -12;-4 -6 16];
b=[-10 32 -16]';
[L U P]=lu(A);
M=L*U
x=inv(U)*inv(L)*b

M =

     4    -2   -10
     2    10   -12
    -4    -6    16


x =

     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);
fprintf('\n Autovectores (Columnas de la matriz)\n')
X(:,1)
fprintf('\n Autovalores (Diagonal)\n')
D

 Autovectores (Columnas de la matriz)

ans =

    0.7071
    0.0000
    0.7071


 Autovalores (Diagonal)

D =

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

Ejercicio 5

%PARA EL CIRCUITO DE LA FIGURA CALCULAR V1,V2 Y LA POTENCIA ENTREGADA POR CADA FUENTE

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+003 *

   1.5114 + 0.9092i
   0.6636 + 0.6342i

Ejercicio 6

%FUNCION RECURSIVA PARA LA TORRE DE HANOI, PROBARLA PARA VALOR 5

% function hanoi(n,i,a,f)
% if n>0
%     hanoi(n-1, i, f, a);
%     fprintf('mover disco %d de %c a %c\n', n, i, f);
%     hanoi(n-1, a, i, f);
% end

hanoi (5,'a','b','c')

mover disco 1 de a a c
mover disco 2 de a a b
mover disco 1 de c a b
mover disco 3 de a a c
mover disco 1 de b a a
mover disco 2 de b a c
mover disco 1 de a a c
mover disco 4 de a a b
mover disco 1 de c a b
mover disco 2 de c a a
mover disco 1 de b a a
mover disco 3 de c a b
mover disco 1 de a a c
mover disco 2 de a a b
mover disco 1 de c a b
mover disco 5 de a a c
mover disco 1 de b a a
mover disco 2 de b a c
mover disco 1 de a a c
mover disco 3 de b a a
mover disco 1 de c a b
mover disco 2 de c a a
mover disco 1 de b a a
mover disco 4 de b a c
mover disco 1 de a a c
mover disco 2 de a a b
mover disco 1 de c a b
mover disco 3 de a a c
mover disco 1 de b a a
mover disco 2 de b a c
mover disco 1 de a a c

Ejercicio 7

%POLINOMIO DE ORDEN 2 QUE MEJOR AJUSTE LOS DATOS DE LA TABLA

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);
plot(x,y,'*',x,yc);
xlabel('x'),ylabel('y'),grid,title('Ajuste polinómico')
legend('Datos','Ajuste polinómico',2)

p =

    4.0233    2.0107    9.6783

Ejercicio 8

%GRAFICAR LAS SIGUIENTES FUNCIONES A PARTIR DE LOS DATOS DADOS:

omegat=0:0.05:3*pi;
v=120*sin(omegat);
i=100*sin(omegat-(pi/4));
subplot(2,2,1)
plot(omegat,v,omegat,i)
title('Gráfica Tensión e Intensidad'),xlabel('\omegat(radianes)')

p=v.*i;
subplot(2,2,2)
plot(omegat,p)
title('Potencia'),xlabel('\omegat (radianes)'),ylabel('watios')
Fm=3.0;
fa=Fm*sin(omegat);
fb=Fm*sin(omegat-2*pi/3);
fc=Fm*sin(omegat-4*pi/3);
subplot(2,2,3)
plot(omegat,fa,omegat,fb,omegat,fc)
title('Fm trifasico'),xlabel('\omegat (radianes)')

subplot(2,2,4)
fr=3.0;
plot(-fr*cos(omegat),fr*sin(omegat))
title('Radio fr')

Ejercicio 9

%GRAFICAR CURVA PARAMÉTRICA

t=0:0.1:16*pi;
X=exp(-0.03*t).*cos(t);
Y=exp(-0.03*t).*sin(t);
Z=t;
subplot(1,1,1)
plot3(X,Y,Z), axis off

Ejercicio 10

%GRAFICAR LA CURVA

[x,y]=meshgrid(-4:0.3:4);
z=sin(x).*cos(y).*exp(-(x.^2+y.^2).^0.5);
mesh(x,y,z)

Ejercicio 11

%HALLAR LAS RAÍCES DEL POLINOMIO

p=[1 0 -35 50 24];
r=roots(p)

r =

   -6.4910
    4.8706
    2.0000
   -0.3796

Ejercicio 12

%RESOLVER LA ECUACIÓN DIFERENCIAL

% function y = HalfSine(t, y, z)
% h = sin(pi*t/5).*(t<=5);
% y = [y(2); -2*z*y(2)-y(1)+h];
Ejemploode
[t, yy] = ode45(@HalfSine, [0 35], [1 0], [ ], 0.15);
plot(t, yy(:,1))

Ejercicio 13

%GRÁFICA DE LA SEÑAL EN EL TIEMPO Y GRÁFICA ESPECTRAL DE AMPLITUDES

figure
k = 5; m = 10; f0 = 10;B0 = 2.5; N = 2^m; T = 2^k/f0;
ts = (0:N-1)*T/N; df = (0:N/2-1)/T;

SampledSignal = B0*sin(2*pi*f0*ts)+B0/2*sin(2*pi*f0*2*ts);
An=abs(fft(SampledSignal,N))/N;
subplot(4,2,1)
plot(ts,SampledSignal)
subplot(4,2,2)
plot(df,2*An(1:N/2))

SampledSignal1 = exp(-2*ts).*sin(2*pi*f0*ts);
An1=abs(fft(SampledSignal1,N))/N;
subplot(4,2,3)
plot(ts,SampledSignal1)
subplot(4,2,4)
plot(df,2*An1(1:N/2))

SampledSignal2 = sin(2*pi*f0*ts+5*sin(2*pi*f0/10*ts));
An2=abs(fft(SampledSignal2,N))/N;
subplot(4,2,5)
plot(ts,SampledSignal2)
subplot(4,2,6)
plot(df,2*An2(1:N/2))

SampledSignal3 = sin(2*pi*f0*ts-5*exp(-2*ts));
An3=abs(fft(SampledSignal3,N))/N;
subplot(4,2,7)
plot(ts,SampledSignal3)
subplot(4,2,8)
plot(df,2*An3(1:N/2))

Ejercicio 14

%GRAFICAR WINDTUNNEL

subplot(1,1,1)
A = imread('WindTunnel.jpg', 'jpeg');
image(A)
hold on
figure
subplot(2,1,1)
r= A(200, :, 1);
plot(r, 'r');
hold on
subplot(2,1,2)
hist(r,0:15:255);

Ejercicio 15

%GRAFICAR EN COORDENADAS POLARES

t = linspace(-pi,pi,300);
r = 2-4*cos(t);
figure
polar(t,r);
title('Gráfico polar de r=2-4*cos(t)')

Published with MATLAB® 7.13