基于MATLAB的幂级数求和与展开（Taylor和Fourier算法）

已于 2023-12-25 15:07:34 修改
阅读量1.4w 收藏 111
点赞数 17
分类专栏： MATLAB 文章标签： matlab Taylor展开 Fourier幂级数 幂级数求和
于 2022-05-06 10:00:00 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forest_LL/article/details/124573090
版权

目录

一. Taylor幂级数展开

1.1 单变量

例题1

例题2

1.2 多变量

例题3

二. Fourier级数展开

例题4

例题5

三. 级数求和

例题6

例题7

例题8

例题9

一. Taylor幂级数展开

1.1 单变量

在x=0点附近的Taylor幂级数如下：

在x=a点展开： 

MATLAB格式：

%x=0进行Taylor幂级数展开
taylor(f,x,k)

%x=a进行Taylor幂级数展开
taylor(f,x,k,a)

例题1

求函数Taylor幂级数展开的前9项；并求关于x=2和x=a的Taylor幂级数展开。

f(x)=\frac{sinx}{x^2+4x+3}

解：

MATLAB代码：

clc;clear;
syms x a;
f=sin(x)/(x^2+4*x+3);
y1=taylor(f,x,'order',9);
pretty(y1)

taylor(f,x,2,'order',9)

taylor(f,x,a,'order',5)

运行结果：


ans =sin(2)/15 - ((131623*cos(2))/35880468750 + (875225059*sin(2))/34445250000000)*(x - 2)^8 + (x - 2)*(cos(2)/15 - (8*sin(2))/225) - ((623*cos(2))/11390625 + (585671*sin(2))/2733750000)*(x - 2)^6 + ((262453*cos(2))/19136250000 + (397361*sin(2))/5125781250)*(x - 2)^7 - (x - 2)^2*((8*cos(2))/225 + (127*sin(2))/6750) + (x - 2)^3*((23*cos(2))/6750 + (628*sin(2))/50625) + (x - 2)^4*((28*cos(2))/50625 - (15697*sin(2))/6075000) + (x - 2)^5*((203*cos(2))/6075000 + (6277*sin(2))/11390625)
 
ans =sin(a)/(a^2 + 4*a + 3) - (cos(a)/(a^2 + 4*a + 3) - (sin(a)*(2*a + 4))/(a^2 + 4*a + 3)^2)*(a - x) - (a - x)^2*(sin(a)/(2*(a^2 + 4*a + 3)) + (cos(a)*(2*a + 4))/(a^2 + 4*a + 3)^2 - (sin(a)*((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3)))/(a^2 + 4*a + 3)) - (a - x)^3*((sin(a)*((2*a + 4)/(a^2 + 4*a + 3)^2 - (((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3))*(2*a + 4))/(a^2 + 4*a + 3)))/(a^2 + 4*a + 3) - cos(a)/(6*(a^2 + 4*a + 3)) + (cos(a)*((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3)))/(a^2 + 4*a + 3) + (sin(a)*(2*a + 4))/(2*(a^2 + 4*a + 3)^2)) + (a - x)^4*(sin(a)/(24*(a^2 + 4*a + 3)) - (sin(a)*(((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3))/(a^2 + 4*a + 3) + ((2*a + 4)*((2*a + 4)/(a^2 + 4*a + 3)^2 - (((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3))*(2*a + 4))/(a^2 + 4*a + 3)))/(a^2 + 4*a + 3)))/(a^2 + 4*a + 3) + (cos(a)*((2*a + 4)/(a^2 + 4*a + 3)^2 - (((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3))*(2*a + 4))/(a^2 + 4*a + 3)))/(a^2 + 4*a + 3) + (cos(a)*(2*a + 4))/(6*(a^2 + 4*a + 3)^2) - (sin(a)*((2*a + 4)^2/(a^2 + 4*a + 3)^2 - 1/(a^2 + 4*a + 3)))/(2*(a^2 + 4*a + 3)))

例题2

对sinx进行Taylor幂级数展开，并观察不同阶次的近似效果。

解：

MATLAB代码：

clc;clear;
x0=-2*pi:0.01:2*pi;
y0=sin(x0);
syms x;
y=sin(x);
plot(x0,y0,'r-'),axis([-2*pi,2*pi,-1.5,1.5]);
hold on 
for n=[8:2:16]
    p=taylor(y,x,'order',n),
    y1=subs(p,x,x0);
    line(x0,y1);
end

运行结果：

p =- x^7/5040 + x^5/120 - x^3/6 + x
 
 
p =x^9/362880 - x^7/5040 + x^5/120 - x^3/6 + x
 
 
p =- x^11/39916800 + x^9/362880 - x^7/5040 + x^5/120 - x^3/6 + x
 
 
p =x^13/6227020800 - x^11/39916800 + x^9/362880 - x^7/5040 + x^5/120 - x^3/6 + x
 
 
p =- x^15/1307674368000 + x^13/6227020800 - x^11/39916800 + x^9/362880 - x^7/5040 + x^5/120 - x^3/6 + x

1.2 多变量

多变量函数f(x_1,x_2,\cdots,x_n)在(a_1,a_2,\cdots,a_n)的Taylor幂级数展开如下：

例题3

求z=f(x,y)的各种Taylor幂级数展开。

z=f(x,y)=(x^2-2x)e^{-x^2-y^2-xy}

解：

MATLAB代码：

clc;clear;
syms x y a;
f=(x^2-2*x)*exp(-x^2-y^2-x*y);
F1=taylor(f,[x,y],[0,0],'order',8)
F2=taylor(f,[x,y],[1,a],'order',3)

运行结果：

F1 =
x^7/3 + x^6*y + x^6/2 + 2*x^5*y^2 + x^5*y - x^5 + (7*x^4*y^3)/3 + (3*x^4*y^2)/2 - 2*x^4*y - x^4 + 2*x^3*y^4 + x^3*y^3 - 3*x^3*y^2 - x^3*y + 2*x^3 + x^2*y^5 + (x^2*y^4)/2 - 2*x^2*y^3 - x^2*y^2 + 2*x^2*y + x^2 + (x*y^6)/3 - x*y^4 + 2*x*y^2 - 2*x
 
 
F2 =
(exp(- a^2 - a - 1) - exp(- a^2 - a - 1)*((a/2 + 1)*(a + 2) - 1))*(x - 1)^2 - exp(- a^2 - a - 1) - exp(- a^2 - a - 1)*(a - y)^2*((2*a + 1)*(a + 1/2) - 1) + exp(- a^2 - a - 1)*(a + 2)*(x - 1) - exp(- a^2 - a - 1)*(2*a + 1)*(a - y) + exp(- a^2 - a - 1)*(a - y)*(x - 1)*((2*a + 1)*(a/2 + 1) + (a + 2)*(a + 1/2) - 1)

二. Fourier级数展开

给定f(x)以及如下条件：

x\in[-L,L],T=2L,f(t)=f(kT+t)

Fourier级数，如下：

 编写Fourier级数的解析函数：

function [A,B,F]=fseries(f,x,n,a,b)
if nargin==3,a=-pi;
    b=pi;
end
L=(b-a)/2;
if a+b, f=subs(f,x,x-L-a);
end %变量区域互换
A=int(f,x,-L,L)/L;
B=[];
F=A/2; %计算a0
for i=1:n
    an=int(f*cos(i*pi*x/L),x,-L,L)/L;
    bn=int(f*cos(i*pi*x/L),x,-L,L)/L;
    A=[A,an];
    B=[B,bn];
    F=F+an*cos(i*pi*x/L)+bn*sin(i*pi*x/L);
end
if a+b,
    F=subs(F,x,x+L+a);
end %换回变量区域
end

例题4

求y的Fourier级数展开。

y=x(x-\pi)(x-2\pi),\quad x\in(0,2\pi)

解：

MATLAB代码：

clc;clear;
syms x;
f=x*(x-pi)*(x-2*pi);
[A,B,F]=fseries(f,x,6,0,2*pi)


function [A,B,F]=fseries(f,x,n,a,b)
if nargin==3,a=-pi;
    b=pi;
end
L=(b-a)/2;
if a+b, f=subs(f,x,x-L-a);
end %变量区域互换
A=int(f,x,-L,L)/L;
B=[];
F=A/2; %计算a0
for i=1:n
    an=int(f*cos(i*pi*x/L),x,-L,L)/L;
    bn=int(f*sin(i*pi*x/L),x,-L,L)/L;
    A=[A,an];
    B=[B,bn];
    F=F+an*cos(i*pi*x/L)+bn*sin(i*pi*x/L);
end
if a+b,
    F=subs(F,x,x+L+a);
end %换回变量区域
end

运行结果：

A =[ -16*pi^3, 24*pi, -6*pi, (8*pi)/3, -(3*pi)/2, (24*pi)/25, -(2*pi)/3]
 
 
B =
 [ -(12*pi - 24*pi^3)/pi, ((3*pi)/2 - 12*pi^3)/pi, -((4*pi)/9 - 8*pi^3)/pi, ((3*pi)/16 - 6*pi^3)/pi, (24*pi^2)/5 - 12/125, (pi/18 - 4*pi^3)/pi]
 
 
F =
(sin(x)*(12*pi - 24*pi^3))/pi - 24*pi*cos(x) - 8*pi^3 - 6*pi*cos(2*x) - (8*pi*cos(3*x))/3 - (3*pi*cos(4*x))/2 - (24*pi*cos(5*x))/25 - (2*pi*cos(6*x))/3 - sin(5*x)*((24*pi^2)/5 - 12/125) + (sin(2*x)*((3*pi)/2 - 12*pi^3))/pi + (sin(3*x)*((4*pi)/9 - 8*pi^3))/pi + (sin(4*x)*((3*pi)/16 - 6*pi^3))/pi + (sin(6*x)*(pi/18 - 4*pi^3))/pi

例题5

对(-\pi,\pi)区间内的方波信号，当x\geq 0时，y=1，否则y=-1。对它进行Fourier级数拟合，并观察用多少项能有较好的拟合效果。

解：

MATLAB代码：

clc;clear;
syms x;
f=abs(x)/x; %定义方波信号
xx=[-pi:pi/200:pi];
xx=xx(xx~=0);
xx=sort([xx,-eps,eps]); %剔除零点
yy=subs(f,x,xx);
plot(xx,yy,'r-.'), %绘制出理论值
hold on  %保持坐标系

for n=2:20
    [a,b,f1]=fseries(f,x,n);
    y1=subs(f1,x,xx);
    plot(xx,y1)
end


function [A,B,F]=fseries(f,x,n,a,b)
if nargin==3,a=-pi;
    b=pi;
end
L=(b-a)/2;
if a+b, f=subs(f,x,x-L-a);
end %变量区域互换
A=int(f,x,-L,L)/L;
B=[];
F=A/2; %计算a0
for i=1:n
    an=int(f*cos(i*pi*x/L),x,-L,L)/L;
    bn=int(f*sin(i*pi*x/L),x,-L,L)/L;
    A=[A,an];
    B=[B,bn];
    F=F+an*cos(i*pi*x/L)+bn*sin(i*pi*x/L);
end
if a+b,
    F=subs(F,x,x+L+a);
end %换回变量区域
end

运行结果：

三. 级数求和

给定：

S=\sum_{k=k_0}^{k_n}f_k

符号工具箱中有提供：

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

例题6

计算求和：

S=2^0+2^1+2^2+2^3+2^4+\cdots+2^{62}+2^{63}=\sum_{i=0}^{63}2^i

解：

MATLAB代码：

clc;clear;
format long;

%数值计算
sum(2.^[0:63])

%把2定义为符号量，计算会更准确
sum(sym(2).^[0:63])

%工具箱
syms k;
symsum(2^k,0,63)

运行结果：

ans =1.844674407370955e+19

 
ans =18446744073709551615
 
 
ans =18446744073709551615

例题7

试求无穷级数的和

S=\frac{1}{1\times 4}+\frac{1}{4\times 7}+\frac{1}{7\times 10}+\cdots+\frac{1}{(3n-2)\times (3n+1)}+\cdots

解：

MATLAB代码：

clc;clear;

%符号工具箱
syms n;
s=symsum(1/((3*n-2)*(3*n+1)),n,1,inf)

%数值计算
format long %以长型方式显示得出的结果
m=1:1000000;
s1=sum(1./((3*m-2).*(3*m+1))) %双精度有效位16

运行结果：

s =1/3
 

s1 =0.333333222222265

例题8

求解：

J=2\sum_{n=0}^\infty\frac{1}{(2n+1)(2x+1)^{2n+1}}

解：

MATLAB代码：

clc;clear;
syms n x;
s1=symsum(2/((2*n+1)*(2*x+1)^(2*n+1)),n,0,inf)

运行结果：

 
s1 =piecewise(1 < abs(2*x + 1), 2*atanh(1/(2*x + 1)))

例题9

求：

\lim_{n\to\infty}[(1+\frac{1}{2}+\frac{1}{3}+\frac{1}{4}+\cdots++\frac{1}{n})-lnn]

解：

MATLAB代码：

clc;clear;
syms m n;
limit(symsum(1/m,m,1,n)-log(n),n,inf);
vpa(ans,70)  %显示70位有效数字

运行结果：

ans =
 
0.5772156649015328606065120900824024310421593359399235988057672348848677

【MATLAB】求级数的和
Fanjufei的博客
02-06 2884
syms ksymsum(1/k^2,1,Inf) %k值为1到无穷大在 MATLAB 中，我们可以通过编写程序来进行级数求和的实验。下面我将详细介绍一下 MATLAB 中进行级数求和实验的流程，并提供一份总结。级数求和实验的主要目的是让用户通过实际计算和观察，深入理解级数求和的方法和原理，提高对级数性质的认识。知识点级数的定义符号级数求和示例演示级数求和是数学中一个重要的概念，它涉及到了数列和无限求和的问题。下面我将详细解释级数求和的方法。首先，我们需要了解级数的定义。
matlab从无到有系列（四）：符号数学基础
左手の明天的博客
03-01 2588
matlab从无到有系列文章汇总： matlab从无到有系列（一）：认识matlab matlab从无到有系列（二）：矩阵运算基础 matlab从无到有系列（三）：数值计算基础
MATLAB泰勒展开
奔流聚海
06-12 3万+
MATLAB函数：taylor 题目：对y=exp(-x)进行4阶泰勒展开，并验证： 函数调用格式： taylor(fcn,x,x0,'Order',6);%对函数fcn在点x0处，进行6阶泰勒展开； MATLAB代码： clc clear all close all syms x y=exp(-x);
matlab计算级数和
BATP2516的博客
07-06 1万+
级数(series)是指将数列的项依次用加号连接起来的函数。典型的级数有正项级数、交错级数、幂级数、傅里叶级数等。 —百度百科 matlab提供了直接计算级数的两个方法，symsum的两种形式（做题刚开始以为是重载，直接导致了上一篇帖子的产生）； 使用help symsum，查看文档，matlab有如图两种求级数的方法： ①F = symsum(f,k,a,b) 官方解释： 简单来说：函数返回级数f关于求和索引k从a到b的和，如果不指定索引k，matlab会调用symvar函数检索第一个值作为索引；特
matlab集训七（微分方程，级数求和）
sunshine博客
05-14 2793
求解微分方程（组） 1.解析解 在 MATLAB 中用 D 表示导数，例如 Dy表示 y’，D2y表示 y’’，Dy(0)=5 表示 y’(0)=5 等，符号常微分方程求解通过函数 dsolve 实现 调用格式为： S=dsolve(eqn) S=dsolve(eqn, cond) S=dsolve(eqn, cond,Name,Value） 例1： syms u(t) eqn=diff(u, ...
利用MATLAB计算级数之和与泰勒级数的应用
Nick的博客
08-19 1万+
文章目录摘要1 级数求和1.1 MATLAB级数求和函数1.2 级数求和应用实例1.3 特例2 泰勒级数2.1 泰勒级数函数的调用2.2 泰勒级数的应用结语 摘要 本文是《科学计算与MATLAB语言》专题七第3小节的学习笔记，如果大家有时间的话，可以去听听课，没有的话，可以看看我的笔记，还是很不错的。本节主要讲了如何利用MATLAB求级数之和以及泰勒级数如何计算，文中每个代码我都跑过一遍，可以直接复制到MATLAB中运行。 1 级数求和 1.1 MATLAB级数求和函数 求无穷级数的和需要符号表达式求和函数
MATLAB 之 级数计算和方程符号求解
最新发布
每天进步一点点！
06-20 4436
【代码】MATLAB 之 级数计算和方程符号求解。
Matlab学习笔记（三）
Enderman_xiaohei的博客
04-02 1920
主要内容（大概） 一、复习积分、级数 二、泰勒级数和傅里叶级数 三、Matlab求解 —— 微分和代数 积分复习 1.quad做无穷（广义）有点问题，int比较精确，用int最好，下面的例子用来对比： 例，计算下列广义积分 ∫−∞∞ecosx−x2dx\int_{-\infty }^{\infty}e^{cosx-x^{2}}dx∫−∞∞​ecosx−x2dx clear format long ...
MATLAB在高等数学教学中的几种应用.pdf
07-03
符号工具箱提供了符号变量的声明方式，可以解决函数的极限、单变量与多变量函数微积分、常微分方程、级数求和、Taylor幂级数展开、Fourier级数展开等高等数学的重要内容。例如，在MATLAB中通过声明符号变量，使用...
[渝粤教育] 东北大学 现代科学运算—MATLAB语言与应用 参考 资料
渝粤题库
02-01 195
教育 -现代科学运算—MATLAB语言与应用-章节资料考试资料-东北大学【】 01-01 本课程的主要内容 1、【判断题】 A、正确 B、错误 参考资料【 】 2、【判断题】 A、正确 B、错误 参考资料【 】 01-02 为什么学习计算机数学语言 1、【判断题】 A、正确 B、错误 参考资料【 】 2、【判断题】 A、正确 B、错误 参考资料【 】 01-03 解析解与数值解 1、【单选题】 A、 B、 C、 D、 参考资料【 】 2、【判断题】 A、正确 B、错误
幂级数的求和方法（杜道渊 caj文件）
11-04
幂级数的算法汇总 常用方法提供 有需要的同学可以参考
【MATLAB符号计算揭秘】：数学表达式处理与符号运算的7大妙招
[【MATLAB符号计算揭秘】：数学表达式处理与符号运算的7大妙招](https://cdn.educba.com/academy/wp-content/uploads/2020/04/Simplify-Matlab.jpg) # 1. MATLAB符号计算概述 MATLAB符号计算是该平台的一个重要组成...
Matlab小笔记随手记——求级数和
qq_43216397的博客
12-19 1160
在MATLAB中，用于级数求和的命令是symsum（）， 命令的应用格式为： symsum（f，v，a，b） 即求通式f在指定自变量v取遍[a,b]中所有整数时的和 其中： f为级数的通项表达式，v是通项中的求和变量，a和b分别为求和变量的起点和终点。 如果a,b缺省，则v从0变到v-1， 如果a,b,v缺省，则系统对f中的默认变量求和 ...
matlab编程级数的和函数,数学实验8 用Matlab软件求级数的和函数的泰勒级数和傅氏级数.pdf...
weixin_36303305的博客
03-16 1467
数学实验8 用Matlab软件求级数的和函数的泰勒级数和傅氏级数数学实验八 用Matlab 软件求级 的和、函 的泰勒级 和傅氏级一、求级 的和Matlab 中，可用symsum 函数求数列或级数的和，其调用格式为：syms nsymsum(f(n), n a ,b)其中，f (n) 为数列或级数的通项，n 为自变量，a 为该数列或级数所求和的起始项数，b 为该数列或级数所求和的结束项数．此格式...
【Matlab】求解级数展开和求和
m0_64548999的博客
05-27 3301
求解级数展开和求和
MATLAB----级数
sqyao_147的博客
08-26 5921
本篇参考于：中国大学慕课，专题六，“7.3级数" 1.级数求和 2.泰勒级数 sum也可以进行级数求和，但是项数必须有限，并且要表示成向量形式，如果级数是无限的且项数非常多，用sum无法解决 1.级数求和 调用格式：symsum(s,v,n,m) 其中，s表示一个级数的通项，是一个符号表达式。v是求和变量，v省略时使用系统的默认变量。n和m是求和变量v的初值和末值。 在级数计算中，因为小数都表示有理分数的形式，容易导致分子或分母出现极大整数从而无法计算的情况，出现了具体的级数，..
【MATLAB求和】
热门推荐
qq_52624396的博客
05-10 10万+
提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档 文章目录MATLAB怎么求和∑一、题目二、代码1.syms函数2.直接用循环补充 MATLAB怎么求和∑ 提示：这里可以添加本文要记录的大概内容： 学习一下用matlab进行简单的函数求和。 提示：以下是本篇文章正文内容，下面案例可供参考 一、题目 二、代码 1.syms函数 代码如下（示例）： syms k symsum((k.^2+7*k).^3,1,100) 2.直接用循环 代码如下（示例）： s=0; for i=1:1.
matlab将函数展开成幂级数,解析函数展开成幂级数的方法分析.doc
weixin_30988079的博客
03-19 2163
解析函数展开成幂级数的方法分析解析函数展开成幂级数的方法分析姓 名：学 号：20专 业：物理指导教师：解析函数展开成幂级数的方法分析摘要关键词一 前 言解析函数边值问题和广义解析函数边值问题在奇异积分方程方面有广泛的应用，它们在弹性力学、流体力学方面也有重要的应用。这些方面的理论及其应用，主要是由苏联学者建立和发展起来的。自20世纪60年代以来，中国的数学工作者在这些方面也做...
Matlab中求和
小陈同学的博客
06-19 2万+
Matlab中求和 一、sum函数 1、如果x是一个向量，则计算向量的元素之和，S=sum（A）如： >> x = randi(10,1,10) x = 7 1 9 10 7 8 8 4 7 2 >> sum(x) ans = 63 %如果x的数组，则是计算数组的列之和，如： >> x = pascal(4,2) x = -1 -1 -1 -1 3
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

唠嗑！

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值