Edgar, prueba definiendo el archivo pentaphi en el editor de matlab, siguiendo la secuencia:
file->new-> mfile
function [x,y]=pentaphi(r,s,phi)
%Primero trazamos una circunferencia de radio r
%donde sera inscrito el pentagono original
x=-r:.1:r; %Se define el dominio de las x
yp=sqrt(r^2-x.^2); %Se define la ecuacion que determina a la circunferencia
yn=-sqrt(r^2-x.^2);
plot(x,yp) %Graficamos la curva
hold on
plot(x,yn)
grid, axis equal
pause(3)
hold on
theta=360/5 %Encontramos el angulo que determina el lado del pentagono
t1=deg2rad(90-theta) %Usamos un angulo auxiliar, en radianes, para las abscisas
t2=deg2rad(theta)
x=[0 r*cos(t1) r*cos(t2/2-pi/2) -r*cos(t2/2-pi/2) -r*cos(t1) 0]
y=[r r*sin(t1) r*sin(t2/2-pi/2) r*sin(t2/2-pi/2) r*sin(t1) r]
line(x,y) %trazamos el pentagono
pause(3)
r=r*s %redimensionamos con el factor s
x=-r:.1:r; %Se redefine el dominio
yp=sqrt(r^2-x.^2); %Se define la ecuacion que determina a la circunferencia modificada
yn=-sqrt(r^2-x.^2);
plot(x,yp) %Graficamos la curva
hold on
plot(x,yn)
pause(3)
theta=360/5 %Encontramos el angulo que determina el lado del nuevo pentagono
t1=deg2rad(90-theta) %Usamos un angulo auxiliar, en radianes, para las abscisas
t2=deg2rad(theta)
x=[0 r*cos(t1) r*cos(t2/2-pi/2) -r*cos(t2/2-pi/2) -r*cos(t1) 0]
y=[r r*sin(t1) r*sin(t2/2-pi/2) r*sin(t2/2-pi/2) r*sin(t1) r]
line(x,y)
pause(3)
R=[cos(phi) sin(phi);sin(-phi) cos(phi)]; %Usamos una matriz de rotacion R en un angulo phi a favor de las manecillas del reloj
x=[0 r*cos(t1) r*cos(t2/2-pi/2) -r*cos(t2/2-pi/2) -r*cos(t1) 0]
y=[r r*sin(t1) r*sin(t2/2-pi/2) r*sin(t2/2-pi/2) r*sin(t1) r]
P=R*[x;y]
plot(P(1,:),P(2,:))
Al final, debes escoger un valor de r, por ejemplo r=5, un valor de s=2 y el valor de phi=35*pi/180 y teclear:
>>pentaphi(5,2,35*pi/180)
para obtener lo deseado.
Cualquier duda que tengas puedes escribirme a mi dirección de correo o a la alterna
[email protected]