ch4 数值计算.docVIP

数值计算 与符号计算相比，数值计算在科研和工程中的应用更为广泛。 和十年、二十年前相比较，在当今计算机软硬件的支持下，人们所能拥有的计算能力已经得到了巨大的提升。 数值微积分 近似数值极限及导数 MATLAB数值计算中，没有专门的求极限和导数的指令。 原因：数值精度有限，不能表示无穷小量，不能准确描述一个数的邻域。 【例4.1-1】设，，试用机器零阈值eps替代理论0计算极限，。 x=eps; L1=(1-cos(2*x))/(x*sin(x)), L2=sin(x)/x, syms t f1=(1-cos(2*t))/(t*sin(t)); f2=sin(t)/t; Ls1=limit(f1,t,0) Ls2=limit(f2,t,0) 理论分析表明：， 。 借助符号计算所求的极限与理论值一致； 用数值法近似计算的极限与理论不一致。 提醒：除非数值近似法求的极限经过理论验证，否则绝不要借助数值法求取极限。 函数在点处的导数 【例4.1-2】已知，求该函数在区间 中的近似导函数。 d=pi/100; t=0:d:2*pi; x=sin(t); dt=5*eps; x_eps=sin(t+dt); % 自变量的增量为5eps。 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) 4.1-1 增量过小引起有效数字严重丢失后的毛刺曲线 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-2 增量适当所得导函数比较光滑 结论：自变量增量的选取一定要避免太小。 MATLAB提供的与“求导”概念有关的指令： dx = diff(X) 求差分 当X是向量时，dx = X(2:n) – X(1:n-1)； 当X是矩阵时，dx = X(2:n,:) – X(1:n-1,:)； dx的长度比X的“长度”短一个元素（或一行）。 x = [1 2 3 4 5]; y = diff(x) FX = gradient(F) 求梯度 [FX,FY] = gradient(F) 二元函数的梯度 当F是向量时，FX(1) = F(2) – F(1), FX(end) = F(end) – F(end-1), FX(2:end-1) = (F(3:end) – F(1:end-2))/2, FX的“长度”与F相同。 当F是矩阵时，FX和FY是与F同样大小的矩阵。 FX的每行给出F相应行元素间的“梯度”； FY的每列给出F相应列元素间的“梯度”。 【例4.1-3】已知，采用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) % 宏观上看，diff和gradient所求得近似导数大体相同 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)