Coursea-吴恩达-machine learning学习笔记（五）【week 2之Octave/Matlab Tutorial】

基础操作：

在$Octave$中的部分符号表示：
不等于：~=$\qquad$等于：==
逻辑与：&&$\qquad$逻辑或：||$\qquad$异或运算：$XOR(A,B)$

变量赋值：$a=3\qquad$如果给一个变量赋值，不希望在屏幕显示结果，可以在命令后面加上分号$(如：a=3；)$

复杂的屏幕显示可以用$disp()$命令来实现；
$format$用来控制输出格式，例如：$format\ long$

建立矩阵：$A=[1\quad 2\ ;3\quad 4\ ;5\quad 6]$
$A$为$3\times 2$矩阵，‘；’作用为矩阵内换行

$V=1:0.1:2\quad$$V$表示一组值，从数值$1$开始，步长为$0.1$，直到增加到2。

$ones(2,3)$表示生成一个$2\times 3$的矩阵，矩阵中的所有元素都为1；
$2*ones(2,3)$生成一个$2\times 3$的矩阵，矩阵中的所有元素都为2；
$zeros(1,3)$生成一个$1\times 3$的矩阵，元素全部是0；
$rand(1,3)$表示$1\times 3$矩阵中的元素均为随机值(取自均值为$0.5$，幅度为$0\sim 1$的均匀分布)；
$W=randn(1,3)$表示$1\times 3$矩阵，元素取自均值为$0$，方差为$1$的高斯分布。

$hist(W)$命令表示绘制直方图；
$eye(4)$表示生成一个$4\times 4$的单位矩阵；
$help$命令用来了解其他命令的功能，例如：$help\quad eye$。

移动数据：

$size(A)$命令返回矩阵$A$的维度，结果形式为一个$1\times 2$的矩阵；
$size(A,1)$返回A的行数，$size(A,2)$返回A的列数；

若存在向量$V=[1,2,3,4]$，$length(V)$返回最大维度的大小；

$pwd$命令显示当前路径；
用$cd$命令切换到所需路径下；
$ls$命令显示路径下所有文件；

在$Octave$中，用$load('文件名')$来读入文件；

$who$命令可以显示出在当前$Octave$空间中的所有变量名；
$whos$命令可以详细显示所有变量的信息；
$clear\ $变量名：用来清除变量；
$clear$：清除所有变量；

$save\ $文件名 变量名：表示把变量存入文件内；
$save\ $文件名 变量名 $-ascii$：把变量以固定编码形式存入文件内；

$A(3,2)$取矩阵$A$的$A_{32}$元素；
$A(2,:)$取矩阵$A$的第$2$行元素；
$A([1\ 3],:)$表示取矩阵的第$1、3$行元素；
$A=[A,[100;101;102]]$表示在矩阵$A$的右侧加一列；
$A(:)$表示把矩阵$A$的所有元素放入一个向量。

计算数据：

在$Octave$中，$*$就是矩阵乘法，设$A$为$m*n$矩阵，$B$为$n*t$矩阵，则$C=A*B$为$m*t$矩阵，$C_{ij}$元素为$A$的第$i$行乘以$B$的第$j$列；

$.*$只应用于同型矩阵，若$A$为$m*n$矩阵，$B$也必须是$m*n$矩阵，$C=A.*B$也为$m*n$矩阵，$C_{ij}=A_{ij}*B_{ij}$；

若$C=A.\land2$，则$C_{ij}=A_{ij}*A_{ij}$；

$V=[1;2;3]$，则：
$C=1./V$，则$C_{ij}=1/V_{ij}$；
$C=log(V)$，则$C_{ij}=log(V_{ij})$；
$C=exp(V)$，则$C_{ij}=exp(V_{ij})$；
$-V=-1*V$；
$V+1$，则$V_{ij}=V_{ij}+1$。

若$a=[1\ 15\ 2\ 0.5]$，$max(a)=15$，表示求向量元素的最大值；
命令$a<3$的结果：$[1\ 0\ 1\ 1]$；
$find(a<3)$的结果：$[1\ 3\ 4]$给出元素的位置；

$A=magic(3)= \left[ \begin{matrix} 8&1&6\\ 3&5&7\\ 4&9&2 \end{matrix} \right]$
$magic$函数返回一个矩阵，它的所有行、列、对角线的元素加起来都等于相同的值。
$[r,c]=find(A>=7)\quad$结果$r=\left[ \begin{matrix} 1\\ 3\\ 2 \end{matrix} \right]\quad c=\left[ \begin{matrix} 1\\ 2\\ 3 \end{matrix} \right]$
意思是找出$A$中大于等于$7$的元素，$r$表示行位置，$c$表示列位置。

$sum(a)$为求$a$向量所有元素之和；
$prod(a)$为求$a$向量所有元素乘积；
$floor(a)$为$a$向量元素均向下取整；
$cell(a)$为$a$向量元素均向上取整。

$C=max(A,B)$，则$C_{ij}=max(A_{ij},B_{ij})$；
$max(A,[],1)$求每一列最大值；
$max(A,[],2)$求每一行最大值；

上述$magic$矩阵$A$中：
$max(A)=[8\ 9\ 7]$；
$max(max(A))=9$。

$sum(A,1)$求得每一列的总和；
$sum(A,2)$求得每一行的总和。

$eye(3)=\left[ \begin{matrix} 1&0&0\\ 0&1&0\\ 0&0&1 \end{matrix} \right]\qquad flipud(eye(3))=\left[ \begin{matrix} 0&0&1\\ 0&1&0\\ 1&0&0 \end{matrix} \right]$

绘制数据：

t=[0:0.01:0.98];#生成一个数据集
y_1=sin(2*pi*4*t);
plot(t,y_1);#正弦函数

y_2=cos(2*pi*4*t);
plot(t,y_2);#余弦函数

plot(t,y_1);
hold on;#将新图跟旧图画在一起
plot(t,y_2,'r');#'r'表示颜色

#设置坐标名
xlabel('time');
ylabel('value');
#设置曲线名
legend('sin','cos');
#设置图片名
title('myplot');
#保存图片
print -dpng 'myplot.png'

#绘制两幅图
figure(1);
plot(t,y_1);
figure(2);
plot(t,y_2);

#将图片分为1*2的格子，用第一个绘图
subplot(1,2,1);

plot(t,y_2);
#设置x,y轴的起始刻度
axis([0.5 1 -1 1]);

A=magic(5);
#将A呈现为一个不同颜色的拼图
imagesc(A);

控制语句：

$for$语句：

for i=1:10,
    V(i)=2^i;
end;

$while$语句：

i=1;
while i<=5,
    V(i)=100;
    i=i+1;
end;

$while+if$语句：

i=1;
while true,
    v(i)=999;
    i=i+1;
    if i==6,
        break;
    end;
end;

函数定义：

function y=squareThisNumber(x)

y=x^2;

上述代码存放在$squareThisNumber.m$文件中，调用时，切换到$squareThisNumber.m$所在路径，执行$squareThisNumber(5)$即可。

将$squareThisNumber.m$路径通过$addpath('XXX')$添加，则无需切换路径，也可执行。

函数返回多个值：

function [y_1,y_2]=SquareAndCubeThisNumber(x)

y_1=x^2;
y_2=x^3;