t=1;alpha=pi/3;%角度0,2*pith=0:pi/1000:pi*2;%半径r=50;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold onth=0:pi/1000:pi*2;%半径r=61;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold on% 画15度50曲线13th=0:pi/1000:pi/12;%半径r=50;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold onth=23*pi/12:pi/1000:2*pi;%半径r=50;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold on% 画30度50曲线2th=pi*11/12:pi/1000:13*pi/12;%半径r=50;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold on% 画右上过度曲线th=pi/12:pi/1000:5*pi/12;%半径r=50;R=61;r01=r+2*(R-r)*(th-pi/12).2/alpha2;r02=2*r-R+4*(R-r)*(th-pi/12)-(th-pi/12).2/(2*alpha)/alpha;xt=r01.*cos(th);yt=r01.*sin(th);xt1=r02.*cos(th);yt1=r02.*sin(th);plot(xt,yt,-)hold onplot(xt1,yt1,-)hold on% 画61,30度线th=5*pi/12:pi/1000:7*pi/12;%半径r=61;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold onth=17*pi/12:pi/1000:19*pi/12;%半径r=61;xt=r*cos(th);yt=r*sin(th);plot(xt,yt,-)hold on% 画左上过度曲线th=pi/12:pi/1000:5*pi/12;%半径r=50;R=61;r01=r+2*(R-r)*(th-pi/12).2/alpha2;r02=2*r-R+4*(R-r)*(th-pi/12)-(th-pi/12).2/(2*alpha)/alpha;xt=r01.*cos(pi-th);yt=r01.*sin(pi-th);xt1=r02.*cos(pi-th);yt1=r02.*sin(pi-th);plot(xt,yt,-)hold onplot(xt1,yt1,-)hold on% 画左下过度曲线th=pi/12:pi/1000:5*pi/12;%半径r=50;R=61;r01=r+2*(R-r)*(th-pi/12).2/alpha2;r02=2*r-R+4*(R-r)*(th-pi/12)-(th-pi/12).2/(2*alpha)/alpha;xt=r01.*cos(pi+th);yt=r01.*sin(pi+th);xt1=r02.*cos(pi+th);yt1=r02.*sin(pi+th);plot(xt,yt,-)hold onplot(xt1,yt1,-)hold on% 画右下过度曲线th=pi/12:pi/1000:5*pi/12;%半径r=50;R=61;r01=r+2*(R-r)*(th-pi/12).2/alpha2;r02=2*r-R+4*(R-r)*(th-pi/12)-(th-pi/12).2/(2*alpha)/alpha;xt=r01.*cos(-th);yt=r01.*sin(-th);xt1=r02.*cos(-th);yt1=r02.*sin(-th);plot(xt,yt,-)hold onplot(xt1,yt1,-)hold onPublished with MATLAB? R2014b