![薛定宇教授大讲堂(卷Ⅱ):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=1738888111-ZhBCkKRquyDb3uJP2Fpzm48LuWQBHNId-0-2362daf0c8bb25a49f57d75e434366fa)
这类序列称为函数项序列(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=1738888111-NFvPrGH9hBisgbZDJN0N9fck9lizLALd-0-7252c89ebced02581bcab8760942fdd4)
解 可以考虑用循环方式计算该序列的前n=40项,序列的各项用火柴杆图的形式显示出来,如图2-15所示。可以看到,序列的值随着n的增加将缓慢减小。
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23117.jpg?sign=1738888111-WQHW1y2qHPzDMOY5UobWiWOSYwimbORd-0-42d65be412fd0b4cc5c9dca0acd68177)
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23118.jpg?sign=1738888111-1u6nQ3Or4QYG9j05XQ5eo0S84MBrQJ5Q-0-8d26746b90f7a3752da7d45db4e1b9fd)
图2-15 序列的演化趋势
可以预见,如果n足够大,这个序列将收敛到某个值上。从理论上可知,这个收敛的值是Euler常数γ。如果计算到第10000项,这个序列的值为0.577265664068165,比较接近于γ=0.5772156649015328606065120900824。
![](https://epubservercos.yuewen.com/30E2CE/16499866805000006/epubprivate/OEBPS/Images/Figure-P37_23121.jpg?sign=1738888111-RL7t0zXeusmn50heqqJy98tzqj8t0QJ5-0-d7467eb8258cc907d651eff287ae28a1)
上面语句中,如果加到10000000000项,则可以得出0.577215665057043。