作者huangkurt (kurt)
标题Re: [问题] 回圈取值
时间Mon Oct 6 15:08:42 2014
这是用原本的程式改写的
不是很懂原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);
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.113.106.128
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1412579324.A.C37.html
※ 编辑: huangkurt (140.113.106.128), 10/06/2014 15:31:48