function [theta rho X]=RosesConfig(a,k,delt,k1,k2)

%
%õ߷̵ıͼ
%rho=a*sin(k*theta)
%rho=a*cos(k*theta)
%kΪΪʱõҶӵĸΪkΪpi
%kΪΪżʱõҶӵĸΪ2*kΪ2*pi
%kΪN/DNĸDΪʱõҶΪNΪD*piNĸDһΪżʱõҶΪ2*NΪ2*D*pi

if k==fix(k)
   if k/2 == fix(k/2)
    % kΪΪżʱõҶӵĸΪ2*kΪ2*piõ߷rho=10*cos(4*theta)Ϊͼ
    theta=delt:delt:2*pi;
    rho=a*cos(k*theta);
   else
    %%%kΪΪʱõҶӵĸΪkΪpiõ߷rho=10*cos(3*theta)Ϊͼ
    theta=delt:delt:pi;
    rho=a*cos(k*theta);
   end
else
    if k1/2 ~= fix(k1/2)
        if k2/2 ~= fix(k2/2)
            theta=delt:delt:k2*pi;
            rho=a*cos(k*theta);
        else
            theta=delt:delt:2*k2*pi;
            rho=a*cos(k*theta);
        end
    else
        %% TODO
            theta=delt:delt:2*k2*pi;
            rho=a*cos(k*theta);
    end
end
%[theta, r]=cart2pol(x,y) %ѿתΪ

[x,y]=pol2cart(theta, rho); %תΪѿ
X = [x' y'];