《计算物理(研究生用)》[第2篇].ppt

§2.1 实验数据的统计评价;正态分布的概率密度曲线如图所示，阴影部分面积表征数据的值落在该范围内的概率。理论可以证明， 在区间（μ-σ，μ+σ）的面积占曲线下总面积的68.27%。也就是说，x的值落在（μ-σ，μ+σ）区间的概率是68.27%。;二、算术平均值和它的标准误差 ;由误差理论可以证明平均值的标准误差;由误差理论可以证明测量值的标准偏差;§2.2 实验数据的插值法;y(x) 经常取x的多项式，作为代数插值。 构造出的代数多项式y(x)，必须满足下列条件： （1）y(x)是一个不超过n次的多项式。这样有n +1个待定系数，对应n +1个测量点。 （2）在插值点xi (i =0,1,2,…n)上， 确保f(x)的插值函数y(x)与f(x)相等，即 yi=y(xi) i =0,1,2,…n;§2.2.1 拉格朗日（Lagrange)插值法;二、 Lagrange抛物线插值（三点）;三、 Lagrange n次多项式插值;MATLAB中实现Lagrange n次多项式插值：;§2.2.2 分段线性插值;注意：Runge现象说明，在大范围内使用高次插值，逼近的效果往往是不理想的。;二、分段插值的概念;四、分段三次插值;MATLAB实现：分段线性插值在MATLAB中不需要自编程序，自身提供了一个一维分段插值函数： interp1(x,y,xi,’method’) 其中 (x,y)：一组节点；xi：插值点； method：插值方法{可选nearest：线性最近项插值；linear：线性插值（默认）；spline：三次条样插值；cubic：三次插值}。;§2.3 实验数据的拟合;§2.3.1 实验数据的一元线性拟合;由拟合直线得 δi = yi - A0- A1xi ;解得 ;三、一元线性拟合的相关系数 ;§2.3.2 多元线性拟合;偏差的平方和为;§2.3.3 多项式曲线拟合;用k次多项式拟合 Y=A0+A1x+ A2x2+…+ Akxk;注意：多项式拟合的正规方程具有唯一解，但多项式的阶次较高时，系数的微小差异有时可能会引起解的巨大变化。为了避免这种情况发生，多项式的阶次不要过高，一般在5次方以内为宜。;在MATLAB中编写程序，例jswlx_4_3_1.m %jswlx_4_3_1.m 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;§2.4 快速傅立叶变换(FFT);§2.4.1 离散傅立叶变换（DFT）;采用周期2π的函数eikx (i为虚数单位,k=0,1,…,N-1 )的线性组合; 上两式都含因子e-i2π/N，令ω= e-i2π/N;§2.4.2 快速傅立叶变换（FFT）;采样点条件：N=2m （m为正整数）。这一条件不难满足，可以大大简化编程和提高计算效率。;具体算法：;N =2m， F(x)可以分解为偶次幂和奇次幂两部分 F(x) ={f0x0 + f2x2 +…+ fN-2xN-2}+ x { f1x0 + f3x2 + …+fN-1xN-2 } 若令 f even(x2)= f0x0 + f2x2 +…+ fN-2xN-2 f odd(x2)= f1x0 + f3x2 + …+fN-1xN-2 则有 F(x) = f even(x2)+ xfodd(x2) F(-x) = f even(x2)-xfodd(x2);§2.4.3 快速傅立叶变换的应用;【例】给出一组采集数据{fl}=[4,3,2,1,0,1,2,3] ，用FFT算法求{fl}的离散频谱{Fk}。其中k,l=0,1,2,3,4,5,6,7。;三、 MATLAB实现快速傅立叶变换（实例二）;MATLAB实现：程序文件名jswlx_5_3_1.m;§2.5 数值积分及其应用;在MATLAB中矩形和梯形求积分，都有专门的函数。它们分别为 (1)矩形求积 用下面的求和函数来实现 cumsum(X) 其中X：向量；cumsum