分形二元器件

N=2;S=3;
n=(2*N-1)^S;
lamda=0.532; 
f=2e5;
k=2*pi/lamda;
nn=1000;
R1=sqrt(lamda*f*n);
% L=sqrt(lamda*f*14*2);
L=R1;
x=linspace(-L,L,nn);
y=linspace(-L,L,nn);  
[xx,yy]=meshgrid(x,y);

a=R1; 
s=((xx.^2+yy.^2))./a^2-0.5;
t=ones(nn);
% p=2/((2*N-1)^m);
 for m=0:S
 p=2/((2*N-1)^m);
t=t.*rectpuls(mod(s+(p-1)/2,p)/p,1).*rectpuls(s,1);%分形公式
 end
figure

imshow(t)