资源描述:
《matlab语言概述cha》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、矩阵、M文件的运用、逻辑运算、多项式运算例分别建立3×3、3×2和与矩阵A同样大小的零矩阵。(1)建立一个3×3零矩阵。zeros(3)(2)建立一个3×2零矩阵。zeros(3,2)(3)设A为2×3矩阵,则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。A=[123;456];%产生一个2×3阶矩阵Azeros(size(A))%产生一个与矩阵A同样大小的零矩阵ROT90(A,K)istheK*90degreerotationofA,K=+-1,+-2,...对角矩阵符号X=diag(v,k):whenvisavectorofncomponents,retu
2、rnsasquarematrixXofordern+abs(k),withtheelementsofvonthekthdiagonal.k=0representsthemaindiagonal,k>0abovethemaindiagonal,andk<0belowthemaindiagonal.K=0K>0K<0K=0K=0K=0K>0K=0K<0K>0K=0例:v=[123];x=diag(v,1);y=diag(v);z=diag(v,-1);v=diag(X,k)formatrixX,returnsacolumnvectorvformedfromtheelementsofthe
3、kthdiagonalofX.x=[123;456;789];v=diag(x,1);w=diag(x,0);u=diag(x,-1);x=[12310;45611;78912];v=diag(x,1);w=diag(x,0);u=diag(x,-1);例先建立5×5矩阵A,然后将A的第一行元素乘以1,第二行乘以2,…,第五行乘以5。A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19];D=diag(1:5);D*A%用D左乘A,对A的每行乘以一个指定常数三角阵三角阵又进一步分为上三角阵和下三角
4、阵,所谓上三角阵,即矩阵的对角线以下的元素全为0的一种矩阵,而下三角阵则是对角线以上的元素全为0的一种矩阵。(1)上三角矩阵求矩阵A的上三角阵的MATLAB函数是triu(A)。triu(A)函数也有另一种形式triu(A,k),其功能是求矩阵A的第k条对角线以上的元素。例如,提取矩阵A的第2条对角线以上的元素,形成新的矩阵B。(2)下三角矩阵在MATLAB中,提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k),其用法与提取上三角矩阵的函数triu(A)和triu(A,k)完全相同。用于专门学科的特殊矩阵(1)魔方矩阵魔方矩阵有一个有趣的性质,其每行、每列及两条对
5、角线上的元素和都相等。对于n阶魔方阵,其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n),其功能是生成一个n阶魔方阵。例将101~125等25个数填入一个5行5列的表格中,使其每行每列及对角线的和均为565。M=100+magic(5)范得蒙矩阵范得蒙(Vandermonde)矩阵最后一列全为1,倒数第二列为一个指定的向量,其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中,函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如,A=vander([1;2;3;5])即可得到上述范得蒙矩
6、阵。伴随矩阵MATLAB生成伴随矩阵的函数是compan(p),其中p是一个多项式的系数向量,高次幂系数排在前,低次幂排在后。例如,为了求多项式的x3-7x+6的伴随矩阵,可使用命令:p=[1,0,-7,6];compan(p)二、在MATLAB中.m文件分为两种:函数(Function),函数接受用户的输入参数,执行函数的具体运算,最后输出结果。1.脚本文件:包括注解和指令两部分。注解部分必须以“%”开头。脚本文件使用与存放的变量和数据都在工作空间中,指向工作空间来操作。例:假设当前目录下有一个命令M文件:%solver.m%usedtosolveA*x=b%whereA=[-1.
7、512;3-11;-135],b=[2.5;5;8].A=[-1.512;3-11;-135];b=[2.5;5;8];x=Ab在命令窗口中执行solver命令,即可得到方程组的解。x=0.7500-0.62502.1250在命令窗口中键入typesolver.m即可在命令窗口中看到该文件。2.函数文件:该文件主要包含以下五部分(1)函数定义行:位于函数的第一行,格式为:function[A1,A2,A3,…]=函数名(B1,B2,B3,…)A1,A2,