matlab从入门到秃头

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

求1²⁺²2+3²⁺⁴2+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

求1²⁺²2+3²⁺⁴2+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