数值分析助教讲得MATLAB课件.pptVIP

* * * * * * * * Hust 4.2 多项式的运算 对多项式的运算主要包括求多项式的值、四则运算和求根等。 多项式的求值 MATLAB 7.0提供了两个函数来对多项式进行求值，即polyval和polyvalm。前者以数组为计算单位，后者以矩阵为计算单位。 Y=polyval(P,X)计算以向量P为系数的多项式在点X的值。如果X是矩阵或向量，那么该命令将对矩阵或向量的每一个元素都进行计算。 Y=polyvalm(P,X)计算以P为系数的多项式在方阵X的值。X必须为方阵。 Hust 求多项式的根 在MATLAB 7.0里，多项式由一个行向量表示，设为p,它的系数按降序排列，使用roots函数可以求出该多项式的根。其使用格式为roots(p) 多项式的四则运算 加法和减法 对多项式 的加法和减法， MATLAB 7.0不提供专用的函数。如果两个多项式的向量阶数相同，标准的数组加法有效。当两个多项式的向量阶数不同时，需要在低阶多项式的前边补0，使得它与相加的高阶多项式有相同的阶数。 乘法 使用conv函数对多项式进行乘法运算。其使用格式为c=conv(a,b)，其中a和b为两个多项式的系数向量，c为相乘所生成的多项式的系数向量。 Hust 除法 使用decon函数来完成除法，具体使用类似于乘法。 求导和积分 分别使用函数polyder和polyint来求多项式的导数与积分。 若P为多项式的系数向量，polyder(P)对该多项式求导并返回求导后的系数向量。 polyder(A,B)相当于polyder(A*B) [Q,D]= polyder(B,A)返回多项式B/A的求导值，并以Q/D的形式表示出来 polyint(P,K)返回多项式P的积分，以K为积分常数 polyint(P)返回多项式P的积分，以K=0为积分常数 Hust 拟合和插值 ⒈ 最小二乘法实现曲线拟合 使用polyfit函数求解最小二乘曲线拟合问题，使用方法如下： polyfit(x,y,n)用最小二乘法对所给数据进行n阶多项式拟合，返回拟合多项式p(x)，使得p(x(i))~=y(i) [p,s]=polyfit(x,y,n)不仅返回多项式的系数p，还返回用函数polyval获得的误差分析报告。如果数据中的错误数据y服从独立正态分布，那么polyval函数对错误数据的预报率为50% [p,s,mu]=polyfit(x,y,n)返回拟合多项式的系数xhat=((x-mu(1))/mu(2))，这里mu(1)=mean(x)，mu(2)=std(x) 当n=length(x)时，或是x被重复赋值，系统将产生出错信息 Hust 2. 曲线插值 Lagrange插值 Hermite插值 这两种插值没有专用的函数，需要编写专门的函数或M文件 三次样条插值 在MATLAB 7.0中，实现基本的三次样条插值的函数有spline、ppval、mkpp和unmkpp。 pp=spline(X,Y)返回精确的三次样条插值多项式 yy=spline(X,Y,xx)相当于yy=ppval(spline(X,Y),xx)，使用三次样条法求出yy，yy是隐函数Y在向量xx所在点的值。X和Y是已知的数值点。如果Y是一个矩阵，那么该命令将对Y的每一列进行计算。 Hust 例：三次样条插值的使用 解：首先使用YY=spline(X,Y,XX)命令如下。 x=0:12; y=tan(pi*x/25); xi=linspace(0,12); yi=spline(x,y,xi); plot(x,y,o,xi,yi),title(Spline fit) 所绘制的三次样条插值曲线如下图所示 Hust 该方法适合于只需要一组内插值的情况。不过，如果需要从相同数据集里获取另一组内插值，再次计算三次样条系数是没有意义的。在这种情况下，可以调用仅带前两个参量的pp=spline(X,Y)。继续在命令窗口中输入如下命令，并按Enter键确认。 x=-4:4; y=[0 .15 1.12 2.36 2.36 1.46 .49 .06 0]; cs=spline(x,[0 y 0]); xx=linspace(-4,4,101); plot(x,y,o,xx,ppval(cs,xx),-); 该程序所绘制的插值曲线如下图所示 Hust 1. 基本的绘图命令 1.1 基本的绘图操作 绘制基本线形图的函数表 图形处理 函数名 功能描述 plot 在x轴和y轴都按线形比例绘制二维图形 plot3 在x轴、y轴和z轴都按线形比例绘制三维图形 loglog 在x轴和y轴按对数比例绘制二维图形 semilogx 在x轴按对数比例，y轴按线形比例绘制二维图形 semil