原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
材料工程
计算机技术应用;第七章、Matlab在材料工程计算中的应用;7.1Matlab语言介绍;7.1.1Matlab语言概述;⑴编程简单使用方便
Matlab语言矩阵和向量运算的基本数据单元既不需要指定维数，也不需要说明数据类型，并且数学表达形式和运算规则与书写的习惯相近。
例7.1.下列语句串实现矩阵乘矩阵的运算：
A＝［1 3 0；2 1 1；4 2 5］；
B＝［6 7 1；1 3 2；1 0 1］；
C＝A*B
其中A和B是3阶方阵，运算的结果是
C＝9 16 7
14 17 5
31 34 13 ;; ; ;7.1.3环境窗口、语言结构和编程方法;;;;对于例7.2，这里也采用两种方法实现其中的作图功能。 方法1：首先建立以的fig．m为名的脚本文件，内容为例7.2中的语句串。然后在命令窗口输人命令fig，即可绘出函数的图形。
方法2：首先建立三个函数文件:
function y= fun1（x） function y= fun2（x）
y=sin2x； y=sinx2；
function y= fun3（x）
y=sinx2；
并建立一个以fig.m为名的脚本文件，其内容为语句串
X＝0：0.0l：6； y3＝fun3（x）；
y1＝fun1（x）； plot（x，y1，x，y2,x,y3）
y2＝fun2（x）；
然后在命令窗口输入命令fig，即可绘出函数的图形。 ;5.1.4主要语法和符号;;⑷常用数学函数;;;;;;;;;;;;;分割图形显示窗口：;;;;7.2数据分析;7.2.1统计数字特征;;;7.2.2剔除异常值; ;;;5.2.3方差分析;;;;;;;5.2.4回归分析;;;;;;7.3数值计算;5.3.1方程求根;;例7.16求方程 x3-x2-4x+4=0在0点附近的一个解和小于0的一个解。
第一步建立方程的函数文件
function y=f1(x)
y=x^3－x^2－4*x＋4;
第二步求0点附近的解，将初值取为0。
x0=0;
x=fzero(‘f1’,x0)
运行程序后得到
x=
1.0000
为了求得小于0的解，取x0为[-10，0]。
x0=[-10，0]; 运行程序后得到
x=fzero(‘f1’,x0) x=
-2
;7.3.2线性方程组求解;例7.17解方程组
解：
A=[1 2 1; 2 –1 1; 1 –1 2]; b=[0; 1; 3]
x=A \ b
运行程序得到
x=
-0.5000
-0.5000
1.5000
;例7.18解方程组
这是一个不定方程，有无穷多组解。
解：
A=[2 1 –1 1; 3 -2 1 -1; 1 4 –3 5]; b=[1; 4; -2]
x=A \ b
运行程序给出警告和一组解：
Warning: Rank deficient, rank=2 tol=5.2545e-015
x= 0.7778
0
0
-0.5556 ;7.3.3非线性方程组求解;例7.19 解方程组
解：(把方程右端非零项移到左端)
f=inline（’[sin(x(1))+(x(2))^2+log(x(3))-7;
3*x(1)+2^(x(2))-x(3)^3+1;
x(1)+x(2)+x
文档评论（0）