数学实验教程实验(级数)

合集下载

最新实验1无穷级数(基础实验)

实验1无穷级数(基础实验)项目四无穷级数与微分方程实验1 无穷级数（基础实验）实验目的观察无穷级数部分和的变化趋势,进一步理解级数的审敛法以及幂级数部分和对函数的逼近. 掌握用Mathematica求无穷级数的和, 求幂级数的收敛域, 展开函数为幂级数以及展开周期函数为傅里叶级数的方法.基本命令1. 求无穷和的命令Sum该命令可用来求无穷和. 例如,输入Sum[1/n^2,{n,l,Infinity}]则输出无穷级数的和为«Skip Record If...»命令Sum与数学中的求和号«Skip Record If...»相当.2. 将函数展开为幂级数的命令Series该命令的基本格式为Series[f[x],{x,x0,n}]它将«Skip Record If...»展开成关于«Skip Record If...»的幂级数. 幂级数的最高次幂为«Skip Record If...»余项用«Skip Record If...»表示. 例如,输入Series[y[x],{x,0,5}]则输出带皮亚诺余项的麦克劳林级数«Skip Record If...»3. 去掉余项的命令Normal在将«Skip Record If...»展开成幂级数后, 有时为了近似计算或作图, 需要把余项去掉. 只要使用Normal命令. 例如,输入Series[Exp[x],{x,0,6}]Normal[%]则输出«Skip Record If...»«Skip Record If...»4. 强制求值的命令Evaluate如果函数是用Normal命令定义的, 则当对它进行作图或数值计算时, 可能会出现问题. 例如,输入fx=Normal[Series[Exp[x],{x,0,3}]]Plot[fx,{x,-3,3}]则只能输出去掉余项后的展开式«Skip Record If...»而得不到函数的图形. 这时要使用强制求值命令Evaluate, 改成输入Plot[Evaluate[fx],{x,-3,3}]则输出上述函数的图形.5. 作散点图的命令ListPlotListPlot [ ]为平面内作散点图的命令, 其对象是数集,例如,输入仅供学习与交流，如有侵权请联系网站删除谢谢127仅供学习与交流，如有侵权请联系网站删除谢谢127ListPlot[Table[j^2,{j,16}],PlotStyle->PointSize[0,012]]则输出坐标为«Skip Record If...»的散点图(图1.1).图1.16. 符号“/；”用于定义某种规则,“/；”后面是条件. 例如,输入Clear[g,gf];g[x_]:=x/;0<=x<1g[x_]:=-x/;-1<=x<0g[x_]:=g[x –2]/;x>=1则得到分段的周期函数«Skip Record If...»再输入gf=Plot[g[x],{x,-1,6}]则输出函数«Skip Record If...»的图形1.2.图1.2注:用Which 命令也可以定义分段函数, 从这个例子中看到用“…(表达式)/; …(条件)”来定义周期性分段函数更方便些. 用Plot 命令可以作出分段函数的图形, 但用Mathematica 命令求分段函数的导数或积分时往往会有问题. 用Which 定义的分段函数可以求导但不能积分. Mathematica 内部函数中有一些也是分段函数. 如:Mod[x,1],Abs[x],Floor[x]和UnitStep[x]. 其中只有单位阶跃函数UnitStep[x]可以用Mathematica 命令来求导和求定积分. 因此在求分段函数的傅里叶系数时, 对分段函数的积分往往要分区来积. 在被积函数可以用单位阶跃函数UnitStep的四则运算和复合运算表达时, 计算傅里叶系数就比较方便了.实验举例数项级数例1.1 (教材例1.1)(1)观察级数«Skip Record If...»的部分和序列的变化趋势.(2) 观察级数«Skip Record If...»的部分和序列的变化趋势.输入s[n_]=Sum[1/k^2,{k,n}];data=Table[s[n],{n,100}];ListPlot[data];N[Sum[1/k^2,{k,Infinity}]]N[Sum[1/k^2,{k,Infinity}],40]则输出(1)级数的近似值为1.64493.输入s[n_]=Sum[1/k,{k,n}];data=Table[s[n],{n,50}];ListPlot[data,PlotStyle->PointSize[0.02]];则输出(2)例1.2 (教材例1.2) 画出级数«Skip Record If...»的部分和分布图.输入命令Clear[sn,g];sn=0;n=1;g={};m=3;While[1/n>10^-m,sn=sn+(-1)^(n-1)/n;g=Append[g,Graphics[{RGBColor[Abs[Sin[n]],0,1/n],仅供学习与交流，如有侵权请联系网站删除谢谢127Line[{{sn,0},{sn,1}}]}]];n++];Show[g,PlotRange->{-0.2,1.3},Axes->True];则输出所给级数部分和的图形（图1.5）,从图中可观察到它收敛于0.693附近的一个数.例1.3 求«Skip Record If...»的值.输入Sum[x^(3k),{k,1,Infinity}]得到和函数«Skip Record If...»例1.4 (教材例1.3)设«Skip Record If...»求«Skip Record If...».输入Clear[a];a[n_]=10^n/(n!);vals=Table[a[n],{n,1,25}];ListPlot[vals,PlotStyle->PointSize[0.012]]则输出«Skip Record If...»的散点图（1.6）,从图中可观察«Skip Record If...»的变化趋势. 输入 Sum[a[n],{n,l,Infinity}]图1.6求幂级数的收敛域例1.5 (教材例1.4)求«Skip Record If...»的收敛域与和函数.仅供学习与交流，如有侵权请联系网站删除谢谢127输入Clear[a];a[n_]=4^(2n)*(x-3)^n/(n+1);stepone=a[n+1]/a[n]//Simplify则输出«Skip Record If...»再输入steptwo=Limit[stepone,n->Infinity]则输出«Skip Record If...»这里对a[n+1]和a[n]都没有加绝对值. 因此上式的绝对值小于1时, 幂级数收敛; 大于1 时发散. 为了求出收敛区间的端点, 输入ydd=Solve[steptwo==1,x]zdd=Solve[steptwo==-1,x]则输出«Skip Record If...»由此可知,当«Skip Record If...»时,级数收敛,当«Skip Record If...»或«Skip Record If...»时,级数发散.为了判断端点的敛散性, 输入Simplify[a[n]/.x->(49/16)]则输出右端点处幂级数的一般项为«Skip Record If...»因此,在端点«Skip Record If...»处,级数发散. 再输入Simplify[a[n]/.x->(47/16)]则输出左端点处幂级数的一般项为«Skip Record If...»因此,在端点«Skip Record If...»处, 级数收敛.也可以在收敛域内求得这个级数的和函数. 输入Sum[4^(2n)*(x-3)^n/(n+1),{n,0,Infinity}]则输出«Skip Record If...»函数的幂级数展开例1.6 (教材例1.5)求«Skip Record If...»的6阶麦克劳林展开式.输入Series[Cos[x],{x,0,6}]则输出«Skip Record If...»注:这是带皮亚诺余项的麦克劳林展开式.例1.6 (教材例1.6)求«Skip Record If...»在«Skip Record If...»处的6阶泰勒展开式.输入Series[Log[x],{x,1,6}]则输出«Skip Record If...»例1.7 (教材例1.7) 求«Skip Record If...»的5阶泰勒展开式.输入仅供学习与交流，如有侵权请联系网站删除谢谢127serl=Series[ArcTan[x],{x,0,5}];Poly=Normal[serl]则输出«Skip Record If...»的近似多项式«Skip Record If...»通过作图把«Skip Record If...»和它的近似多项式进行比较. 输入Plot[Evaluate[{ArcTan[x],Poly}],{x,-3/2,3/2},PlotStyle->{Dashing[{0.01}],GrayLevel[0]},AspectRatio->l]则输出所作图形（图1.7.例1.9 求«Skip Record If...»在«Skip Record If...»处的8阶泰勒展开, 并通过作图比较函数和它的近似多项式.输入Clear[f];f[x_]=Exp[-(x-1)^2*(x+1)^2];poly2=Normal[Series[f[x],{x,1,8}]]Plot[Evaluate[{f[x],poly2}],{x,-1.75,1.75},PlotRange->{-2,3/2},PlotStyle->{Dashing[{0.01}],GrayLevel[0]}]则得到近似多项式和它们的图1.8.«Skip Record If...»«Skip Record If...»仅供学习与交流，如有侵权请联系网站删除谢谢127图1.8例1.10 求函数«Skip Record If...»在«Skip Record If...»处的«Skip Record If...»阶泰勒展开, 通过作图比较函数和它的近似多项式, 并形成动画进一步观察.因为«Skip Record If...»所以输入Do[Plot[{Sum[(-1)^j*x^(2j+1)/(2j+1)!,{j,0,k}],Sin[x]},{x,-40,40},PlotStyle->{RGBColor[1,0,0],RGBColor[0,0,1]}],{k,1,45}]则输出为«Skip Record If...»的3阶和91阶泰勒展开的图形. 选中其中一幅图形,双击后形成动画. 图1.9是最后一幅图.例1.11 利用幂级数展开式计算«Skip Record If...»(精确到«Skip Record If...»).因为«Skip Record If...»根据«Skip Record If...»在«Skip Record If...»处的展开式有«Skip Record If...»故前«Skip Record If...»项部分和为«Skip Record If...»仅供学习与交流，如有侵权请联系网站删除谢谢127输入命令s[n_]=3(1-1/(5*3^4)-Sum[Product[5i-1,{i,1,k-1}]/(5^k k!3^(4k)),{k,2,n-1}]);r[n_]=Product[5i-1,{i,1,n-1}]/5^n/n!3^(4n-5)/80;delta=10^(-10);n0=100;Do[Print["n=",n,",","s[n]=",N[s[n],20]];If[r[n]<delta,Break[]];If[n==n0,Print["failed"]],{n,n0}]则输出结果为«Skip Record If...»傅里叶级数例1.12 (教材例1.8) 设«Skip Record If...»是以«Skip Record If...»为周期的周期函数,它在«Skip Record If...»的表达式是«Skip Record If...»将«Skip Record If...»展开成傅里叶级数.输入Clear[g];g[x_]:=-1/;-Pi<=x<0g[x_]:=1/;0<=x<Pig[x_]:=g[x-2Pi]/;Pi<=xPlot[g[x],{x,-Pi,5 Pi},PlotStyle->{RGBColor[0,1,0]}];则输出«Skip Record If...»的图形 (图1.10).因为«Skip Record If...»是奇函数, 所以它的傅里叶展开式中只含正弦项. 输入b2[n_]:=b2[n]=2 Integrate[1*Sin[n*x],{x,0,Pi}]/Pi;fourier2[n_,x_]:=Sum[b2[k]*Sin[k*x],{k,1,n}];tu[n_]:=Plot[{g[x],Evaluate[fourier2[n,x]]}, {x,-Pi,5 Pi},PlotStyle->{RGBColor[0,1,0],RGBColor[1,0.3,0.5]},DisplayFunction->Identity];(*tu[n]是以n为参数的作图命令*)tu2=Table[tu[n],{n,1,30,5}];(*tu2是用Table命令作出的6个图形的集合*)toshow=Partition[tu2,2];(*Partition是对集合tu2作分割, 2为分割的参数*)Show[GraphicsArray[toshow]](*GraphicsArray是把图形排列的命令*)则输出6个排列着的图形（图1.11）,每两个图形排成一行. 可以看到«Skip Record If...»越大, «Skip Record If...»的傅里叶级数的前«Skip Record If...»项和与«Skip Record If...»越接近.仅供学习与交流，如有侵权请联系网站删除谢谢127实验习题1.求下列级数的和:(1)«Skip Record If...» (2)«Skip Record If...» (3)«Skip Record If...» (4)«SkipRecord If...»2. 求幂级数«Skip Record If...»的收敛域与和函数.3. 求函数«Skip Record If...»的6阶麦克劳林多项式.4. 求«Skip Record If...»的6阶麦克劳林多项式.5. 设«Skip Record If...»,求«Skip Record If...»的5阶和10阶麦克劳林多项式,把两个近似多项式和函数的图形作在一个坐标系内.6. 设«Skip Record If...»在一个周期内的表达式为«Skip Record If...», 将它展开为傅里叶级数(取6项), 并作图.7. 设«Skip Record If...»在一个周期内的表达式为«Skip Record If...», 将它展开为傅里叶级数(取8项), 并作图.8. 求级数«Skip Record If...»的和的近似值.仅供学习与交流，如有侵权请联系网站删除谢谢127。

数学实验-数列与级数

数学实验-数列与级数实验四数列与级数⼀．实验名称: 数列与级数。

⼆．实验⽬的: 学习使⽤Mathematica 发现数列与级数的极限与规律以及极限状态的性质,通过作出Fibonacci 数列的折线图，考察Fibonacci 数列的极限与规律。

三．实验环境: Mathematica 系统，Word ⽂档，课本。

四．实验的基本理论和⽅法：⽤计算机计算出Fibonacci 数列每⼀项的值，并在⼆维平⾯上画出顺次连接点()n F n ,，观察其单调性以及增加速度；猜测通项公式满⾜n n cr F =，并进⾏尝试，带⼊差分⽅程，从中解出特征根。

3n+1问题只需对奇数进⾏分析，如果对每个n ，数列中有某⼀项⼩于n ，从n 开始产⽣的数列最后都落于124→→中。

五．实验内容和步骤：1. Fabonacci 数列的极限规律（1）递推关系式为n n n F F F +=++12 ， .,2,1 =n .2,121==F F的数列被称为Fibonacci 数列。

（2）画Fibonacci 数列折线图输⼊程序：FibShow n _Integer :Modulet ,i ,For i 1,i n,i ,AppendTo t , i ,Fibonacci i ListPlot t ,PlotJoined True FibShow 20运⾏结果得：510151000200030004000（3）⽤直线去拟合（i,)log(i F ），i=1，2，….的函数。

输⼊程序：FibFit n _Integer :Modulet ,i ,For i 1,i n,i ,AppendTo t , i ,Log Fibonacci i Fit t , 1,x ,x FibFit 1000运⾏结果得：0.8039010.481（4）演奏Fibonacci 数列的函数输⼊程序：FibPlay n _Integer :Modulet ,i ,For i 1,i n,i ,AppendTo t ,Mod Fibonacci i ,n ;ListPlay t ,PlayRange 0,n ,SampleRate FibPlay 50运⾏结果得：（5）显⽰点列()()i i sin ,，n i ,,2,1 =的函数输⼊程序：PlotList n _Integer :Modulet ,i ,For i 1,i n,i ,AppendTo t , i ,Sin i ListPlot t ,PlotStyle P ointSize 0.005 PlotList 100运⾏结果得： 20406080-1-0.50.512. 3n+1问题3n+1问题起源于20世纪50年代，⼜称为Syracuse 猜想，⾓⾕猜想，Collatz 问题，Hasse 算法问题，Ulamw 问题，Thwaites 猜想等等，⽬前有⼈验证50237.1?≤n ，猜想仍然成⽴。

实验八、级数及运算

实验八级数及运算【实验类型】验证性【实验学时】1学时【实验目的】1．掌握用MA TLAB 判定常数项级数的敛散性的方法。

2．掌握用MA TLAB 进行幂级数求和的方法。

3．掌握用MA TLAB 将函数展开成幂级数的方法；【实验内容】1．熟悉有关级数收敛、发散的判定方法和级数求和； 2．利用MATLAB 判断常数项级数的敛散性； 3．熟悉有关幂级数的各种运算；4．利用MATLAB 进行幂级数的求和运算； 5．利用MATLAB 进行函数的幂级数展开；【实验方法与步骤】一、实验的基本理论与方法 1、常数项级数的审敛法：（1）级数收敛的必要条件：若级数∑∞=1n nu收敛，则必有0lim =∞→n n u 。

（2）比较审敛法的极限形式：设有正项级数∑∞=1n n u ，∑∞=1n n v ，若ρ=∞→nnn v u lim，)0(+∞<<ρ，则级数∑∞=1n n u ，∑∞=1n n v 同时收敛或发散。

（3）比值审敛法：设有正项级数∑∞=1n n u ，若ρ=+∞→nn n u u 1lim，当1<ρ时级数收敛；当1>ρ时级数发散。

（4）条件收敛与绝对收敛：若级数∑∞=1n nu收敛（级数∑∞=1n nu必收敛），则称级数绝对收敛；若级数∑∞=1n nu发散，而级数∑∞=1n nu收敛，则称级数条件收敛。

2、幂级数展开的唯一性：若函数)(x f 在含点0x 的某一区间内能展开为幂级数，则必为Taylor 级数+-++-''+-'+=n n x x x f n x x x f x x x f x f x f ))((!1))((!21))(()()(00)(200000 二、实验使用的MATLAB 函数1、 taylor(f,x,k) 将)(x f 按x=0进行Taylor 幂级数展开taylor(f,x,k,a) 将)(x f 按x=a 进行Taylor 幂级数展开其中，f 为函数的符号表达式，x 为自变量，若函数只有一个自变量，则x 可以省略。

实验8 无穷级数与数值计算

西安理工大学应用数学系
（2）相应的Matlab代码为 syms x taylor((1+x)^0.5, x, 7, 1)
运行结果为 ans =2^(1/2) + (2^(1/2)*(x - 1))/4 - (2^(1/2)*(x - 1)^2)/32 +
(2^(1/2)*(x - 1)^3)/128 - (5*2^(1/2)*(x - 1)^4)/2048 + (7*2^(1/2)*(x - 1)^5)/8192 - (21*2^(1/2)*(x 1)^6)/65536
(n
3
1)!

105,取
n8
即可。
西安理工大学应用数学系
同理取 x 1
2
,令误差
2n
1
3 (n

1)!

105
,取 n 6 即可。
计算e的Matlab代码为
syms x y y=taylor(exp(x),x,9)
a=subs(y,x,'1') %将符号表达式y中的变量x换成变量1 e=vpa(a,10) % 将符号表达式a转化成10位精度的近似解
1 x 1 1 x 1 x2 1 3 x3 (1)n1 1 35 (2n 3) xn （6.6）
2 24 246 取
2 46 2n
于是有
, 代入(6.6), 有
1 x 1 1 x 1 x2 13 x3 (1)n1 1 3 5 (2n 3) xn
西安理工大学应用数学系
hold on for i=1:length(x)
y(i)=x(i)-x(i)^3/6+x(i)^5/120; plot(x(i),y(i),'.'),pause(.005) end hold on for i=1:length(x) y(i)=x(i)-x(i)^3/6+x(i)^5/120-

数学实验一

数学实验报告册姓名：马会兰学号：200771010423班级：07级数4实验一：（微积分基础）一．实验目的：学会使用Mathematica 的一些基本功能，验证或观察得出微积分的几个基本结论。

二．实验环境：在Mathematica 环境下结合教材进行实验。

三．实验的基本理论和方法：Mathematica 能够进行初等数学和高等数学的数值计算、符号计算、画图等各种事情。

四．实验的内容和步骤：练习1：泰勒（Taylor ）级数⑴在同一坐标系里作出函数36x y x =-及其导数'sin y x =，0.8y x =，y x =与1.2y x =的图像。

Mathematica 语句如下：321321 （图1-1）结果分析：从上图中可以发现，在具有不同斜率k 的过原点的直线y kx =中，k=1时的直线y x =与正弦曲线sin y x =在原点附近最接近，如上图所示。

观察发现：从原点出发沿直线y x =前进与沿正弦曲线sin y x =前进的方向时一致的，在原点的附近的很小一段旅程中两条路线几乎一样，但继续下去，就分开了，因此能不能用越来越高次的多项式函数去逼近sin y x =呢？请看下面。

⑵在同一坐标系里作出区间[,]x ππ∈-上正弦函数s i n y x =及多项式函数36x y x =-，356120x x y x =-+，3573!5!7!x x x y x =-+-的图像。

3211.00.5Mathematica 语句如下：运行的结果：n a ，n A 的值为：结果分析：可以看出n a 的值与n A 的值越来越接近，最后而这达到相等的地步。

⑵在同一坐标系中画出下面三个函数的图象：101(1)10x x y =+，1011(1)10x x y +=+，y e = 观察当x 增大时图像的走向。

Ⅰ.函数在区间[1，4]内的图象 Mathematica 语句如下：图像如下：（图2-1）Ⅱ. 函数在区间[3，5]内的图象Mathematica 语句如下：图像如下：（图2-2）Ⅲ. 函数在区间[5，6]内的图象 Mathematica 语句如下：图像如下：（图2-3）结果分析：通过观察可以看出，当n 增大时1(1)n n an =+递增，11(1)n n A n+=+递减。

数学物理实验第三节(泰勒级数展开)

复变函数的泰勒级数和实变函数的运算法则一样，但要注意复数运算和实数运算的异同，在计算的时候，考虑全面！
13
展开公式
2 3 n z z z ez 1 z 2! 3! n! z3 z5 z 2 n1 n sin z z 1 3! 5! 2n 1!
m m
可求得收敛半径为1，由此可得
m m
m(m 1) 2 m(m 1)(m 2) 3 m (1 z ) 1 1 z z z ... 2! 3! 1! z 1
9
m(m 1) 2 m(m 1)(m 2) 3 m (1 z ) 1 1 z z z ... 2! 3! 1! z 1
由 R lim | ak | k a k 1
可知泰勒级数的收敛半径为无限大，只要z
是有限的，则泰勒级数就是收敛的！
例2
在z0＝0的邻域上把 f1 ( z) sin z, f 2 ( z) cos z 展开
( z) sin z 解： f1 ( z) sin z 的前四阶导数是 f1( z) cos z, f1
iz iz
例
1 在 z 1上有一个奇点 z 1, 而它在 z 1 内处处解析 . 解 2 1 z
它在 z 1 内可展开成 z的幂级数 .
1 把函数展开成z的幂级数. 2 1 z
1 z 1 1 z z 2 1n z n , z
依次进行下去，可得到与前完全一样的展开式，这样就证明了解析函数可以展开为唯一的泰勒级数，泰勒级数与解析函数有密切的关系。
4
二、解析函数展为泰勒级数举例：例1 在z0＝0的邻域上把 f ( z) e z 展开

实验数列与级数

实验3 数列与级数级数是微积分乃至整个数学分析最重要的基本内容之一。

远在公元前三世纪，古希腊人Archimedes 就采用了数列极限的思想来计算曲边三角形的面积。

本实验的目的是通过计算机发现数列的规律、极限状态的性质。

所谓一个无穷数列是指按一定顺序排列的一串数字1a ,2a ,... ,n a , (1)而一个无穷级数则是用无穷项数字构成的和式∑∞=1n n a= 1a +2a + (2)数列与级数有密不可分的关系。

给定一个无穷级数（2），它唯一地确定了一个无穷数列 1S , 2S ,…其中n S = 1a +2a +…+n a , n = 1,2 ,… .反过来，给定一个无穷数列（1），它也唯一地确定了一个无穷级数∑∞=1n n b这里1b = 1a ,1--=n n n a a b ,n = 2 ,3 ,… 。

并且，无穷级数的和就是相应的无穷数列的极限。

因此，无穷数列与无穷级数是可以相互转化的。

给定的数列{n a } ,人们最关心的问题是：1．数列n a 有什么规律与性质？2．当n →∞时，数列n a 的极限是什么？3．极限是否是一个有限的数字？还是无穷大？抑或根本不存在？4．如果极限是无穷大，那么它趋于无穷大的阶是什么？5．如果数列的极限根本不存在，那么在无穷大的极限状态又怎么样？对于给定的一个无穷级数，也可以提出上述类似的问题。

本实验将通过计算机图示的方法来帮助我们发现数列的规律及其极限行为。

我们以Fibonacci 数列为例来探讨上述问题。

3.1 Fibonacci 数列给定如下的数列1，1，2，3，5，8，13，21，34，55，89，……其递推关系式由n n n F F F +=++12, 1=n ，2，…， 11=F ,12=F (3)给出，该数列被称为Fibonacci 数列。

Fibonacci 数列经常以著名的养兔问题提出来。

某人养了一对兔子（公母各一只）。

一月后，这对兔子生了一对小兔。

实验三级数

实验三级数【实验目的】 1．了解级数的有关理论。

2．了解函数的Taylor 展开式3．学习掌握MATLAB 软件有关的命令。

【实验内容】1．求函数x y sin =的Taylor 级数，并考察它的Taylor 展开式的前几项构成的多项式函数向x y sin =的图形的逼近的情况2．计算级数∑∞=121n n的值3．验证Euler 公式 5771.0)ln 131211(lim =-++++=∞→n nC n【实验准备】1．级数的基本概念数项级数：称用加号将数列n a 的项连成的式子+++++n a a a a 321为（常数项）无穷级数，简记为∑∞=1n na。

称级数∑∞=1n na前n 项构成的和∑==++++=nk k n n a a a a a S 1321为级数的部分和。

若S S n n =∞→lim ，则称级数∑∞=1n na收敛，其和为S 。

Taylor 级数：设函数)(x f 在包含a x =的区域内具有各阶导数，则称幂级数+-++-+-+=-∑∞=n n n n n a x n a f a x a f a x a f a f a x n a f )(!)()(!2)())((')()(!)()(2)2(0)(为函数)(x f 在a x =的Taylor 级数，当0=a 时称为Maclaurin(麦克劳林)级数。

2．级数的MATLAB 命令MATLAB 中主要用symsum,taylor 求级数的和及进行Taylor 展开。

【实验方法与步骤】练习1求函数x y sin =的Taylor 级数，并考察它的Taylor 展开式的前几项构成的多项式函数向x y sin =的图形的逼近的情况clearsyms x;taylor(sin(x),0,1)ans = 0clear syms x;taylor(sin(x),0,2)ans = x clear syms x;taylor(sin(x),0,3)ans = xclear syms x;taylor(sin(x),0,4)ans =x-1/6*x^3clear syms x;taylor(sin(x),0,5)ans =x-1/6*x^3clear syms x;taylor(sin(x),0,6) ans =x^5/120 - x^3/6 + x在区间[0,]π上做出函数sin()y x =和其泰勒展开式的前几项构成的多项式335,,3!3!5!x x x y x y x y x ==-=-+的图形clearx=0:0.01:pi; y1=sin(x); y2=x;y3=x-x.^3/6;y4=x-x.^3/6+x.^5/120;plot(x,y1,x,y2,'r:',x,y3,'k-.',x,y4,'c--')练习2 利用幂级数计算指数函数2312!3!!nxx x x e x n =++++++ 建立M 脚本文件ex0302x=input('x='); n=input('n='); y=1;for i=1:ny=y+x^i/prod(1:i); endvpa(y,10)取x=1;n=10,对比vpa(exp(1),10) 再取x=2;n=10,对比vpa(exp(2),10) 发现计算精度不高。

数学物理实验第五节(洛朗级数展开)

例4: 在z0=0的邻域上把 e1/ z 展开为洛朗级数
解: 我们知道ex在原点邻域上的展开式为
ez 1 zk 1 1 z 1 z2 1 z3 ...(| z | )
k0 k!
1! 2! 3!
把z全换成1/z,可得到以下结果:
1
ez

k 0
1 1 k
CR1
的积分,展开如下:
CR1 CR1

1
z

k 0
(z z0 )k
( z0 )k1
而对于沿C R2的积分,考虑到
z z0

z0
用以下方法将其展开
6

1
z

(

1 z0 ) (z

z0 )

z
1 z0
1
1
z0
z z0
1
上任一点z, f (z)可展开成为幂级数

k
f z ak z z0 ,
其中ak

1
2i
C

fk z0 k1
d
k

1,2,

并且展开式唯一
积分路径C为位于环域内按逆时针方向饶内圆一周的任一闭合
曲线.
5
证明: 为了避免讨论在圆周上函数的
三. 收敛环的确定
设正幂部分的收敛半径为R1；而负幂部分在变换 =1/(z-z0)下的级数的收敛半径为1/R2 ，则其在
|z-z0|>R2外收敛。如果R2<R1，那么双边幂级数就在环状域 R2<|z-z0|<R1 内收敛，所以 R2<|z-z0|<R1 给出了双边幂级数的环状收敛域，称为收敛环。

高等数学数学实验报告无穷级数

高等数学数学实验报告实验人员：院（系） _电子科学与工程学院_ 学号_06211623_ 姓名_吴晓锋_ 实验地点：计算机中心机房实验一一、实验题目观察∑∞=1!n n n n 的部分和序列的变化趋势，并求和二、实验目的和意义学会如何利用幂级数的部分和对函数进行逼近以及函数值的近似计算。

三、计算公式∑∞=1!n n n n四、程序设计（1）逼近（2）求和五、程序运行结果N[Sum[n!/n n,{n,Infinity}],50]Output= 1.87985386217525853349六、结果的讨论和分析通过利用mathematics可以直观的看出逼近图像，利用Table命令可以生成部分和的序列的数据点，同时控制点的疏密程度以利于观测。

利于软件求部分和十分快速，精确，不失为一种求和的好方法。

实验二一、实验题目观察函数,0()1,0x xf xxππ--≤<⎧=⎨≤<⎩展成的Fourier级数的部分和逼近()f x的情况。

二、实验目的和意义本实验的目的是用Mathematica显示级数部分和的变化趋势；学会如何利用幂级数的部分和对函数进行逼近以及函数值的近似计算；展示Fourier级数对周期函数的逼近情况。

三、计算公式⎰=ππ-f(x )dx π1a ⎰=ππ-nx dx x )cos (f π1n a ⎰=ππ-nx dx x )sin (f π1n b四、程序设计五、程序运行结果六、结果的讨论和分析如初值对结果的影响；不同方法的比较；该方法的特点和改进；整个实验过程中（包括程序编写，上机调试等）出现的问题及其处理等广泛的问题，以此扩大知识面和对实验环节的认识。

实验十一级数与方程符号求解

实验十一级数与方程符号求解一、实验目的1、掌握级数求和的方法2、掌握将函数展开为泰勒级数的方法3、掌握微分方程符号求解的方法4、掌握代数方程符号求解的方法二、实验内容 1、级数符号求和（1）计算101121n S n ==-∑>> n=sym('n');s=symsum(1/(2*n-1),1,10) s =31037876/14549535 （2）求级数211n n n x∞-=∑的和函数，并求215n n n ∞=∑之和。

>> n=sym('n'); x=sym('x');s=symsum(n*n*x^(n-1),n,1,inf)s =(-x-1)/(x-1)^3>> s1=symsum(n*n/5^n,n,1,inf)s1 =15/322、将ln x在1x 处按5次多项式展开为泰勒级数。

>> x=sym('x');taylor(log(x),6,1)ans =x-1-1/2*(x-1)^2+1/3*(x-1)^3-1/4*(x-1)^4+1/5*(x-1)^5 3、求下列方程的符号解。

()()()()()251ln 121sin 210335sin 78.50100043580x x xx xe x x y +-=++-=+-=-=+-=⎪⎩(1) >> x=solve('log(1+x)-5/(1+sin(x)-2)','x') x =-3.9222501966216401268562035062927-1.5308457578423990317770669453365*i(2) >> x=solve('x*x+9*sqrt(x+1)-1','x') x =-1+(1/12*(972+12*6465^(1/2))^(1/3)+2/(972+12*6465^(1/2))^(1/3)+1/2*i*3^(1/2)*(-1/6*(972+12*6465^(1/2))^(1/3)+4/(972+12*6465^(1/2))^(1/3)))^2-1+(1/12*(972+12*6465^(1/2))^(1/3)+2/(972+12*6465^(1/2))^(1/3)-1/2*i*3^(1/2)*(-1/6*(972+12*6465^(1/2))^(1/3)+4/(972+12*6465^(1/2))^(1/3)))^2-1(3) >> x=solve('3*x*exp(x)+5*sin(x)-78.5','x')x =1.1438132465904248768705131138349-3.5140576459458125 696291745801210*i(4) >> [x y]=solve('sqrt(x*x+y*y)-100','3*x+5*y-8','x,y')x =12/17-10/17*21246^(1/2)12/17+10/17*21246^(1/2)y =20/17+6/17*21246^(1/2)20/17-6/17*21246^(1/2)4、求微分方程初值问题的符号解，并与数值解进行比较。

数列与级数(mathematica数学实验报告)

可以看出， Fibonacci数列的变化速度非常快，且单调递增趋于无穷；从图象中也可明显看出n 取值越大，图像越陡，即递增越快。

事实上，由Fibonacci 数列的递推关系式2112,1,2,...,1,1n n n F F F n F F ++=+===，（1）容易得到12113/22,n n n n n F F F F F ++++<=+< （2）因此，n F 的阶应该在()3/2n与2n 之间。

为进一步研究Fibonacci 数列n F 的特性，我们将n F 取对数，在直角坐标系中画出顺次连接点()(),log ,1,2,...n n F n N =的折线图。

此时的折线图近乎于一条直线。

因此，我们猜测()log n F 是n 的线性函数。

取1000N =，对上述数据进行拟合可得()log 0.8039030.481211n F n ≈-=，（3）故0.447567 1.61803n n F ≈⨯. （4）2.下面，我们分别取50,100,500,1000n =，利用Mathematica 编程，用直线去拟合上述数据()(),log ,1,2,...n n F n N =，由此来求数列n F 的近似表示。

过程如下：可以看出，给定的n值越大，线性拟合的结果便趋于稳定，而且，对每一组拟合的线性方程，其系数与黄金分割数有着紧密的联系。

由计算机观察得到的上述结果我们似乎可F的通项具有形式以猜测数列nn n F cr = （5）将上式代入递推公式（1）得21r r =+ （6）从而()15/2r =+.因为数列趋于无穷，故取()15/2r =+。

于是152n n F c ⎛⎫+= ⎪ ⎪⎝⎭ （7）然而，公式（7）并不满足121F F ==，即并非数列n F 的通项公式.不过，它仍然是数列n F 的主项.3.取一组整数50,100,500,1000,5000,10000n =，将Fibonacci 数列模n 得到一周期数列，将该周期数列的值作为高音，编程演奏它.运行结果如下：根据运行结果，明显可以看出，n的取值越大，图像上的点越稠密.实验结果和结果分析：附录：。

数学实验 Mathematic实验九无穷级数讲解

天水师范学院数学与统计学院实验报告实验项目名称无穷级数所属课程名称数学实验实验类型微积分实验实验日期2011.11.16班级学号姓名成绩fx=Normal[Series[Exp[x]，{x，0，3}]]Plot[fx，{x，-3，3}]则只能得到去掉余项后的展开式，得不到函数的图形，这时要使用强制求值命令Evaluate，改成输入Plot[Evaluate[fx]，{x，-3，3}]便可以得到函数的图形5.作散点图命令ListPlot.ListPlot[Table[j^2，{j，16}]，PlotStyle PointSize[0.012]] 6.用符号“／；”定义分段函数.符号“／；”用于定义某种规则，“／；”后面是条件，例如输入Clear[g，gf]；g[x_]:=x/；0x<1g[x_]:=-x/；-1x<0g[x_]:=g[x-2]/；x 1gf=Plot[g[x]，{x，-1，6}]用which命令也可以定义分段函数，从这个例子中看到，用“…(表达式)/；…(条件)”来定义周期性分段函数更方便些.用Plot命令可以作出分段函数的图形，但用Mathematica命令求分段函数的导数或积分时往往会有问题.用which定义的分段函数可以求导，但不能积分.Mathematica内部函数中有一些也是分段函数，如：Mod[x，1]，Abs[x]，Floor[x]和Unitstep[x].其中只有单位阶跃函数Uniltstep[x]可以用Mathematica命令来求导和求定积分,因此在求ListPlot[vals，PlotStyle PointSize[0.012]] Sum[a[n]，{n，1，Infinity}]2.求幂级数的收敛域.例9.4 求24(3)1n nnxn∞=-+∑收敛域与和函数.Clear[a]；a[n_]=4^(2n)*(x-3)^n/(n+1)；stepone=a[n+1]/a[n]//Simplifysteptwo=Limit[stepone，n Infinity]ydd=Solve[steptwo1，x]zdd=Solve[steptwo-1，x]Simplify[a[n]/.x(49/16)]Simplify[a[n]/.x(47/16)]Sum[4^(2n)*(x-3)^n/(n+1)，{n，0，Infinity}] 3.函数的幂级数展开.例9.5 求cos x的6阶麦克劳林展开式.Series[Cos[x]，{x，0，6}]例9.6 求ln x在1x=处的6阶泰勒展开式.Series[Log[x]，{x，1，6}]例9.7 求arctan x的5阶麦克劳林展开式.ser1=Series[ArcTan[x]，{x，0，5}]；poly=Normal[ser1]Plot[Evaluate[{ArcTan[x]，poly}]，{x ，-3/2，3/2}，PlotStyle {Dashing[{0.01}]，GrayLevel[0]}，AspectRatio 1]例9.8 求22(1)(1)x x e --+在1x =处的8阶泰勒展开，并通过作图比较函数和它的近似多项式.Clear[f]；f[x_]=Exp[-(x-1)^2*(x+1)^2]； poly2=Normal[Series[f[x]，{x ，1，8}]] Plot[Evaluate[{f[x]，poly2}]，{x ，-1.5，1.5}，PlotRange {-2，3/2}，PlotStyle {Dashing[{0.01}]，GrayLevel[0]}]例9.9 求函数x sin 在0=x 处的3，5，7，…，9l 阶泰勒展开，通过作图比较函数和它的近似多项式，并形成动画进一步观察.Do[Plot[{Sum[(-1)^j*x^(2j+1)/(2j+1)!，{j ，0，k}]，Sin[x]}，{x ，-40，40}，PlotStyle {RGBColor[1，0，0]，RGBColor[0，0，1]}]，{k ，1，45}] 4.傅里叶级数.例9.10 设()f x 是周期为2的周期函数它在一个周期内的表达式为1,01(),10x f x x x ≤<⎧=⎨--≤<⎩求它的傅立叶级数展开式的前5项和前8项，作出()f x 和它的近似三角级数的图形.Clear[f ，a ，b ，fs ，L]； f[x_]:=1/；0x<1 f[x_]:=-x/；-1x<0 f[x_]:=f[x-2]/；1x gf=Plot[f[x]，{x ，-1，5}] Clear[L ，a ，b ，fs ，f1，f2]； L=1；a[n_]:=(Integrate[-x*Cos[n*Pi*x/L]，{x，-L，0}]+Integrate[Cos[n*Pi*x/L]，{x ，0，L}])/Lb[n_]:=(Integrate[-x*Sin[n*Pi*x/L]，{x，-L，0}]+Integrate[Sin[n*Pi*x/L]，{x ，0，L}])/Lfs[k_，x_]:=a[k]*Cos[k*Pi*x/L]+b[k]*Sin[k*Pi*x/L] fourier[n_，x_]:=a[0]/2+Sum[fs[k ，x]，{k ，1，n}] f1=fourier[5，x]//N f2=fourier[10，x]//NPlot[Evaluate[{f[x]，f1}]，{x ，-1，5}，PlotStyle {GrayLevel[0]，GrayLevel[0.4]}]Plot[Evaluate[{f[x]，f2}]，{x ，-1，5}，PlotStyle {GrayLevel[0]，GrayLevel[0.4]}]设)(x g 是以2Pi 为周期的周期函数，它在],[ππ-的表达式是1,0()1,0x g x x ππ--≤<⎧=⎨≤<⎩，将)(x g 展开成傅里叶级数. Clear[g]；g[x_]:=-1/；-Pi x<0 g[x_]:=1/；0x<Pi g[x_]:=g[x-2Pi]/；Pi xPlot[g[x]，{x ，-Pi ，5Pi}，PlotStyle {RGBColor[0，1，0]}]； Clear[b2，fourier2，tu ，tu2，toshow]；b2[n_]:=b2[n]=2Integrate[1*Sin[n*x]，{x ，0，Pi}]/Pi ； fourier2[n_，x_]:=Sum[b2[k]*Sin[k*x]，{k ，1，n}]； tu[n_]:=Plot[{g[x]，Evaluate[fourier2[n ，x]]}，{x ，-Pi ，5Pi}， PlotStyle{RGBColor[0，1，0]，RGBColor[1，0.3，0.5]}，DisplayFunctionIdentity]；tu2=Table[tu[n]，{n ，1，30，5}]； toshow=Partition[tu2，2]； Show[GraphicsArray[toshow]]【实验结论】（结果）1.用Mathematica 求无穷级数的和；2.求幂级数的收敛域；3.展开函数为幂级数以及展开周期函数为傅里叶级数.附录1：源程序1Sum k2^k,k,1,Infinity2Sum12k1^2,k,1,Infinity 283Sum12k^2,k,1,Infinity 2244Sum1^k1k,k,1,Infinity Log2Clear a;a n_x1^2n15^n; stepone a n1a n Simplify11x25steptwo Limit stepone,n Infinity11x25ydd Solve steptwo1,xzdd Solve steptwo1,xx15,x15x15,x15x15,x15Simplify a n.x1Sqrt5Sin kkSimplify a n.x1Sqrt5Sin kkSum x1^2n15^n,n,0,Infinity 51x62x x2Series1x Log1x,x,0,6Log2Log x Log x2O x7Series ArcSin x,x,0,6x x363x540O x7Clear f;f x_x x^21;Series f x,x,0,5Series f x,x,0,10p1Normal Series f x,x,0,5p2Normal Series f x,x,0,10p3Plot Evaluate f x,p1,p2,x,3,3,PlotRange2,32, PlotStyle Dashing0.01,GrayLevel0x x3x5O x6x x3x5x7x9O x11x x3x5x x3x5x7x9GraphicsClear f,a,b,fs,L;f x_:1x^2;12x12 f x_:f x1;x12gf Plot f x,x,1,5GraphicsL12;a n_:Integrate x Cos n Pi x L,x,L,0Integrate Cos n Pi x l,x,0,LLb n_:Integrate x Sin n Pi x L,x,L,0Integrate Sin n Pi x l,x,0,L L fs k_,x_:a k Cos k Pi x L b k Sin k Pi x Lfourier n_,x_:a02Sum fs k,x,k,1,6f1fourier5,x Nf2fourier10,x N0.625 2.Cos 6.28319x0.05066060.31831l Sin 1.5708l0.31831l Cos12.5664x Sin 3.14159l2.Cos18.8496x0.005628950.106103l Sin 4.71239l0.159155l Cos25.1327x Sin 6.28319l2.Cos31.4159x0.002026420.063662l Sin 7.85398l0.106103l Cos37.6991x Sin 9.42478l2.0.07957750.31831l0.31831l Cos 1.5708lSin 6.28319x2.0.03978870.159155l0.159155l Cos3.14159lSin12.5664x2.0.02652580.106103l0.106103l Cos 4.71239lSin18.8496x2.0.01989440.0795775l0.0795775l Cos 6.28319lSin25.1327x2.0.01591550.063662l0.063662l Cos 7.85398lSin31.4159x2.0.01326290.0530516l0.0530516l Cos 9.42478lSin37.6991x0.625 2.Cos 6.28319x0.05066060.31831l Sin 1.5708l0.31831l Cos12.5664x Sin 3.14159l2.Cos18.8496x0.005628950.106103l Sin 4.71239l0.159155l Cos25.1327x Sin 6.28319l2.Cos31.4159x0.002026420.063662l Sin 7.85398l0.106103l Cos37.6991x Sin 9.42478l2.0.07957750.31831l0.31831l Cos 1.5708lSin 6.28319x2.0.03978870.159155l0.159155l Cos3.14159lSin12.5664x2.0.02652580.106103l0.106103l Cos 4.71239lSin18.8496x2.0.01989440.0795775l0.0795775l Cos 6.28319lSin25.1327x2.0.01591550.063662l0.063662l Cos 7.85398lSin31.4159x2.0.01326290.0530516l0.0530516l Cos 9.42478lSin37.6991xPlot Evaluate f x,f1,x,1,5,PlotStyle GrayLevel0,GrayLevel0.4Plot Evaluate f x,f2,x,1,5,PlotStyle GrayLevel0,GrayLevel0.4Plot::plnr:f x is not a machine size real number at x 1..Plot::plnr:f x is not a machine size real number at x0.756598.Plot::plnr:f x is not a machine size real number at x0.629911.General::stop:Further output of Plot::plnr will be suppressed during this calculation.GraphicsPlot::plnr:f x is not a machine size real number at x 1..Plot::plnr:f x is not a machine size real number at x0.756598.Plot::plnr:f x is not a machine size real number at x0.629911.General::stop:Further output of Plot::plnr will be suppressed during this calculation.GraphicsClear f,a,b,fs,L;f x_:1;0x1f x_:2x;1x2 f x_:f x2;x2 gf Plot f x,x,1,5GraphicsL1;a n_:Integrate x Cos n Pi x L,x,L,0Integrate Cos n Pi x l,x,0,LLb n_:Integrate x Sin n Pi x L,x,L,0Integrate Sin n Pi x l,x,0,L L fs k_,x_:a k Cos k Pi x L b k Sin k Pi x Lfourier n_,x_:a02Sum fs k,x,k,1,8f1fourier5,x Nf2fourier10,x N0.75Cos 3.14159x0.2026420.31831l Sin 3.14159l0.159155l Cos 6.28319x Sin 6.28319lCos9.42478x0.02251580.106103l Sin 9.42478l0.0795775l Cos12.5664x Sin 12.5664lCos15.708x0.008105690.063662l Sin 15.708l0.0530516l Cos18.8496x Sin 18.8496lCos21.9911x0.004135560.0454728l Sin 21.9911l0.0397887l Cos25.1327x Sin 25.1327l0.318310.31831l0.31831l Cos 3.14159lSin 3.14159x0.1591550.159155l0.159155l Cos 6.28319lSin 6.28319x0.1061030.106103l0.106103l Cos 9.42478lSin9.42478x0.07957750.0795775l0.0795775l Cos 12.5664lSin12.5664x0.0636620.063662l0.063662l Cos 15.708lSin15.708x0.05305160.0530516l0.0530516l Cos 18.8496lSin18.8496x0.04547280.0454728l0.0454728l Cos 21.9911lSin21.9911x0.03978870.0397887l0.0397887l Cos 25.1327lSin25.1327x0.75Cos 3.14159x0.2026420.31831l Sin 3.14159l0.159155l Cos 6.28319x Sin 6.28319lCos9.42478x0.02251580.106103l Sin 9.42478l0.0795775l Cos12.5664x Sin 12.5664lCos15.708x0.008105690.063662l Sin 15.708l0.0530516l Cos18.8496x Sin 18.8496lCos21.9911x0.004135560.0454728l Sin 21.9911l0.0397887l Cos25.1327x Sin 25.1327l0.318310.31831l0.31831l Cos 3.14159lSin 3.14159x0.1591550.159155l0.159155l Cos 6.28319lSin 6.28319x0.1061030.106103l0.106103l Cos 9.42478lSin9.42478x0.07957750.0795775l0.0795775l Cos 12.5664lSin12.5664x0.0636620.063662l0.063662l Cos 15.708lSin15.708x0.05305160.0530516l0.0530516l Cos 18.8496lSin18.8496x0.04547280.0454728l0.0454728l Cos 21.9911lSin21.9911x0.03978870.0397887l0.0397887l Cos 25.1327lSin25.1327xPlot Evaluate f x,f1,x,1,5,PlotStyle GrayLevel0,GrayLevel0.4Plot Evaluate f x,f2,x,1,5,PlotStyle GrayLevel0,GrayLevel0.4Plot::plnr:f x is not a machine size real number at x 1..Plot::plnr:f x is not a machine size real number at x0.756598. Plot::plnr:f x is not a machine size real number at x0.491147. General::stop:Further output of Plot::plnr will be suppressed during this calculation.GraphicsPlot::plnr:f x is not a machine size real number at x 1..Plot::plnr:f x is not a machine size real number at x0.756598. Plot::plnr:f x is not a machine size real number at x0.491147. General::stop:Further output of Plot::plnr will be suppressed during thisGraphicsClear a;a n_Sin k k;vals Table a k,k,1,50;ListPlot vals,PlotStyle PointSize0.015Graphics附录2：实验报告填写说明1.实验项目名称：要求与实验教学大纲一致。

数学实验教程实验9(级数

实验9 级数实验目的1．理解幂级数的概念，并会用软件将函数展开成幂级数 2．理解Fourier 级数的概念，并将函数展开成Fourier 级数实验准备1．数项级数、幂级数的收敛性判断； 2．幂级数的展开、级数求和； 3．Fourier 级数的概念、展开方法；实验内容1．函数的幂级数展开 2．收敛级数的和 3．Fourier 级数展开软件命令表9-1 Matlab 级数操作命令实验示例【例9.1】级数观察观察下列级数的部分和序列的变化趋势，并求和。

1. 11n n ∞=∑； 2. 11(1)n n n ∞=-∑。

【步骤】：Step1：计算部分和n S ； Step2：描点观察。

【程序】：clearclc clffor n=1:100 for k=1:np1(k)=1/k;p2(k)=(-1)^k/k; ends1(n)=sum(p1); s2(n)=sum(p2); endplot(s1) plot(s2) syms i;symsum(1/i,i,1,inf))symsum((-1)^i/i,i,1,inf))【输出】：图 9-1 部分和序列收敛性观察级数（1）发散；调和级数（2）收敛，收敛于ln2。

【例9.2】调和级数实验—欧拉常数记11()ni H n i ==∑，()()ln C n H n n =-，研究C(n)的极限值是否存在。

【程序】：%图形观察h(1)=1;for i=2:10^5h(i)=h(i-1)+double(1/i); c(i)=h(i)-log(i); end plot(c)% 求极限syms k nlimit(symsum(1/k,k,1,n)-log(n),n,inf) 【例9.3】函数的幂级数展开将下列函数在指定点处展开成幂级数，并计算近似值，至少保留三位小数。

1．330()1,1,9f x x x += 2．011()arctan,1,arctan 12x f x x x -==+； 3．0()sin(1),0,sin1f x x x =+=。

级数实验

大学数学实验报告
无穷级数与函数逼近
班级＿＿＿＿＿＿＿
姓名＿＿＿＿＿＿＿
学号＿＿＿＿＿＿＿
指导教师＿＿＿＿＿＿＿
实验时间＿＿＿＿＿＿＿
实验所用软件及版本： 4.1
实验目的：
（1）掌握用软件进行级数运算、求傅立叶级数的语句和方法。
（2）用软件显示级数的部分和的变化趋势。
（3）学会如何利用幂级数的部分和对函数的逼近以及进行函数值的近似计算。
（2）将在点处展开到的4次幂。
3.利用幂级数展开式计算
4.设周期为的周期函数在一个周期内的表达式为
试生成的傅里叶级数，并从图上观察该级数的部分和逼近的情况。
实验步骤与实验结论
选作实验：
设是周期为的周期函数，它在上的表达式为，将分解成简谐振动的叠加形式。
实验步骤与实验结论
创新实验：
试在同一屏幕上作出及它的阶傅里叶多项式的图形，观察逼近的情况。
实验步骤与实验结论
教师评：
（4）展示傅里叶级数对周期函数的逼近情况。
实验涉及的语句：
语句一：
－求级数的和（和函数），其中为级数的通项，为的起始值，为终止值，当从1开始时，可以省略。
语句二：
－把函数在点处展开到 - 的次幂。
语句三：
语句四：
必作实验：
1.观察的部分和序列的变化趋势，并求和。
2.（1）分别将，，在点处展开到的5次幂。

实验十一级数与方程符号求解.doc

实验十一级数与方程符号求解一、实验目的1、掌握级数求和的方法2、掌握将函数展开为泰勒级数的方法3、掌握微分方程符号求解的方法4、掌握代数方程符号求解的方法二、实验内容1、级数符号求和10（1）计算S = Z>> n=sym(’ n’)；s=symsum(l/ (2*n~l),1,10)31037876/14549535oo(2)求级数的和函数，并求艺二n=l"=1 5 >> n=sym(，n，):x=sym(’ x’)；s二symsum(n本n本x" (n一1)，n, 1，inf)>> sl=symsum(n*n/5^n, n, 1, inf)15/322、将1M在;v = l处按5次多项式展开为泰勒级数。

>> x=sym(’ x’)； taylor (log(x), 6, 1) ans =x-l-l/2*(x-l)"2+l/3*(x-l)"3-l/4*(x-l)"4+1/5*(x_l厂53、求下列方程的符号解。

(2) r + 9A / x +1 —1—0 (3) 3x^ v +5 sin A :-78.5 = 0(4) |V %2+夕2 -100 = 0[3^ + 5^-8 = 0 (1)〉〉x=solve(’log(l+x)-5/(l+sin(x)-2)’，’ x’)-3. 9222501966216401268562035062927-1. 5308457578423990317770 669453365*i(2)〉〉x=solve (’x*x+9*sqrt (x+l)-l’，’ x’)-1+(1/12*(972+12*6465"(1/2)厂(1/3)+2/(972+12*6465" (1/2))^(1 /3)+l/2*i*3" (1/2)*(-1/6*(972+12*6465"(1/2)厂(1/3)+4/(972+12 *6465" (1/2))"(1/3))厂2-1+(1/12*(972+12*6465" (1/2)V(1/3)+2/(972+12*6465"(1/2)V(1 /3)-l/2*i*3" (1/2)*(-1/6*(972+12*6465" (1/2)厂(1/3)+4/(972+12 ⑴ ln(l + x) - 5 1 + sinx=2*646571/2))71/3))厂2-1(3)>> x=solve (’ 3*x*exp (x) +5*sin (x) -78. 5’，’ x’)1.1438132465904248768705131138349-3. 51405764594581256962917 45801210*i(4)>> [x y] =solve (’ sqrt (x*x+y*y) -100’，’ 3*x+5*y-8’，’ x，y’)12/17-10/17*2124671/2)12/17+10/17*21246^(1/2)20/17+6/17*21246"(1/2)20/17-6/17*2124671/2)4、求微分方程初值问题的符号解，并与数值解进行比较。

实验九常微分方程与级数

例:求函数y=sin(x)在x=0点处的5阶taylor展开式及在处的6阶taylor展开式. 解：输入命令： >>syms x; >>taylor(sin(x)) 结果：ans = x-1/6*x^3+1/120*x^5 >>taylor(sin(x),pi/3,7) 结果：ans =1/2*3^(1/2)+1/2*x-1/6*pi-1/4*3^(1/2)*(x-1/3*pi)^2-1/12*(x-1/3*pi)^3+1/48*3^(1/2)*(x-1/3*pi)^4+1/240*(x-1/3*pi)^5-1/1440*3^(1/2)*(x-1/3*pi)^6
入命令: >>dsolve('D2y=x+Dy', 'y(0)=1', 'Dy(0)=0', 'x') 结果: ans = -1/2*x^2+exp(x)-x
用Matlab求已知函数的泰勒展开式求已知函数的泰勒展开式用命令taylor,其几种调用格式为: taylor(函数f(x)) 求f(x)在x=0点的5次taylor多项式. taylor(函数f(x),n) 求f(x)在x=0点的n-1次taylor多项式. taylor(函数f(x),a,n) 求f(x)在x=a点的n-1次taylor多项式.
实验九常微分方程与级数
CLICK HERE TO ADD A TITLE
实验目的: 1.掌握用Matlab命令dsolve求解微分方程. 2.学习Matlab泰勒级数展开命令taylor. 3.加强对幂级数和函数概念的理解. 4.如何用Matlab判断数项级数的敛散性.
一.用Matlab命令dsolve求解微分方程一般地,含有未知函数及未知函数的导数或微分的方程称为微分方程.微分方程在Matlab中输入时应注意: 应输入Dy, 应输入D2y, 应输入D3y等等. 1.Matlab求解微分方程的命令是dsolve,调用格式: dsolve(‘微分方程’)给出微分方程的解析解,表示为t的函数,把y 当作t的函数求解,这是系统默认的. 2.dsolve(‘微分方程’,’x’) 给出微分方程的解析解,表示为x的函数.