《计算物理(本科)》[实验例题].pptVIP

实验一 实验数据的统计处理法 一、统计直方图 合肥工业大学电子科学与应用物理学院 * 直方图：以区间频数和该区间作为边长画成长方形，这些长方形组成的图就是统计直方图。 MATLAB中实现统计直方图的函数是hist(x,y)，x：向量，测量的一组数据；y：分成区间的个数，或自定义区间的向量，缺省默认10个区间。 MATLAB中实现算术平均值的函数是mean(x) x=1:100; mean(x) ans = 50.5000 x=rand(1,100); mean(x) ans = 0.5000 MATLAB中实现计算测量值标准偏差的函数是std(x)，例如： x=1:100; std(x) ans = 29.0115 实验二 实验数据的插值法 x=[-5:1:5]; y=1./(1+x.^2); x0=[-5:0.1:5]; y0=lagrange(x,y,x0); y1=1./(1+x0.^2); plot(x0,y0,--r) hold on plot(x0,y1,-b) 【实例】取n=10,用Lagrange插值进行计算 实验三 实验数据的拟合 在MATLAB中编写程序: x=[0.5 1.0 1.5 2.0 2.5 3.0]; y=[1.75 2.45 3.81 4.80 8.00 8.60]; a=polyfit(x,y,2) x1=[0.5:0.05:3.0]; y1=a(3)+a(2)*x1+a(1)*x1.^2; plot(x, y,*) hold on plot(x1, y1, -r) 运行结果： a = 0.4900 1.2501 0.8560 【例】设y=span{1,x, x2 },用多项式曲线拟合如表所示的数据 x 0.5 1.0 1.5 2.0 2.5 3.0 y 1.75 2.45 3.81 4.80 8.00 8.60 【例】课本中P53例子： 2x1+3x2+x3=9 x1+x2+x3=4 x1-2x2-x3=-4 MATLAB程序： a=[2.0 3.0 1.0 1.0 1.0 1.0 1.0 -2.0 -1.0]; b=[9 4 -4]; x=a\b 运算结果： x = 1.0000 2.0000 1.0000 实验四 线性代数方程组的解法 【例】设初值为0，精度10-6，用Jacobi迭代法求方程组解 10 x1 - x2 =9 - x1+10x2 - 2x3=7 -2x2+10x3=6 MALAB实现程序: a=[10 -1 0; -1 10 -2;0 -2 10]; b=[9;7;6]; x0=[0; 0; 0]; jacobi(a,b,x0) 实验五 数值积分及其应用 解： d=pi/1000; x=0:d:3*pi; y=exp(-0.5*x).*sin(x+pi/6); II=trapz(x,y) 【例】用梯形法求下列积分 结果： II= 0.9008 解：如图所示建立坐标系，那么(x0,y0)处的电势为 【例1】设长为l=0.02(m)的杆上，均匀带有电荷q=16*10-6(库仑)，可认为是线电荷，求在杆所在平面内任意一点的电势。 x y (x0,y0) dx o -l/2 l/2 r 在MATLAB中用Simpson求积法对某一点(x0,y0)处的电势进行计算。 1、建立(x0,y0)的积分函数 %jswlx_7_3_1f.m function v=jswlx_7_3_1f(x) x0=1;y0=1; v=16e-6/(0.02*4*pi*8.85e-12)./sqrt((x-x0).*(x-x0)+y0*y0); 2、在MATLAB窗口中输入 v=quad(‘jswlx_7_3_1f, -0.01,0.01,1e-10) 计算结果： v = 1.0173e+005 实验六 常微分方程数值解法及其物理中的应用 【例3】阻尼振动：求解阻尼振动方程 已知质量m=10，倔强系数k=10，阻尼系数c=2，初速度v0=0，初位置x0=10。 解：化一阶方程组 x’=v , x(0)=10 v’=-0.2v-x , v(0)=0 1、建立微分方程组m文件jswlx_8_4_3f.m function y=jswlx_8_4_3f(t,x) y=[x(2);-0.2*x(2)-x(1)] 2