- 1、本文档共21页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MatLab4研讨
MatLab 数学建模
第四讲 数值计算
符号数学工具箱
符号表达式的运算 numeric 符号到数值的转换 pretty 显示悦目的符号输出 subs 替代子表达式 sym 建立符号矩阵或表达式 symadd 符号加法 symdiv 符号除法 symmul 符号乘法 symop 符号运算 sympow 符号表达式的幂运算 symrat 有理近似 symsub 符号减法 symvar 求符号变量
符号表达式的简化 collect 合并同类项 expand 展开 factor 因式 simple 求解最简形式 simplify 简化 symsum 和级数
符号多项式 charpoly 特征多项式 horner 嵌套多项式表示 numden 分子或分母的提取 poly2sym 多项式向量到符号的转换 sym2poly 符号到多项式向量的转换
符号微积分 diff 微分 int 积分 jordan 约当标准形 taylor 泰勒级数展开
符号可变精度算术 digits 设置可变精度 vpa 可变精度计算
求解符号方程 compose 函数的复合 dsolve 微分方程的求解 finverse 函数逆 linsolve 齐次线性方程组的求解 solve 代数方程的求解
符号线性代数 charploy 特征多项式 determ 矩阵行列式的值 eigensys 特征值和特征向量 inverse 矩阵逆 jordan 约当标准形 linsolve 齐次线性方程组的解 transpose 矩阵的转置
方程求解
求解单个代数方程
MATLAB具有求解符号表达式的工具,如果表达式不是一个方程式(不含等号),则在求解之前函数solve将表达式置成等于0。
solve( a*x^2+b*x+c ) % solve for the roots of the eqution
ans=
[1/2/a*(-b+(b^2-4*a*c)^1/2)]
[1/2/a*(-b-(b^2-4*a*c)^1/2)]
结果是符号向量,其元素是方程的2个解。如果想对非缺省x变量求解,solve必须指定变量。
solve( a*x^2+b*x+c , b ) % solve for b
ans=
-(a*x^2+c)/x
带有等号的符号方程也可以求解。
f=solve( cos(x)=sin(x) ) % solve for x
f=
1/4*pi
t=solve( tan(2*x)=sin(x) )
t=
[ 0]
[acos(1/2+1/2*3^(1/2))]
[acos(1/2=1/2*3^(1/2))]
并得到数值解。
numeric(f)
ans=
0.7854
numeric(t)
ans=
0
0 + 0.8314i
1.9455
注意在求解周期函数方程时,有无穷多的解。在这种情况下,solve对解的搜索范围限制在接近于零的有限范围,并返回非唯一的解的子集。
如果不能求得符号解,就计算可变精度解。
x=solve( exp(x)=tan(x) )
x=
1.306326940423079
代数方程组求解
可以同时求解若干代数方程,语句solve(s1,s2,.....,sn)对缺省变量求解n个方程,语句solve(s1,s2,...,sn, v1,v2,...,vn )对n个 v1,v2,...vn 的未知数求解n个方程。
solve(f) 解符号方程式f。
solve(f1,…,fn) 解由f1,…,fn组成的联立方程式。
我们先定义以下的方程式:
eq1 = x-3=4; % 注意也可写成eq1=x-7
eq2 = x*2-x-6=0; % 注意也可写成eq2=x*2-x-6
eq3 = x2+2*x+4=0;
eq4 = 3*x+2*y-z=10;
eq5 = -x+3*y+2*z=5;
eq6 = x-y-z=-1;
solve(eq1)
ans=
7
solve(eq2)
ans=
[[3],[-2]] % 原方程式有二个根3, -2
solve(eq3)
ans=
[[-1+*3^(1/2)],[-1-i*3^(1/2)]] % 注意实根和虚根的表示式
solve(eq4,eq5,eq6) % 解三个联立方程式
ans=
您可能关注的文档
- chapter_02__化学电源研讨.ppt
- chp3.2复杂电力系统潮流计算机计算11v研讨.ppt
- CIP应用研讨.ppt
- CO2气体保护焊及其工艺问题浅谈研讨.doc
- CPU和散热器1研讨.ppt
- CPU的结构和功能研讨.ppt
- CTP印刷中常见故障分析及处理研讨.ppt
- ct原理及应用研讨.pptx
- CT成像-基本概念研讨.doc
- C语言-地址与指针研讨.ppt
- 原电池电动势的测定实验报告.pdf
- 与业主、设计、总包、监理和他承包人的配合措施.pdf
- 公司管理流程.pptx
- 2024_2025学年新教材高中地理第1章地球的运动素养综合训练新人教版选择性必修1.doc
- 2024_2025学年新教材高中地理第3章大气的运动第1节常见天气系统第1课时锋与天气分层作业新人教版选择性必修1.doc
- 2024_2025学年新教材高中地理第1章地球的运动第2节地球运动的地理意义第4课时正午太阳高度的变化四季更替和五带划分分层作业课件新人教版选择性必修1.pptx
- 2024_2025学年新教材高中地理第2章地表形态的塑造第2节构造地貌的形成第1课时地质构造与地貌课件新人教版选择性必修1.pptx
- 2024_2025学年新教材高中地理第1章地球的运动问题研究人类是否需要人造月亮课件新人教版选择性必修1.pptx
- 五片小雪花课件.pdf
- 2024_2025学年新教材高中地理第3章大气的运动第2节气压带和风带第1课时气压带和风带的形成分层作业课件新人教版选择性必修1.pptx
文档评论(0)