ejercicio 1

Contents

apartado d

Están en este fichero las funciones de los apartados a b c porque son necesarias para ejecutar el d. Las podría haber puesto fuera del fichero pero para evitar problemas con los directorios las reescribo dentro del apartado.

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)

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


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

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

Published with MATLAB® R2017a