基本语法
变量命名规则
- 区分大小写
- 长度不超过63位
- 字母开头,可以有字母、下划线和数字组成,但不能使用标点
- 应该简洁明了
命令行窗口
>>>clc
清楚命令窗口
>>> claer all
清理工作区内容
注释
%%
注释符
数据类型
1.数字
1+1
2*2
1/1
10-7
2.字符串
matlab
s = 'a'
abs(s) %%Ascll码
char(65) %% 65代表 'a'
num2str(65)%% 代表65
str = 'I Love MATLAB'
length(str) %%计算字符串长度
3.矩阵
matlab
A = [1 2 3;4 5 2;3 2 7;4 6 8] %% 3*3矩阵
B = A' %% 矩阵转置
C = A(:) %% 把A矩阵拉长(竖着拉长)
D = inv(A) %% 矩阵的求逆【矩阵必须是方阵】
A * D
E = zeros(10,5,3) %% 生成三维 10*5的0矩阵
E(:,:,1) = rand(10,5)
E(:,:,2) = randi(5,10,5)
E(:,:,1) = randn(10,5)
1.rand
生成均匀分布的伪随机数。分布在(0~1)之间
主要语法:rand(m,n)
生成m行n列的均匀分布的伪随机数
rand(m,n,'double')
生成指定精度的均匀分布的伪随机数,参数可以 是'single'
rand(RandStream,m,n)
利用RandStram(类似于随机种子)生成伪随机数
2.randn
生成标准正态分布的伪随机数(均值为0,方差为1)
主要语法:与rand
一样
3.randi
生成均匀分布的伪随机数
主要语法:randi(iMax)在开区间(0,iMax)生成均匀分布的伪随机整数
randi(iMax,m,n)在开区间(0,iMax)生成mXn型随机矩阵
r = randi([iMin,iMax],m,n)在开区间(iMin,iMax)生成mXn型随机矩阵
4.元胞数组
定义
matlab
A = cell(1,6) %%
A{2} = eye(3) %% 生成3*3对角线为1单位矩阵
A{5} = magic(5)%% 生成n阶的幻方矩阵
B = A{5}
5.结构体
matlab
books = struct('name',{{'Machine','Data Mining'}},'price,[30,40]')
MATLAB矩阵操作
矩阵定义与构造
matlab
A = [1 2 3 5 8 5 4 6]
B = 1:3:9 %%[1,9]步长为3
C = repmat(B,3,1) %% 让B重复三行一列(横着重复3次 竖着重复1次)
D = ones(2,4)%% 生成2*4全1的矩阵
矩阵运算
matlab
A = [1 2 3 4;5 6 7 8]
B = [1 1 2 2;2 2 1 1]
C = A + B
D = A - B
E = A * B'
F = A. *B %% 对应项相乘
G = A / B %%相当于A乘以B的逆
H = A./B %% 对应项相除
矩阵下标
注意:下标从1开始
matlab
A = magic(5);
B = A(2,3)
C = A(:,4)
D = A(:,4)
[m,n] = find(A > 20)
循环结构
for
结构
for 循环变量 = 初值 : 步长 :终值
[执行语句]
end
for 循环变量 = 初值 : 终值 (步长 1)
[执行语句]
end
matlab
%% 求和1^2 + 2^2 + 3^3 + 4^2 + 5^2
sum = 0
for n = 1 : 5
sum = sum + n^2;
end
sum
while
while 条件表达式
[执行语句]
end
matlab
s = 0;
n = 1;
while n<=10
s = s+n;
n = n+1;
end
s;
分支结构
if ... else结构
if 条件表达式
[语句体]
end
if...else...end
if 表达式
[语句体1]
else
[语句体2]
end
switch ... case ... end结构
switch 表达式(数值或字符串)
case 数值或字符串
语句体;
...
otherwise
语句体n;
end
二维平面画图
matlab
%% 画图
x = 0:0.01:2*pi;%% 步长0.01
y = sin(x);
%%figure 建立幕布
figure
plot(x,y); %%绘制当前x,y图
title('y = sin(x)');
xlabel('x');
ylabel('sin(x)');
%%画图
xlim([0 2*pi]) %% 规定区间
x = 0:0.01:20;
y1 = 200*exp(-0.05*x).*sin(x);
y2 = 0.08*exp(-0.5*x).*sin(10*x);
figure
[AX,H1,H2] = plotyy(x,y1,x,y2,'plot');%%共用一个坐标系
set(get(AX(1),'Ylabel'),'String','Slow Decay');
set(get(AX(2),'Ylabel'),'String','Fast Decay');
xlabel('Time (\musec)');
title('Multiple Decay Rates');
set(H1,'LineStyle','--');
set(H2,'LineStyle',':');