CP010-计算物理实验模拟.pptVIP

简单物理实验的模拟 简谐振动实验的模拟--原理 物体受到合力： 物体的运动方程： 物体运动轨迹： 简谐振动实验的模拟--流程图 简谐振动实验的模拟--源代码 function m=cp011 A=12;phai=0; t=(0:200)/200*2*pi; y=A*cos(t+phai); n=length(y); m=moviein(n); for i=1:n plot([0 0], [14 y(i)]); axis([-1 1 -15 15]); axis off; hold on; % 定义 线 色、 点 型（点）、点的大小（ 40 ）、擦除方式（ xor) line(ydata,y(i),xdata,0,Color ,[1 0 0], Marker ,. , MarkerSize ,60, EraseMode , xor); m(:,i)=getframe; drawnow; % 刷新屏幕 21 hold off end movie(m,8,300); 简谐振动实验的模拟--效果 振动合成的模拟--原理 振动方向相同的两个简谐振动： 则它们的合振动为： 当 时： 但一般条件下，合振动的解析式是很难求出，但可以利用计算机模拟合振动。 振动合成的模拟--流程图 振动合成的模拟--源代码 function m=cp012 A1=8;A2=12;w1=2*pi/5;w2=2*pi/7; wy1=2*(A1+A2)+4+A2;wy2=wy1+4+A2+A1;wy3=wy2+A1+4; phai1=pi/3;phai2=pi/2; T=70; t=0:0.1:T; y1=A1*cos(w1*t+phai1); y2=A2*cos(w2*t+phai2); y3=y1+y2; y1=y1+wy2; y2=y2+wy1; y3=y3+(A1+A2); n=length(t); for i=1:n-1 plot([t(i),t(i+1)],[y1(i),y1(i+1)]); axis([0 T 0 wy3]); hold on plot([t(i),t(i+1)],[y2(i),y2(i+1)]); plot([t(i),t(i+1)],[y3(i),y3(i+1)]); pause(0.01); end 振动合成的模拟--效果 驻波的模拟--原理 正向波： 反向波： 驻波： 驻波的模拟--流程图 驻波的模拟--源代码 function m=cp013 A=10;w=2*pi/5;v=2;phai=pi/3; T=70;X=100; t=0:0.1:T; x=0:0.01:X; n=length(t); m=length(x); for i=1:n y1=A*cos(w*(t(i)-x/v)+phai)+54; y2=A*cos(w*(t(i)+x/v)+phai)+54; y=y1+y2-108+20; plot(x,y1); axis([0 X 0 65]); hold on; plot(x,y2); plot(x,y); pause(0.0001); hold off; end 驻波的模拟--效果 光的多缝衍射模拟--原理 单缝衍射在屏上的光强分布： 多缝干涉在屏上的光强分布： 多缝衍射在屏上的光强分布： 光的多缝衍射模拟--流程图 光的多缝衍射模拟--源代码 function m=cp014 a=1.3*(1e-6);d=3.38*(1e-6);n=8;k=pi/(632*(1e-9)); NN=2; shida=[-pi/NN:pi/180/10:pi/NN]; u=k*a*sin(shida); c=k*d*sin(shida); i=(sin(u)./u).^2; ii=((sin(n*c)).^2)./(sin(c).^2)/n/n; iii=i.*ii; subplot(3,1,1) plot(shida,i); title(单缝衍射); axis off; subplot(3,1,2) plot(shida,ii); title(多缝干涉); axis off; subplot(3,1,3) plot(shida,iii); title(多缝衍射