


下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验二MATLAB的符号运算符号计算的特点:一,运算以推理解析的方式进行,因此不受计 算误差积累问题困扰;二,符号计算,或给出完全正确的封闭解,或 给出任意精度的数值解当封闭解不存在时;三,符号计算指令的 调用比拟简单;四,计算所需时间较长。在 MATLAB,符号计算虽以 数值计算的补充身份出现,但涉及符号计算的指令使用、运算符操 作、计算结果可视化、程序编制以及在线帮助系统都是十分完整、便 捷的。一、实验目的1、熟悉符号对象和符号表达式;2、掌握符号表达式和符号函数的操作;3、学习利用符号代数方程的求解 。二、实验根本知识1、根本运算指令;2. 掌握符号运算的应用 三、上机练习1. 1通过符
2、号计算求y(t) sint的导数dy。 2然后根据此结dt果,求兽和等。dt 11 0dt t 2>> clear allsyms ty=abs(si n( t)df=diff(y)df1=limit(df,t,0,'left')df2=subs(df,'t',sym(pi/2) y =abs(si n( t) df =sig n( si n( t)*cos(t)df1 =-1 df2 =02 x2222.计算二重积分1 1 (x y )dydx>>clear allsyms x yf=xA2+yA2Fin t=(i nt(i nt(f,
3、y,1,xA2),x,1,2) f = xA2 + yA2Fint =1006/1053. f x2y2(1)求今;求11 fdx. dy1>> clear all syms x y f=xA2+yA2; dfdy2=diff(f,y,2) dfdy2 =>> in tf=i nt(f,x,-1,1) intf =2*yA2 + 2/34. 解方程:3x24x 60>> clear allsyms xeq=solve(3*xA2+4*x-6)eq =-2/3+1/3*22八(1/2)-2/3-1/3*22八(1/2)5. 方程:d3y o dyr 3 2y
4、6dx dx(1)求通解;(2)求满足y (0) 1,y (0) 2,y(0) 3的特解.>>clear all>> y=dsolve('D3y+3*Dy+2*y-6=0','x')y =3+C1*exp(1+2A(1/2)F(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+C2*exp(- 1/2*(1+2A(1/2)A(1/3)*(-1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*si 门(1/2*3八(1/2)*(1+2A(1
5、/2)A(1/3)*(1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+C3*exp(- 1/2*(1+2A(1/2)A(1/3)*(-1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*cos(1/2*3A(1/2)*(1+2A(1/2)A(1/3)*(1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)>> y二dsolve('D3y+3*Dy+2*y-6=0','y(0)=3,D2y(0)=1,Dy(0)=2','
6、;x')3+(-1/12*(1+2八(1/2)八(1/3)*2八(1/2)+1/6*(1+2八(1/2)八(2/3)+1/3*(1+2 八(1/2)八(1/3)-1/4*(1+2A(1/2)A(2/3)*2A(1/2)*exp(1+2A(1/2)A(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+(1/3*3A(1/2)* (1+2A(1/2)A(1/3)-1/12*3A(1/2)*(1+2A(1/2)A(1/3)*2A(1/2)- 1/6*3A(1/2)*(1+2A(1/2)A(2/3)+1/4*3A(1/2)*(1+2A
7、(1/2)A(2/3)*2 A(1/2)*exp(-1/2*(1+2A(1/2)F(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*sin(1/2*3A(1/ 2)*(1+2A(1/2)A(1/3)*(1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+(1/12*(1+2A(1 /2)A(1/3)*2A(1/2)-1/6*(1+2A(1/2)A(2/3)- 1/3*(1+2A(1/2)A(1/3)+1/4*(1+2A(1/2)A(2/3)*2A(1/2)*exp(- 1/2*(1+2A(
8、1/2)A(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*cos(1/2*3A(1/ 2)*(1+2A(1/2)A(1/3)*(1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)226. 求方程x2 y2 1, xy 2 的解。>>clear all>> syms x y>> S=solve('xA2+yA2-1=0','x*y-2=0','x','y')S =x: 4x1 symy: 4
9、x1 sym>> disp('S.x'),disp(S.x),disp('S.y'),disp(S.y)S.x-1/2*(1/2*5A(1/2)+1/2*i*3A(1/2)A3+1/4*5A(1/2)+1/4*i*3A(1/2) -1/2*(1/2*5A(1/2)-1/2*i*3A(1/2)A3+1/4*5A(1/2)-1/4*i*3A(1/2) -1/2*(-1/2*5A(1/2)+1/2*i*3A(1/2)A3-1/4*5A(1/2)+1/4*i*3A(1/2) -1/2*(-1/2*5A(1/2)-1/2*i*3A(1/2)A3-1/4*5A(1/2)-1/4*i*3A(1/2)S.y1/2*5八(1/2)+1/2*i*3八(1/2)1/2*5八(1/2)-1/2*i*3八(1/2)-1/2*5A(1/2)+1/2*i*3A(1/2)-1/2*5A(1/2)-1/2*i*3A(1/2)实验小结:本次试验中,是为了让我们更好的学习MATLAB符号运算,让我们熟练掌握符号计算的解题理念和模式,经过本次实验,能够使学生 感受到MATLA优缺点,对以后的解题很有帮助。第一题说明了极限和导数的符号运算,学会应用 diff 的不同情况 和 limit 的求取极限作用。第二题是对积分
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 新解读《HG-T 3089-2001燃油用O形橡胶密封圈材料》新解读
- 级配碎石底基层施工方案
- 专题梳理 动词不规则变化自测表 课件
- 低能负离子在锥形玻璃管中的传输行为研究
- 合作股份公司管理制度
- 物理中考一轮复习教案 第二十九讲 电功、电功率
- 仓储店线下活动策划方案
- 仓库拍卖活动方案
- 仓鼠食品活动方案
- 代理记账推广活动方案
- 【MOOC】政府审计学-南京审计大学 中国大学慕课MOOC答案
- 《基督教概论》课件
- 虚拟现实技术导论 习题答案或解题思路 梁晓辉
- 计算机应用技术专业调研报告(高职)
- 2024NEA水性气硅涂膏隔热保温墙体构造
- 山西省太原市(2024年-2025年小学四年级语文)部编版期末考试((上下)学期)试卷及答案
- BPC10完整版本.0技术培训V1.0
- 2024年新高考II卷高考历史试卷(真题+答案)
- 2024年黑龙江医疗卫生事业单位招聘(药学)备考试题库(含答案)
- 2024年新高考1卷数学真题试卷及答案
- 湖北省武汉市洪山区2023-2024学年七年级下学期期末考试语文试卷
评论
0/150
提交评论