%%*Tarea Matlab*


%%Ejercicio 1

x = [10 20 40 60 80];
y = [x; log(x); log10(x)];
fileID = fopen('ej1.txt','w');
fprintf(fileID,'%5s %12s %12s\r\n','x','log(x)','log10(x)');
fprintf(fileID,'%5d %12.4f %12.6f\r\n\r\n',y);
fclose(fileID);


%%Ejercicio 2

function IMC
i = 1;
while i > 0
p = input('Introduzca su peso en kg: ');
a = input('Introduzca su altura en m: ');
imc = p/(a^2);
disp(['Su IMC es: ' num2str(imc,5)])
if imc < 16
    disp('Infrapeso: Delgadez Severa')
elseif 16 <= imc && imc <= 16.99
    disp('Infrapeso: Delgadez moderada')
elseif 17 <= imc && imc <= 18.49
    disp('Infrapeso: Delgadez aceptable')
elseif 18.5 <= imc && imc <= 24.99
     disp('Peso normal')
elseif 25 <= imc && imc <= 29.99
     disp('Sobrepeso')
elseif 30 <= imc && imc <= 34.99
     disp('Obeso: Tipo I')
elseif 35 <= imc && imc <= 40
     disp('Obeso: Tipo II')
else
     disp('Obeso: Tipo III')
end
i = i+1;
cont = input('¿Desea continuar? Y/N ', 's');
if cont == 'n'
    break
elseif cont == 'N'
    break
end
end


%%Ejercicio 3

function base
dec = 0:20;
mat = [dec; str2num(dec2bin(dec))];
fileID = fopen('ej3.txt','w');
fprintf(fileID,'%5s %5s\r\n','Decimal','Binario');
fprintf(fileID,'%5.3f %5.3f\r\n',mat);
fclose(fileID);


%%Ejercicio 5

%Matriz de resistencias
R = [100 0 0 0 -1 1 0; 0 200 0 0 0 -1 1; 0 0 50 0 -1 0 0; 0 0 0 100 0 -1 0; 0 300 0 0 0 0 -1; 1 0 1 0 0 0 0; 1 -1 0 -1 0 0 0];
%Matriz de temperaturas
T = [0; 0; -50; -50; -50; 25 ;0];
Rinv = inv(R);
%Resultados
Q = Rinv*T


%%Ejercicio 6

R = [100 0 0 0 -1 1 0; 0 200 0 0 0 -1 1; 0 0 50 0 -1 0 0; 0 0 0 100 0 -1 0; 0 300 0 0 0 0 -1; 1 0 1 0 0 0 0; 1 -1 0 -1 0 0 0];
%Número de condición
NumC = rcond(R)
%Normas 1 y 2
Norm1 = norm(R,1)
Norm2 = norm(R
%Inf
%Inf = Inf(R)
%Frobenius
Fro = norm(R,'fro')
%LU
[L,U] = lu(R)
%Resolución
X = inv(U)*inv(L)*[0;0;-50;-50;-50;25;0]


%%Ejercicio 7

A = [0 1 -1; -6 -11 6; -6 -11 5];
%Autovectores y autovalores
[AuVec AuVal] = eig(A)


%%Ejercicio 8

A11=complex(1.5,-2);
A12=complex(-0.35,1.2);
A21=complex(-0.35,1.2);
A22=complex(0.9,1.6);
A=[A11 A12;A21 A22];
B11=complex(30,40);
B21=complex(20,15);
Ainv = inv(A);
B = [B11; B21];
X = Ainv*B


%%Ejercicio 11

[X,Y]=meshgrid(-4:0.3:4,-4:0.3:4);
Z=sin(X).*cos(Y).*(-exp((X.^2)+(Y.^2)).^0.5);
surf(X,Y,Z)
title('Ejercicio 11')
xlabel('x')
ylabel('y')


%%Ejercicio 12

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


%%Ejercicio 13

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;
g1 = Bo*sin(2*pi*fo*ts)+Bo/2*sin(2*pi*fo*2*ts);
g2 = exp(-2*ts).*sin(2*pi*fo*ts);
g3=sin(2*pi*fo*ts+5*sin(2*pi*(fo/10)*ts));
g4=sin(2*pi*fo*ts-5*exp(-2*ts));
An1 = abs(fft(g1, N))/N;
subplot(2,2,1); plot(df, 2*An1(1:N/2))
title('G1')
An2 = abs(fft(g2, N))/N;
subplot(2,2,2); plot(df, 2*An2(1:N/2))
title('G2')
An3 = abs(fft(g3, N))/N;
subplot(2,2,3); plot(df, 2*An3(1:N/2))
title('G3')
An4 = abs(fft(g4, N))/N;
subplot(2,2,4); plot(df, 2*An4(1:N/2))
title('G4')


%%Ejercicio 14

alpha=-pi:0.1;pi;
r=2-4*cos(alpha);
figure
ezpolar('r')

Error: File: C:\Users\Fernando\Desktop\Untitled.m Line: 17 Column: 1
Function definitions are not permitted in this context.

Published with MATLAB® R2013b