这是用原本的程式改写的
不是很懂原PO要不同的判断值(0.2 or 0.5)
但如果要加的话只需要在程式内用一些判断式应该就没问题了
clear all;close all;clc;
n=11;
g1=[3 0];
g2=[0 3];
x=linspace(-g1(1)/2,g1(1)/2,n+1);
y=linspace(-g2(2)/2,g2(2)/2,n+1);
xy=fliplr(x)'*ones(1,n+1)+1i*ones(n+1,1)*x;
r1=0.1;
r2=0.2;
c=reshape(ones(3,1)*(-1:1)+1i*(-1:1)'*ones(1,3),1,[]);
diff=ones([size(xy),length(c)]);
for nc=1:length(c)
diff(:,:,nc)=abs(xy-c(nc));
end
eta=zeros([size(xy),length(c)]);
eta(diff<sqrt(r1))=0.2;
eta=sum(eta,3);
contourf(x,y,eta);