Matlab从入门到秃头
基本知识
变量命名规则
- 变量名区分大小写
- 变量名长度不超过63位
- 变量名以字母开头,可以由字母、数字、下划线组成但不能使用标点
- 变量名应当简洁明了,通过变量名可以直观看出变量所表示的物理意义
常见命令
清空环境变量及命令
- clear all 清除workspace所有变量
- clc 清除command window所有命令
注释
注释整行:%%+注释内容
注释文字:%+注释内容
数据类型
数字
指基本数字的运算,加减乘除
1+2
1-2
3*8
8/2
字符与字符串
-
字符串
%% 字符与字符串 s='a' abs(s)%求字符的ascii码 char(97)%把ascii码转化成字符 num2str(65)%把数字转化为字符串 str='hello world' length(str) % 计算长度 doc num2str
矩阵
%%
% 矩阵
A=[1 2 3;4 5 6;3 2 7;]% 构造一个矩阵,行与行之间用分号隔开
B=A' % B矩阵是A矩阵的转秩
C=A(:) % C矩阵是A矩阵把每一列依次拼起来
D=inv(A) % D矩阵是A矩阵的逆矩阵
A*D % 求A矩阵和D矩阵的积
E=zeros(10,5,3) %创造10行5列全是0的三维矩阵
E(:,:,1)=rand(10,5) % 生成10行5列均匀分布在(0,1)之间的随机数
E(:,:,2)=randi([1,6],10,5) % 生成10行5列在开区间(1,6)之间的随机矩阵
E(:,:,3)=randn(10,5) % 生成10行5列服从标准正态分布的伪随机数
元胞数组
%%
% 元胞数组
A=cell(1,5)% 建立1行5列的元胞数组,其中每个子数组的数据类型可以不一样
A{2}=eye(3)% A{2}为一个3阶的单位矩阵,注意,下标从1开始
B=A{2}
A{5}=magic(5)%A{5}为一个5阶的幻方(也就是每行每列以及对角线的和相等)
C=A{5}
结构体
-
其实就是键值对
%% % 结构体 books=struct('name',{{'hello world','data machine'}},'price',[30,40]) books.name %{'hello world' 'data machine'} books.name(1) %取出的是cell {'hello world'} books.name{1} %取出来的是cell里面的值 'hello world'
四则运算
矩阵的定义和构造
%%
% 矩阵的定义和构造
A=[1 2 3; 5 8 5 ;4 6 7] %分号换行
B=1:2:9 %表示从1开始,每次增加2,最大值为9
C=repmat(B,3,2) %把B矩阵重复3行2列
D=ones(2,4) %生成2行4列全为1的矩阵
矩阵的四则运算
%%
% 矩阵的四则运算
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' % A和B的转秩序相乘
F=A.*B % A和B的对应项相乘
G=A/B %A*B的逆
H=A./B %对应项相除
矩阵的下标
%%
% 矩阵的下标
A=magic(5)
B=A(2,3) %取第二行,第三列
C=A(3,:) %取第三行
D=A(:,4) %取第四列
[m,n]=find(A>20)%找A中大于20的序号值,m存储行,n存储列
流程控制
for循环
for 循环变量=初值:步长:终值
循环语句1
循环语句2
end
求12+22+32+42+5^2的和
%%
sum=0;
for n=1:5
sum=sum+n^2;
end
sum
99乘法表
%%九九乘法表
for i =1:9
for j=1:9
a(i,j)=i*j
end
end
while循环
while 条件表达式
执行语句1
执行语句2
执行语句n
end
求12+22+32+42+5^2的和
n=1;
s=0;
while n <=5
s=s+n*n;
n=n+1;
end
s
if else end
%%
a=100;
b=20;
c=0;
if a>b
c=1;
else
c=2;
end
c
switch语句
switch 表达式
case数值或字符串1
语句体1;
case数值或字符串2
语句体2;
otherwise
语句体n
end
n=2;
res=0;
switch(n)
case 0
res=-1;
case 1
res=0;
otherwise
res=1;
end
res
绘图
二维平面绘图
demo1
x=0:0.01:2*pi; %定义x,范围为0到2pai,步长为0.01
y=sin(x); %定义对应的y为sin(x)
figure %建立一个幕布
plot(x,y) %绘制当前的二维的平面的图
title('y=sin(x)')%设置标题
xlabel('x') %x轴上标的标签
ylabel('sin(x)') %y轴上标的标签
xlim([0 2*pi]) %显示的范围是x从0到2pai
demo2
%%
%二维平面绘图
x=0:0.01:20;
y1=100*exp(-0.05*x).*sin(x); %exp返回指数 ex
y2=0.8*exp(-0.5*x).*sin(10*x);
figure
plotyy(x,y1,x,y2,'plot')
[Ax,H1,H2]=plotyy(x,y1,x,y2,'plot');%两个y共享一个x,调用函数plot来绘制图像
set(get(Ax(1),'Ylabel'),'String','Slow Decay'); %第一个y轴的名字叫Slow Decay
set(get(Ax(2),'Ylabel'),'String','Fast Decay'); %第二个y轴的名字叫Fast Decay
xlabel('Time(\musec)');
title('multiple decay rates');
set(H1,'LineStyle','--');% 把h1设置为虚线
set(H2,'LineStyle',':');%把h2设置为点线
三维绘图
%三维平面绘图
t=0:pi/50:10*pi;
plot3(sin(t),cos(t),t);% xyz轴对应的点分别取sin(t),cos(t),t
xlabel('sin(t)')
ylabel('cos(t)')
zlabel('t')
grid on %在原图之中加上网格线
axis square%让各个方向变成正方形,起到限制作用
完结撒花感谢陪伴
%%双峰函数
[x,y,z]=peaks(30);
mesh(x,y,z)
grid