实验三源代码
内容一:
x=linspace(-0.003,0.003,1000);
y=linspace(-0.003,0.003,1000);
R=0.855;
la=589.3*10.^(-9);
[X,Y]=meshgrid(x,y);
r2=X.^2+Y.^2;
I=abs(sin(pi*r2/R/la)).^2;
imshow(I);
内容二:
x=linspace(-0.003,0.003,1000);
y=linspace(-0.003,0.003,1000);
R=0.855;
la=589.3*10.^(-9);
[X,Y]=meshgrid(x,y);
r2=X.^2+Y.^2;
F=moviein(50);
for delta=1:51;
I=abs(sin(pi*r2/R/la+2*pi*(del
ta)*0.04*10.^(-6)/la)).^2;
imshow(I);
F(delta)=getframe;
end
movie(F,10)
movie2avi(F,'dynamic');
思考题2:
x=linspace(-0.003,0.003,1000);
y=linspace(-0.003,0.003,1000);
R=0.855;
la=589.3*10.^(-9);
[X,Y]=meshgrid(x,y);
r2=X.^2+Y.^2;
F=moviein(50);
for delta=1:51;
I=abs(sin(pi*r2/R/la+2*pi*(51-
delta)*0.04*10.^(-6)/la)).^2;
imshow(I);
F(delta)=getframe;
end
movie(F,10)
movie2avi(F,'dynamic');
评论0