Untitled1.zip

r0=4*10^-3;%内环半径 a=1/8; %锯齿比值宽和基圆半径的比值 M=72; %锯齿数 y=[]; i=1; o=[]; r=[]; for m=0:M-1 for theta=(2*m-1)/M*pi:pi/M/10:(2*m+1)/M*pi o(i)=theta; y(i)=rect((theta-(2*m*pi)/M)*M/pi); r(i)=r0.*(1+a.*y(i));%方齿齿形分 i=i+1; end end % h=polar([0 2*pi],[0 5*10^-3]); % hold on % cla [x,y]=pol2cart(o,r); % polar(o,r); plot(x,y) % hold on % % % end axis equal axis off