matlab实2：微积分基础.pptVIP

数学实验 微积分基础 实验二：微积分基础 * * 主要内容 函数及其图像 级数及其图像 连乘积及其图像 1.1 函数及其图像 例1 正弦函数的逐步逼近 过原点的直线 y=0.8x 和 y=0.2x。问哪条直线与y=sin(x)过原点的切线最接近。 第一步：画出这三条曲线 x=-pi:0.1:pi; y=sin(x); y1=0.8*x; y2=0.2*x; plot(x,y,r); hold on; plot(x,y1,b); plot(x,y2,g); hold off; 第二步：调整直线 y=0.8x x=-pi:0.1:pi; y=sin(x); y1=0.8*x; %对照“组” y2=0.95*x; %调整直线 斜率 分别为： 0.9 0.95 1 plot(x,y,r); hold on; plot(x,y1,b); plot(x,y2,g); hold off; 第三步：试图找高次逼近 %分析：y=sin(x) 为奇函数，x的幂次取为奇数 x=-pi:0.1:pi; y=sin(x); y1=x-x.^3/5 ; %对照“组” y2=x-x.^3/6 ; %调整x的幂 的系数 plot(x,y,r); hold on; plot(x,y1,b); plot(x,y2,g); legend(y=sin(x),y=x-x.^3/5,y=x-x.^3/6); %加上图例 hold off; %分析：y=sin(x) 为奇函数，x的幂次取为奇数 x=-pi:0.1:pi; y=sin(x); plot(x,y,r,linewidth,2); axis([-4,3,-4,3]); y1=x ; y2=x-x.^3/6 ; y3=x-x.^3/6+x.^5/120 ; y4=x-x.^3/6+x.^5/120-x.^7/(120*6*7); %这就是泰勒展开的头4项 hold on; plot(x,y1,b); plot(x,y2,g); plot(x,y3,m); plot(x,y4,k); legend(y=sin(x),y=x,y=x-x.^3/6, ... y=x-x.^3/6+x.^5/120,... y=x-x.^3/6+x.^5/120-x.^7/(120*6*7)?); %图例 hold off; % 这里可以看出泰勒级数是如何逼近正弦函数的。。。试一试余弦。 第四步：再进一步 比较 例2 正弦函数的叠加 y=sin(x)+sin3x/3+ sin5x/5+sin7x/7+ sin9x/9+… + sinmx/m m=1,3,5,7,9,… 第一步：分别画出前1、2、3 4项之和，对应的曲线 x=-pi:0.1:pi; y=sin(x); plot(x,y,r); hold on; m=9; for k=3:2:m y=y+sin(k*x)/k; end plot(x,y,b); hold off; 第二步：让m逐步变大，如，m=5,50,150,500…,观察相应曲线的变化， 分析它的波形趋势。。。 % 第四步：一般的 1,sin(x),cos(x)都以2pi为周期，它们的组合也是。 试一试， f(x)=b0+a1 sin(x)+b1cos(x)+…+ ak sin(kx)+bkcos(kx)+… 。 例2 正弦函数的叠加 例3 无极限的函数列：y=sin(1/x)在0附近的情况 % 这是在0附近震荡的一条曲线。 % 第一步：画出函数在[-1,1]之间的图像 a=1 x=-a:0.000001:a; % 让区间逐步缩小：a=0.1，0.01，0.001 y=sin(1./x); plot(x,y); title(a=1;x=-a:0.000001:a;y=sin(1./x)); % 第二步：逐步缩小区间，看看震荡的内部。。。 % 第三步：取一些离散点,让 xn - 0 k=1000:8000; x=1./k; % 取一列趋向0的点列 y=sin(1./x); plot(x,y,.); title(k=1000:8000; x=1./k; y=sin(1./x); ); % 看看混沌中的有序和美丽 例3 无极限的函数列 % 第四步：变变步长看看 m=2; %变步长：取另一些离散点m=2,3,5,8, … k=1000:m:8000; x=1./k; % 取一列趋向0的点列 y=sin(1./x); plot(x,y,.); title(m=2; k=1000:m:8000; x=1./k; y=sin(1./x); ); %