基于MATLAB的数值分析.pptVIP

下载本文档

33
0
约1.48万字
约 55页
2017-02-17 发布于重庆
举报
版权申诉

基于MATLAB的数值分析.ppt

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的数值分析

* A=[1 2 3 4; 5 6 7 8; 9 10 11 12] B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6] * A=[1 0 1 1 2;0 1 -1 2 3;3 0 5 1 0;2 3 1 2 1], vr=[1,3];vc=[1,3];A1=A(vr,vc) * A11=A(1:2,1:2),A12=A(1:2,3:5), A21=A(3:4,1:2),A22=A(3:4,3:5) B=[A11 A12;A21 A22] * A=[1 0 1 1 2;0 1 -1 2 3;3 0 5 1 0;2 3 1 2 1]; A(1,:)=[0 0 0 0 0]; A --------------------------------------- A([1,3],:)=[ ] * y1=rand(1,5), y2=rand(1,5), rand(seed,3), x1=rand(1,5), rand(seed,3), x2=rand(1,5) * C=A1+B1 D=A1-B1 syms c, cA=c*A1 A2=A1(:,1:3), B1 G=A2*B1 * rand(seed,1);A=rand(3,6), Asort=sort(A), Amax=max(A), Asum=sum(A) 4.矩阵的修改和提取修改矩阵A，将它的第1行变为0。程序： A=[1 0 1 1 2;0 1 -1 2 3; 3 0 5 1 0;2 3 1 2 1] A(1,:)=[0 0 0 0 0]; A 删除上面矩阵A的第1、3行。程序： A([1,3],:)=[ ] 结果 5.生成特殊矩阵全1阵 ones(n), ones(m,n), ones(size(A)) 全零阵： zeros(n) ，zeros(m,n), zeros(size(A)) 常常用于对某个矩阵或向量赋0初值单位阵： eye(n)，eye(m,n) 随机阵： rand(m,n)， rand(n)=rand(n,n)用于随机模拟，常和rand(seed,k)配合使用。将rand指令运行多次，观察结果。程序： y1=rand(1,5), y2=rand(1,5), rand(seed,3), x1=rand(1,5), rand(seed,3), x2=rand(1,5) 结果 a=-4:4 A=reshape(a,3,3) a = -4 -3 -2 -1 0 1 4 A = -4 -1 2 -3 0 3 -2 1 4 a1=diag(A,1) a1 = -1 3 A1=diag(a1,-1) A1 = 0 0 0 -1 0 0 0 3 0 A = -4 -1 2 -3 0 3 -2 1 4 flipud(A) ans = -2 1 4 -3 0 3 -4 -1 2 fliplr(A) ans = 2 -1 -4 3 0 -3 4 1 -2 rot90(A) ans = 2 3 4 -1 0 1 -4 -3 -2 tril(A) ans = -4 0 0 -3 0 0 -2 1 4 triu(A) ans= 0 -1 2 0 0 3 0 0 0 6.矩阵的数组运算 + - .* ./ .\ .^ C=A+B, C=A-B C=A.*B C=A./B或B.\A C=A.^p, C=p.^A C=s+B C=s*A C=s./B 或 B.\s C=A./B ? for i=1:3 for j=1:3 c(i,j)=a(i,j)/b(i,j); end end sin(A) l