最近一直在查找矩阵函数之类的问题,上午正好有机会和大家分享一下.
机器学习- Matlab 编程常用命令速览
--总结自Ng-ML-class Octave/Matlab Tutorial Coursera
A、Basic operations and Moving data around
1 在命令行模式用shift + 回车便可附加下一行输出
2 length命令apply到矩阵时返回较高的一维的dimension
3 help + 命令是显示命令的简要帮助信息
doc + 命令 是显示命令的详细帮助文档
4 who 命令 显示 以后全部创立的变量
whos 命令 显示以后全部创立变量的详细信息
5 保存变量到.mat 文件save hello.mat b 以二进制压缩保存数据
save hello.txt v -ascii 以可读形式文件保存 即文本格式
6 ：means every elements in this col
7 A([1 3], :) 获取第 1、3两行全部列的数据
8 A = [A, [100; 101; 102]] 在A矩阵后面加一列col vector [100,101,102]
9 size(A) 返回一个1行2列矩阵 表明第1和第2个dimensional 的巨细
10 C = [A B]等价于C = [A, B] []为向后面的列添加，连接两个矩阵 [] 为concat 连接矩阵或者字符串
11 C= [A; B] ;号表现向上面行添加，因此会增长响应行数，列数不变
B、Computing on data
12 A.*B是矩阵/向量点乘 A*B是矩阵相乘
13 log(v) 和exp(v)求以e为底的对数和指数
14 abs()求绝对值
15 A‘ 求A的转置矩阵
16 max函数返回矩阵中最大元素的值和索引 [val, ind] = max(a)
17 A < 3 会判断A当中的每个是不是小于3，若小于3，对应位置返回true，否则对于位置返回false
18 find(A<3) 返回矩阵中全部值小于3的索引
19 magic Magic square.
magic(N) is an N-by-N matrix constructed from the integers
1 through N^2 with equal row, column, and diagonal sums.
Produces valid magic squares for all N > 0 except N = 2.
20 [r, c] = find(A >= 7) 返回值大于等于7的element的row及col的索引
21 prod(a) 求矩阵a里面全部元素的乘积
floor(a) 对矩阵a中元素向下取整
ceil(a)对矩阵a中元素向上取整
rand(3) 生成3X3的随机方阵
max（A,[],1） 求矩阵A的每一列的最大值（最后一维是1表明为dimension 1）
max（A,[],2） 求矩阵A的每一行的最大值
sum(A, 1) 对矩阵A第一维度（即每列）求和（注意matlab中第一维默许是列，然后是行，再然后依次类推。。。）
sum(A, 2) 对矩阵A第二维度（即每行）求和
sum(sum(A.*eye(9))) 求矩阵A的对角线元素之和
22 矩阵翻转操纵
flipud Flip matrix in up/down direction. 将矩阵上下翻转，
类似还有阁下翻转 fliplr, rot90, flipdim.
flipud(X) returns X with columns preserved and rows flipped
in the up/down direction. For example,
X = 1 4 becomes 3 6
2 5 2 5
3 6 1 4
23 pinv(A) 及inv(A) 求矩阵A 的逆矩阵
C、Plotting data
24 t = [0.1 ： 0.01 ： 0.98]
y = sin(t)
plot(t, y) 画正弦曲线
25 hold on; 保存以后曲线，画下一条曲线
26 xlabel 标定x轴说明
27 legend('sin','cos') 添加图例
28 title('my plot') 添加图片标题
29 print -dpng 'myPlot.png' 保存图片
30 线条颜色标注控制
b blue . point - solid
g green o circle : dotted
r red x x-mark -. dashdot
c cyan + plus -- dashed
m magenta * star (none) no line
y yellow s square
k black d diamond
w white v triangle (down)
^ triangle (up)
< triangle (left)
> triangle (right)
p pentagram
h hexagram
31 subplot 画子图，即将多个图合成一个图
>> subplot(1,2,1) % Divides plot a 1X2 grid, access the 1st element
>> plot(t,y1)
>> subplot(1,2,2)
>> plot(t,y2)
32 grid 加上网格
33 figure(1)
>> plot(t, y1)
>> figure(2)
>> plot(t,y2)
将y1与y2的曲线画到文件名为figure1和figure2的两个文件中
34 axis([0.5 1 -1 1]) 设定x轴范围为0.5～1，y轴范围为-1～1
35 clf Clear current figure.
36 imagesc(A)
imagesc Scale data and display as image.
把矩阵A画成黑色小方格
37 imagesc(A), colorbar, colormap gray;
colorbar 显示颜色渐变条，表明颜色含义
colormap 设置colormap性子 即RGB三色
A color map matrix may have any number of rows, but it must have
exactly 3 columns. Each row is interpreted as a color, with the
first element specifying the intensity of red light, the second
green, and the third blue. Color intensity can be specified on the
interval 0.0 to 1.0.
38 a = 1; b = 2; c=3; 不会打出a b c的值
a = 1, b = 2, c=3 会打出a b c 的值
D、Control statements: for, while, if statements
39 for 循环
for i = 1:10,
v(i) = 2 ^i;
end;
>> v
v =
2
4
8
16
32
64
128
256
512
1024
>> indices = 1:10
for i = indices,
disp(i)
end;
40 while 循环
>> while i <= 5,
v(i) = 100;
i = i+1;
end;
>> v
v =
100
100
100
100
100
64
128
256
512
1024
>> i = 1;
>> while true,
v(i) = 999;
i = i + 1;
if i == 6,
break;
end;
end;
>> exit
v =
999
999
999
999
999
64
128
256
512
1024
41 if elseif 判断分支语句
v(1)=2
v =
2
999
999
999
999
64
128
256
512
1024
>> if v(1) == 1,
disp('The value is one');
elseif v(1) == 2,
disp('The value is two');
else
disp('The value is not one or two.');
end;
The value is two
42 函数的定义和应用
定义函数
squareThisNumber.m文件
function
y = squareThisNumber(x)
每日一道理
闷热的天，蝉儿耐不住寂寞地不停在鸣叫，我孤单一人，寂静的身旁没有一个知音，想疯狂地听摇滚乐，听歇斯底里的歌声，那只为逃避无人的世界里那浓烈的孤单气息。一个人是清冷，两个人便是精彩，于是，莫名的冲动让我格外想念旧日的好友，怀念过去的日子，尽管不够现实的遐想追回不了曾经一切，但却希望思绪可以飞扬于闭上双目后的世界中，印有微笑，印有舞动的身姿，翩翩起舞……
y = x^2;
调用
squareThisNumber(5)
ans =
25
43 addpath 增长matlab搜索函数的路径
44 matlab里面定义的函数可以返回多于一个值， 这是其与C C++等编程语言
的不同之处，C\C++里面的函数有独一返回值
可以答应多个返回值可以带来编程上的方便
function
[a,b] = squareAndCubeThisNumber(x)
a = x^2;
b = x^3;
调用
>> [x1,x2] = squareAndCubeThisNumber(5)
x1 =
25
x2 =
125
45 cost function J 函数示例
function
J = costFunctionJ(X, y, theta)
% X is the "design matrix" containing our training examples.
% Y is the class labels
m = size(X,1);
% number of training examples
predictions = X*theta;
% predictions of hypothesis on all m examples
sqrErrors = (predictions - y).^2;
% squared errors
J = 1/(2*m) * sum(sqrErrors);
调用
X = [1 1; 1 2; 1 3];
y = [1; 2; 3];
theta = [0;1];
>> j = costFunctionJ(X, y, theta)
j =
0
% squared errors is 0 in this example
>> theta = [0;0]
theta =
0
0
>> j = costFunctionJ(X, y, theta)
j =
2.3333
% squared errors is 2.3333
in this example
% which is (1^2 + 2^2 + 3^2) / (2 * 3) = 2.3333
E、Vectorization
Matlab中“向量化”的实现算法，可以减少不必要的循环，对向量全部值的计算
操纵一致，这是更倏地更高效的算法实现思路，要注意体会matlab中“向量化”
“矩阵化”操纵变量实现算法与其他编程语言Java , C, C++的不同
文章结束给大家分享下程序员的一些笑话语录：
大家喝的是啤酒，这时你入座了。
你给自己倒了杯可乐，这叫低配置。
你给自已倒了杯啤酒，这叫标准配置。
你给自己倒了杯茶水，这茶的颜色还跟啤酒一样，这叫木马。
你给自己倒了杯可乐，还滴了几滴醋，不仅颜色跟啤酒一样，而且不冒热气还有泡泡，这叫超级木马。
你的同事给你倒了杯白酒，这叫推荐配置。
菜过三巡，你就不跟他们客气了。
你向对面的人敬酒，这叫p2p。
你向对面的人敬酒，他回敬你，你又再敬他……，这叫tcp。
你向一桌人挨个敬酒，这叫令牌环。
你说只要是兄弟就干了这杯，这叫广播。
有一个人过来向这桌敬酒，你说不行你先过了我这关，这叫防火墙。
你的小弟们过来敬你酒，这叫一对多。
你是boss，所有人过来敬你酒，这叫服务器。
酒是一样的，可是喝酒的人是不同的。
你越喝脸越红，这叫频繁分配释放资源。
你越喝脸越白，这叫资源不释放。
你已经醉了，却说我还能喝，叫做资源额度不足。
你明明能喝，却说我已经醉了，叫做资源保留。
喝酒喝到最后的结果都一样
你突然跑向厕所，这叫捕获异常。
你在厕所吐了，反而觉得状态不错，这叫清空内存。
你在台面上吐了，觉得很惭愧，这叫程序异常。
你在boss面前吐了，觉得很害怕，这叫系统崩溃。
你吐到了boss身上，只能索性晕倒了，这叫硬件休克。
---------------------------------
原创文章 By
矩阵和函数
---------------------------------