Re: [問題] 迴圈取值

看板MATLAB作者 (kurt)時間9年前 (2014/10/06 15:08), 9年前編輯推噓0(000)
留言0則, 0人參與, 最新討論串3/3 (看更多)
這是用原本的程式改寫的 不是很懂原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://www.ptt.cc/bbs/MATLAB/M.1412579324.A.C37.html ※ 編輯: huangkurt (140.113.106.128), 10/06/2014 15:31:48
文章代碼(AID): #1KCZ_ymt (MATLAB)
文章代碼(AID): #1KCZ_ymt (MATLAB)