pcf_circle.m

来自「Photonic crystal fiber section area Grap」· M 代码 · 共 40 行

function p
l1=.160
l2=.265
T=l1+l2

for r=19:-1:1
    if mod(r,2)==0
        circle(0,0,r-.3)
    else
        circle1(0,0,r)
    end
end
circle2(0,0,3)
end


function circle(x,y,r)
fai=0:pi/100:2*pi
x0=r*sin(fai)+x
y0=r*cos(fai)+y
plot(x0,y0)

end

function circle1(x,y,r)
fai=0:pi/100:2*pi
x0=r*sin(fai)+x
y0=r*cos(fai)+y
plot(x0,y0)
fill(x0,y0,'g')
hold on
end

function circle2(x,y,r)
fai=0:pi/100:2*pi
x0=r*sin(fai)+x
y0=r*cos(fai)+y
plot(x0,y0)

end