5.3.6微分方程的数值积分-Read.ppt

　　从以上的结果可以看出，傅里叶级数在函数的间断点附近误差是比较大的，由此也可以理解为什么狄利克雷条件与函数的间断点数量有关。对于连续平滑的曲线，不必用很高阶次的傅里叶级数就能较好地近似描述它。 　　【例 5-4-7】方波可用相应频率的基波及其奇次谐波合成， 这也是将它展开为正弦级数的出发点。现要求用MATLAB来演示这一结论。 　　解： 　　◆ 建模 　　一个以原点为奇对称中心的方波y(t)可以用奇次正弦波的叠加来逼近， 即 　　方波宽度为π， 周期为2π， 可以用MATLAB程序来检验这种逼近的程度和特征。 　　◆ MATLAB程序 t=linspace(0，pi，501); %设定时间数组，它有501个点 y=sin(t); plot(t，y)，pause %频率为w=1(f=1/2π)的正弦基波 y=sin(t) + sin(3*t)/3; plot(t，y)，pause %叠加三次谐波用1， 3， 5， 7， 9次谐波叠加 y=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9; plot(t，y)，pause %为了绘制三维曲面， 要把各次波形数据存为一个三维数组， 所以必须重新定义y， 重新编程。由于打算求到19次谐波， 因此开始就把点取得较密 y=zeros(10，max(size(t))); x=zeros(size(t)); for k=1∶2∶19 x=x + sin(k*t)/k; y((k+1)/2,: )=x; %叠加各次谐波 end %将各波形叠合绘出 pause，plot(t，y(1∶9,∶))， %将各波形绘成三维网格图， 看出增加谐波阶次对方波逼近程度的影响 pause， mesh(y，k，t)，pause 　　◆ 程序运行结果 　　程序运行中将出现多幅画面，这里只给出最后得出的三维曲面图， 如图5-41所示。取的阶次愈高， 愈接近于方波， 但总是消除不了边缘上的尖峰， 这称为吉布斯效应。可以在命令窗内再键入命令 rotate3d 然后用鼠标拖动三维曲面旋转， 就容易看清吉布斯效应。 　　通过本例， 说明傅里叶级数的阶数愈高，其合成的波形愈接近于原函数波形。该例题还解释了在原函数的间断点附近出现的吉布斯效应。 图 5-41 不同次数谐波叠合的三维曲面 本节习题 　　1. 对下列各序列，计算并画出其前25项，判断它们是否收敛。若收敛，极限L是多少？ 　　(a)　　　　　 (b) 　　(c)　　 (d)  　　2. 画出级数　　　　　　的部分和分布图。 　　3. 对下列级数，(i) 求前25项的和; (ii) 求N→∞时级数的和。 　　(a)  　　(b)  　　　　　(c)  　　4. 设序列 　　 。 　　(i) 求出并画出它前400个点的图形; 　　(ii) 求出其前400个点范围内的部分和 　　　　　　 ，并画出(k，sk)曲线; 　　(iii) k=355时，sk与数355/113有什么关系？ k达到多少时又会出现类似的现象？ 　　5. 求幂级数　　　　　　　的收敛域及和函数。 　　6. 设 ，(b) f(x)=(1+x)ln(1+x)，求f(x)的5阶和 10阶麦克劳林多项式，把两个多项式和函数的图形画在同一坐标系中。 　　7. 对下列函数 　　(a) 　　(b) f(x)=(1+x)3/2， ≤x≤2 　　(c) f(x)=cos x sin 2x，|x|≤2 (d) f(x)=e-x cos2x，|x|≤1 　　(i) 在x=0处求出一次、二次和三次泰勒多项式p1(x)，p2(x)，p3(x); 　　(ii) 在给定的区间把函数和三个近似多项式的图形画在同一幅图上; 　　(iii) 讨论在给定的区间中，用三种近似式替代原函数的最大误差是多少？ 　　8. 对下列周期性函数（周期为［-π，π］)，求出其0~9阶傅里叶级数的系数。 (a) (b) (c) (d) f=sint cost (e) f=tcost 　　9. 对上题中的傅里叶系数取0~3阶， 组成傅里叶级数，画出其图形， 与上题的原始图形作比较， 试