先附上程式碼
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)的地方有圓