2013年最新Matlab语言及其在工程领域中的应用课件_1.ppt

第二章 MATLAB在数值分析中的应用;t = dat(:,1); CA1 = dat(:,2); CA2 = dat(:,3); plot(t,CA1,o,t,CA2,^) axis([0 t(end)+0.5 min(min(dat(:,2:3)))-0.1 max(max(dat(:,2:3)))+0.05]) xlabel(Time (min)) ylabel(C_A_1, C_A_2 (kmol/m^3)) legend(C_A_1,C_A_2) title(Concentration profiles); 在平面上给定n个点(xk,yk),可以唯一确定一个最多n-1次的多项式通过这些点，这个多项式叫插值多项式;插值多项式例子;Polyinterp (Lagrange 插值形式);Lagrange插值多项式基函数;x0 =0:3; y0 = [-5 -6 -1 16]; x = -.25:.001:3.25; %-------------------------------- tic A = vander(x0); X= A\y0; y1=polyval(X,x); t1=toc %-------------------------------- tic y2 = polyinterp(x0,y0,x); t2=toc %-------------------------------- plot(x0,y0,o,x,y1,-,x,y2,b) legend(data,vand method,polyinterp);;x00=data(:,1); y00=data(:,2); x0 =da00(:,1); y0 =da00(:,2); x =0:.1:9; %-------------------------------- tic A = vander(x0); X= A\y0; y1=polyval(X,x); t1=toc %-------------------------------- tic y2 = polyinterp(x0,y0,x); t2=toc %-------------------------------- plot(x00,y00,*,x0,y0,o,x,y1,-,x,y2,b) legend(da00,data,vand method,polyinterp);Polyinterp(symbol);0;x0=[-5:1:5]; y0=1./(1+x0.^2); x=[-5:0.1:5]; y=polyinterp(x0,y0,x); y1=1./(1+x.^2); plot(x0,y0,r--) hold on plot(x,y1,b-);Hermite插值;function y=hermite(x0,y0,y1,x);分段插值;0;k(x = x0(j)) = j;;x0=[-5:1:5]; y0=1./(1+x0.^2); x=[-5:0.1:5]; y=interp1(x0,y0,x,’linear’); y1=1./(1+x.^2); y2= piecelin(x0,y0,x) plot(x0,y0,‘^) hold on plot(x,y,’o’,x,y1, r,x,y2,’-’);0;三次样条;在三次样条中，要寻找三次多项式，以逼近每对数据点间的曲线。在样条术语中，这些数据点称之为断点。因为，两点只能决定一条直线，而在两点间的曲线可用无限多的三次多项式近似。因此，为使结果具有唯一性。在三次样条中，增加了三次多项式的约束条件。通过限定每个三次多项式的一阶和二阶导数，使其在断点处相等，就可以较好地确定所有内部三次多项式。此外，近似多项式通过这些断点的斜率和曲率是连续的。然而??第一个和最后一个三次多项式在第一个和最后一个断点以外，没有伴随多项式。因此必须通过其它方法确定其余的约束。最常用的方法，也是函数spline所采用的方法，就是采用非扭结(not-a-knot)条件。这个条件强迫第一个和第二个三次多项式的三阶导数相等。对最后一个和倒数第二个三次多项式也做同样地处理。 基于上述描述，人们可能猜想到，寻找三次样条多项式需要求解大量的线性方程。实际上，给定N个断点，就要寻找N-1个三次多项式，每个多项式有4个未知系数。这样，所求解的方程组包含有4*(N-1)个未知数。把每个三次多项式列成特殊形式，并且运用各种约束，通过求解N个具有N个未知系数的方程组，就能确定三次多项式。这样，如果有50个断点，就有50个具有50个未知系数的方程组。幸好，用稀疏矩阵，这些方程式能够简明地列出并求解，这就是函数spline所使用的计算未知系数的方法。; x = 1:6; y = [16 18