高等数学实验指导书8

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验八 无穷级数

8.1 实验目的

掌握利用Mathematica 软件进行级数的有关计算的方法; 通过实验进一步熟悉无

穷级数的一些基本概念。

8.2 实验内容

一、 常数项级数 例1 计算∑∞

=+1

)1(1

n n n 。 [实验]输入:

得结果:1

例2 设,,2)1(21

∑==-+=n i i n n n

n u s u 计算1s 、2s 、…、30s ,及∑∞

=1n n u 。

[实验]输入:

得结果:

{0.50000000,1.2500000 ,1.3750000,1.5625000,1.5937500,1.6406250,1.6484375,1.6601563 ,1.6621094 ,1.6650391 ,1.6655273 ,1.6662598 ,1.6663818 ,1.6665649,1.6665955,1.6666412,1.6666489 ,1.6666603 ,1.6666622 ,1.6666651 ,1.6666656 ,1.6666663,1.6666664,1.6666666,1.6666666 ,1.6666666 ,1.6666666 ,1.6666667 ,1.6666667 ,1.6666667} 3

例3 设,,1

1

∑===n

i i n n u s n u 计算10s 、100s 、1000s 、10000s 、100000s 、200000s 、

10ln 10-s 、100ln 100-s 、1000ln 1000-s 、10000ln 10000-s 、100000ln 100000-s 、

200000

ln 200000-s ,及∑∞

=1

n n u 、∑=∞

→-k

n n k k u 1

)ln (lim 。 [实验]输入:

得结果:

2.92897

5.18738 7.48547 9.78761

12.0901 12.7833 0.626383 0.582207 0.577716 0.577266 0.577221 0.577218

再输入:

得结果:∞

0.577215664901532860606512090083

例4 设,,1

)

1(1

1

∑=-=-=n

i i n n n u s n u 计算10s 、100s 、1000s 、10000s 、100000s 、200000s ,及∑∞

=1n n u 。

[实验]输入:

得结果:

0.645635

0.688172 0.692647

0.693097

0.693142

0.693145

Log[2]

0.693147180559945309417232121458

二、幂级数

例5 求幂级数∑∞

=-

1

)1 (

n

n

n

n

x的和函数。

[实验]输入:

得结果:-Log[1+x]

例6 求幂级数∑∞

=+

01

n

n

n

x的和函数。

[实验

]输入:

得结果:

再输入:Limit [f , x→0]

得:1

例7 求sin x在x0=0处的19阶幂级数展开式;再分别取n=1,3,5,…,19,作出sin x在x0=0处的n阶幂级数展开式及sin x的图象,观察sin x的幂级数逼近sin x的效果。

[实验]输入:Series[Sin[x],{x,0,19}]

得结果:

再输入:

得结果:

双击其中任一幅图象,便可观察到sin x 的幂级数逼近sin x 的效果。

说明:Series[f , x , x 0, n

] generates a power series expansion for f .

三、 傅里叶级数

例8 设)(x f 是周期为π2的周期函数,它在),[ππ-上的表达式为⎩⎨

⎧<≤<≤--=.

0,

1,

0,1)(ππx x x f 可将)(x f 展开成傅里叶级数,请作出它们

的有关图像,并观察该傅里叶级数逼近)(x f 的效果。

[实验]输入:

f[x_]:=Which[x<-2 Pi,-1,x<-Pi,1,x<0,-1,x=2Pi,1]; b[n_]:=(1-(-1)^n)*2/n/Pi; For[i=1,i<=31,i+=5,

s[x_]:=Sum[b[n]*Sin[n*x],{n,1,i}]; Plot[{f[x],s[x]},{x,-3Pi,3Pi},PlotStyle → {RGBColor[0,1,0],RGBColor[1,0.3,0.5]}]]; 得结果:

双击其中任一幅图象,便可观察到该傅里叶级数逼近)

f的效

(x

果。

相关文档
最新文档