导航:首页 > 编程知识 > 如何用matlab编程矩阵的初等变换

如何用matlab编程矩阵的初等变换

发布时间:2023-01-16 22:00:07

1. 用Matlab怎么将一个矩阵化为行最简矩阵

对矩阵实行初等行变换,化为行最简形矩阵,其特点是:非零行的第一个非零元为1,且这些非零元所在的列的其他元素都为0
x=[2 -1 -1 1 2;1 1 -2 1 4;4 -6 2 -2 4;3 6 -9 7 9]
object=x;
[m,n]=size(object);
for i=1:m
temp(1:m,1)=object(1:m,i);
if( all(abs(temp(i:m,1))<eps))
continue;
end
N_zero=find(temp);
clear temp;
[p,q]=size(N_zero);
clear Molus;
if( all(abs(object(i,1:n))<eps))
continue;
end

z=1;
for j=1:p
T_M=object(N_zero(j,1),i)/object(N_zero(i,1),i);
if z==1
Molus=T_M;
z=z+1;
continue;
end
Molus=[Molus,T_M];
z=z+1;
end
%初等行变换:把某一行所有元素的k倍加到另一行对应的元素上去(第j行的k倍加到第i行上,记做ri+k*rj)
for j=1:p
if j==i
continue;
end
object(N_zero(j,1),1:n)=object(N_zero(j,1),1:n)-object(N_zero(i,1),1:n)*Molus(j);
end
%初等行变换:对调两行(对调i,j两行,记做ri<->rj)
temp(1:m,1)=object(1:m,i);
N_zero=find(temp);
clear temp;
if N_zero(1)==i
object(i,1:n)=object(i,1:n)/object(i,i);
continue;
end
temp(1,1:n)=object(i,1:n);
object(i,1:n)=object(N_zero(1,1),1:n);
object(N_zero(1,1),1:n)=temp(1,1:n);
%初等行变换:以数k~=0乘某一行中的所有元素(第i行乘k,记作ri*k);
object(i,1:n)=object(i,1:n)/object(i,i);
end
L=abs(object)<eps;
object(L)=0;
a=object;

2. 如何用matlab编写矩阵运算程序

一、矩阵的表示方法

1、矩阵元素必须在”[]”内;

2、矩阵的同行元素之间用空格(或”,”)隔开;

3、矩阵的行与行之间用”;”(或回车符)隔开;

4、矩阵的元素可以是数值、变量、表达式或函数;

5、矩阵的尺寸不必预先定义。

二、矩阵的运算

1、算术运算

MATLAB的基本算术运算有:+(加)、-(减)、*(乘)、/(右除)、(左除)、^(乘方)、’(转置)。运算是在矩阵意义下进行的,单个数据的算术运算只是一种特例。

(1)矩阵加减运算假定有两个矩阵A和B,则可以由A+B和A-B实现矩阵的加减运算。运算规则是:若A和B矩阵的维数相同,则可以执行矩阵的加减运算,A和B矩阵的相应元素相加减。如果A与B的维数不相同,则MATLAB将给出错误信息,提示用户两个矩阵的维数不匹配。

(2)矩阵乘法假定有两个矩阵A和B,若A为m*n矩阵,B为n*p矩阵,则C=A*B为m*p矩阵。

(3)矩阵除法在MATLAB中,有两种矩阵除法运算:和/,分别表示左除和右除。如果A矩阵是非奇异方阵,则AB和B/A运算可以实现。

AB等效于A的逆左乘B矩阵,也就是inv(A)*B,而B/A等效于A矩阵的逆右乘B矩阵,也就是B*inv(A)。对于含有标量的运算,两种除法运算的结果相同。对于矩阵来说,左除和右除表示两种不同的除数矩阵和被除数矩阵的关系,一般AB≠B/A。

(4)矩阵的乘方一个矩阵的乘方运算可以表示成A^x,要求A为方阵,x为标量。

(5)矩阵的转置对实数矩阵进行行列互换,对复数矩阵,共轭转置,特殊的,操作符.’共轭不转置(见点运算);

(6)点运算在MATLAB中,有一种特殊的运算,因为其运算符是在有关算术运算符前面加点,所以叫点运算。点运算符有.*、./、.和.^。两矩阵进行点运算是指它们的对应元素进行相关运算,要求两矩阵的维参数相同。

2、关系运算

MATLAB提供了6种关系运算符:&lt;(小于)、&lt;=(小于或等于)、&gt;(大于)、&gt;=(大于或等于)、==(等于)、~=(不等于)。关系运算符的运算法则为:

(1)当两个比较量是标量时,直接比较两数的大小。若关系成立,关系表达式结果为1,否则为0;

(2)当参与比较的量是两个维数相同的矩阵时,比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行,并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵,它的元素由0或1组成;

(3)当参与比较的一个是标量,而另一个是矩阵时,则把标量与矩阵的每一个元素按标量关系运算规则逐个比较,并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵,它的元素由0或1组成。

3、逻辑运算

MATLAB提供了3种逻辑运算符:&(与)、|(或)和~(非)。逻辑运算的运算法则为:

(1)在逻辑运算中,确认非零元素为真,用1表示,零元素为假,用0表示;

(2)设参与逻辑运算的是两个标量a和b,那么,a&b a,b全为非零时,运算结果为1,否则为0。a|b a,b中只要有一个非零,运算结果为1。~a当a是零时,运算结果为1;当a非零时,运算结果为0。

(3)若参与逻辑运算的是两个同维矩阵,那么运算将对矩阵相同位置上的元素按标量规则逐个进行。最终运算结果是一个与原矩阵同维的矩阵,其元素由1或0组成;

(4)若参与逻辑运算的一个是标量,一个是矩阵,那么运算将在标量与矩阵中的每个元素之间按标量规则逐个进行。最终运算结果是一个与矩阵同维的矩阵,其元素由1或0组成;

(5)逻辑非是单目运算符,也服从矩阵运算规则;

(6)在算术、关系、逻辑运算中,算术运算优先级最高,逻辑运算优先级最低。

(2)如何用matlab编程矩阵的初等变换扩展阅读:

1.获取矩阵元素

可以通过下标(行列索引)引用矩阵的元素,如Matrix(m,n)。

也可以采用矩阵元素的序号来引用矩阵元素。

矩阵元素的序号就是相应元素在内存中的排列顺序。

在MATLAB中,矩阵元素按列存储。

序号(Index)与下标(Subscript)是一一对应的,以m*n矩阵A为例,矩阵元素A(i,j)的序号为(j-1)*m+i。

其相互转换关系也可利用sub2ind和ind2sub函数求得。

2.矩阵拆分

利用冒号表达式获得子矩阵:

(1)A(:,j)表示取A矩阵的第j列全部元素;A(i,:)表示A矩阵第i行的全部元素;A(i,j)表示取A矩阵第i行、第j列的元素。

(2)A(i:i+m,:)表示取A矩阵第i~i+m行的全部元素;A(:,k:k+m)表示取A矩阵第k~k+m列的全部元素,A(i:i+m,k:k+m)表示取A矩阵第i~i+m行内,并在第k~k+m列中的所有元素。

此外,还可利用一般向量和end运算符来表示矩阵下标,从而获得子矩阵。end表示某一维的末尾元素下标。

利用空矩阵删除矩阵的元素:

在MATLAB中,定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。注意,X=[]与clear X不同,clear是将X从工作空间中删除,而空矩阵则存在于工作空间中,只是维数为0。

3、特殊矩阵

(1)魔方矩阵魔方矩阵有一个有趣的性质,其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵,其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n),其功能是生成一个n阶魔方阵。

(2)范得蒙矩阵范得蒙(Vandermonde)矩阵最后一列全为1,倒数第二列为一个指定的向量,其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中,函数vander(V)生成以向量V为基础向量的范得蒙矩阵。

(3)希尔伯特矩阵在MATLAB中,生成希尔伯特矩阵的函数是hilb(n)。使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中,有一个专门求希尔伯特矩阵的逆的函数invhilb(n),其功能是求n阶的希尔伯特矩阵的逆矩阵。

(4)托普利兹矩阵托普利兹(Toeplitz)矩阵除第一行第一列外,其他每个元素都与左上角的元素相同。生成托普利兹矩阵的函数是toeplitz(x,y),它生成一个以x为第一列,y为第一行的托普利兹矩阵。这里x,y均为向量,两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。

(5)伴随矩阵MATLAB生成伴随矩阵的函数是compan(p),其中p是一个多项式的系数向量,高次幂系数排在前,低次幂排在后。

(6)帕斯卡矩阵我们知道,二次项(x+y)n展开后的系数随n的增大组成一个三角形表,称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。函数pascal(n)生成一个n阶帕斯卡矩阵。

网络——MATLAB

3. matlab化矩阵标准型

在MATLAB中,计算矩阵A的特征值和特征向量的函数是eig(A),常用的调用格式有3种:
(1) E=eig(A):求矩阵A的全部特征值,构成向量E。
(2) [V,D]=eig(A):求矩阵A的全部特征值,构成对角阵D,并求A的特征向量构成V的列向量。
(3) [V,D]=eig(A,‘nobalance’):与第2种格式类似,但第2种格式中先对A作相似变换后求矩阵A的特征值和特征向量,而格式3直接求矩阵A的特征值和特征向量。

其他的很好解决吧

4. MATLAB怎么求解两个矩阵之间的变换矩阵

1:A=rand(4,3);A=[A(1,:),A(2,:);A(3,:),A(4,:)];2:用reshape函数b=reshape(a,[m n]);b=reshape(a,[2 6]);注意转换后和转换前的元素个数要相等也就是2x6要等于3x41 矩阵的建立1.直接输入法最简单的建立矩阵的方法是从键盘直接输入矩阵的元素。具体方法如下:将矩阵的元素用方括号括起来,按 矩阵行的顺序输入各元素,同一行的各元素之间用空格或逗号分隔,不同行的元素之间用分号分隔。2.利用M文件建立矩阵对于比较大且比较复杂的矩阵, 可以为它专门建立一个M文件。下面通过一个简单例子来说明如何利用M文件创建矩阵。例2-2 利用M文件建立MYMAT矩阵。(1) 启动有关编辑程序或MATLAB文本编辑器,并输入待建矩阵:(2) 把输入的内容以纯文本方式存盘(设文件名为mymatrix.m)。(3) 在MATLAB命令窗口中输入mymatrix,即运行该M文件,就会自动建立一个名为MYMAT的矩阵,可供以后使用。3.利用冒号表达式建立一个向量冒号表达式可以产生一个行向量,一般格式是:e1:e2:e3其中e1为初始值,e2为步长,e3为终止值。在MATLAB中,还可以用linspace函数产生行向量。其调用格式为:linspace(a,b,n)其中a和b是生成向量的第一个和最后一个元素,n是元素总数。显然,linspace(a,b,n)与a:(b-a)/(n-1):b等价。4.建立大矩阵大矩阵可由方括号中的小矩阵或向量建立起来。

5. 如何使用Matlab进行矩阵变换

如何根据题主的要求,使用Matlab进行矩阵变换?我们可以按下列步骤来考虑:

第一步:A和B矩阵组合成C矩阵。即

C=[A;B];

第二步:使用size()和max()函数,使用zeros()函数创建满足题主要求大小的D零矩阵,即

[m,n]=size(C);

D=zeros(m,max(max(C)));

第三步:使用for循环语句,将满足A和B矩阵中的数值,放入D矩阵中。即

for i=1:length(A)

D1=A(i);

D(1,D1)=A(i);

end

for j=1:length(B)

D2=B(j);

D(2,D2)=B(j);

end

第四步:显示变换后的D矩阵

D

阅读全文

与如何用matlab编程矩阵的初等变换相关的资料

热点内容
苹果手表固件描述文件下载 浏览:330
android获取sd卡视频文件 浏览:949
苹果手机设置通用网络设置 浏览:83
md298zpa是什么版本 浏览:317
srslog文件在哪个目录 浏览:948
无法找到文件中可删除的图片 浏览:739
dnf90版本副职业 浏览:848
c只读打开文件 浏览:575
如何在电脑文件上添加图片 浏览:297
xslist网站怎么登录 浏览:735
ftp不能直接打开文件 浏览:145
ps调整后的xps文件 浏览:572
小米如何取消wifi和数据同时使用 浏览:347
微信数据6个g怎么清理 浏览:533
找厂房去哪个app 浏览:881
linuxmini 浏览:997
如何找编程类的工作 浏览:286
jsp从mysql读取时间 浏览:680
有什么app可以存app 浏览:603
游戏编程从哪里学的 浏览:738

友情链接