Bienvenido a la web de Vicente Canencia

%Ejercicio 2
A=[4 -2 -10; 2 10 -12; -4 -6 16];
B=[-10; 32; -16];
q=A/B;


%Ejercicio 4
C=[4 -2 -10; 2 10 -12; -4 -6 16];
[X,D]=eig(C);

%Ejercicio 5
Y=[1.5-20j -3.5+1.2j;-0.35-1.2j 0.9-1.6j];
I=[30+40j;20+35j];
disp('Solucion');
V=Y/I;
S=V.*conj(Y);

%Ejercicio 6
%Resuelve el problema de Hanoi utilizando una funcion recursiva
%para valor de 5 discos
%function hanoi(n, k, a, f)
  %  if n > 0
  %  hanoi(n-1, k, f, a);
  %  fprintf('mover disco %d de %c a %c\n', n, i, f);
  %  hanoi(n-1, a, k, f);
  %  end
%end
hanoi(5,'a','b','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)
plot(x,y,'x',x,yc)
xlabel('x'),ylabel('y'),grid,title('Ajuste polinomico')
legend('Datos',('Ajuste polinomico'),4)


%Ejercicio 8

x=[0:0.05:3*pi];
figure;
% Graficar v = 120 seno wt e i = 100 seno(wt - p/4 ) en función de wt en la parte superior izquierda
subplot(2,2,1);
plot(x,120*sin(x));
hold on;
plot(x, 100*sin(x-pi*.25));
% Graficar p = vi en la parte superior derecha
subplot(2,2,2);
plot(x,120*sin(x)*100.*sin(x-(pi*.25)));
% Para Fm = 3.0, graficar fa = Fm seno wt, fb = Fm seno(wt – 2 p/3) y fc =
% Fm seno(wt– 4 p/3) en función de wt en la parte inferior izquierda
subplot(2,2,3);
plot(x,3*sin(x));
hold on;
plot(x,3*sin(x-2*(pi/3)));
hold on;
plot(x,3*sin(x-4*(pi/3)));
% Para fR = 3.0, construir un círculo de radio fR en la parte inferior derecha
subplot(2,2,4);
ang=0:0.01:2*pi;
xp=3*cos(ang);
yp=3*sin(ang);
plot(xp,yp);


%Ejercicio 11
%f=x^4-35*x^2+50*x+24
f=[1 0 -35 50 24]
sol=roots(f)

%Ejercicio 12
[t, h] = ode45(@HalfSine, [0 35], [1 0], [], 0.15);
figure;
plot(t, h(:,1))

%Ejercicio 13
% Base ejemplo de la transformada de Fourier
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;

% Primer caso
SampledSignal = Bo*sin(2*pi*fo*ts)+Bo/2*sin(2*pi*fo*2*ts);
An = abs(fft(SampledSignal, N))/N;
plot(df, 2*An(1:N/2));
hold on

% Segundo caso
SampledSignal = exp(-2*ts).*sin(2*pi*fo*ts);
An = abs(fft(SampledSignal, N))/N;
plot(df, 2*An(1:N/2))
hold on
    % Tercer caso
SampledSignal = sin(2*pi*fo*ts+sin(2*pi*fo*.1*ts));
An = abs(fft(SampledSignal, N))/N;
plot(df, 2*An(1:N/2));
hold on

% Cuarto caso
SampledSignal = sin(2*pi*fo*ts-5*exp(-2*ts));
An = abs(fft(SampledSignal, N))/N;
plot(df, 2*An(1:N/2))
hold on


%Ejercicio 14
v = imread('WindTunnel.jpg');
image(v)
figure;
row = 2;
red = v(row, :, 1);
gr  = v(row, :, 2);
bl  = v(row, :, 3);
plot(red, 'r');
hold on
plot(gr, 'g');
plot(bl, 'b');

Error using /
Matrix dimensions must agree.

Error in tarea_matlab (line 4)
q=A/B;

Published with MATLAB® R2015a