jdg.m

matlab没有提供画椭圆的函数

function jdg(xc,yc,a,b,k)
% 画个椭园)
% (xc,yc) is the position of center
% a is long radio
% b is short radio
% k is the position of the angle
% Example:
% jdg(0,0,3,5,pi/6)
t1=0:.02:pi;
t2=pi:.02:2*pi;
z1=exp(i*t1);
z2=exp(i*t2);
z1=(a*real(z1)+i*b*imag(z1))*exp(i*k);
z2=(a*real(z2)+i*b*imag(z2))*exp(i*k);
z1=z1+xc+yc*i;
z2=z2+xc+yc*i;
plot(z1,'r')
hold on
plot(z2,'r')

hold off
grid on