由(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