先附上程式码
clear
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);
r1=0.1;
r2=0.2;
cx=[-1 0 1];
cy=[-1 0 1];
for nx=1:n+1
for ny=1:n+1
ra(nx,ny)=(x(nx)-cx(2))^2+(y(ny)-cy(2))^2;
rb(nx,ny)=(x(nx)-cx(1))^2+(y(ny)-cy(1))^2;
if (r1>ra(nx,ny))
eta(nx,ny)=0.2;
elseif(r1<ra(nx,ny) & r2>rb(nx,ny))
eta(nx,ny)=0.5;
else
eta(nx,ny)=1;
end
end
end
contourf(x,y,eta);
这样出现的图会是以(x,y)=(0,0)以及(-1-1)为圆心画两个圆
我想问的是说我有没有办法用个循环
让他分别在
(x,y)=(-1,-1) (-1,0) (-1,1) (0,-1) (0,0) (0,1) (1,-1) (1,0) (1,1)
以不同的值画圆呢
有试着把在cx和cy写循环
但是出来的图形都只有在(x,y)=(1,1)的地方有圆