计算物理学pptchp4.ppt

计算物理 (Computational Physics) 非线性曲线拟合 目录 实验数据的拟合 用 k (k n) 次多项式拟合, 拟合函数为： y(x) = A0 + A1 x + A2 x2 + ··· + Ak xk 最小二乘原则： 测量值： x x1 x2 … xn y y1 y2 … yn 计算物理 (Computational Physics) 非线性曲线拟合 目录 实验数据的拟合 矩阵形式： 解得： 令： 计算物理 (Computational Physics) 目录 实验数据的拟合 function f=ls(x,y,k) % 最小二乘曲线拟合,生成k次多项式 sx=size(x); %size( ) 显示矩阵大小 if sx(1)==1 x=x; end sy=size(y); if sy(1)==1 y=y; end n=max(sx); if length(y)~=n error(Vectors must be the same lengths); elseif (n-1)k error(k is large); %error() 程序到此结束,显示引号中内容 end X1=ones(n,1); X2=repmat(x,1,k); %repmat(), 由向量生成矩阵 X=cumprod([X1 X2],2); %cumprod(), 生成连积矩阵 A=(X*X)\(X*y); dx=(x(n)-x(1))/(n*32); xx=x(1):dx:x(n); m=length(xx); XX1=ones(m,1); XX2=repmat(xx,1,k); XX=cumprod([XX1 XX2],2); f=XX*A; plot(xx,f,g); hold on plot(x,y,.,Markersize,14) x=0:10; y=cos(x); ls(x,y,5); 计算物理 (Computational Physics) 目录 实验数据的拟合 size( ) 显示矩阵大小 A=[1 2 3;4 5 6]; size(A) ans = 2 3 error() 程序到此结束,显示引号中内容 x=0:10; y=[1 2 3]; ls(x,y,5); ??? Error using == ls Vectors must be the same lengths repmat(), 由向量生成矩阵 x=[1 2 3]; repmat(x,2,1) ans = 1 2 3 1 2 3 cumprod( ), 生成连积矩阵 B = cumprod(A, dim) 此调用格式返回A中由标量 dim所指定的维数的累加积: A=[1 2 3;4 5 6] A = 1 2 3 4 5 6 cumprod(A) %沿着第一维（列）的累加积 ans = 1 2 3 4 10 18 cumprod(A,2) %沿着第二维（行）的累加积 ans = 1 2 6 4 20 120 计算物理 (Computational Physics) 非线性曲线拟合 目录 实验数据的拟合 Polyfit 多项式拟合命令 格式： p=polyfit(x,y,n) 其中： p为返回的多项式向量，从最高次幂到最低次幂的系数x为准备拟合的自变量，y为应变量n为拟合的次数 计算物理 (Computational Physics) 非线性曲线拟合 目录 实验数据的拟合 例： cos(x)的曲线拟合 x=0:10; y=cos(x); a=polyfit(x,y,5) ; x1=0:0.05:10; y1=a(6)+a(5)*x1+a(4)*x1.^2+a(3)*x1.^3+a(2)*x1.^4+a(1)*x1.^5; plot(x,y,*) hold on plot (x1,y1,-r) 计算物理 (Computational Physics) 非线性曲线拟合 目录 实验数据的拟合 多项式拟合的异常现象-龙格现象 利用多项式对某一函数的近似逼近，这样，利用多项式就可以计算相应的函数值。一般情况下，多项式的次数越多，需要的数据就越多，而预测也就越准确。但是当多项式的方次较高时，系数的微小差异有时会引起解的巨大变化。 函数 在区间[-5,5]上做多次多项式插值插值。 计算物理 (Computational Physics) 非线性