![薛定宇教授大讲堂(卷Ⅱ):MATLAB微积分运算](https://wfqqreader-1252317822.image.myqcloud.com/cover/148/29977148/b_29977148.jpg)
上QQ阅读APP看书,第一时间看更新
2.5 序列与函数项序列
定义2-15 按某种规律排列的数称为序列(sequence),又称数列,本书将统称为序列。
在日常生活与科学研究中经常会遇到序列,例如
1,2,…,n,…
20,21,22,…,2n−1,2n,…
在这两个序列中,n与2n−1分别称为序列的通项。在MATLAB下可以用一个变量名描述序列的通项,有了通项,就相当于将序列直接描述出来了。
定义2-16 序列的通项还可以为自变量x的函数,例如
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P36_23111.jpg?sign=1738943165-BSG2rQwq6djZGTcIQYZlofiKIrItcue7-0-82584f7584d127b7092244e0269b54a3)
这类序列称为函数项序列(functional sequence)。
当然,如果用n的函数形式描述序列,将使得其取值更方便。
例2-23 试用MATLAB描述式(2-5-1)中给出的函数项序列,并写出序列的第1128项。
解 其实在MATLAB下要描述函数项序列,只需描述其通项就足够了,该序列的第1128项为4 sin(2255x)/(3383π)。
>> syms n x; f(x,n)=4*sin((2*n-1)*x)/(3*n-1)/pi, f(x,1128)
例2-24 试观察下面的序列,考虑如果n的值增大这个序列将如何变化。
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23115.jpg?sign=1738943165-3zoZ3IHQu8RKZ6n0xYRcRQc9lJlxIVWN-0-913dd49f3fb81a5312f3b7d6dd6087b7)
解 可以考虑用循环方式计算该序列的前n=40项,序列的各项用火柴杆图的形式显示出来,如图2-15所示。可以看到,序列的值随着n的增加将缓慢减小。
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23117.jpg?sign=1738943165-uTYA4yWbUfTkLhavSn96hgbwe0tVSDtj-0-d12276cf3e5d543b1f2f7673168c3a35)
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23118.jpg?sign=1738943165-HbHHhuMWoo4BvdQLH23nzccvqgFHBMTS-0-b42c09d78f677156ef363531c9c0b195)
图2-15 序列的演化趋势
可以预见,如果n足够大,这个序列将收敛到某个值上。从理论上可知,这个收敛的值是Euler常数γ。如果计算到第10000项,这个序列的值为0.577265664068165,比较接近于γ=0.5772156649015328606065120900824。
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23121.jpg?sign=1738943165-hmFpBlp1naS0vtfQZyAhDeE0KHa9Ue4J-0-c3e0f86496865e078b2ebfd27fd46e2e)
上面语句中,如果加到10000000000项,则可以得出0.577215665057043。