MATLAB论文12010245326张金锭.docVIP

《MATLAB语言》课程论文 Matlab在基础物理学中的应用 姓名：张金锭 学号：12010245326 专业：电子信息工程 班级：2010级 指导老师：汤全武 学院：物理电气信息学院 完成日期：2011年12月10日 MATLAB在基础物理学中的应用 (姓名：张金锭 学号12010245326 2010级) [摘要]基于MATLAB强大的绘图和仿真功能,对基础物理学中的力学中机械波的相干叠加与驻波的形成，光学中的双缝干涉、光栅光强分布，电磁学中的等量同号点电荷的电场线的绘制和带电粒子在均匀电磁场中的运动等问题进行仿真, 来简便、直观、高效分析物理问题。 [关键字]MATLAB语言 基础物理学 仿真 图形绘制 一、问题的提出 基础物理学是理工科大学生的一门基础课，其内容描述了自然界的现象与规律，让学生在理解物理规律的基础上去解释现实生活中的一些物理现象，M atlab是美国Mathw orks公司于80年代推出的大型数学软件, 通过多年的升级换代,现在已发展成为集数值计算、符号计算、可视化功能以及诸多的工具箱为一体的大型科学计算软件, 它在美国等西方发达国家已广泛应用于科研院所、工程技术等各个部门, 并成为大学生、研究生必备的工具软件，应用MATLAB强大的数学功能通过实例的建模、编程和仿真过程,来简便、直观、高效分析物理问题。 二、在力学中应用举例 （一）驻波形成—火柴杆图 本例用Matlab来研究机械波的相干叠加与驻波的形成。考察分别沿x轴正向和负向传播的两列相干横波，它们的方程为： (1) (2) 设时间从t=0开始到t=10结束，考察区间为[0,4]，令k=π，则λ=2π/k=2，在考察区间上恰好能观察到两个完整波形。令v=1，A1=-0.4，A2=0.4，则方程为： (3) (4) 现在就可以编程演示机械波的叠加和驻波的形成了。 程序如下： %横波驻波的形成 t=0:0.1:10;x=0:1/15:4; %时间数组和位置坐标数组 for i=1:100 %设置循环 %由右向左传播的行波 x1=x(x=(4-t(i))); %挑出坐标位置大于4-t(i)的质元，将其位置坐标赋给新的变量x1 y11=0.4*sin(pi*(x1+t(i))); %这些位置上的质元在t(i)时刻的位移，赋给变量y11 x2=x(x(4-t(i))); %挑出坐标位置小于4-t(i)的质元，将其位置坐标赋给新的变量x2 y12=x2-x2; %这些位置上的质元在t(i)时刻的位移，赋给变量y12 y1=[y12 y11]; %将所有质元的位移组合到同一个列矢量y1中 %由左向右传播的行波 x3=x(x=t(i)); y21=-0.4*sin(pi*(x3-t(i))); x4=x(xt(i)); y22=x4-x4; y2=[y21 y22]; y3=y1+y2; %质元同时参加两个振动，实现了行波的合成 y=[y1;y2;y3]; %将t(i)时刻向左、向右以及合成的波的质 %元位移存放在矢量y中，以便作图 %分区作图 for j=1:3 subplot(3,1,j) %分区作图 stem(x,y(j,:),b:.) ; %火柴杆图 axis([0,4,-1,1]); %指定坐标轴范围 grid on %开启网格线 end pause(0); %暂停的技巧，否则屏幕由于刷新 %过快而导致不显示任何图象 end 运行结果如图1所示。 图1某时刻行波及驻波的波形 程序在10s内每隔0.1s取一个时刻，在考察区间上每隔1/15取一个质元。 在t4s之前，波前还未走过整个区间，只有部分质元参与了振动，其余质元仍处于静止状态，位移为0，因此需要根据当前时刻t(i)找出已经运动和尚未运动的质元，分别对其位移y11、y12赋值。注意y12=x2-x2和y12=0两种赋值是不等效的，前者是长度与x2相同的矢量，后者只是一个数值（标量）。 为了有明显的对比，程序采用了分区作图，分别绘制了当前时刻的两列行波波形以及