Re: [問題] 點光源模擬已回收
close
clear
wavelength=633e-6;
k=2*pi/wavelength;
L=200;N=128;
z=L/sqrt(N*wavelength);
ymax=L/2;y=linspace(-ymax, ymax, N);
xmax=ymax;x=linspace(-xmax, xmax, N);
for ii=1:N
for jj=1:N
r(ii,jj)=sqrt(x(ii)^2+y(jj)^2+z^2);
O(ii,jj)=1/r(ii,jj)*exp(i*k*r(ii,jj));
end
end
lgray=zeros(256,3);
for ii=0:255
lgray(ii+1,:)=(255-ii)/255;
end
imagesc(abs(O)); colorbar;
colormap(lgray);
axis equal tight
% 試試看唄
% 應該可以跑吧
※ 引述《lexjiang (Lex)》之銘言：
: wavelength=633e-6;
: k=2*pi/wavelength;
: L=200;N=128;
: z=L/sqrt(N*wavelength);
: ymax=L/2;y=linspace(-ymax, ymax, N);
: xmax=ymax;x=linspace(-xmax, xmax, N);
: for i=1:N
: for j=1:N
: r(i,j)=sqrt(x(i)^2+y(j)^2+z^2);
: O(i,j)=1/r(i,j)*exp(i*k*r(i,j));
: end
: end
: lgray=zeros(256,3);
: for i=0:255
: lgray(i+1,:)=(255-i)/255;
: end
: imagesc(O);
: colormap(lgray);
: 只有黑底而已 ..該怎麼修改呢
--
拙僧の肉棒を試させろ?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.116.204.213
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
