四数值计算.PPTVIP

第四章 数值计算 数值计算内容: 4.1.1 多项式构造 4.3.2 差分和累计指令 沿第n维对X元素按模增大排列 [XS,KK]=sort(X,n) 用梯形法沿第n维求函数Y关于x的累计积分 cumtrapz(x,Y,n) 沿第n维求累计和 cumsum(X, n) 沿第n维求累计乘积 cumprod(X,n) 用梯形法沿第n维求函数关于自变量x的积分 trapz(x,Y,n) 沿第n维求和 sum(X, n) 沿第n维求乘积 prod(X, n) 对Z求x,y方向梯度 [DZx,DZy]=gradient(z,dx,dy) 沿第n维求m次差分 (二维：n=1; 列,n=2,行) diff(X, m, n) 五点Laplacian del2(U,hx,hy) 　　　　　　　　　含　　　义 函　　数 diff, trapz, cumtrapz指令的演示 a=[1,2,3,4,5]; b=[a, a+1, a+2,a+3] diff(a)=diff(a,1)=[1,1,1,1] diff(b)=[1,1,1,1 ; 1,1,1,1 ; 1,1,1,1 ; 1,1,1,1] trapz(b)=[12 16 20 24] cumtrapz(b) 0 0 0 0 1.50 2.50 3.50 4.50 4.00 6.00 8.00 10.00 7.50 10.50 13.50 16.50 12.00 16.00 20.00 24.00 1 2 3 4 2 3 4 5 3 4 5 6 4 5 6 7 5 6 7 8 例4.3.2-1　求1+2+3+…+100以及50! x=1:1:100; sum_x=sum(x); sum_x = 5050 a=1:1:50; prod_a=prod(a) prod_a = 3.0414e+064 cumsum, cumprod的用法 上两个函数分别是求向量的累计和和累计乘积.如果a=1:1:n, 则 cumsum(a)=[1, 1+2, 1+2+3, ..., 1+2+3+…+n] cumprod(a)=[1, 1*2, 1*2*3, …, 1*2*3*…*n] 例4.3.2-1 求f(x)=3x2在区间[0,2]的积分 dt=0.001; t=(0:dt:2); y=3*t.^2; s1=dt*sum(y) s2=dt*trapz(y) s=dt*cumsum(y); s3=s(end) s=dt*cumtrapz(y); s4=s(end) matlab还有更精良的积分指令: quad, quda8 s1 = 8.0060 s2 = 8.000001 s3 = 8.0060 s4 = 8.000001 4.4 MATLAB泛函指令 在MATLAB中,凡以函数为输入宗量的指令,都被统称为泛函指令. 最常见的有: z=fzero(fun,x0)---求一元函数零点指令的最简单格式． x=fsolve(fun,x0)---解非线性方程组的最简单格式 x=fminbnd(fun,x1,x2)---求函数在区间(x1,x2)中极小值的指令 最简格式． x=fminsearch(fun,x0)---单纯形法求多元函数极值点指令的最简格式. x=fminunc(fun,x0)---拟牛顿法求多元函数极值点指令的最简格式 a=lsqnonlin(fun,a0)---解非线性最小二乘问题指令的最简格式. q=quad(fun,a,b)---采用递推自适应simpson法计算积分． q=quadl(fun,a,b)---采用递推自适应Lobatto法求数值积分． SS=dblquad(fun,inmin,inmax,outmin,outmax)---二重闭环数值积分 [t,YY]=ode45(fun,tspan,Y0)---采用4,5阶Runge-Kutta方程解算ODE初值问题 指令中被处理的函数fun,可以取:字符串表达式,内联函数, M函数文件的函数句柄 4.4.1 求函数零点 对于任意函数f(x)=0,零点情况复杂,没有一个通用解法.一般来说,零点的数值计算过程是:先猜测一个初始零点或该零点所在的区间,然