Matlab note

%解释

edit 创建一个MATLAB文件，function要与文件名一致

seq = randseq(45);

MySample =fix( rand(12,4)*60) 随机生成一个12*4的矩阵，同时，数值在60以内。

dim1 = MySample(:,1);  定义维：第一列为第一维

std(dim2)^2 

cov(MySample) 求协方差

B=repmat(A,m,n) B是由A平铺成m*n的高维数组

mean(B) 求每一维的均值,也就是每一列的均值

mean(B,2) 求每一行的均值

inv(A)  求逆矩阵，A must square.

矩阵中行的分割用“，”，列的分割用“；”

A.’    求转置

c=a+b    求加法

d=a-b   求减法

e=a*b    a的列数等于b的行数

a.*b=[a1*b1,a2*b2,…an*bn]

a.^b=[a1^b1,a2^b2,…an^bn]

a./b=[a1/b1,a2/b2,…an/bn]

a.\b=[b1/a1,b2/a2,…bn/an]

若 a,b为一般向量， 如a=[1,2,3]   b=[3,4,5] ，dot(a,b) 求点积；cross(a,b) 求叉积；conv(a,b)  求卷积

size(A)  求A的行数和列数; [r,c]=size(A) r,c分别返回行列；size(A,n)

zscore(X)  标准化：使平均值为0，标准差为1；

消除单位影响和变异大小因素的影响。

m(1:3,2) 输出第2列，的1~3行

ones(3,5) 建立都是1的3*5矩阵

zeros(2,3) 建立都是0的2*3矩阵

eye(4) 建立4*4的单位矩阵; eye(2,3) 某个部分是单位矩阵，其他是零.因为MATLAB的元素是按列存储

b([1 3 5])=[]   去除1，3，5 

a=rand(n,m); b=a(:) 把a 变成(n+m)*1的矩阵

a(i1:i3,j4:j5)取矩阵a的i1~i3行,j4~j5列构成新的矩阵.

A(i2:-1:i1,:) 以逆序提取矩阵A的第i1~i2行，构成新矩阵 ??

B(:,j2:-1:j1) 以逆序提取矩阵的第j1~j2列，构成新矩阵 ??

A(i1:i2,:)=[]删除i1~i2行，构成矩阵A

C(:,j1:j2)=[]删除j1~j2列 构成矩阵C

 n=m(1:end-1, 2:end) 提取m的1~倒数第一行，2~最后一列 作为矩阵n

P=fix(rand(3)*10) 形成一个3*3矩阵，其中的元素小于10,也就是1~9随机排列，不能形成n*m格式。

u=o(:,[2,3]) 与y=o(:,2:3) 形成的矩阵一样，语句有什么不同？？？

[A B],[A;B]拼接矩阵,前一个是行，后一个是列

det(a) 计算方阵a的行列式

inv(A)   计算方阵A的逆矩阵

[V,D]=eig(I) 返回方阵的特征值 和 特征向量

e=eig(A)返回列向量—特征值

trace(A) 迹，方阵的对角线的和

与或非&|~

 roots(p)求多项式的根，p是多项式向量 求x^3-x^2+x-1=0的根roots([1,-1,1,-1])

solve(fun) 求fun=0的符号解，若不能求得精确的符号解，可以计算变精度的数值解 x^9+x^8+1=0的根 solve('x^9+x^8+1’)给出方程数值解（32位有效数字的符号量）

solve(fun,var) 对指定变量var 求代数方程fun=0的符号解。ax^3+bX^2+c=0  syms a b c x; f=a*x^3+b*X^2+c; solve(f)  如果不知名变量，系统默认为x，也可指定自变量，比如指定b为自变量： syms a b c x;f=a*x^3+b*X^2+c; solve(f,b)     ???没有定义function或者变量'x'

fsolve(fun,x0) 求非线性方程fun=0在估计值x0附近的近似解 x=e^(-x) 在0附近的根：fsolve(‘x-exp(-x)’,0)

fzero(fun,x0) 求函数fun在x0附近的零点。x-10^x+2=0在x0=0.5附近的根为fzero(‘x-10^x+2’,0.5)