Matlab绘制双波干涉

最新推荐文章于 2022-04-11 19:07:34 发布

Reckonist

最新推荐文章于 2022-04-11 19:07:34 发布

阅读量3.8k

点赞数 11

分类专栏： matlab 文章标签： matlab 物理学

转载或引用须注明出处，并标明页面地址；访问者可将页面提供的内容用于个人学习、研究或欣赏，以及其他非商业性或非盈利性用途。

本文链接：https://blog.csdn.net/qq_37077097/article/details/116918837

版权

matlab 专栏收录该内容

4 篇文章 1 订阅

订阅专栏

直接上图

在这里插入图片描述

源波

在这里插入图片描述

双波叠加

在这里插入图片描述

干涉条纹

代码

%波的干涉 
clc
clear all
f=10;%波源频率(两波同频)
A=2;%幅值
x=-5:0.01:5;
y=-5:0.01:5;
L_x=length(x);
z=zeros(1,L_x);

[X1,Y1,Z1]=griddata(x,y,z,linspace(min(x),max(x),150)',...
    linspace(min(y),max(y),150),'cubic');%插值
[X2,Y2,Z2]=griddata(x,y,z,linspace(min(x),max(x),150)',...
    linspace(min(y),max(y),150),'cubic');%插值
[X3,Y3,Z3]=griddata(x,y,z,linspace(min(x),max(x),150)',...
    linspace(min(y),max(y),150),'cubic');%插值
N=length(X1);

for i=1:N
    for j=1:N
        r1=sqrt((X1(i,j)+2)^2+Y1(i,j)^2);
        r2=sqrt((X2(i,j)-2)^2+Y2(i,j)^2);
        Z1(i,j)=A*sin(r1*f);
        Z2(i,j)=A*sin(r2*f);
        Z3(i,j)=Z1(i,j)+Z2(i,j);
    end
end

figure
surf(X1,Y1,Z1);%三维曲面1
hold on
surf(X2,Y2,Z2);%三维曲面
colorbar('vert')
view(-8,87)

figure
surf(X3,Y3,Z3);%三维曲面
colorbar('vert')
view(-8,87)

figure
contourf(X3,Y3,Z3),shading interp;
colorbar('vert')