MATLAB 基本命令 ...

最新推荐文章于 2023-10-18 22:15:43 发布

weixin_30566063

最新推荐文章于 2023-10-18 22:15:43 发布

阅读量454

点赞数

文章标签： matlab

原文链接：http://www.cnblogs.com/tolic/p/7142211.html

版权

注：每条命令后面跟不跟分号的区别就是：不加分号执行命令后立即显示出来，加上分号不显示。
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:2

v =
    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(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 %M的第一列都除以3，第二列都除以2

>> M ./ Y

ans =
    0.3333    1.0000
    1.0000    2.0000
    1.6667    3.0000

M .* Y %M的第一列都乘以3，第二列都乘以2

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.*B % A和B的对应元素相乘（所有要就AB的维度要相同）

A.^2 %得到A的每一个元素的平方

A ./ 2 %A的每一个元素除以2

>> A ./ 2

ans =
    0.5000    1.0000
    1.5000    2.0000
    2.5000    3.0000

2 ./ A %2除以A的每一个元素

>> 2 ./ A

ans =
    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 < 3

ans =
  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
     1
colunm =
     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 figure

imagesc(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 20

disp(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 priceY

load 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

转载于:https://www.cnblogs.com/tolic/p/7142211.html