用计算机模拟实验的创新点,大学物理计算机模拟实验的MATLAB开发(1)

由(1)-(4)式,并考虑到几何关系cos

θ=(x 0-vt )/(x 0-vt )2+y 02

编程如下,能模拟声源发出的信号和听者接

收到的信号。

clear

x0=500;y0=30;v =50;

c =340;w1=1000;t =0:0.001:30;u1=sin (w13t )+sin (1.13w13t );

sita =aco s ((x0-v 3t )/sqrt ((x0-v 3t ).^2+y0^2));

w2=c 3w1/(c -v 3cos (sita ));t2=t +sqrt ((x0-v 3t ).^2+y0^2)/c ;u2=sin (w2.3t2)+sin (1.13w2.3t2);sound (u1);　%模拟声波信号pause (5);%暂停5秒

sound (u2)

%模拟听者接收到的信号

以此为平台,可以引申出声源和听者相对运动的多种情况。2　磁场中电子运动的模拟

以一定速度射入磁场中的电子,通过理论计算可以求出它在磁场中的运动轨迹,而计算机模拟实验能动态地展示电子在磁场中的运动情况。

设一初位置为(0,0,0),初速度为v 0=2i +5j +105k (m/s )的电子进入的匀强磁场B =0.1k

(T )中,忽略重力作用,电子在运动过程中所受的

洛仑兹力为

F =-e v ×B =-ev y B i +ev x B j

(5)

式中e =1.6×10-19C 为电子所带电量的数值。依据牛顿第二定律,有

-ev y B i +ev x B j =m d v

d t =m d v x d t i +m d v y d t

j

(6)

式中m =9.1×10-31Kg 为电子质量,v x =d x d t ,v y =d y d t ,d v x d t =d 2x d t 2,d v y d t =d 2y d t 2,于是(6)式可写成分量方程形式

d 2x d t

2=-eB m d y d t d 2y

d t 2=

eB m d x d t d 2z d t 2

=0(7)

由(7)式和电子运动的初始条件,先解微分方程,然后利用动画功能展示电子在磁场中的螺旋运动。编程如下,其运动轨迹如图2所示.

图2　电子在磁场中的运动轨迹

[x ,y ]=dsolve (′D2x =-(e 3B/m )

3

Dy ′,

′D2y =(e 3B/m )3Dx ′,...

′x (0)=0,Dx (0)=2′,′y (0)=0,Dy (0)=5′

);　%解微分方程组X =char (x );Y =char (y );

X1=vectorize (X );Y1=vectorize (Y );

e =1.6e -19;m =9.1e -31;B =0.1;Vz =1e5;t =linspace (0,1.15e -9,100);X2=eval (X1);Y2=eval (Y1);Z =Vz 3t ;comet3(X2,Y2,Z ) %动画展示电子的

运动3　小球弹跳运动的模拟

弹跳球是指小球被抛出后,运动一段时间落地,在与地面接触的过程中,损失一定的能量,再向上弹起,如此循环至小球最终静止在地面上。通过计算机作图显示小球的高度和速度的变化情况,研究小球的运动演化过程。

设小球在离地面10m 处以15m/s 的初速度竖直向上抛。设落地时小球的弹性系数为0.8.打开SIMUL IN K 工具箱,选取两个积分模块、一个增益模块、一个常数模块、一个示波器模块和一个用来设置起始值的模块构成如图3所示的框图,并设置相应的参数,双击Scope 图标,然后点击Start 开始模拟,呈现图4和图5的图形,分别直观地展示小球整个运动过程中高度和速度的变化情况。

6