%水波干涉的動(dòng)畫演示

a=4;w=4;A0=0.3;v=1;%水波源的位置及方程各參量

[x,y]=meshgrid([-4*pi:pi/20:4*pi],[-4*pi:pi/20:4*pi]);%產(chǎn)生網(wǎng)絡(luò)坐標(biāo)

r1=sqrt((x-a).^2+y.^2)+5;%空間點(diǎn)（x，y）到波源1（a，0）的距離

c1=1./r1;%波源1傳播到各點(diǎn)的振幅

r2=sqrt((x+a).^2+y.^2)+5;%空間點(diǎn)（x，y）到波源2（-a，0）的距離

c2=1./r2;%波源2傳播到各點(diǎn)的振幅

t=0:0.1:10;%設(shè)定運(yùn)行時(shí)間

for j=1:41%此循環(huán)產(chǎn)生動(dòng)畫

??? z=A0*c1.*cos(w*(t(j)-r1/v))+0.3*c2.*cos(w*(t(j)-r2/v));%波動(dòng)疊加

surf(x,y,z)%繪制波動(dòng)疊加曲面

colormap(gray);%創(chuàng)建線性灰階色圖給圖形表面著色

axis equal;%使每個(gè)坐標(biāo)軸的刻度增量相同

shading interp;%用細(xì)致添加描影法渲染圖形

view([180 90])%在方位角180°，仰角90°觀察圖形

axis([-15 15 -15 15 -0.1 0.1])%設(shè)定x、y、z、軸的范圍

axis off%去掉坐標(biāo)軸

m(j)=getframe;%保存繪制好的圖形到變量數(shù)組m（j）

end

movie(m,10)%放映儲(chǔ)存在m矩陣中電影動(dòng)畫10遍

%水波干涉的動(dòng)畫演示

a=4;w=4;A0=0.3;v=1;%水波源的位置及方程各參量

[x,y]=meshgrid([-4*pi:pi/20:4*pi],[-4*pi:pi/20:4*pi]);%產(chǎn)生網(wǎng)絡(luò)坐標(biāo)

r1=sqrt((x-a).^2+(y-a).^2)+5;%空間點(diǎn)（x，y）到波源1（a，a）的距離

c1=1./r1;%波源1傳播到各點(diǎn)的振幅

r2=sqrt((x-a).^2+(y+a).^2)+5;%空間點(diǎn)（x，y）到波源2（a，-a）的距離

c2=1./r2;%波源2傳播到各點(diǎn)的振幅

r3=sqrt((x+a).^2+(y+a).^2)+5;%空間點(diǎn)（x，y）到波源3（-a，-a）的距離

c3=1./r3;%波源3傳播到各點(diǎn)的振幅

r4=sqrt((x+a).^2+(y-a).^2)+5;%空間點(diǎn)（x，y）到波源3（-a，a）的距離

c4=1./r4;%波源4傳播到各點(diǎn)的振幅

t=0:0.1:10;%設(shè)定運(yùn)行時(shí)間

for j=1:41%此循環(huán)產(chǎn)生動(dòng)畫

??? z=A0*c1.*cos(w*(t(j)-r1/v))+0.3*c2.*cos(w*(t(j)-r2/v))+0.3*c3.*cos(w*(t(j)-r3/v))+0.3*c4.*cos(w*(t(j)-r4/v));%波動(dòng)疊加

surf(x,y,z)%繪制波動(dòng)疊加曲面

colormap(gray);%創(chuàng)建線性灰階色圖給圖形表面著色

axis equal;%使每個(gè)坐標(biāo)軸的刻度增量相同

shading interp;%用細(xì)致添加描影法渲染圖形

view([180 90])%在方位角180°，仰角90°觀察圖形

axis([-15 15 -15 15 -15 15])%設(shè)定x、y、z、軸的范圍

axis off%去掉坐標(biāo)軸

m(j)=getframe;%保存繪制好的圖形到變量數(shù)組m（j）

end

movie(m,10)%放映儲(chǔ)存在m矩陣中電影動(dòng)畫10遍