




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、符号表达式及其应用符号表达式及其应用微积分基本运算微积分基本运算taylor级数展开级数展开常微分方程求解常微分方程求解第第4章章 高等数学计算高等数学计算matlab的的符号计算符号计算符号计算相对于数值计算而存在符号计算相对于数值计算而存在,数值计算的特点是数值计算的特点是“近似近似”, 符号计算的特点是符号计算的特点是“准确准确”. 早期的早期的matlab不能实现符号计算不能实现符号计算,1993年年mathworks公司购买软件公司购买软件maple的使用权的使用权,并开发实现符号计算的工具箱并开发实现符号计算的工具箱symbolic math toolbox. matlab符号计算
2、能实现微积分计算、线性代数计算、符号计算能实现微积分计算、线性代数计算、方程求解、积分变换等。在符号计算中方程求解、积分变换等。在符号计算中, 符号表达式是符号表达式是主要操作对象主要操作对象.符号表达式符号表达式符号变量、运算符、函数、数字组成符号变量、运算符、函数、数字组成表达式中的符号变量必须事先创建表达式中的符号变量必须事先创建.创建的常方法是用创建的常方法是用关键词关键词syms 空格后跟空格后跟符号变量符号变量。2/16例例1. 用符号表达式定义用符号表达式定义 f = e 0.2x sin (0.5x)并绘图并绘图.syms x ; f = exp(-0.2*x)*sin(0.5
3、*x); ezplot(f,0,2*pi)3/16syms 符号变量符号变量1 符号变量符号变量2 f = exp(-1/5*x) * sin(1/2*x)ezplot( f )绘表达式绘表达式f=f(x)的图的图绘图区域绘图区域 -2*pix2*piezplot(f, a,b)绘图区域绘图区域:axa时的极限时的极限symsum(f,k,m,n) 求级数求级数7/16 nmkkf)(例例5.绘函数绘函数 的图形的图形,并分别计算并分别计算 一阶和二阶导数一阶和二阶导数syms xf=3*x2+6*x-1;g=x2+x-3;r=f/g;ezplot(r,-10,10)r1=diff(r,x);
4、simplify(r1)f1,g1=numden(r1)r2=diff(r,x,2);simplify(r2)f2,g2=numden(r2)8/163163)(22 xxxxxr222)3()17163()( xxxxxr3223)3()4151243(2)( xxxxxxrx=double(solve(f1,x);x = -1.4648 -3.8685例例6. 求例求例5中函数的极值点与拐点中函数的极值点与拐点x1=x(1);x2=x(2);y1=subs(r,x,x1);y2=subs(r,x,x2); ezplot(r,-10,10),hold onplot(x1,y1,ro,x2,y
5、2,ro)9/16极值点分布极值点分布xx=double(solve(f2)xx = -5.2635 -1.3682 - 0.8511i -1.3682 + 0.8511ix0=xx(1)y0=subs(r,x,x0)ezplot(r,-10,10),hold onplot(x0,y0,ro)求例求例5 5中函数的拐点中函数的拐点10/16拐点分布拐点分布例例7. 计算不定积分计算不定积分 dxbxeax)sin(syms a b xf=exp(a*x)*sin(b*x); int(f,x);pretty(ans)11/16 b exp(a x) cos(b x) a exp(a x) sin
6、(b x) - - + - 2 2 2 2 a + b a + bg=simplify(ans) g = exp(a*x)*(-b*cos(b*x)+a*sin(b*x)/(a2+b2)s2=symsum(k2,k,1,n); factor(s2) ans = 1/6*n*(n+1)*(2*n+1) nkks11例例8.计算级数计算级数 nkks122syms k ns1=symsum(k,k,1,n);factor(s1) ans = 1/2*n*(n+1)s3=symsum(1/k2,k,1,inf) s3= 1/6*pi2 1231kks12/16taylor(f,n,x) n-1次麦克
7、劳林多项式展开次麦克劳林多项式展开taylor(f,n,x,a) a点的点的n-1次泰勒多项式展开次泰勒多项式展开.)cos(451)(xxf 例例9. 麦克劳林展开函数麦克劳林展开函数syms xf=1/(5+4*cos(x)t=taylor(f,8)pretty(t) 2 4 49 6 1/9 + 2/81 x + 5/1458 x + - x 13122013/16命令格式命令格式:dsolve(eq1,con1,x)y的一阶导数的一阶导数 dy, y的二阶导数的二阶导数 d2y例例10.解微分方程解微分方程14/16y = dsolve(dy=1/(1+x2)-2*y2,y(0) =
8、0,x)y = 2*x/(2*x2+2)22211yxy 符号解符号解: y(x)= x / (1 + x 2) 练习题与思考题练习题与思考题15/161.用用matlab求极限求极限xxxx30sinsintanlim 2.用用matlab求导数求导数xxxfcos1sin1)( 3.用用matlab求定积分求定积分 022dxexx4.用用matlab求级数和求级数和 12)12(1nnn5.计算曲线段计算曲线段 f(x)=exp(a x)sin(b x), 绕绕x轴旋转的旋转曲面面积轴旋转的旋转曲面面积 20 x6. 编写一个函数文件,用于计算椭圆周长,要求输编写一个函数文件,用于计算椭圆周长,要求输入变
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年江西中考数学试卷真题答案解读及备考指导
- 高中化学课程中校园植物化学成分分析与应用研究论文
- 小学生网络互动游戏对认知发展影响分析论文
- 高中语文课程思政教育中的历史教育价值挖掘与传承论文
- 中国医药级酮咯酸氨丁三醇行业市场前景预测及投资价值评估分析报告
- 节电方案与管理制度
- 英文版公司管理制度
- 电工学试题集和试题集及答案
- 小学语文《夜色》课件
- 财务管理学自考历年真题
- NB/T 11123-2023煤矿安全双重预防机制规范
- 内墙抹灰安全技术交底
- 中学美术校本教材《素描》
- 上海市2023年初中毕业物理课程终结性评价指南
- 中国教师发展基金会国家教师科研基金规划全国重点课题实施指南
- (通桥【2018】8370)《铁路桥梁快速更换型伸缩缝安装图》
- “双减”背景下的初中数学课堂教学设计与思考 论文
- 义务教育语文课程标准(2022)测试题带答案(20套)
- (医院护理论文)护理本科生学习适应现状及影响分析
- 保护性约束完整版
- 明源地产erp3.04-费用管理操作手册
评论
0/150
提交评论