实验六MATLAB在普通物理中的应用.docVIP

实验六：MATLAB在普通物理中的应用 一、实验目的 掌握用MATLAB解决普通物理中的难题。 二、实验仪器 1、计算机 2、MATLAB 软件环境 三、实验内容 1、一质点在xy平面内运动，运动方程x=x(t),y=y(t),求质点的轨迹方程并画出轨迹图。 (分别2组运动方程x=3t,y=2t^2+3；x=tcos(t)，y=tsin（t）运算) 2、设电荷均匀分布在从x=-L到x=L，通过原点的线段上，其密度为q(单位为C/m)， 试求出在z-y平面上的电位分布。 3、用毕奥-萨伐定律计算位于x-z平面上的电流环在x-y平面上产生的磁场分布。 4、振动的合成及拍频现象 5、用MATLAB程序来计算演示光的单缝衍射现象 四、实验过程 1．函数文件：zzz clear,close all x=input(x(t)=,s); y=input(y(t)=,s); tf=input(tf=); Ns=100; t=linspace(0,tf,Ns); xPlot=eval(x); yPlot=eval(y); end, plot(xPlot, yPlot) zzz x=3*t y=2*(t.^2)+3 tf=20 2．建立函数文件：zzz： E0=8.85e-12; C0=1/4/pi/E0 ; q=input(q=);L=input(L=);N=input(N=);Nr=input(Nr=); L0=linspace(-L, L, N+1); L1=L0(1:N); L2=L0(2:N+1); Lm=(L1+L2)/2; dL=2*L/N; R=linspace(0, 10, Nr+1); for k=1:Nr+1 Rk=sqrt(Lm.^2+R(k)^2); Vk=C0*dL*q./Rk; V(k)=sum(Vk); end [max(V), min(V)] plot(R, V), grid zzz q=1 L=5 N=50 Nr=50 ans = 1.0e+010 * 9.3199 0.8654 3．建立函数文件：zzz: clear all; mu0=4*pi*1e-7; I0=1;Rh=1; C0=mu0/(4*pi) * I0; NGx=20; NGy=20; x=linspace(-3, 3, 20);y=x; Nh=20; theta0=linspace(0, 2*pi, Nh+1); theta1=theta0(1: Nh); x1=Rh*cos(theta1); z1=Rh*sin(theta1); theta2=theta0(2: Nh+1); x2=Rh*cos(theta2); z2=Rh*sin(theta2); dly=0;dlx=x2-x1;dlz=z2-z1; yc=0;xc=(x2+x1)/2;zc=(z2+z1)/2; for i=1:NGy for j=1:NGx rx=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=dlz.*rx-dlx.*rz; Bx(i, j)=sum(C0*dlXr_x./r3); By(i, j)=sum(C0*dlXr_y./r3); end end clf;quiver(x, y, Bx, By); gtext(x),gtext(y) 4.建立函数文件：zzz： clear all 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(2), sound(y2); pause(2), sound(y), pause subplot(1, 1, 1) zzz 振幅1=1.2 频率1=300 振幅2=1.8