第3章 MATLAB数值运算培训讲解.pptVIP

MATLAB基础及其应用教程 第3章: MATLAB数值运算;第3章 提纲; ;3.1.2 多项式的四则运算;1、多项式相加、减（polyadd）;2、多项式乘法（conv）;例5：求 的“商”及“余”多项式。;3.1.3 多项式求值、求根运算和多项式构造;例10：用多项式的根构造多项式;多项式函数;3.2 插值和拟合;1、多项式插值（interp1）;; x=0:10;y=cos(x); xi=0:.25:10; y0=cos(xi); %精确值 y1=interp1(x,y,xi); %线性插值结果 y2=interp1(x,y,xi,cubic); %三次方程式插值结果 y3=interp1(x,y,xi,spline); %三次样条式插值结果 plot(xi,y0,o,xi,y1,xi,y2,-.,xi,y3) %画出函数计算值和3种插值比较图;3种插值方法比较图;2、多项式拟合（polyfit）;3～5阶多项式拟合曲线图;插值拟合函数;3.3.1 微分和差分（diff）; S1=6*x^3-4*x^2+b*x-5; S2=sin(a); S3=(1-t^3)/(1+t^4); diff(S1); diff(S1,2); diff(S1,’b’); diff(S2); diff(S3);;3.3.2 数值积分（cumsum、trapz、quad、quadl）;梯形法数值积分;辛普森数值积分;科茨数值积分:采用高阶自适应 Lobatto 方法求积分;MATLAB 中二重积分和三重积分分别由函数 dblquad() 和函数 triplequad() 来实现。首先介绍函数 dblquad()，该函数的基本格式如下： q = dblquad(fun,xmin,xmax,ymin,ymax)，函数的参数分别为函数句柄、两个自变量的积分限，返回积分结果。 q = dblquad(fun,xmin,xmax,ymin,ymax,tol)，指定积分结果的精度。 triplequad() 函数的调用格式和 dblquad() 基本相同，在调用 triplequad() 函数时，需要六个参数指定积分限。;3.4 线性方程组的数值解;例15：求解线性方程组;消去法;2、迭代法;function tx=jacobi(A,b,imax,x0,tol) del=10^-10; tx=[x0];n=length(x0); for i=1:n dg=A(i,i); if abs(dg)del disp(diagonal element is too small); return end end ;例18：求解线性方程组;Gauss-Seidel迭代法;function tx=gseidel(A,b,imax,x0,tol) del=10^-10; tx=[x0];n=length(x0); for i=1:n dg=A(i,i); if abs(dg)del disp(diagonal element is too small); return end end ;例18：求解线性方程组;SOR迭代法;function tx=gseidel(A,b,imax,x0,tol) del=10^-10; tx=[x0];n=length(x0); for i=1:n dg=A(i,i); if abs(dg)del disp(diagonal element is too small); return end end ;例19：求解线性方程组;3.5 稀疏矩阵;例20：创建4×5阶矩阵的稀疏矩阵;例21：创建矩阵A的对角稀疏矩阵;2、稀疏矩阵的存储;3.6 常微分方程的数值解;1、泰勒（Taylor）公式;例23：用前向欧拉法和改进的欧拉法解初值问题;3、Runge-Kutta法;（3）四阶Runge-Kutta公式;function f=exam2fun(x,y) f=-y-x*y.^2; f=f(:);;思考题：某厂房容积为 。经测定，空气中含有0.2%的二氧化碳。开动通风设备，以 的速度输入含有0.05%二氧化碳的新鲜空气，同时又排出同等数量的室内空气。求30