TAREA MATLAB

Contents

EJERCICIO 1

apartado a)

type factr

function x = factr(n)
    if n==0
        x=1;
    else
        x = n * factr(n-1);
    
end

apartado b)

type combina

function x=combina(n,i)
    x=factr(n)/(factr(i)*factr(n-i));
end

apartado c)

type bernstein

function x=bernstein(n,i,t)
    x=combina(n,i)*t.^i.*(1-t).^(n-i);
end

apartado d)

t=linspace(0,1);
v=[1 2 4 4.6;1 3 -1 1.5];
plot(v(1,:),v(2,:),'-o');
n=size(v);
n=n(2);
s=size(t);
x=zeros(n,s(2));
y=zeros(n,s(2));
for i=1:n
    x(i,:)=bernstein(n-1,i-1,t)*v(1,i);
    y(i,:)=bernstein(n-1,i-1,t)*v(2,i);
end
a=sum(x);
b=sum(y);
hold on
plot(a,b)

% grafico
t=linspace(0,1);
n=3;
for i=0:n
    b=bernstein(n,i,t);
    plot(t,b);
    hold on;
end
xlabel('t');
ylabel('Polinomio de Bernstein')
title('Polinomio de Bernstein grado 3')
legend('B0,0','B3,1','B3,2','B3,3');

EJERCICIO 2

Apartado a)

[Y]=xlsread('sotaventogaliciaanual.xlsx','B:B');
hist(Y)
title('Histograma de frecuencias')

Apartado b)

Y=xlsread('sotaventogaliciaanual.xlsx','B:B');
x=0.5:1:max(Y);
viento=hist(Y,x);
frec=viento/sum(viento);
f=@(a,x) (a(1)/a(2))*((x/a(2)).^(a(1)-1)).*exp(-(x/a(2)).^a(1));
k=std(Y);
c=mean(Y);
a0=[k c];  %valor inicial
af=nlinfit(x,frec,f,a0)
hold on
%diagrama de frecuencias
bar(x,frec,'y');
%curva de ajuste
x=linspace(0,max(Y),100);
y=f(af,x);
plot(x,y,'b')
title('Ajuste a la función Weibull')
xlabel('Velocidad')
ylabel('Frecuencia')
hold off

af =

    2.3849    6.0208

Apartado c)

X=xlsread('sotavento_curva potencia.xlsx','A:A');
Y=xlsread('sotavento_curva potencia.xlsx','B:B');
title('Curva de potencia interpolada')
xq=linspace(0,25);
vq1 = interp1(X,Y,xq,'pchip');
plot(X,Y,'o',xq,vq1,':.');
title('Interpolación curva de potencia');
legend('Valores de la curva de potencia','Interpolación','Location','southeast');
grid on

EJERCICIO 3

%mD2x(t)+cDx(t)+kx
m=20;
k=20;
c=5;c1=40;c2=200;
x0=[1,0];%1 es el desplazamiento inicial y cero la velocidad inicial
tf=40; %Tiempo final de 40
f=@(t,x) [x(2);-c*x(2)/m-k*x(1)/m];
[t,x]=ode45(f,[0,tf],x0);
plot(t,x(:,1))
grid on
xlabel('t')
ylabel('x')
title('sistema masa-resorte-amortiguador')
hold on
f1=@(t,x) [x(2);-c1*x(2)/m-k*x(1)/m];
[t,x]=ode45(f1,[0,tf],x0);
plot(t,x(:,1))
f2=@(t,x) [x(2);-c2*x(2)/m-k*x(1)/m];
[t,x]=ode45(f2,[0,tf],x0);
plot(t,x(:,1))

Published with MATLAB® R2017a