数值计算：近似数值极限与导数示例分析.pdfVIP

第4章数值计算

4.1.1近似数值极限及导数

【例4.1-1】

x=eps;

L1=(1-cos(2*x))/(x*sin(x)),%

L2=sin(x)/x,%

L1

0

L2

1

symst

f1=(1-cos(2*t))/(t*sin(t));

f2=sin(t)/t;

Ls1=limit(f1,t,0)

Ls2=limit(f2,t,0)

Ls1

2

Ls2

1

【例4.1‑2】

(1)dpi/100;t=0:d:2*pi;x

sin(t);dt=5*eps;%x_epssin(t+

dt);dxdt_eps(x_eps‑x)/dt;%

plot(t,x,LineWidth,5)holdon

plot(t,dxdt_eps)holdoff

legend(x(t),dx/dt)xlabel(t)

第4章数值计算

4.1.1近似数值极限及导数

【例4.1-1】