matlab中矩阵和向量的创建
1、向量的创建
1)直接输⼊:
⾏向量:a=[1,2,3,4,5]
列向量:a=[1;2;3;4;5]
2)⽤“:”⽣成向量
a=J:K ⽣成的⾏向量是a=[J,J+1,…,K]
a=J:D:K ⽣成⾏向量a=[J,J+D,…,J+m*D],m=fix((K-J)/D)
3)函数linspace ⽤来⽣成数据按等差形式排列的⾏向量
x=linspace(X1,X2):在X1和X2间⽣成100个线性分布的数据,相邻的两个数据的差保持不变。构成等差数列。
x=linspace(X1,X2,n): 在X1和X2间⽣成n个线性分布的数据,相邻的两个数据的差保持不变。构成等差数列。
4)函数logspace⽤来⽣成等⽐形式排列的⾏向量
X=logspace(x1,x2) 在x1和x2之间⽣成50个对数等分数据的⾏向量。构成等⽐数列,数列的第⼀项x(1)=10x1,x(50)=10x2
X=logspace(x1,x2,n) 在x1和x2之间⽣成n个对数等分数据的⾏向量。构成等⽐数列,数列的第⼀项x(1)=10x1,x(n)=10x2
注:向量的的转置:x=(0,5)’
2、矩阵的创建
1)直接输⼊:
将数据括在[]中,同⼀⾏的元素⽤空格或逗号隔开,每⼀⾏可以⽤回车或是分号结束。
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
2)函数eye,⽣成单位矩阵
eye(n) :⽣成n*n阶单位E
eye(m,n):⽣成m*n的矩阵E,对⾓线元素为1,其他为0
eye(size(A)):⽣成⼀个矩阵A⼤⼩相同的单位矩阵
eye(m,n,classname):对⾓线上⽣成的元素是1,数据类型⽤classname指定。其数据类型可以是:duoble、single、int8、uint8、int16、uint16、int32、uint32 。
3)函数ones ⽤ones⽣成全1的矩阵
ones(n) : ⽣成n*n的全1矩阵
ones(m,n) : ⽣成m*n的全1矩阵
ones(size(A)) : ⽣成与矩阵A⼤⼩相同的全1矩阵
ones(m,n,p,…)⽣成m*n*p*….的全1的多维矩阵
ones(m,n,…,classname)制定数据类型为classname
4)函数zeros 函数zeros⽣成全0矩阵
zeros(n):⽣成n*n的全0矩阵
zeros(m,n:)⽣成m*n的全0矩阵
zeros(size(A)): ⽣成与矩阵A⼤⼩相同的全0矩阵
zeros (m,n,p,…)⽣成m*n*p*….的全0的多维矩阵
zeros (m,n,…,classname)指定数据类型为classname
linspace函数调用的格式为5)函数rand 函数rand⽤来⽣成[0,1]之间均匀分布的随机函数,其调⽤格式是:Y=rand:⽣成⼀个随机数
Y=rand(n):⽣成n*n的随机矩阵
Y=rand(m,n):⽣成m*n的随机矩阵
Y=rand(size(A)):⽣成与矩阵A⼤⼩相同的随机矩阵
Y=rand(m,n,p,…):⽣成m*n*p*…的随机数多维数组
6)函数randn 函数rand⽤来⽣成服从正态分布的随机函数,其调⽤格式是:
Y=randn:⽣成⼀个服从标准正态分布的随机数
Y=randn(n):⽣成n*n的服从标准正态分布的随机矩阵
Y=randn(m,n):⽣成m*n的服从标准正态分布的随机矩阵
Y=randn(size(A)):⽣成与矩阵A⼤⼩相同的服从标准正态分布的随机矩阵
Y=randn(m,n,p,…):⽣成m*n*p*…的服从标准正态分布的随机数多维数组
3、矩阵元素的提取与替换
1)单个元素的提取
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊b=a(1,2)
b =
2
2)提取矩阵中某⼀⾏的元素,
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊b=a(1,:)
b =
1 2 3
3)提取矩阵中某⼀列:
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊b=a(:,1)
b =
1
3
4)提取矩阵中的多⾏元素
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊b=a([1,2],:)
b =
1 2 3
3 4 5
5)提取矩阵中的多列元素
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊b=a(:,[1,3])
b =
1 3
3 5
6)提取矩阵中多⾏多列交叉点上的元素如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊b=a([1,2],[1,3])
b =
1 3
3 5
7)单个元素的替换:
如:a=[1,2,3;3,4,5],运⾏后:
a =
1 2 3
3 4 5
输⼊:a(2,3)=-1
a =
1 2 3
3 4 -1
4、矩阵元素的重排和复制排列
1)矩阵元素的重排
B=reshape(A,m,n):返回的是⼀个m*n矩阵B,矩阵B的元素就是矩阵A的元素,若矩阵A的元素不是m*n个则提⽰错误。B=reshape(A,m,n,p):返回的是⼀个多维的数组B,数组B中的元素个数和矩阵A中的元素个数相等
B=reshape(A,…,[],…):可以默认其中的⼀个维数
B=reshape(A,siz) : 由向量siz指定数组B的维数,要求siz的各元素之积等于矩阵A的元素个数
2)矩阵的复制排列函数是repmat
B=repmat(A,n):返回B是⼀个n*n块⼤⼩的矩阵,每⼀块矩阵都是A
B=repmat(A,m,n):返回值是由m*n个块组成的⼤矩阵,每⼀个块都是矩阵A。
B=repmat(A,[m,n,p,…]):返回值B是⼀个多维数组形式的块,每⼀个块都是矩阵A
5、矩阵的翻转和旋转
1)矩阵的左右翻转左右翻转函数是fliplr,调⽤格式:
B=fliplr(A):将矩阵A左右翻转成矩阵B。
输⼊:A=[1,2,3;3,4,2]
A =
1 2 3
3 4 2
输⼊:B=fliplr(A)
B =
3 2 1
2 4 3
2)矩阵上下翻转函数:flipud,调⽤格式:
B=flipud(A):把矩阵A上下翻转成矩阵B
3)多维数组翻转函数:flipdim,调⽤格式:
B=flipdim(A,dim):把矩阵或多维数组A沿指定维数翻转成B
4)矩阵的旋转函数:rot90,调⽤格式:
B=rot90(A):矩阵B是矩阵A沿逆时针⽅向旋转90。得到的
B=rot90(A,k):矩阵B是矩阵A沿逆时针⽅向旋转k*90。得到的(要想顺时针旋转,k取-1)
6、矩阵的⽣成与提取函数
1)对⾓线函数对⾓线函数diag既可以⽤来⽣成矩阵,⼜可以来提取矩阵的对⾓线元素,其调⽤格式:
a) A=diag(v,k):当v是有n个元素的向量,返回矩阵A是⾏列数为n+|k|的⽅阵。向量v的元素位于A的第k条对⾓线上。K=0 对应主对⾓线,k>0对应主对⾓线以上,k<0对应主对⾓线以下。
b) A=diag(v):将向量v的元素放在⽅阵A的主对⾓线上,等同于A=diag(v,k)中k=0的情况。
c) v=diag(A,k):提取矩阵A的第k条对⾓线上的元素于列向量v中。
d) v=diag(A):提取矩阵A的主对⾓线元素于v中,这种调⽤等同于v=diag(A,k)中k=0的情况。
2)下三⾓阵的提取⽤函数tril,调⽤格式:
a) L=tril(A): 提取矩阵A的下三⾓部分
b) L=tril(A,k):提取矩阵A的第k条对⾓线以下部分。K=0 对应主对⾓线,k>0对应主对⾓线以上,k<0对应主对⾓线以下。
3)上三⾓阵的提取函数triu,调⽤格式:
a) U=triu(A): 提取矩阵A的上三⾓部分元素
b) U=triu(A,k): 提取矩阵A的第k条对⾓线以上的元素。K=0 对应主对⾓线,k>0对应主对⾓线以上,k<0对应主对⾓线以下。
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论