母函数

本文介绍了母函数的一般形式及其在求解线性常系数齐次递推关系中的应用,包括普通型和指数型母函数,并通过实例展示了如何利用母函数解决具体的数学问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

母函数的一般形式

1.普通型:

G(x)=1+a1x1+a2x2+...+anxn+...G(x)=1+a1x1+a2x2+...+anxn+...

2.指数型:

G(x)=1+x11!+x22!+...+xnn!+...G(x)=1+x11!+x22!+...+xnn!+...

一些常用的Taylor展开式

1 普通型:

11x=k=0xk11−x=∑k=0∞xk

(1ax)1=k=0akxk(1−ax)−1=∑k=0∞akxk

2 指数型:
1)

ex=k=0xkk!ex=∑k=0∞xkk!

2)

ex=k=0(1)xxkk!e−x=∑k=0∞(−1)x∗xkk!

3)

12(exex)=k=0x2k+1(2k+1)!12(ex−e−x)=∑k=0∞x2k+1(2k+1)!

4)

12(ex+ex)=k=0x2k(2k)!12(ex+e−x)=∑k=0∞x2k(2k)!


例:求n位十进制数出现偶数个5的数的个数

an为出现偶数个5的数的个数,bnbn为出现奇数个5的数的个数
an=9an1+bn1an=9an−1+bn−1
bn=9bn1=an1bn=9bn−1=an−1
a1=8,b1=1a1=8,b1=1

设母函数
Ga(x)=a1+a2x+a3x2+...Ga(x)=a1+a2x+a3x2+...
Gb(x)=b1+b2x+b3x2+...Gb(x)=b1+b2x+b3x2+...

x:a2=9a1+b1x:a2=9a1+b1
x2:a3=9a2+b2x2:a3=9a2+b2
x3:a4=9a3+b3x3:a4=9a3+b3
......
求和得
Ga(x)a1=9xGa(x)+xGb(x)Ga(x)−a1=9xGa(x)+xGb(x)
整理得
(19x)Ga(x)xGb(x)=a1=8(1−9x)Ga(x)−xGb(x)=a1=8
同理可得
(19x)Gb(x)xGa(x)=b1=1(1−9x)Gb(x)−xGa(x)=b1=1
解方程可得
Ga(x)=71x+8(18x)(110x)Ga(x)=−71x+8(1−8x)(1−10x)
分解成部分分式,
A,BA,B使得A18x+B110x=Ga(x)A1−8x+B1−10x=Ga(x)成立
解得A=72,b=92A=72,b=92

Ga(x)=12(718x+9110x)∴Ga(x)=12(71−8x+91−10x)
taylor展开

(1ax)1=k=0akxk(1−ax)−1=∑k=0∞akxk

Ga(x)=12k=0(78k+910k)xkGa(x)=12∑k=0∞(7∗8k+9∗10k)xk
an=728n+9210k∴an=72∗8n+92∗10k

线性常系数齐次递推关系

递推关系:

an+c1an1+c2an2+...+ckank=0an+c1an−1+c2an−2+...+ckan−k=0

设母函数为

G(x)=a0x0+a1x1+a2x2+...+anxn+...G(x)=a0x0+a1x1+a2x2+...+anxn+...

可得

xk(ak+c1ak1+...+cka0)=0xk(ak+c1ak−1+...+cka0)=0

其中第 i 项的和为

cixi(G(x)j=0ki1aixi)cixi(G(x)−∑j=0k−i−1aixi)

总和为

i=0k(cixi(G(x)j=0ki1aixi))=0∑i=0k(cixi(G(x)−∑j=0k−i−1aixi))=0

整理得

G(x)=ki=0(cixiki1j=0aixi)ki=0cixiG(x)=∑i=0k(cixi∑j=0k−i−1aixi)∑i=0kcixi

设分子ki=0(cixiki1j=0aixi)=P(x)∑i=0k(cixi∑j=0k−i−1aixi)=P(x)

G(x)=P(x)ki=0cixiG(x)=P(x)∑i=0kcixi

单独看分母,提取xkxk,设m=x1m=x−1,得

xk(mk+c1mk1+...+ck1m1+ck)xk(mk+c1mk−1+...+ck−1m1+ck)

特征多项式(可直接通过递推关系写出):

C(x)=xk+c1xk1+...+ck1x+ckC(x)=xk+c1xk−1+...+ck−1x+ck

特征多项式可进行因式分解

C(x)=(xa1)(xa2)...(xak)C(x)=(x−a1)(x−a2)...(x−ak)

分母得

xk(1a1x1)(1a2x1)...(1akx1)xk(1−a1x−1)(1−a2x−1)...(1−akx−1)

原式得

G(x)=P(x)(1a1x)(1a2x)...(1akx)G(x)=P(x)(1−a1x)(1−a2x)...(1−akx)

分类讨论:
1)无重根

待定系数A1,A2...AkA1,A2...Ak

G(x)=i=1kAi1aixG(x)=∑i=1kAi1−aix

利用下式taylor展开可得

(1ax)1=k=0akxk(1−ax)−1=∑k=0∞akxk

an=i=1kAianian=∑i=1kAiain

A1,A2...AkA1,A2...Ak 可通过初始值列线性方程组求解

2)有重根
ββC(x)的重根,可得

G(x)=i=1kAi(1βx)iG(x)=∑i=1kAi(1−βx)i

利用二项式定理可得

(a+b)n=C0nan+C1nan1b+...+Cnnbn(a+b)n=Cn0an+Cn1abn−1+...+Cnnbn

an=i=1kAjCnn+i1βnan=∑i=1kAjCn+i−1nβn

an=(i=1kAini1)βnan=(∑i=1kAini−1)βn


例:求Sn=nk=0k2Sn=∑k=0nk2

解:
Sn=12+22+...+n2
可得
1) SnSn1=n2Sn−Sn−1=n2
2) Sn1Sn2=(n1)2Sn−1−Sn−2=(n−1)2
3) Sn2Sn3=(n2)2Sn−2−Sn−3=(n−2)2
4)Sn3Sn4=(n3)2Sn−3−Sn−4=(n−3)2

1) - 2)得
Sn2Sn1+Sn2=2n1Sn−2Sn−1+Sn−2=2n−1
2) - 3)得
Sn12Sn2+Sn3=2n3Sn−1−2Sn−2+Sn−3=2n−3

可得
Sn3Sn1+3Sn2Sn3=2Sn−3Sn−1+3Sn−2−Sn−3=2
同理可得
Sn13Sn2+3Sn3Sn4=2Sn−1−3Sn−2+3Sn−3−Sn−4=2

Sn4Sn1+6Sn24Sn3+Sn4=0∴Sn−4Sn−1+6Sn−2−4Sn−3+Sn−4=0
S0=0,S1=1,S2=5,S4=14S0=0,S1=1,S2=5,S4=14
至此,可用母函数特征多项式求出各项。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值