matlab物理案例分析
Matlab在物理上的应用举例
1. 单列波%%单列波t=0:0.001:10;A=input('振幅A=');w=input('频率w=');a=input('a=');y=A.*sin(w.*t+a);plot(t,y);pause(1),sound(y);ylabel('y'),xlabel('t')2. %%光的单缝衍射现象Lambda=500e-9; %a=input('a='); % 可取0.2e-3,1e-3,2e-3三种情况z=1 %ymax=3*Lambda*z/a; %Ny=51; %ys=linspace(-ymax,ymax,Ny); %NPoints=51; %yPoint=linspace(-a/2,a/2,NPoints); %for j=1:Ny %L=sqrt((ys(j)-yPoint).^2+z^2); %Phi=2*pi.*(L-z)./Lambda; %SumCos=sum(cos(Phi)); %SumSin=sum(sin(Phi)); %B(j)=(SumCos^2+SumSin^2)/NPoints^2; %endclf,plot(ys,B,'*',ys,B);grid; %3. %%用毕奥－沙伐尔定律计算电流环产生的磁场mu0=4*pi*1e-7;I0=5.0;Rh=1;C0=mu0/(4*pi)*I0;NGx=21;NGy=21;x=linspace(-Rh,Rh,NGx);y=linspace(-3,3,20);y=x;Nh=20;theta0=linspace(0,2*pi,Nh+1);theta1=theta0(1:Nh);y1=Rh*cos(theta1);z1=Rh*sin(theta1);theta2=theta0(2:Nh+1);y2=Rh*cos(theta2);z2=Rh*sin(theta2);xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;for i=1:NGyfor j=1:NGxrx=x(j)-xc;ry=y(i)-yc;rz=0-zc;r3=sqrt(rx.^2+ry.^2+rz.^2).^3;dlXr_x=dly.*rz-dlz.*ry;dlXr_y=dly.*rx-dlx.*rz;Bx(i,j)=sum(C0*dlXr_x./r3);By(i,j)=sum(C0*dlXr_y./r3); endendclf;quiver(x,y,Bx,By);4. %%多普勒效应x0=500;v=50;y0=20;c=330;w=1000;t=0:0.001:30;r=sqrt((x0-v*t).^2+y0.^2);t1=t-r/c;u=sin(w*t)+sin(1.1*w*t);u1=sin(w*t1)+sin(1.1*w*t1); sound(u);pause(5);sound(u1);5.亥姆霍兹线圈clear allmu0=4*pi*1e-7;I0=5.0;Rh=1;C0=mu0/(4*pi)*I0;NGx=21;NGy=21;x=linspace(-Rh,Rh,NGx);y=linspace(-Rh,Rh,NGy);Nh=20;theta0=linspace(0,2*pi,Nh+1); theta1=theta0(1:Nh);y1=Rh*cos(theta1);z1=Rh*sin(theta1);theta2=theta0(2:Nh+1);y2=Rh*cos(theta2);z2=Rh*sin(theta2);xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;for i=1:NGyfor j=1:NGxrx=x(j)-xc;ry=y(i)-yc;rz=0-zc;r3=sqrt(rx.^2+ry.^2+rz.^2).^3;dlXr_x=dly.*rz-dlz.*ry;dlXr_y=dly.*rx-dlx.*rz;Bx(i,j)=sum(C0*dlXr_x./r3);By(i,j)=sum(C0*dlXr_y./r3);endendBax=Bx(:,11:21)+Bx(:,1:11);Bay=By(:,11:21)+By(:,1:11);subplot(1,2,1)mesh(x(11:21),y,Bax);xlabel('x');ylabel('y'); subplot(1,2,2),plot(y,Bax),grid,xlabel('y');ylabel('Bx'); ６．库仑引力clear all;N=input('电荷数目N：');for ic=1:N %fprintf('-----\n 对电荷#%g\n',ic);rc=input('电荷位置[x y]（米）：');x(ic)=rc(1); %y(ic)=rc(2); %q(ic)=input('输入电荷量（库仑）'); endE0=8.85e-12; %C0=1/(4*pi*E0); %for ic=1:NFx=0.0; Fy=0:0; %for jc=1:N %if(ic~=jc) %xij=x(ic)-x(jc); yij=y(ic)-y(jc);Rij=sqrt(xij^2+yij^2); %Fx=Fx+C0*q(ic)*xij/Rij^3;Fy=Fy+C0*q(ic)*yij/Rij^3;endendfprintf('其它电荷作用在电荷#%g上的合力为：\n',ic); %fprintf('x－分量：%gN\n',Fx);fprintf('y－分量：%gN\n',Fy)end７．李萨如图形% lisaru.msyms t a1 a2 w1 w2x=cos(w1.*t+a1);y=sin(w2.*t+a2);a1=input('a1=');a2=input('a2=');w1=input('w1=');w2=input('w2=');tf=10;Ns=1000;t=linspace(0,tf,Ns);dt=tf/(Ns-1); %分Ns个点,求出时间增量dtxplot=eval(x);yplot=eval(y); %计算Ns个点的位置x(t),y(t)figure(gcf);subplot(1,2,1),for i=1:750plot(yplot(1:i),xplot(1:i)); %画点的轨迹图axis('equal'); grid ; %使两轴比例相同pause(0.01)end８．耦合振子m1=2;m2=2;K1=16;K2=4; % x0=[1;0];xd0=[6;6];tf=10; %M=[m1,0;0,m2];K=[K1+K2,-K2;-K2,K1+K2] ; %u(:,s)�[u,L]=eig(K,M) ; %t=linspace(0,tf,101);x=zeros(2,101); % for s=1:2 %alfa=sqrt(u(:,s)'*M*u(:,s)) ; %u(:,s)= u(:,s)/alfa; %w(j)=sqrt(L(j,j)); % xt=u(:,j)*(u(:,j)'*M*x0*cos(w(j)*t)+u(:,j)'*M*xd0/w(j)*sin(w(j)*t));x=x+xt; %endfor r=1:2 %x1,x1?? subplot(2,1,r)plot(t.x(r,:)),grid;xlabel('xxx');ylabel(['yyy',num2str(r)]);end９．拍频%%t=0:0.001:10;a1=input('??1=');w1=input('??1=');a2=input('??2=');w2=input('??2=');y1=a1*sin(w1*t);y2=a2*sin(w2*t);y=y1+y2;subplot(3,1,1),plot(t,y1),ylabel('y1')subplot(3,1,2),plot(t,y2),ylabel('y2')subplot(3,1,3),plot(t,y),ylabel('y'),xlabel('t')pause,sound(y1);pause(5),sound(y2);pause(5),sound(y),pausesubplot(1,1,1)１０。
Matlab语言在物理实验数据处理中的应用
计算周期的平方值
户一 ｐ ｏ ｌ ｙ ｆ ｉ ｔ （ Ｌ， Ｔ丁， １ ） ；
单摆周期公式 ， 测得的 ｇ 值 完全符合所在地重力 加 速度 的大 小 。
用 Ｍａ ｔ ｌ ａ ｂ语 言不 仅 可以拟合 直线 ， 还 可 以拟 合 曲线 。例 如 ， 在利 用 霍 尔 效应 测 量 螺 线 管 内轴 线 上磁 感应 强度 的分 布实 验 中 ， 磁 场 随 位置 的分 布图线 是一 条 曲线 ， 要 在 坐 标 纸上 手 工 画 出这 条
参 考文 献 ：
［ １ ］ 张志勇． 掌 握 和精通 ＭＡＴ Ｌ ＡＢ Ｅ Ｍ］ ． 北京： 北 京航
空航天 大学 出版社 ， １ ９ ９ ７ ．
Ｅ ２ ］ 钟 季康 ， 鲍 鸿 吉． 大 学 物 理 习题 计 算 机 解 法——
摆长 Ｌ ／ ｃ ｍ
图 １ 单摆 测 重 力加 速 度 实 验 与 拟 合 曲线
应为一条直线 ， 其斜率为 ｋ 一 ， 因此， 重力加速
收 稿 日期 ： ２ ０ １ １ — ０ ７ １ ０
２ ． Ｏ １ １ ］ ； 输 入周 期 的实验数 据
万方数据
４ ４
Ｍａ ｔ ｌ ａ ｂ语言在物理实验数据处理 中的应用
Ｉ Ｔ 一 ＊Ｔ ：
方 成正 比 ， 用 线性 回归 函数 比较合 理 ， 从 而验 证 了
对 Ｘ坐标轴 加标 注
ｙ ｌ ａ ｂ ｅ ｌ （ ‘ Ｔ ＇ ２ ／ ｓ ２ ’ ） ：
曲线拟合效果越好 ， Ｍａ ｔ ｌ ａ ｂ 语言提供 了各种多项 式 拟合 方案 ， 可根 据 具 体情 况 选 取 拟合 效 果 最 好
的那 条 曲线 。
０ ５ ４ ３ ３ ０ ２ ５ ２ Ｏ ｌ ５ １ ０ ０ ５ ．Ｏ
适合用matlab解决的经典物理例题
适合用Matlab解决的经典物理例题在物理学领域，经典物理例题一直是学习和研究的重要内容。
而Matlab作为一种强大的数学软件，非常适合解决各种物理问题。
本文将从力学、电磁学和热力学等多个方面，选取一些经典的物理例题，通过Matlab进行分析和求解，展示Matlab在解决物理问题时的强大用途。
1. 简谐振动问题简谐振动是物理学中一个重要的模型，涉及到弹簧振子、单摆等问题。
通过Matlab可以很方便地求解简谐振动的运动规律。
对于弹簧振子的运动方程，可以通过Matlab进行数值模拟，得到振动的周期、频率、位移等参数，从而更好地理解简谐振动的特性。
2. 电场问题在电磁学中，电场是一个重要的研究对象。
通过Matlab可以很容易地分析不同形状的电荷分布所产生的电场分布。
可以通过Matlab计算出点电荷、均匀带电细棒等情况下的电场分布，并绘制出电场线图，直观地展现电场的分布规律。
这样的分析对于理解电场的性质和相互作用具有重要意义。
3. 热传导问题热传导是热力学研究的一个重要方面，涉及到导热方程的求解和热量分布的分析。
通过Matlab可以对不同材料和形状的热传导问题进行数值模拟和求解。
可以通过Matlab计算出棒状材料中的温度分布随时间的演化，从而得到材料的热传导性能。
这样的分析对于工程实践中的热设计和材料选型具有重要指导意义。
4. 万有引力问题在力学中，万有引力是一个经典的例题，涉及到行星轨道、卫星运动等问题。
通过Matlab可以很方便地进行万有引力场下的物体运动模拟。
可以通过Matlab计算地球和月球的引力作用下的月球轨道，从而揭示天体运动的规律和特性。
这样的模拟对于探索宇宙中天体运动规律具有重要帮助。
总结回顾：通过以上例题的分析，我们不仅了解了Matlab在经典物理例题中的应用，也可以发现Matlab在解决物理问题时的便捷和高效。
当然，实际物理问题可能具有更多的复杂性和多样性，需要结合理论分析和实验数据进行综合研究。
Matlab在基础物理学中的应用举例4
Matlab在基础物理学中的应用举例(四)（光学：牛顿环的模拟）2011-02-13 13:30:18| 分类：Matlab Sky | 标签：matlab 光学|字号大中小订阅下面是模拟牛顿环的图案，给出了一个用户交换操作界面的实现程序：close all;clear;figure('Position',[90 164 873 483]); % 设置图形窗口位置L=632.8;R=5;H=0.005; % 设置默认参数值：波长L、半径R和厚度Ha1=axes('Position',[0.83,0.3,0.15,0.4]); % 设置坐标轴位置，在其上绘制简易的原理图hold on;axis([0,1,0,1]); % 设置坐标轴范围plot([0.25,0.25],[0.5,0.8],'r'); % 画竖直线，入射光线plot([0.5,0.5],[0.5,0.8],'r'); % 画竖直线，入射光线plot([0.75,0.75],[0.5,0.8],'r'); % 画竖直线，入射光线fill([0.22,0.25,0.27],[0.5,0.44,0.5],'r'); % 画出箭头fill([0.47,0.5,0.52],[0.5,0.44,0.5],'r'); % 画出箭头fill([0.72,0.75,0.77],[0.5,0.44,0.5],'r'); % 画出箭头z=1.8*exp(i*(linspace(-0.21,0.21,30)-pi/2))+2.1*i+0.5; % 生成平凸透镜的轮廓fill(real(z),imag(z),'w'); % 利用fill函数画平凸透镜的轮廓rectangle('Position',[0.1,0.18,0.8,0.12]); % 画玻璃板底面set(gca,'xtick',[],'ytick',[],'box','on'); % 删去坐标轴刻度a2=axes('Position',[0.4,0.16,0.4,0.7]); % 生成第二个坐标轴，在其上画牛顿环[x,y]=meshgrid(linspace(-0.005,0.005,200)); % 生成坐标矩阵r2=(x.^2+y.^2); % 坐标矩阵各点到中心的距离Di=[2*H+2*(R-sqrt(R^2-r2))*1e9]/L; % 计算光程差In=abs(cos(Di*pi*2)); % 得到牛顿环上的光强数值cr=abs(L-560)/200; % 红色的分量值cg=1-cr; % 绿色的分量值cb=abs(L-600)/240; % 蓝色的分量值Ik(:,:,1)=In*cr; % 转化为RGB三基色表示颜色矩阵数值（红色）Ik(:,:,2)=In*cg; % 转化为RGB三基色表示颜色矩阵数值（绿色）Ik(:,:,3)=In*cb; % 转化为RGB三基色表示颜色矩阵数值（蓝色）Pc=imshow(Ik,[]); % 绘制缺省参数下的数值title('the pattern of Newton''s rings','fontsize',12); % 添加图题% 实时显示波长的文本框Lt=uicontrol(gcf,'style','text',...'unit','normalized','position',[0.06,0.86,0.21,0.06],...'BackgroundColor',0.7*[1,1,1],'ForegroundColor',[0.8,0.1,0.9],...'string','wavelength: 632.8nm','fontsize',12,'fontname','times new roman');s1=uicontrol(gcf,'style','slider',...'unit','normalized','position',[0.06,0.76,0.21,0.04],...'BackgroundColor',0.7*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'SliderStep',[0.01,0.01],'value',(632.8-360)/400,...'callback',['L=get(s1,''value'')*400+360;',...'set(Lt,''string'',[''wavelength: '',num2str(L/10),''nm'']);',...'Di=[2*H+2*(R-sqrt(R^2-r2))*1e9]/L;',...'In=abs(cos(Di*pi*2));cr=abs(L-560)/200;cg=1-cr;',...'cb=abs(L-600)/240;Ik(:,:,1)=In*cr;Ik(:,:,2)=In*cg;',...'Ik(:,:,3)=In*cb;set(Pc,''CData'',Ik);']); % 利用滑动条改变波长数值uicontrol(gcf,'style','text',...'unit','normalized','position',[0.04,0.81,0.08,0.04],...'BackgroundColor',0.8*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'string','360','fontsize',12,'fontname','times new roman'); % 显示最小波长数值uicontrol(gcf,'style','text',...'unit','normalized','position',[0.22,0.81,0.08,0.04],...'BackgroundColor',0.8*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'string','760','fontsize',12,'fontname','times new roman'); % 显示最大波长数值Rt=uicontrol(gcf,'style','text',...'unit','normalized','position',[0.06,0.66,0.23,0.06],...'BackgroundColor',0.7*[1,1,1],'ForegroundColor',[0.8,0.1,0.9],...'string','radii:','fontsize',12,'fontname','times new roman'); %实时显示半径数值的文本框s2=uicontrol(gcf,'style','slider',...'unit','normalized','position',[0.06,0.56,0.21,0.04],...'BackgroundColor',0.7*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'SliderStep',[0.01,0.01],...'callback',['R=get(s2,''value'')*7+5;',...'set(Rt,''string'',[''radii: 5m'',num2str(R),''m'']);',...'Di=[2*H+2*(R-sqrt(R^2-r2))*1e9]/L;',...'In=abs(cos(Di*pi*2));cr=abs(L-560)/200;cg=1-cr;',...'cb=abs(L-600)/240;Ik(:,:,1)=In*cr;Ik(:,:,2)=In*cg;',...'Ik(:,:,3)=In*cb;set(Pc,''CData'',Ik);']); % 通过滑动条改变半径数值uicontrol(gcf,'style','text',...'unit','normalized','position',[0.04,0.61,0.08,0.04],...'BackgroundColor',0.8*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'string','5','fontsize',12,'fontname','times new roman'); % 显示最小半径uicontrol(gcf,'style','text',...'unit','normalized','position',[0.22,0.61,0.08,0.04],...'BackgroundColor',0.8*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'string','12','fontsize',12,'fontname','times new roman'); % 显示最大半径% 实时显示厚度的文本框Ht=uicontrol(gcf,'style','text',...'unit','normalized','position',[0.06,0.46,0.23,0.06],...'BackgroundColor',0.7*[1,1,1],'ForegroundColor',[0.8,0.1,0.9],...'string','thickness: 5nm','fontsize',12,'fontname','times new roman');s3=uicontrol(gcf,'style','slider',...'unit','normalized','position',[0.06,0.36,0.21,0.04],...'BackgroundColor',0.7*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'SliderStep',[0.01,0.01],'value',0.05,...'callback',['H=get(s3,''value'')*0.01;',...'set(Ht,''string'',[''thickness: '',num2str(H),''nm'']);',...'Di=[2*H+2*(R-sqrt(R^2-r2))*1e9]/L;',...'In=abs(cos(Di*pi*2));cr=abs(L-560)/200;cg=1-cr;',...'cb=abs(L-600)/240;Ik(:,:,1)=In*cr;Ik(:,:,2)=In*cg;',...'Ik(:,:,3)=In*cb;set(Pc,''CData'',Ik);']); % 通过滑动条改变厚度数值uicontrol(gcf,'style','text',...'unit','normalized','position',[0.04,0.41,0.08,0.04],...'BackgroundColor',0.8*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'string','0','fontsize',12,'fontname','times new roman'); % 显示最小厚度uicontrol(gcf,'style','text',...'unit','normalized','position',[0.22,0.41,0.08,0.04],...'BackgroundColor',0.8*[1,1,1],'ForegroundColor',[0.1,0.1,0.9],...'string','0.01','fontsize',12,'fontname','times new roman'); % 显示最大厚度运行结果：。
MATLAB在力学机械中的应用举例(课件PPT)
wx=diff(vx)./dt(1:Ldt-1);wy=diff(vy)./dt(1:Ldt-1); %二次导数
[t(2:Ldt),x(2:Ldt),y(2:Ldt),wx,wy]
%显示数据
10
子程序ex713f
函数程序应另存成一个文件ex713f.m function zprime=ex713f(t,z) global vt vm zprime=[0;0]; % 给出t0之前zprime初值 zprime(1)=-vt-vm*z(1)/sqrt(z(1)^2+z(2)^2); zprime(2)=-vm*z(2)/sqrt(z(1)^2+z(2)^2); %上面两句可换成一个矩阵语句： zprime=-vt*[1;0]-vm*z/sqrt(z(1)^2+z(2)^2);
5
线性数学模型
对杆件1：ΣX=0 Nax + Ncx = 0 ΣY=0 Nay + Ncy - G1 = 0; ΣM=0 Ncy*L1*cos(theta1)-Ncx*L1*sin(theta1)-…
G1*L1/2*cos(theta1)=0; 对杆件2: ΣX=0 Nbx - Ncx = 0; ΣY=0 Nby - Ncy - G2 = 0; ΣM=0 Ncy*L2*cos(theta2)+ …
从而 w3 = L1w1cos(π/2-θ1+θ2)/ (L3cos(θ3-π/2-θ2)) 由杆2两端点a和b的速度沿杆长垂直方向的分量之差,可以求
出杆2的角速度. w2 = (-(L3sin(θ3-π/2-θ2))- L1w1sin(π/2-θ1+θ2))/L2 2. 求运动全过程的角位置,角速度,角加速度曲线,这只有借助 于计算工具才能做到,因为用手工算一个点就不胜其烦, 算 几十个点是很难想象的.而由MATLAB编程调用fzero函数时, 要求给出一个近似猜测值,若连续算几十点,前一个解就可 作为后一个解的猜测值,所以反而带来了方便. 这样,本书将提供两个程序ex714a.m和ex714b.m来表述这两种 方法,它们所要调用的函数程序命名为ex714f.m.
MATLAB应用实例分析例分析
MATLAB应用实例分析例分析Matlab应用例题选讲仅举一些运用MATLAB的例子，这些问题在数学建模中时常遇到，希望能帮助同学们在短时间内方便、快捷的使用MATLAB 解决数学建模中的问题，并善用这一工具。
常用控制命令:clc:%清屏; clear:%清变量; save:%保存变量; load:%导入变量一、利用公式直接进行赋值计算本金P以每年n次，每次i%的增值率(n与i的乘积为每年增值额的百分比)增加，当增加到r×P 时所花费的时间T为:(利用复利计息公式可得到下式) lnrnT() r，P,P(1，0.01i),T,r,2,i,0.5,n,12nln(1，0.01i)MATLAB 的表达形式及结果如下:>> r=2;i=0.5;n=12; %变量赋值>> T=log(r)/(n*log(1+0.01*i)) 计算结果显示为:T = 11.5813即所花费的时间为T=11.5813 年。
分析:上面的问题是一个利用公式直接进行赋值计算问题，实际中若变量在某个范围变化取很多值时，使用MATLAB，将倍感方便，轻松得到结果，其绘图功能还能将结果轻松的显示出来，变量之间的变化规律将一目了然。
若r在[1,9]变化，i在[0.5,3.5]变化;我们将MATLAB的表达式作如下改动，结果如图1。
r=1:0.5:9;i=0.5:0.5:3.5;n=12;p=1./(n*log(1+0.01*i));T=log(r')*p;plot(r,T)xlabel('r') %给x轴加标题ylabel('T') %给y轴加标题q=ones(1,length(i));text(7*q-0.2,[T(14,1:5)+0.5,T(14,6)-0.1,T(14,7)-0.9],num2str(i'))40350.5302520T 1151.510 22.55 33.50123456789r图11从图1中既可以看到T随r的变化规律，而且还能看到i的不同取值对T—r 曲线的影响(图中的六条曲线分别代表i的不同取值)。
matlab 物理光学实验matlab 物理光学实验
1.工程光学系列之一：杨氏双缝干涉matlab1.基本原理杨氏干涉实验是两点光源干涉实验的典型代表。
杨氏干涉实验以极简单的装置和巧妙构思实现了普通光源干涉。
无论从经典光学还是从现代光学的角度来看，杨氏实验都具有十分重要的意义。
杨氏双缝实验的装置如图2-18所示，按照惠更斯-菲涅耳原理，线光源S上的点将作为次波源向前发射次波（球面波），形成交叠的波场。
在较远的地方放置一观察屏，屏上可以观测到一组几乎是平行的直线条纹。
图杨氏干涉实验原理图2.matlab源代码clearlam=500e-9;a=2e-3;D=1;ym=5*lam*D/a;xs=ym;n=101;ys=linspace(-ym,ym,n);for i=1:nr1=sqrt((ys(i)-a/2).^2+D^2); r2=sqrt((ys(i)+a/2).^2+D^2); phi=2*pi*(r2-r1)./lam;B(i,:)=sum(4*cos(phi/2).^2); endN=255;Br=(B/4.0)*N;subplot(1,2,1)image(xs,ys,Br); colormap(gray(N)); subplot(1,2,2)plot(B,ys)3.实验现象2. 工程光学系列之二：等倾干涉matlab2.1 基本原理等倾干涉是薄膜干涉的一种。
薄膜此时是均匀的，光线以倾角i入射，上下两条反射光线经过透镜作用会汇聚一起，形成干涉。
化学教案权倾一时化学教案内外无不造门者化学教案唯景仁不至试卷试题年三十化学教案方为著作佐郎试卷试题桓玄诛元图等倾干涉薄膜由于入射角相同的光经薄膜两表面反射形成的反射光在相遇点有相同的光程差，也就是说，凡入射角相同的就形成同一条纹，故这些倾斜度不同的光束经薄膜反射所形成的干涉花样是一些明暗相间的同心圆环．这种干涉称为等倾干涉。
倾角i相同时，干涉情况一样（因此叫做“等倾干涉”）2.2 matlab源代码%等倾干涉clear allclose allclc%%k=2000;s=500;D=0.2;bochang=s*10^(-9);theta=0.15;d=k*bochang/4;rMax=D*tan(theta/2);N=501;for i=1:Nx(i)=(i-1)*2*rMax/(N-1)-rMax; for j=1:Ny(j)=(j-1)*2*rMax/(N-1)-rMax; r(i,j)=sqrt(x(i)^2+y(j)^2);delta(i,j)=2*d/sqrt(1+r(i,j)^2/D^2); Phi(i,j)=2*pi*delta(i,j)/bochang; B(i,j)=4*cos(Phi(i,j)/2)^2;endendNCLevels=255;Br=(B/4.0)*NCLevels;figure(1);image(x,y,Br);colormap(gray(NCLevels));2.3 实验现象3. 工程光学系列之三：夫琅禾费矩孔衍射matlab3.1 实验原理衍射的定义：光波在传播过程遇到障碍物时，光束偏离直线传播，强度发生重新分布的现象。
Matlab在物理教学中的应用
%绘 制 F1 和 F2 axis([0 20 0 20]); text(M1,N1,'F1');
%标 注力的名称 text(M2,N2,'F2'); hold on; y1=0:N1/100:N1;
x2=0:M2/100:M2;
Yf2=N2;
plot(x1,Yf1,'k:',x2,Yf2,'k:');
text(0,Yf1,'Yf1');
%标注 Y 上的分力名称
text(0,Yf2,'Yf2');
hold on;
Xh=M1+M2;
%求合力
Yh=N1+N2;
Fh=sqrt(Xh.^2+Yh.^2)
The application Matlab in physics teaching
An YuYu (College of Physics and Information Science,Tianshui Normal University, Tianshui
Gansu 741001) Abstract：With a brief introduction of matlab by concrete examples of matlab in physics teaching on the specific applications. software will be introduced in physics teaching and learning, the contents more image, to simplify the theory and application of reasoning, is computer aid physics teaching produced a better way. Key words: Matlab；Physics teaching ; Application