TAREA MATLAB

Contents

EJERCICIO 1

apartado a)

type factr

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

apartado b)

type combina

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

apartado c)

type bernstein

function ber=bernstein(n,i,t)
    ber=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));
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('Polinomios de Bernstein de grado 3')
legend('B_3_,_0','B_3_,_1','B_3_,_2','B_3_,_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,'c');
%curva de ajuste
x=linspace(0,max(Y),100);
y=f(af,x);
plot(x,y,'r')
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® R2015b