- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验3 方程和方程组的求解
一、实验目的
熟悉MATLAB中关于求解方程和方程组的各种命令,掌握利用MATLAB进行线性方程组、非线性方程、非线性方程组的求解。
二、相关知识
在MATLAB中,由函数solve()、null()、fsolve(),fzero等来解决线性方程(组)和非线性方程(组)的求解问题,其具体格式如下:
X=solve(‘eqn1’,’eqn2’,…,’eqnN’,’var1’,’var2’,…,’varN’)
X=fsolve(fun,x0,options)
函数solve用来解符号方程、方程组,以及超越方程,如三角函数方程等非线性方程。参数’eqnN’为方程组中的第N个方程,’varN’则是第N个变量。
函数null(A)则用来解线性方程组AX=O的基础解系,实际是求系数矩阵A的零空间,在null函数中可加入参数’r’,表示有理基。通过求系数矩阵的秩和增广矩阵的秩,可以判定方程组是否有解,以及是否需要求基础解系。
另外,还可以用函数fzero求解非线性方程。用法与fsolve类似,请大家自己查看帮助系统。
例1:求解方程 的MATLAB程序为:X=solve(‘x^2-x-6=0’,’x’)
结果为:X=3, -2
例2:求解方程组 的程序为:
[X,Y]=solve(x^2+y-6=0,y^2+x-6=0,x,y)
结果为:X =2, -3, 1/2-1/2*21^(1/2), 1/2+1/2*21^(1/2)
Y =2, -3, 1/2+1/2*21^(1/2), 1/2-1/2*21^(1/2)
例3:求解方程组 的程序为:
clear
format rat
A=[5, 0, 4, 2;1, -1, 2, 1;4, 1, 2,0;1,1,1,1]
B=[3;1;1;0]
X=A\B
结果请大家自己运行。
例4:求方程组 的通解的程序为:
clear
format rat
A=[1,2,2,1;2,1,-2,-2;1,-1,-4,-3]
C=null(A,r) %求出矩阵A的解空间的有理基。
结果如下:C =
2 5/3
-2 -4/3
1 0
0 1
接着,用命令:
syms k1 k2
X=k1*C(:,1)+k2*C(:,2)
求出的通解为:X =
[ 2*k1+5/3*k2]
[ -2*k1-4/3*k2]
[ k1]
[ k2]
例5:求方程组 的通解的程序为:
clear
format rat
A=sym([1,2,2,1;2,1,-2,-2;1,-1,-4,-3])
b=sym([1;2;2])
B=[A,b]
n=length(A(1,:))
RA=rank(eval(A))
RB=rank(eval(B))
if(RA==RBRA==n)
X=eval(A\B) %在方程组满秩时,求出唯一解
elseif(RA==RBRAn)
C=eval(A\b) %在方程组不满秩时,求出特解
D=null(eval(A),‘r’) %求出矩阵A的零空间的基,即方程组的基础解系
syms k1 k2
X=k1*D(:,1)+k2*D(:,2)+C %求出方程组的全部解
else
fprintf(No Solution for the Equations)
end
结果请大家自己运行。
非线性方程组的求解——函数fsolve
例6:求解非线性方程组
先建立存放函数的m文件,文件名必须与函数名一致,这里为sy6_6.m,内容如下:
function y=sy6_6(x)
y(1)=x(1)-0.5*sin(x(1))-0.3*cos(x(2))
y(2)=x(2)-0.5*cos(x(1))+0.3*sin(x(2))
接着,我们建立另一个m文件sy6_6_1.m,其内容为:
clear
format short
x0=[0.1,0.1]
fsolve(@sy6_6,x0,optimset(fsolve))
%这里的optimse
您可能关注的文档
- 【2017年整理】实务复习要点(3至5章).doc
- 【2017年整理】实务概念.doc
- 【2017年整理】实心球教学设计.doc
- 【2017年整理】实施例3 制造陶瓷材料整体式安全插座.doc
- 【2017年整理】实施型钢拱架工厂化加工 提升隧道施工安全质量管控能力(中铁22局).doc
- 【2017年整理】实施强制性产品认证的产品目录.doc
- 【2017年整理】实施性施组(6-13).doc
- 【2017年整理】实施细则(土建).doc
- 【2017年整理】实时数据库开发注意事项.doc
- 【2017年整理】实现一个简单的linux字符设备驱动.doc
- 广东省深圳市龙岗区启英学校2023-2024学年七年级上学期期中数学试题.docx
- 广东省深圳市龙岗区宏扬学校2023-2024学年九年级上学期期中数学试题.docx
- 广东省深圳市南山区桃源中学2022-2023学年九年级数学上学期第一次月考(21.1—23.3)数学测试题 .docx
- 广东省深圳市南山区深圳南实集团麒麟中学2022-2023学年九年级上学期期末数学试题.docx
- 广东省阳江市阳春市2023-2024学年九年级上学期期末数学试题.docx
- 广东省深圳市深圳实验学校初中部2022-2023学年九年级上学期月考数学试题.docx
- 广东省深圳市新华中学2023-2024学年九年级上学期月考数学试题.docx
- 广东省深圳市南山外国语集团和南山第二外国语集团联考2023-2024学年七年级上学期期中数学试题.docx
- 广东省深圳市龙华区高峰学校2022-2023学年七年级上学期期中考试数学试卷.docx
- 广东省深圳市南山区2022-2023学年九年级上学期期末考试数学试卷.docx
文档评论(0)