Matlab_级数

级数

01 求级数的和函数

调用格式为：

fn为级数的通项，k为级数项数，k0与kn为级数的开始项和终止项。
（fn需要以符号表达式给出，若fn只有一个自变量，则k可以省略。）

S=symsum(fn,k,k0,kn)

实践

1. 求等比数列$a{q}^n$的和函数

syms n q a;
fn=a*q^n;
s=symsum(fn,n,0,Inf)

运行

s =
 
piecewise(1 <= q, Inf*a, abs(q) < 1, -a/(q - 1))

piecewise表示分段，即：$当q\ge 1时，s=\infty ，当|q|<1时，s=\frac{-a}{q-1}$

02 求和函数可用于检验敛散性

实践

1. $判断级数的敛散性：1、{\sum }_{n=1}^{\infty }\frac{1}{n(n+1)};2、{\sum }_{n=1}^{\infty }n(\frac{4}{3}{)}^n$

求解

%题1
syms n;
fn=1/(n*(n+1));
s=symsum(fn,n,1,Inf)

%题2
syms n;
s=symsum(n*(4/3)^n,1,Inf)

运行

题1:级数收敛
s =
 
1

题2：级数发散
s =
 
Inf

注意

有时候，对于一些复杂级数，matlab无法直接给出的解析解。例如：

1. $判断级数的敛散性：1、{\sum }_{n=1}^{\infty }(-1{)}^n\frac{1}{\sqrt{n}}$

syms n;
s=symsum((-1)^n/n^(1/2),n,1,Inf)

运行

>> syms n;
s=symsum((-1)^n/n^(1/2),n,1,Inf)
 
s =
 
symsum((-1)^n/n^(1/2), n, 1, Inf)

显然，这个答案不是我们想要的结果。此时我们可以尝试使用vpa命令给出答案的精度。

syms n;
s=symsum((-1)^n/n^(1/2),n,1,Inf);
vpa(s,8)        % 8表示s的数值解的精度

运行

>> syms n;
s=symsum((-1)^n/n^(1/2),n,1,Inf);
vpa(s,8)
 
ans =
 
-0.60489864

说明级数收敛

除了上述方法，我们还会遇到其他一些和函数无法给出解析解的情况。可考虑以下思路。

$判断{\sum }_{n=1}^{\infty }sin\frac{1}{n(n+1)}的敛散性$

syms n;
s=symsum(sin(1/(n*(n+1))),n,1,Inf);

运行

>> syms n;
s=symsum(sin(1/(n*(n+1))),n,1,Inf)
 
s =
 
symsum(sin(1/(n*(n + 1))), n, 1, Inf)

显然Matlab无法给出解析解，此时我们可以通过求有限项之和来求级数和。

syms n;
s=symsum(sin(1/(n*(n+1))),n,1,1000);     % n=1000时的和
vpa(s,6)
s=symsum(sin(1/(n*(n+1))),n,1,10000);     % n=10000时的和
vpa(s,6)
s=symsum(sin(1/(n*(n+1))),n,1,100000);     % n=100000时的和
vpa(s,6)

运行

ans =
 
0.977528        % n=1000时的和
 
ans =
 
0.978427        % n=10000时的和
 
ans =
 
0.978517        % n=100000时的和

此时从数学实验的角度分析，可以得出级数收敛于某个数。
当然，更严谨的讲，有时候需要运用比较审敛法等审敛方法来确定级数是否收敛。

03 求收敛半径与收敛域

实践

1. $求幂级数{\sum }_{n=1}^{\infty }\frac{x^n}{3^{n}n}的收敛域与和函数$

求解

A. 第一步：求收敛半径 (利用后项系数与前项系数之比，再求倒数即为收敛半径。)

syms n;
a1=1/(3^n*n);        
a2=1/(3^(n+1)*(n+1));
r=limit(a1/a2,n,Inf)

r = 3       %得出收敛半径为3

B. 第二步：考察端点处级数是否收敛

syms n;
s1=symsum((-3)^n/(3^n*n),n,1,Inf)
s2=symsum((3)^n/(3^n*n),n,1,Inf)

s1 =
 
-log(2)     % 说明级数在x=-3处收敛
 
 
s2 =
 
Inf         % 说明级数在x=3处发散

$则幂级数{\sum }_{n=1}^{\infty }\frac{x^n}{3^{n}n}收敛域为[-3,3)$

C. 第三步：求和函数

syms n x;
s1=symsum(x^n/(3^n*n),n,1,Inf)

运行结果

s1 =
 
piecewise(x == 3, Inf, abs(x) <= 3 & x ~= 3, -log(1 - x/3))     

04 泰勒展开

调用格式为：
f为要展开成级数的函数，v为自变量，a表示在x=a处将函数展开，默认为0；

Value为展开的级数的项数，默认为6，此时为n-1次泰勒式展开；

如果要改变Value的值，则’Order’不可省略。（类似于命令？键值？）

f需要以符号表达式给出，若f只有一个自变量，则v可以省略。

taylor(f,v,a,'Order',Value)

实践

1. $将y=sinx在x=0处展开成5阶泰勒展开式$

syms x;
t51=taylor(sin(x))
t52=taylor(sin(x),x)
t53=taylor(sin(x),x,0)
t54=taylor(sin(x),x,0,'Order',6)

上面4个式子的结果一致：

t5 =
 
x^5/120 - x^3/6 + x

我们还可以尝试着绘制sinx的6阶展开式的图形，直观观察：

syms x;
f=sin(x);
F=ezplot(f,[-4+pi/4,4+pi/4]);
set(F,'LineWidth',2)
hold on
t6=taylor(sin(x),x,0,'Order',7);    %注意6阶展开式，项数应该为7
T6=ezplot(t6,[-4+pi/4,4+pi/4])
set(T6,'LineStyle','--')
hold off

运行