Matlab操作矩阵的相关方法
下面这篇文章主要是对吴恩达老师机器学习中matlab操作的一个整理和归纳
一、基本操作
1.生成矩阵（ones、zeros）
A = [1 2;3 4;5 6] #生成3行4列的矩阵
B = [1 2 3] #B就是一个行向量
C = [1;2;3] #定义c为一个列向量
D = 1:0.1:2 #定义开始值为1，步长为0.1，结束值为2的一个行向量
E = 1:6 #定义开始值为1，步长默认为1，结束值为6的行向量
ones(2,3) #矩阵中所有元素都为1 定义一个2行3列的矩阵
zeros(2,3) #矩阵中所有的元素都为0 定义一个2行3列的矩阵
2.生成随机矩阵（rand、randn）
rand(1,3) #生成1行3列的随机矩阵
randn(2,3) #生成高斯随机矩阵，高斯随机矩阵即为标准差或方差为1
3.生成单位矩阵（eye(n)）
eye(n) #生成n行n列的单位矩阵
4.帮助命令（help）
help 变量名 #可查看函数的API详解
二、移动数据
1.操作.txt文件（load）
1.1 加载.txt文件并且拆分文件的行和列的值
data = load('文件路径') #加载文件
获取多列的数据（获取多行的数据和多列类似，只需要修改第一个参数即可）
data(:,1) #拿到所有行第一列的数据
data(:,1:2) #拿到所有行第一列和第二列的数据
data(:,1:3) #拿到所有行第一列、第二列和第三列的数据
data(:,[1,3]) #拿到所有行第一列和第三列的数据
将矩阵所有的数据扁平化为一列
data(:)
将矩阵所有的数据扁平化为一行
data(:)'
1.2 将数据保存为.txt文件
v = data(:,1) #拿到第一列的数据
save test.txt v -ascii #将数据保存到test.txt文件中
2.矩阵的操作
2.1 获得矩阵的行数和列数（size()）
size(A) #返回一个1行2列的矩阵 分别是矩阵的行数和列数
size(A,1) #返回矩阵的行数
size(A,2) #返回矩阵的列数
2.2 拿到矩阵的最大维度（length()）
length(A) #获得矩阵的行数和列数中维度较大的一个
2.3 通过矩阵索引获取某一个值
A(m,n) #索引到矩阵m行n列的位置
2.4 修改矩阵的某一行或者某一列
A(:,2) = [10;11;12] #修改矩阵第二列的数据
2.5 在矩阵中添加一行新的数据
A = [A,[10;11;12]] #向矩阵中添加一行新的数据
C=[A B]
2.6 矩阵的结合
横向结合：
A = [1 2;3 4;5 6]
B = [11 12;13 14;15 16]
C = [A B]
纵向结合：
C= [A;B]
三、计算数据
1.A.*B(矩阵之间的乘积)
A .*B # A中对应位置元素和B中对应位置元素的乘积
2.A.^2 (矩阵自身的平方)
A.^2 #矩阵A的平方（A矩阵中的每个元素都平方）
3.1./A(矩阵中每个元素的倒数)
1./A 矩阵A中每个元素分别求倒数
4.log(A) (对矩阵中每个元素求对数) ,exp(A)(对A中的每个元素以e的底数)
5.abs(A)(对矩阵中的每个元素求绝对值)
6.-A(对矩阵中的每个元素求相反数)
7.A+1(对矩阵中每个对应的元素+1)
8.A’(A的转置)
9.一些有用的函数
求矩阵中最大的一个值：
max(max(A))
或者
max(A(:)) 首先扁平化A成为一个列向量，然后求最大值
max是默认求每列的最大值：
max(A) #求矩阵A的最大值（如果A是矩阵，会拿到每一列的最大值）
max(A,[],1) #拿到矩阵A中每一列的最大值
max(A,[],2) #拿到矩阵A中每一行的最大值
[val, ind] = max(a) #返回矩阵A中的最大值和索引
A<3 (对应元素的比较 如果小于3返回1，如果大于3返回0)
find(A<3) #找到A中所有小于3的元素，并且返回他们的索引
A=magic(3) #任意行、列、对角的元素相加的和等于相同的值
[r,c] = find(A>=7) #拿到所有大于等于7的元素的所在行和列
sum(A) #获得矩阵中所有元素的和
sum(A,1) #获得矩阵中每一列相加的和
sum(A,2) #获得矩阵中每一行相加的和
sum(sum(A)) #获得所有元素的值
prod(A) #获得矩阵中所有元素的乘积
floor(A) #对矩阵中所有元素向下取整
ceil(A) #对矩阵中所有元素向上取整
10.逆矩阵
pinv(A) #求A得逆矩阵
pinv(A)*A #就会拿到单位矩阵
四、数据绘制
1.绘制正弦函数
t = [0:0.01:0.98];
y1 = sin(2*pi*4*t);
plot(t,y1);
2.绘制余弦函数
t = [0:0.01:0.98];
y2 = cos(2*pi*4*t);
plot(t,y2);
3.同时绘制正弦函数和余弦函数
t = [0:0.01:0.98];
y1 = sin(2*pi*4*t);
y2 = cos(2*pi*4*t);
plot(t,y1);
hold on; #hold on 的作用是在旧的图像上绘制新的图像
plot(t,y2,'r')
xlabel('times'); #添加横轴的label
ylabel('values'); #添加纵轴的label
legend('sin','cos') #将图例添加到右上角
title('my plot') #给图像一个title
close; #关闭图像
figure(1);plot(t,y1);
figure(2);plot(t,y2); #给不同的图像命名
4.将图像分为一个1*2的格子
subplot(1,2,1) #前两个参数的意思是分为1*2的格子，后面一个参数的意思是当前使用第一个格子
5.改变轴的刻度
axis([0.5 1 -1 1]) #(xmin xmax ymin ymax)
6.清除一副图像（clf）
7.可视化一个矩阵
A = magic(5)
imagesc(A);
imagesc(A), colorbar, colormap gray; #一个灰度分布图