注:每条命令后面跟不跟分号的区别就是:不加分号执行命令后立即显示出来,加上分号不显示。
1.矩阵的建立A = [1 2; 3 4; 5 6]; % 建立一个 3x2的矩阵、
>> A = [1 2; 3 4; 5 6]A = 1 2 3 4 5 6
v = [1 2 3]; %建立一个1x3的行向量
>> v = [1 2 3] v = 1 2 3
v = [1;2;3]; %建立一个3维列向量
>> v = [1;2;3]v = 1 2 3
v = 1:0.2:2;%建立一个行向量,从1开始,步长为0.2,到2结束
>> v = 1:0.2:2v = 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000
v = 1:6; %建立一个行向量,从1开始,步长为1,到6结束
v = 1:2:6; %建立一个行向量,从1开始,步长为2,到6结束 ones(2,3); %建立一个2x3的矩阵,元素全是1>> ones(2,3)ans = 1 1 1 1 1 1
zeros(2,3);%建立一个2x3的矩阵,元素全是0
c = 2*ones(2,3);%建立一个2x3的矩阵c,元素全是2>> c = 2*ones(2,3)c = 2 2 2 2 2 2
rand(2,3);%建立一个2x3的随机矩阵
>> rand(2,3)ans = 0.8147 0.1270 0.6324 0.9058 0.9134 0.0975
randn(2,3); %同上
eye(4) % 4 x 4 单位矩阵A = [1 2; 3 4; 5 6];
size(A); %返回A的行数和列数>> A = [1 2; 3 4; 5 6];>> size(A)ans = 3 2
size(A,1); % 返回A的行数
size(A,2);% 返回A的列数 length(A);// 返回max(row,column)2.数据操作
A = [1 2; 3 4; 5 6];
A = 1 2 3 4 5 6
A(3,2);% 得到第三行第二列的元素
>> A(3,2)ans = 6
A(2,:); % 得到第二行的所有元素
>> A(2,:)ans = 3 4
A(:,1) ;%得到第一列的所有元素。‘:’表示沿着每一行/列的所有元素
>> A(:,1) ans = 1 3 5
A([1 3],:)%得到所有行号为1、3的元素,即第一行和第三行
A([1:2:7],:) %得到第一、三、五、七行>> A([1 3],:)ans 1 2 5 6
A(:,[1 2])%得到第一、二列
A(:,2) = [10;11;12] %用10,11,12来替换A的第二列
>> A(:,2) = [10;11;12]A = 1 10 3 11 5 12
A(2,:) = [23 24] %用23 24来替换第二行
>> A(2,:) = [23 24]A = 1 10 23 24 5 12>>
A = [A,[100;101;102]] %在A的右边加上一列
>> A = [A,[100;101;102]] A = 1 10 100 23 24 101 5 12 102
注:如果若干元素在同一行,用逗号或空格隔开;若在同一列用分号隔开
A = [A;[7,8,9]] %在A的下方添加一行
>> A = [A;[7,8,9]]A = 1 10 100 23 24 101 5 12 102 7 8 9
V=A(:) %让A的所有元素变成一列,并赋给V。
>> V = A(:)V = 1 23 5 7 10 24 12 8 100 101 102 9
A = [1 2; 3 4]
B = [6 7;8 9] C = [A B] ;%将A(左边) B(右边)组合成一个矩阵>> C = [A B]C = 1 2 6 7 3 4 8 9
C = [A;B] %将A(上) B(下)组合成一个矩阵,使用的时候注意维度的匹配
>> C = [A;B]C = 1 2 3 4 6 7 8 9
3.矩阵计算
M = [1 2;3 4;5 6] Y = [3 2]
M - Y %M的第一列都减3,第二列都减2>> M - Y ans = -2 0 0 2 2 4
M ./ Y %M的第一列都除以3,第二列都除以2
>> M ./ Yans = 0.3333 1.0000 1.0000 2.0000 1.6667 3.0000
M .* Y %M的第一列都乘以3,第二列都乘以2
>> M .* Y ans = 3 4 9 8 15 12
A = [1 2; 3 4; 5 6]
B = [11 12; 13 14; 15 16] C = [1 1; 2 2] A*C % 3x2 * 2x2 = 3x2 的矩阵>> A*Cans = 5 5 11 11 17 17
A.*B % A和B的对应元素相乘(所有要就AB的维度要相同)
>> A.*Bans = 11 24 39 56 75 96
A.^2 %得到A的每一个元素的平方
>> A.^2ans = 1 4 9 16 25 36
A ./ 2 %A的每一个元素除以2
>> A ./ 2ans = 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000
2 ./ A %2除以A的每一个元素
>> 2 ./ Aans = 2.0000 1.0000 0.6667 0.5000 0.4000 0.3333
log(A)
exp(A) abs(A) %对A的每一个元素进行相应的运算A’ %得到A的转置
>> A'ans = 1 3 5 2 4 6
max(A) %得到每一列的最大元素,如果A是一个1xn的向量,则返回最大值
>> max(A)ans = 5 6
max(max(A))% 返回的就是5,6中的最大值
>> max(max(A))ans = 6
max(A,[],1) // 等同于 max(A)
max(A,[],2) //得到每一行的最大值>> max(A,[],2)ans = 2 4 6
A < 3 %返回元素是否小于3的逻辑判断值
>> A < 3ans = 3×2 logical 数组 1 1 0 0 0 0
find(A<3) %返回元素小于3的元素的编号
>> find(A<3) ans = 1 4
[row,colunm] = find(A<3)%返回元素小于3的元素的行号和列号
>> [row,colunm] = find(A<3)row = 1 1colunm = 1 2 %即A中的(1,1),(1,2)都小于3
sum(A) %计算A中每一列的和,如果A为1xn的向量,则返回所有元素的和
>> sum(A) ans = 9 12
sum(sum(A)) %返回9 12 的和,也就是A中所以元素的和,等同于sum(A(:))
>> sum(sum(A))ans = 21
sum(A,1) %等同于 sum(A),返回每一列的和
sum(A,2) %返回每一行的和>> sum(A,2)ans = 3 7 11
prod(A) %返回每一列各个元素的乘积,等同于prod(A,1)
>> prod(A)ans = 15 48
prod(A,2) %返回每一行各个元素的乘积
**floor(r) %所有元素向下取整
ceil(r) % 所有元素向上取整**>> r=rand(2)r = 0.9649 0.9706 0.1576 0.9572>> floor(r)ans = 0 0>> ceil(r)ans = 1 1 1 1
flipud(eye(3))
>> flipud(eye(3))ans = 0 0 1 0 1 0 1 0 0
pinv(A) %A的逆运算
>> pinv(A)ans = -1.3333 -0.3333 0.6667 1.0833 0.3333 -0.4167
4.作图
t = [0:0.01:0.98]
y1 = sin(2*pi*4*t) plot(t,y1)y2=cos(2*pi*4*t);
hold on; %表示画下一个图的时候,上一个图不消失,可以使两幅图在一个窗口中 plot(t,y2,’r’)%’r’ 表示曲线颜色 xlabel(‘time’) ylabel(‘value’) %添加横纵坐标标签 legend(‘sin’,’cos’) %标识曲线 title(‘my plot’)%添加名称plot(t,y1,’rx’,’MarkerSize’,10) % ‘rx’ 代表用红叉作图,大小为10
print -dpng ‘myPlot’ %导出图片到当前路径
figure(1); plot(t,y1);%打开窗口一,并画图(可以两个命令同时执行,用分号隔开)
figure(2); plot(t,y2); %打开另外一个窗口二,作图subplot(1,2,1) // divides plot in a one-by-two grid, access first element
subplot(1,2,2) plot(t,y2,’r’) axis([0.5 1 -1 1])%设置X的范围 clf //clear figureimagesc(A) //视觉化矩阵A
imagesc(magic(5)), colorbar, colormap gray %用灰度图视觉化5.其他
mean(M,1) % 计算每一列的平均值
mean(M,2) % 计算每一行的平均值 std(M) % 计算每一列的标准差,等同于(M,0,1) std(M,0,2) % 计算每一行的平均值 size(M,1) % 一共多少行 size(M,2) % 一共多少列 numel(A) %返回一共有多少个元素 plot(J,’-b’,’Linewidth’,2) % plot a line in blue with width equal to 2 plot(X,’rx’,’MarkerSize’,10,’LineWidth’,20) % plot figure with red cross and width equal to 20disp(var) // display variable
disp(sprintf(‘2 decimals: %0.2f’,a)) //format with two decimals format long // format short // hist(w) //画出w值的直方分布图 pwd //reture the current path ls //display load filename //loading file load(‘filename’) //loading file who //show the all varable in current memory whos //show the all varable of detailed view // save and load data load priceY //load a data set (a vector) named priceY v = priceY(1:10)// v will be a vector with first 10 elements in priceYload dataX //load a data set (a matrix) named dataX
v = dataX(1:10) // v will be a 1 by 10 matrix with first 10 elements in first column of dataX save hello.mat v; // save the variable V into a file called hello.mat at current direction clear // deletest all of the variables in your workspace