2014年下学期数学实验与数学建模作业习题6
数学建模作业、微分方程实验、北京工业大学

2微分方程实验1、微分方程稳定性分析绘出下列自治系统相应的轨线,并标出随 t 增加的运动方向,确定平■衡点, 并按稳定的、渐近稳定的、或不稳定的进行分类:解:(1)由 f (x ) =x=0, f (y ) =y=0;可得平衡点为(0,0),___ 1 0系数矩阵A,求得特征值入1=1,入2=1;0 1p=-(入1+入2)=-2<0 , q=入1入2=1>0;对照稳定性的情况表,可知平■衡点(0, 0) 是不稳定的。
图形如下:(2)如上题可求得平衡点为(0,0 ),特征值入1=-1,入2=2;p=-(入1+入2)=-1<0 , q-入1入2=-2<0;对照稳定性的情况表,可知平■衡点(0, 0) 是不稳定的。
其图形如下:dx⑴dt dtx, y;dxdtdydt dx x, ⑶尸 2y ;晋 dx y, (4) ? 2x;也 dtx+1, 2y.(3) 如上题可求得平■衡点为(0,0 ),特征值入1=0 + 1.4142i,入2=0 -1.4142i; p=-(入1+入2)= 0, q-入1入2=1.4142>0;对照稳定性的情况表,可知平■衡点(0, 0)是不稳定的。
其图形如下:(4) 如上题可求得平衡点为(1,0 ),特征值入1=-1,入2=-2;p=-(入1+入2)= 3>0, q=入1入2=2>0;对照稳定性的情况表,可知平■衡点(1, 0) 是稳定的。
其图形如下:2、种群增长模型一个片子上的一群病菌趋向丁繁殖成一个圆菌落.设病菌的数目为N,单位成员的增长率为r1,则由Malthus生长律有竺r1 N,但是,处丁周界表面的dt那些病菌由丁寒冷而受到损伤,它们死亡的数量与N2成比例,其比例系数为r2, 求N满足的微分方程.不用求解,图示其解族.方程是否有平衡解,如果有,是否为稳定的?解:由题意很容易列出N满足的微分方程:坐r1N r2N; f(N)dt令f(N)=O,可求得方程的两个平■衡点N1=0,N2=「22/r i21 1d2N 1 5 52 (r1 r2N 2) (r1N r2N 2)dt 2进而求得A d2N 令r dt2 2 0可求得N=r2 /4r〔则N=N1 N=N2 N=r22/4r i2可以把第一象限划为三部分,且从下到上三部分中分0,冬dt2.2 2 c dN cdN c dN cdN 0, ;—0, —r 0; —0, ―rdt dt dt dt则可以画出N (t) 的图形,即微分方程的解族,如下图所示:由图形也可以看出,对丁方程的两个平■衡点,其中N1=0是不稳定的;N2=^2 /「;是稳定的o3、有限资源竞争模型1926年Volterra 提出了两个物种为共同的、有限的食物来源而竞争的模型当[b MX h 2X 2)]x dt dX2 电 2(h i X i h 2X 2)]X 2dt假设也 坦,称垣为物种i 对食物不足的敏感度,(1) 证明当x1(t0)>0时,物种2最终要灭亡; (2) 用图形分析方法来说明物种 2最终要灭亡.解:(1)由上述方程组 f (x1) =[b 1〔S' h 2x 2)]x 1=0,f (x2)=电2 (h 1X 1h 2X 2)]X 2=0,可得方程的平■衡点为R (0,0), P 1 (E,0),P 2 (0, M).2 h 2对平衡点P 。
《数学建模与实验》习题库a

高于 M,因为环境无法支持,数量将会下降。在下面模型中, an 表示 n 年后的鲸鱼数量;
试讨论模型: an+1 − an = k(M − an )(an − m)
6.假设存在某种药物,当其浓度大于 100 毫克/升时,可以治疗疾病。药物的初始浓度是 640 毫克/升,从实验知道药物是以每小时现有量的 20%的比率衰减
参考文献 数学建模(英文版), 机械工业出版社, 北京, 2003. 5. (经典原版书库, 原书名: A First
Course in Mathematical Modeling (Third Edition), by Frank R. Giordano, Maurice D. Weir, William P. Fox.)
Hyundai 车价$12,400 预付$500 月利率 6.5%直到 48 个月
你每个月为买车最多能付 475 美元。利用动力系统模型来决定你应该买哪家公司的车。
1.2 习题
1.从引进到 Tasmania 岛的新环境里的养群数量的增长得到下面的数据。
年
1814 1824 1834 1844 1854 1864
(a)构建一个表示每小时浓度的模型
2
(b)建立一张浓度值表并确定何时浓度到达 100 毫克/升 7.利用习题 6,研制的模型开一个初始剂量处方,以及一个能把浓度保持在高出有效水平 500ppm(即百万分之五百,或万分之五)但低于安全水平 1000ppm 的维持剂量处方。用不 同的值来做实验,直到结果满意为止。 8.在一处古篝火遗址附近发现了一个猿人头骨。考古学家确信该头骨和古篝火是同时代的。 实验测试确定取自篝火的灰烬中,仅留存原来的碳 14 量的 1%,已知碳 14 以与其剩余量成 比例的比率衰减而且碳 14 在 5700 年里衰减掉 50%。构建一个碳 14 测定年代的模型。 9.附表的数据展示了一辆汽车的数据 n(以 5 英里/小时的增量计)以及从刹车到停止的(滑
数学实验与数学建模(校本教材)

第一章数学软件简介第一节 MatLab简介及基本运算一、MatLab简介MATLAB名字由MATrix和 LABoratory 两词的前三个字母组合而成。
那是20世纪七十年代后期的事:时任美国新墨西哥大学计算机科学系主任的Cleve Moler教授出于减轻学生编程负担的动机,为学生设计了一组调用LINPACK和EISPACK库程序的“通俗易用”的接口,此即用FORTRAN编写的萌芽状态的MATLAB。
经几年的校际流传,在Little的推动下,由Little、Moler、Steve Bangert合作,于1984年成立了MathWorks公司,并把MATLAB正式推向市场。
从这时起,MATLAB的内核采用C语言编写,而且除原有的数值计算能力外,还新增了数据图视功能。
MATLAB以商品形式出现后,仅短短几年,就以其良好的开放性和运行的可靠性,使原先控制领域里的封闭式软件包(如英国的UMIST,瑞典的LUND和SIMNON,德国的KEDDC)纷纷淘汰,而改以MATLAB为平台加以重建。
在时间进入20世纪九十年代的时候,MATLAB 已经成为国际控制界公认的标准计算软件。
在欧美大学里,诸如应用代数、数理统计、自动控制、数字信号处理、模拟与数字通信、时间序列分析、动态系统仿真等课程的教科书都把MATLAB作为内容。
这几乎成了九十年代教科书与旧版书籍的区别性标志。
在那里,MATLAB是攻读学位的大学生、硕士生、博士生必须掌握的基本工具。
在国际学术界,MATLAB已经被确认为准确、可靠的科学计算标准软件。
在许多国际一流学术刊物上,(尤其是信息科学刊物),都可以看到MATLAB的应用。
在设计研究单位和工业部门,MATLAB被认作进行高效研究、开发的首选软件工具。
如美国National Instruments公司信号测量、分析软件LabVIEW,Cadence公司信号和通信分析设计软件SPW等,或者直接建筑在MATLAB之上,或者以MATLAB为主要支撑。
2013年下学期数学实验作业

数学实验与数学建模实验报告学院:专业班级:姓名:学号:完成时间:2014 年1 月6日实验一 图形的画法1. 做出下列函数的图像:(1))2sin()(22--=x x x x y ,22≤≤-x (分别用plot 、fplot ) (2)22/9/251x y +=(用参数方程)(3) 在同一图形窗口中,画出四幅不同图形(用subplot 命令):1cos()y x =,2sin(/2)y x pi =-,23cos()y x x pi =-,sin()4x y e =(]2,0[π∈x )2 作出极坐标方程为)cos 1(2t r -=的曲线的图形.3 作出极坐标方程为10/t e r =的对数螺线的图形.4 绘制螺旋线⎪⎩⎪⎨⎧===t z t y t x ,sin 4,cos 4在区间[0,π4]上的图形.在上实验中,显示坐标轴名称。
5 作出函数22y x xye z ---=的图形.6 作出椭球面1194222=++z y x 的图形.(该曲面的参数方程为,cos ,sin sin 3,cos sin 2u z v u y v u x === (ππ20,0≤≤≤≤v u ).)7 作双叶双曲面13.14.15.1222222-=-+z y x 的图形.(曲面的参数方程是,csc 3.1,sin cot 4.1,cos cot 5.1u z v u y v u x ===其中参数πππ<<-≤<v u ,20时对应双叶双曲面的一叶, 参数πππ<<-<≤-v u ,02时对应双叶双曲面的另一叶.)8 作出圆环v z u v y u v x sin 7,sin )cos 38(,cos )cos 38(=+=+=,(πππ22/,2/30≤≤≤≤v u )的图形.9 作出球面22222=++z y x 和柱面1)1(22=+-y x 相交的图形.10 作出锥面222z y x =+和柱面1)1(22=+-y x 相交的图形.11用动画演示由曲线],0[,sin π∈=z z y 绕z 轴旋转产生旋转曲面的过程. (该曲线绕z 轴旋转所得旋转曲面的方程为,sin 222z y x =+ 其参数方程为])2,0[],,0[(,,sin sin ,cos sin ππ∈∈===u z z z u z y u z x ) 12. 画出变上限函数⎰xdt t t 02sin 及其导函数的图形.13.迪卡尔曲线)03(13,1333222=-++=+=axy y x tat y t at x 14.蔓叶线)(1,1322322x a x y tat y t at x -=+=+= 15.摆线)cos 1(),sin (t b y t t a x -=-=16.内摆线(星形线))(sin ,cos 32323233a y x t a y t a x =+==17.圆的渐伸线(渐开线))cos (sin ),sin (cos t t t a y t t t a x -=+=18.空间螺线ct z t b y t a x ===,sin ,cos 19.阿基米德线0,≥=r a r ϕ。
2014年全国数学建模a题解析

承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》(以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛网站下载)。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为,我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。
(论文纸质版与电子版中的以上信息必须一致,只是电子版中无需签名。
以上内容请仔细核对,提交后将不再允许做任何修改。
如填写错误,论文可能被取消评奖资格。
)赛区评阅编号(由赛区组委会评阅前进行编号):编号专用页赛区评阅编号(由赛区组委会评阅前进行编号):全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):嫦娥三号软着陆轨道设计与控制策略摘要嫦娥三号卫星着陆器实现了我国首次地外天体软着陆任务。
要保证准确的在月球预定区域内实现软着陆轨道与控制策略的设计。
问题一运用活力公式[1]来建立速度模型,利用matlab软件代入数值计算出。
所求速度33⨯⨯(=1.692210m/s,=1.613910m/s)v v远近采用轨道六根数[2]来建立近月点,远月点位置的模型。
轨道根数是六个确定椭圆轨道的物理量,也是联系赤道直角坐标与轨道极坐标重要夹角的关系。
通过着陆点的位置求出轨道根数各个值的数据,从而确定近月点,远月点的位置,坐标分别为(19.51W 27.88N 15KM),(160.49 27.885S 100KM)E。
2014年武汉理工大学数学建模训练题目2

2014年武汉理工大学数学建模训练题目
第2题:实验安排问题
(陈建业老师提供)
2014年某专业职业技能考试共有技能考题5道,要求每位考生必须完成其中三道,除第一题为必考项目外,余下考题四选二。
考生总人数79人,每位考生所选选项均已知,见附件1。
受考场数量限制,只能提供A、B两个技能考场,其中A考场能容纳24个人同时开考,B考场能容纳32人同时开考。
为便于管理每个技能考题只能安排在同一个考场,每个技能考题同时开考的最大数量没有限制,但是一旦确定中途就不能更改。
请分别在以下条件下合理设计考场分布和安排实验顺序,使得本次考试所用时间最短。
(1)假定技能考题1考试时间为45分钟,而其他4个技能考题考试时间为30分钟,考生不能提前离开考场。
(2)假定增加C考场能容纳20人,而3个考场不在同一个地方,3个考场之间往返的时间见附件2,请结合(1)的条件给出合理的实验安排。
(3)假定(1)的考试时间均为最长考试时间,考生做完实验可以提前离开考场,这样可以提前安排下一位考生进场,请考虑这种情况下的实验安排。
2014南理工数学建模题目A题
A 题
“十字线”交叉点定位法
在实际工作中,为了确定目标的位置,通常采用“十字线”交叉点定位法(如图1)确定目标的距离和方向,它的关键问题为:如何准确地确定“十字线”的交叉点位置?解决这个问题的难点是因为有以下2种情况的出现:(1)从目标返回的图像,由原来的水平和竖直的直线变成了不水平或者不垂直的(如图2)状态,甚至变成弯曲状态(如图3)。
(2)由于噪声、图像设备缺陷、光线影响等很多因素导致实际图像存在模糊、光强不均匀、线条变形等等(如图4)。
本课题要求建立数学模型,解决以下问题:
(1)建立“十字线”交叉点求解模型,确定图1、图2、图3这三种类型的直线或者曲线交叉点位置(即横坐标和纵坐标)。
(2)建立亚像素级的“十字线”交叉点定位模型(即低于1个像素精度的“十字线”交叉点位置求解模型),确定图1、图2、图3这三种类型的直线或者曲线交叉点位置。
可以直接解决问题(2),而不必一定解决问题(1),但是解决问题(1)可能有助于问题(2)的解决。
(3)针对实际捕捉到的图像(见附件images.rar)建立相应的“十字线”交叉点位置求解模型。
注:图1、图2、图3的原文件见相应的附图。
图1 横竖四直线图.bmp
图2 四直线图.bmp
图3 四曲线图.bmp
图4 images\160803.bmp。
2014版更新高等数学作业题参考答案20140410
高等数学作业题参考答案(2014更新版)一、单项选择题1. D2. B3. B4. A5. B6. B7. A8. B9. B 10. C 11. B 12. B 13. B 14. B 15. C 16. B 17. D 18. B 19. B 20. A 21. B 22. C 23. D 24. A 25. C二、填空题1.[)(]2,11,2 -2. 3-=x3. 可导4. 下5. 母线为z 轴,2240x y z ⎧+=⎨=⎩为准线的圆柱面6. 无限增大 (或∞→)7. )0,1(-;),0(+∞8.(){}x y x y x <<-,9.2e三、计算题1. 解:13121lim -→⎪⎭⎫ ⎝⎛-xx x⎪⎭⎫⎝⎛-⎪⎭⎫ ⎝⎛-⋅-→⎪⎭⎫ ⎝⎛-=13122021lim x x x x x ⎪⎭⎫ ⎝⎛+-⋅-→⎪⎭⎫ ⎝⎛-=2612021lim x x x x 61-=e2. 解:x dx dy x22ln 2+= 2)2(ln 2222+=x dx y d3. 解:b ax x y ++='232,a x y 26+=''因为函数有拐点)1,1(-,所以⎩⎨⎧-==''1)1(0)1(y y ,即⎩⎨⎧-=+++=+11026c b a a 因为在0=x 处有极大值1,所以0)0(='y ,即0=b ,带入上式得⎪⎩⎪⎨⎧==-=103c b a4. 解:dx xe x⎰∞+-22|2e e +∞+∞==-=⎰5. 23323,3xy x y zy y x x z -=∂∂-=∂∂6.⎰⎰---=22111),(y y dxy x f dy7. 解:分离变量得xdx ydy cot tan -=两边积分得⎰⎰-=xdx ydy cot tan从而)sin arccos(x C y =8. 解:4586lim 221+-+-→x x x x x 12lim 1--=→x x x ∞=9. 解:dx xx dy x)5ln 551(254-=10. 解:x y 452--=',无驻点,y '不存在的点为45=x ,但]1,1[45-∉=x1)1(,3)1(==-y y所以最大值是3)1(=-y ,最小值是1)1(=y11. 解:dx xe x⎰∞+-22|2e e +∞+∞==-=⎰12. yx x z 332-=∂∂ ,x y y z 32--=∂∂13.⎰⎰=xxdyy x f dx 2),(1014. 解:分离变量得x dx y y dy sin ln =,两边积分得⎰⎰=x dxy y dy sin ln 两边积分得⎰⎰=x dx y y dy sin ln ,从而原方程的特解为2tan xe y =。
数学建模与数学实验第五版课后答案4
数学建模与数学实验第五版课后答案4.41、27.下列各函数中,奇函数的是()[单选题] *A. y=x^(-4)B. y=x^(-3)(正确答案)C .y=x^4D. y=x^(2/3)2、4.点(-3,-5)关于x 轴的对称点的坐标为()[单选题] *A(-3,5)(正确答案)B(-3,-5)C(3,5)D(3,-5)3、1.(必修1P5B1改编)若集合P={x∈N|x≤2 022},a=45,则( ) [单选题] * A.a∈PB.{a}∈PC.{a}?PD.a?P(正确答案)4、2.在+3,﹣4,﹣8,﹣,0,90中,分数共有()[单选题] *A.1个B.2个C.3个(正确答案)D.4个5、13.在海上,一座灯塔位于一艘船的北偏东40°方向,那么这艘船位于灯塔()[单选题] *A.南偏西50°方向B.南偏西40°方向(正确答案)C.北偏东50°方向D.北偏东40°方向6、5.已知集合A={x|x=3k+1,k∈Z},则下列表示不正确的是( ) [单选题] *A.-2∈AB.2 022?AC.3k2+1?A(正确答案)D.-35∈A7、45.下列运算正确的是()[单选题] *A.(5﹣m)(5+m)=m2﹣25B.(1﹣3m)(1+3m)=1﹣3m2C.(﹣4﹣3n)(﹣4+3n)=﹣9n2+16(正确答案)D.(2ab﹣n)(2ab+n)=4ab2﹣n28、10.如图是丁丁画的一张脸的示意图,如果用表示左眼,用表示右眼,那么嘴的位置可以表示成().[单选题] *A.(1,0)B(-1,0)(正确答案)C(-1,1)D(1,-1)9、6.方程x2=3x的根是()[单选题] *A、x = 3B、x = 0C、x1 =-3, x2 =0D、x1 =3, x2 = 0(正确答案)10、13.在数轴上,下列四个数中离原点最近的数是()[单选题] *A.﹣4(正确答案)B.3C.﹣2D.611、13.设x∈R,则“x3(x的立方)>8”是“|x|>2”的( ) [单选题] * A.充分而不必要条件(正确答案)B.必要而不充分条件C.充要条件D.既不充分也不必要条件12、35.若代数式x2﹣16x+k2是完全平方式,则k等于()[单选题] * A.6B.64C.±64D.±8(正确答案)13、f(x)=-2x+5在x=1处的函数值为()[单选题] *A、-3B、-4C、5D、3(正确答案)14、13.不等式x+3>5的解集为()[单选题] *A. x>1B. x>2(正确答案)C. x>3D. x>415、掷三枚硬币可出现种不同的结果()[单选题] *A、6B、7C、8(正确答案)D、2716、6、已知点A的坐标是,如果且,那么点A在()[单选题] *x轴上y轴上x轴上,但不能包括原点(正确答案)y轴上,但不能包括原点17、49.若(x+2)(x﹣3)=7,(x+2)2+(x﹣3)2的值为()[单选题] * A.11B.15C.39(正确答案)D.5318、1.计算| - 5 + 3|的结果是[单选题] *A. - 2B.2(正确答案)C. - 8D.819、4.一个数是25,另一个数比25的相反数大- 7,则这两个数的和为[单选题] *A.7B. - 7(正确答案)C.57D. - 5720、30.圆的方程+=4,则圆心到直线x-y-4=0的距离是()[单选题] *A.√2(正确答案)B.√2/2C.2√2D.221、用角度制表示为()[单选题] *30°(正确答案)60°120°-30°22、19.对于实数a、b、c,“a>b”是“ac2(c平方)>bc2(c平方) ; ”的()[单选题] * A.充分不必要条件B.必要不充分条件(正确答案)C.充要条件D.既不充分也不必要条件23、10.下列四个数中,属于负数的是().[单选题] *A-3(正确答案)B 3C πD 024、6.有15张大小、形状及背面完全相同的卡片,卡片正面分别画有正三角形、正方形、圆,从这15张卡片中任意抽取一张正面的图形既是轴对称图形,又是中心对称图形的概率是1/3?,则正面画有正三角形的卡片张数为()[单选题] *A.3B.5C.10(正确答案)D.1525、若m·23=2?,则m等于[单选题] *A. 2B. 4C. 6D. 8(正确答案)26、下列各式计算正确的是( ) [单选题] *A. (x3)3=x?B. a?·a?=a2?C. [(-x)3]3=(-x)?(正确答案)D. -(a2)?=a1?27、13.下列说法中,正确的为().[单选题] * A.一个数不是正数就是负数B. 0是最小的数C正数都比0大(正确答案)D. -a是负数28、函数式?的化简结果是()[单选题] *A.sinα-cosαB.±(sinα-cosα)(正确答案)C.sinα·cosαD.cosα-sinα29、6.下列说法正确的是().[单选题] * A.不属于任何象限的点不在坐标轴上就在原点B.横坐标为负数的点在第二、三象限C.横坐标和纵坐标互换后就表示另一个点D.纵坐标为负数的点一定在x轴下方(正确答案)30、下列函数是奇函数的是()[单选题] *A、f(x)=3x(正确答案)B、f(x)=4xC、f(x)= +2x-1D、f(x)=。
数学建模与数学实验ppt课件
02
通过数学实验,可以发现和解决数学理论中的问题,推动数学
理论的发展和完善。
数学实验在科学、工程、经济等领域有广泛应用,为解决实际
03
问题提供有效的工具和方法。
数学实验的常用工具
MATLAB
一种常用的数学计算软件,具有强大的数值 计算、矩阵运算和图形绘制等功能。
Python
一种通用编程语言,广泛用于科学计算、数 据分析和机器学习等领域。
02
03
相互促进
两者都是为了解决实际问题或探 究数学问题而进行的方法和工具。
数学建模为数学实验提供理论指 导,而数学实验可以验证数学建 模的正确性和有效性。
区别
目的
数学建模的主要目的是建立数学模型,描述实际问题中变 量之间的关系;而数学实验则是通过实验手段来探究数学 规律或验证数学结论。
应用领域
数学建模广泛应用于各个领域,如物理、工程、经济等; 而数学实验则更多应用于数学教育和研究领域。
简化模型
在保证模型精度的基础上,对模型进行必要 的简化。
求解模型
求解方法选择
根据模型的特点选择合适的数值计算方法或解 析解法。
编程实现
利用编程语言实现模型的求解过程。
误差分析和收敛性判断
对求解过程进行误差分析,判断求解方法的收敛性和稳定性。
模型验证与优化
数据拟合与检验
将模型结果与实际数据进行对比,检验模型的准确性和适用性。
问题分析
明确问题定义
对问题进行深入理解,明确问题的目标、约束条件和 相关参数。
收集数据和信息
收集与问题相关的数据和背景信息,为建立模型提供 依据。
确定主要影响因素
分析问题中起决定性作用的关键因素,忽略次要因素。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2014年下学期数学实验与数学建模作业习题6
1.求下列级数的和: (1) ∑∞=-1212n n n (2) ∑∞=+1)
12(1n n n (3) 112n n n n ∞=+⋅∑ (4) 313n n n ∞=∑
(5) ∑∞=12
sin n n x
【1】求解代码:
【2】运行结果: syms n x
s1=symsum((2*n-1)/2^n,n,1,inf)
s2=symsum(1/n/(2*n+1),n,1,inf)
s3=symsum((n+1)/n/2^n,n,1,inf)
s4=symsum(n^3/3^n,n,1,inf)
s5=symsum(sin(x)/n^2,n,1,inf)
2.求级数∑∞=021n n 的前n 项和与∑∞=121
n n 的级数和。
【1】求和代码:
【
2】运行结果: syms n x
s1=symsum(1/2^n,n,0,n-1)
s2=symsum(1/n^2,n,1,inf)
s=s1+s2
3.将函数sin x 展开为x 的幂级数,分别展开至5次和20次。
4.将函数(1)m x +展开为x 的幂级数,m 为任意常数。
展开至4次幂。
5.将函数21()53
f x x x =+-展开为(2)x -的幂级数。
6.将函数cos x 展开成()3
x π-
的幂级数,取前10项。
【1】syms n x m t pi
%3题 taylor(sin(x),x,6)
taylor(sin(x),x,21)
%4题
taylor((1+x)^m,x,5)
%5题
taylor(1/(x*x+5*x-3),x,6,2)
%新版本用taylor(1/(x*x+5*x-3),x,'ExpansionPoint',2)
%6题
taylor(cos(x),x,10,pi/3)
【2】运行结果:
3-(1)
3-(2)
4题
5题
6题
ans =
(3^(1/2)*(pi/3 - x)^5)/240 - (3^(1/2)*(pi/3 - x)^3)/12 - (3^(1/2)*(pi/3 - x)^7)/10080 + (3^(1/2)*(pi/3 - x)^9)/725760 - (pi/3 - x)^2/4 + (pi/3 - x)^4/48 - (pi/3 - x)^6/1440 + (pi/3 - x)^8/80640 + (3^(1/2)*(pi/3 - x))/2 + 1/2
7.求函数2()f x x =在[,]ππ-上的傅立叶级数。
8.求出函数32()f x x x =+在区间[,]ππ-上的前11个傅立叶系数,即n =5。
【1】建立afourier 函数和safourier 函数
function [a0,an,bn] = afourier(f,l)
%求解傅里叶级数系数,l 半周期
syms x n l a0
a0=int(f,x,-l,l)/l;
an=int(f*cos(n*pi*x/l),x,-l,l)/l;
bn=int(f*sin(n*pi*x/l),x,-l,l)/l;
function S = safourier(f,a0,an,bn,l,m)
%求解傅里叶级数,l为半周期,m为展开到n=m
syms x n l pi
sn=an*cos((n*pi*x)/l)+bn*sin((n*pi*x)/l);
ssn=symsum(sn,n,1,m);
S=a0/2+ssn;
End
调用命令:syms x a0 m n
[a0,an,bn]=afourier(x*x,pi)
safourier(x*x,a0,an,bn,pi,5)
【2】运行结果:
a0 =(2*l^2)/3
an =(2*(pi^2*l^3*n^2*sin(pi*n) - 2*l^3*sin(pi*n) + 2*pi*l^3*n*cos(pi*n)))/(pi^3*l*n^3)
bn =0
ans =
l^2/3 - (4*l^2*cos((pi*x)/l))/pi^2 + (l^2*cos((2*pi*x)/l))/pi^2 - (4*l^2*cos((3*pi*x)/l))/(9*pi^2) + (l^2*cos((4*pi*x)/l))/(4*pi^2) - (4*l^2*cos((5*pi*x)/l))/(25*pi^2)。