数值计算-Read.DOCVIP

数值计算 数值微积分 近似数值极限及导数 MATLAB数值计算中，没有专门求极限和导数的指令。原因是：有限精度。但这不意味：数值计算不能应用于与“导数”等函数邻域概念有关的问题。事实上，数值计算是解各类微分方程的最主要途径。 dx=diff(X) 求差分 FX = gradient(F) 求一元（函数）梯度 [FX,FY] = gradient(F) 求二元（函数）梯度 【例4.1-1】用数值计算法求极限。 L0_p=sin(eps)/eps L0_n=sin(-eps)/(-eps) 【例4.1-2】已知，求。 dfdx_p=abs(eps)/eps dfdx_n=abs(-eps)/(-eps) 【例4.1-3】已知，求该函数在区间 中的近似导函数。 （1）自变量的增量取eps数量级 d=pi/100; t=0:d:2*pi; x=sin(t); dt=5*eps; %增量为eps数量级 x_eps=sin(t+dt); dxdt_eps=(x_eps-x)/dt; %数值导数 plot(t,x,LineWidth,5) hold on plot(t,dxdt_eps) hold off legend(x(t),dx/dt) xlabel(t) （2）自变量的增量取得适当 x_d=sin(t+d); dxdt_d=(x_d-x)/d; %以d=pi/100为增量算得的数值导数 plot(t,x,LineWidth,5) hold on plot(t,dxdt_d) hold off legend(x(t),dx/dt) xlabel(t) 【例4.1-4】已知，采用diff和gradient计算该函数在区间 中的近似导函数。 clf d=pi/100; %自变量增量 t=0:d:2*pi; x=sin(t); dxdt_diff=diff(x)/d; dxdt_grad=gradient(x)/d; subplot(1,2,1) plot(t,x,b) hold on plot(t,dxdt_grad,m,LineWidth,8) plot(t(1:end-1),dxdt_diff,.k,MarkerSize,8) axis([0,2*pi,-1.1,1.1]) title([0, 2\pi]) legend(x(t),dxdt_{grad},dxdt_{diff},Location,North) xlabel(t),box off hold off subplot(1,2,2) kk=(length(t)-10):length(t); hold on plot(t(kk),dxdt_grad(kk),om,MarkerSize,8) plot(t(kk-1),dxdt_diff(kk-1),.k,MarkerSize,8) title([end-10, end]) legend(dxdt_{grad},dxdt_{diff},Location,SouthEast) xlabel(t),box off hold off 数值求和与近似数值积分 Sx=sum(X) 沿列方向求和 Scs=cumsum(X) 沿列方向求累计和 St=trapz(x, y) 采用梯形法沿列方向求函数y关于自变量x 的积分 Sct=cumtrapz(x, y) 采用梯形法沿列方向求函数y关于自变量x 的累计积分 【例 4.1-5】求积分，其中。 clear d=pi/8; %子区间间隔 t=0:d:pi/2; %包含5个采样点的一维数组 y=0.2+sin(t); %5个点处的函数值数组 s=sum(y); s_sa=d*s; s_ta=d*trapz(y); disp([sum求得积分,blanks(3),trapz求得积分]) disp([s_sa, s_ta]) t2=[t,t(end)+d]; y2=[y,nan]; stairs(t2,y2,:k) hold on plot(t,y,r,LineWidth,3) h=stem(t,y,LineWidth,2); set(h(1),MarkerSize,10) axis([0,pi/2+d,0,1.5]) hold off shg 计算精度可控的数值积分 S1=quad(fun,a,b,tol) 采用递推自适应Simpson法计算积分 S1=quadl(fun,a,b,tol) 采用递推自适应 Lobatto法求数值积分 S2=dblquad(fun,xmin,xmax,ymin,ymax,tol) 二重（闭型）数值积分指令 S3=tr