MATLAB数值计算

简单介绍MATLAB数值计算的命令。

---

%矩阵操作

u=[A;B]
u=[A,B]     %矩阵A,B构成新的矩阵
rand(m,n)   %随机矩阵，元素取值范围为（0,1）
eye(m,n)    %产生单位矩阵
ones(m,n)   %产生全为1的矩阵
rank(A)     %矩阵的秩
det(A)      %矩阵的行列式
inv(A)      %矩阵求逆
A=compan(u) %求多项式的伴随阵，和eig(A)一起使用可以得到多项式的解
[a,b]=eig(A)%求矩阵特征值和特征向量
sparse(i,j) %稀疏矩阵
magic(n)    %魔方矩阵
diag(v)     %对角矩阵
A(*,*)=*    %修改矩阵
str=mat2str(A); %矩阵转为字符串
B/A         %B*inv(A),XA=B的解
A\B         %inv(A)*B，AX=B的解
A^p         %A是方阵，P是大于1的矩阵
C=kron(A,B) %A是m*n矩阵，B是p*q的矩阵，C是mp*nq矩阵，A与B的张量积
[V,D]=eig(X) %矩阵特征值和特征向量，X*V=V*D，D是对角线上元素为特征值的矩阵
b=reshape(a,3,4) %矩阵的变维
rot90(A)    %矩阵旋转
fliplr(A)   %矩阵上翻
flipud(A)   %矩阵下翻
diag(A)     %抽取矩阵对角线
triu(A)     %抽取矩阵上三角
tril(A)     %抽取矩阵下三角
<,<=,>,>=,==,~=    %关系运算，产生一个0/1矩阵

---

%数组操作

a+b         %矩阵对应元素相加
a-b         %矩阵对应元素相减
a.*b        %矩阵对应元素相乘
a*b         %矩阵乘法
a./b=b.\a   %矩阵对应元素相除,a被b除
a.\b=b./a   %矩阵对应元素相除,b被a除
a.^b        %数组乘方,a11^b11,a12^b12,……

---

%集合运算

[b,i,j]=unique(a)   %取向量a的不重复元素构成的集合向量,i,j体现元素在原来向量的位置，如果b是矩阵，则加上参数‘rows’表示返回同行元素组成的矩阵
b=ismember(a,A)     %当a中元素属于A时，bi=1,否则bi=0,对于a中所有元素都进行此操作，得到一个0\1向量
b=ismember(A,B,'rows')  %A,B按行比较,行相同取，否则取0，但是要求列数相同
[c,ai,bi]=intersect(a,b)    %返回向量a,b的交集，ai表在a中的位置，bi表示在b中位置。当a,b为矩阵时，返回相同行，同样要求列数相同
[c,ai,bi]=union(a,b)    %返回a,b的并集，元素唯一。与交类似
[c,i]=setdiff(a,b)      %返回属于a但是不属于b的不同元素的集合,i表示c中元素在A中的位置,a,b也可以是矩阵，矩阵则是按行进行比较

---

%多项式运算

p=poly(a)       %产生根为a的多项式的系数向量,特征多项式是n+1维的,第一个元素为1,如果a是矩阵，生成的是特征多项式
p1=poly2str(p,'x')      %显示多项式形式
c=conv(a,b)     %多项式乘法
[d,r]=deconv(a,b)   %多项式除法,d是a除b后的整数,r是余数
[p,q]=polyder(a,b)  %多项式微分,返回多项式a,b乘积的微分系数p,以及多项式除法的a/b的微分系数