grant's Blog

记录生活中的点点滴滴

Matlab快速入门笔记（配套代码）

发表于 更新于 分类于   阅读次数： Valine：
本文字数： 10k 阅读时长 ≈ 10 分钟

本文参考自B站西电通院科协glanny的1小时matlab速成视频！

视频链接

第1部分：变量定义和基本运算

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
%%
% 建议有C语言或其他编程基础，了解线性代数和矩阵相关知识
% https://ww2.mathworks.cn/help/matlab/
% 加*为了解内容
% 生成矩阵
    % 直接法
        a = [1,2,3;4,5,6;7,8,9];
    % 冒号一维矩阵 a = 开始：步长：结束，步长为1可省略
        b = 1:1:10;  % 1,2,...10
        b = 1:10;  %与上一个等价
    % 函数生成
        % linspace(开始，结束，元素个数)，等差生成指定元素数的一维矩阵，省略个数则生成100个
            c = linspace(0,10,5);
        % 特殊矩阵
            e = eye(4);  % eye(维数)单位阵
            z = zeros(1,4);  % zeros(维数)全零阵
            o = ones(4,1);  % ones(维数)全1阵
            r = rand(4);  % rand(维数)0~1分布随机阵
            rn = randn(4);  % randn(维数)0均值Gaussian分布随机阵
                      
%%
% 矩阵运算
    diag_a = diag(a,1);  % diag(行向量，主对角线上方第k条斜线)用行向量生成对角阵
    tril_a = tril(a,1);  % tril(矩阵，主对角线上方第k条斜线)生成矩阵的下三角阵，triu上三角阵
    % 加、减、乘、乘方
        a*a
    % 点运算
        % a.*b , a./b , a.\b , a.^b  对应元素的*,/,\,^运算
        a.*a
    % 逆矩阵
        pinv(a)  % 伪逆矩阵，当a不是方阵，求广义逆矩阵；当a是可逆方阵，结果与逆矩阵相同
    % 特征值，特征向量
        [v,D] = eig(a);  % 输出v为特征向量，D为特征值对角阵
    % *行列式
        det(a)
    % *秩
        rank(a)
    % *伴随
        compan(b)
        
%%
% 矩阵的修改
    %部分替换
        chg_a = a;
        chg_a(2,3) = 4;  % (行，列)元素替换
        chg_a(1,:) = [2,2,2];  % (行,:)替换行，为[]删除该行
        chg_a(:,1) = [];  % (:,列)替换列，为[]删除该列
    % 转置
        T_a = a';
    % 指定维数拼接
        c1_a = cat(1,a,a);  % 垂直拼接
        c2_a = cat(2,a,a);  % 水平拼接
    % *变维
        rs_a = reshape(a,1,9);  % 元素个数不变，矩阵变为m*n

%%
% 信息获取
    % 矩阵的行列数
        [row_a, col_a] = size(a);  % [行数，列数]
    % 行列中最大的
        len_a = length(a);
        
%%
% 多维数组
    % 创建
        % 直接法
            mul_1(:,:,1) = [1,2,3;2,3,4];
            mul_1(:,:,2) = [3,4,5;4,5,6];
        % *扩展法
            mul_2 = [1,2,3;2,3,4];
            mul_2(:,:,2) = [3,4,5;4,5,6];  % 若不赋值第一页，第一页全为0
        % cat法
            mul_31 = [1,2,3;2,3,4];
            mul_32 = [3,4,5;4,5,6];
            mul_3 = cat(3,mul_31,mul_32);  % 把a1a2按照“3”维连接

%%
% *字符串
    % 创建
        str0 = 'hello world';  % 单引号引起
        str1 = 'I''m a student';  % 字符串中单引号写两遍
        str3 = ['I''m' 'a' 'student'];  % 方括号链接多字符串
        str4 = strcat(str0, str1);  % strcat连接字符串函数
        str5 = strvcat(str0, str1);  % strvcat连接产生多行字符串
        str6 = double(str0);  % 取str0的ASCII值，也可用abs函数
        str7 = char(str6);  % 把ASCII转为字符串
    % 操作
        % 比较
            strcmp(str0, str1);  % 相等为1，不等为0
            strncmp(str0, str1, 3);  % 比较前3个是否相等(n)
            strcmpi(str0, str1);  % 忽略大小写比较(i)
            strncmpi(str0, str1, 3);  % 忽略大小写比较前3个是否相等
        % 查找替换
            strfind(str0, str1);  % 在str0找到str1的位置
            strmatch(str1, str0);  % 在str0字符串数组中找到str1开头的行数
            strtok(str0);  % 截取str0第一个分隔符（空格，tab，回车）前的部分
            strrep(str0, str1, str2);  % 在str0中用str2替换str1
        % 其他
            upper(str0);  % 转大写，lower转小写
            strjust(str0, 'right');  % 将str0右对齐，left左对齐，center中间对齐
            strtrim(str0);  % 删除str0开头结尾空格
            eval(str0);  % 将str0作为代码执行
            
%%
%转换
    % ___2___  -->  如num2str，将数字转字符串； dec2hex，将十进制转十六进制
    str_b = num2str(b);
    % abs，double取ASCII码；char把ASCII转字符串
    abs_str = abs('aAaA');

第2部分：程序结构

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
%%
a = 5;
x = [1, 2]; y =[3, 4];

%%
%选择结构
    %if-elseif-else-end
        if a>0
            disp(x);
        elseif a==0
            disp(a);
        else
            disp(a-1);
        end
        
    %switch-case-otherwise-end
        switch a
            case 0
                disp(a);
            case 1
                disp(a+1);
            otherwise
                disp('aaa');
        end
        
    %try-catch
        try
            z = x*y;
        catch
            z = x.*y;  % 若try出错，则执行
        end
        disp(z);

%%
% 循环结构
    % for 循环变量=初值:步长:终值 - end
    for i=0:1:10  % 步长为负，则初值大于终值
        disp(i);  % 循环体内不可对循环变量做修改
    end
    
    % while-end
    while a>2
        disp(a);
        a = a-1;
    end
    
%%
%程序控制
    %continue 跳过当次循环剩下语句，进入下一循环
    %break 跳出当前循环
    %return 跳出程序并返回
    
%%
%m文件
    %脚本文件：没有输入输出参数，执行后变量结果返回工作空间，可直接运行
        %以下是脚本文件，文件名假设为exp.m
        %**********************************************
            clear
            r = 5;
            s = pi*r*r;
            p = 2*pi*r;
            disp(s)
            disp(p)
        %**********************************************
        %以下是调用
        %**********************************************
            exp
        %**********************************************
    %函数文件：以function开头，有输入输出，变量为局部变量不返回工作空间，需要调用
        %以下是函数文件
        %**********************************************
            function [s, p] = circ(r)  % 文件命名应与函数名一致，系统找文件名circ.m
            %CIRC 计算圆面积和周长  % 简单说明
            %参数：输入参数r:圆半径；输出参数s:圆面积，p:周长  % 详细说明
                s = pi*r*r;
                p = 2*pi*r;
            end
        %**********************************************
        %以下是调用
        %**********************************************
            [a, b] = circ(5);  % 返回为多个参数时，若写a = circ(5)则保留第一个返回值
        %**********************************************
        %*以下是带子函数的函数文件
        %**********************************************
        function y = key(w)  % 主函数放第一个，函数名为key
            if w==0
                y = type0(w);  % 调用子函数type0
            else
                y = type1(w);
            end
        end
        function y0 = type0(a)  % 子函数，各子函数之间顺序无所谓
            y0 = a+1;
        end
        function y1 = type1(a)
            y1 = a+4;
        end
        %**********************************************
        %*函数输入输出参数可以不定
        %nargin：输入参数个数，nargout：输出参数个数
        %varargin：输入参数内容的元胞数组，varargout：输出参数
        %以下是函数文件
        %**********************************************
        function varargout = idk(varargin)
            x = length(varargin);
            varargout{1} = x;
            varargout{2} = x+1;
        end
        %**********************************************

第3部分：图像绘制

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
%%
x = 0:0.1:2*pi;
y1 = sin(x);
y2 = cos(x);

%%
% 二维曲线绘制
    % 基本函数
        % plot(y)
            % y为向量
                plot(y1);  % 纵坐标为y的值；横坐标自动为元素序号(角标+1)，此处为1~9
            % y为矩阵
                figure;  % 开启新绘图窗口，下一次绘图在新窗口
                y = [y1', y2'];
                plot(y);  % 当y为矩阵，按每一列画出曲线，颜色自动区分
        % plot(x, y)
            % xy为向量
                plot(x, y1);  % 先绘制曲线
        % plot(x1, y1, x2, y2...)
            plot(x, y1, x, y2);  % 在同一个窗口同一坐标轴绘制多条曲线
    % 线性图形格式设置
        % 线形颜色数据点
            plot(x, y1, 'b:o');  % 蓝色 点线 圆圈
            % b蓝 g绿 r红 c青 m紫 y黄 k黑 w白
            % -实线 :点线 --虚线 -.点画线
            % .实点 o圆圈 x叉 +十字 *星号 s方块 d钻石 v下三角 ^上三角 <左三角 >右三角 p五角星 h六角星
        % 坐标轴
            plot(x, y1);
            axis([-1*pi, 3*pi, -1.5, 1.5]);  % 规定横纵坐标范围
    % 图形修饰
        % 标题标签
            title('a title');  % 图像标题
            xlabel('this is x');  % x轴标记，同理还有ylabel，zlabel
        %图例设置
            legend('hahaha', 'location', 'best');  % str的顺序与绘图顺序一致; 'best'指图例位置最佳化，还有其他位置
        %图形保持
            plot(x, y1);
            hold on;  % 在原有窗口y1曲线上增加绘制下一个图形
            plot(x, y2);  % y2在同一窗口内被绘制
            hold off;
        %分割绘制
            subplot(2, 2, 1);  % 分割成2x2区域，在第一块区域绘制下一个图形
            plot(x, y1);  % y1被绘制在4块区域的第一块
            subplot(2, 2, 2);  % 分割方法相同，区域改变
            plot(x, y2);  % y2在第二块区域

%%
%*二维特殊图形绘制
    %柱状图
        bar(x, y, width, '参数');  % x横坐标向量，m个元素; y为向量时，每个x画一竖条共m条，矩阵mxn时，每个x画n条;
                                   % width宽度默认0.8，超过1各条会重叠;
                                   % 参数有grouped分组式，stacked堆栈式; 默认grouped
                                   % bar垂直柱状图,barh水平柱状图,bar3三维柱状图,barh3水平三维柱状图(三维多一个参数detached, 且为默认)
    %饼形图
        pie(x, explode, 'lable');  % x为向量显示每个元素占总和百分比, 为矩阵显示每个元素占所有总和百分比
                                   % explode向量与x同长度，为1表示该元素被分离突出显示，默认全0不分离
                                   % pie3绘制三维饼图
    %直方图
        hist(y, n);  % y为向量，把横坐标分为n段绘制
        hist(y, x);  % x为向量，用于指定每段中间值, 若取N = hist(y, x), N为每段元素个数
    %离散数据图
        stairs(x, y, 'b-o');  % 阶梯图，参数同plot
        stem(x, y, 'fill');  % 火柴杆图，参数fill是填充火柴杆，或定义线形
        candle(HI, LO, CL, OP);  % 蜡烛图:HI为最高价格向量,LO为最低价格向量,CL为收盘价格向量,OP为开盘价格向量
    %向量图
        compass(u, v, 'b-o');  % 罗盘图横坐标u纵坐标v
        compass(Z, 'b-o');  % 罗盘图复向量Z
        feather(u, v, 'b-o');  % 羽毛图横坐标u纵坐标v
        feather(Z, 'b-o');  % 羽毛图复向量Z
        quiver(x, y, u, v);  % 以(x, y)为起点(u, v)为终点向量场图
    %极坐标图
        % polar(theta, rho, 'b-o');  % 极角theta, 半径rho
        theta = -pi:0.01:pi;
        rho = sin(theta);
        polar(theta, rho, 'b')
    %对数坐标图
        semilogx(x1, y1, 'b-o');  % 把x轴对数刻度表示, semilogy是y轴对数刻度表示，loglog是两个坐标都用对数表示
    %双纵坐标
        plotyy(x1, y1, x2, y2, 'fun1', 'fun2');  % fun规定了两条条线的绘制方式，如plot,semilogx,semilogy,loglog,stem等
    %函数绘图
        f = 'sin(2*x)';
        ezplot(f, [0, 2*pi]);  % 绘制f并规定横坐标范围，也有[xmin, xmax, ymin, ymax]
        x = '2*cos(t)';
        y = '4*sin(t)';
        ezplot(x, y);  % 绘制x(t),y(t)在[0, 2*pi]图像, 也可以在最后用[tmin, tmax]规定t的范围

%%
%三维曲线曲面绘制
    %三维曲线
        x = 0:0.1:2*pi;
        y = sin(x); z = cos(x);
        plot3(x, y, z, 'b-*');
    %三维曲面
        %三维网格
            x = -5:0.1:5;  % 规定了x轴采样点，也规定了x轴范围
            y = -4:0.1:4;  % 规定了y轴采样点，也规定了y轴范围
            [X, Y] = meshgrid(x, y);  % 得到了xoy面网格点
            Z = X.^2+Y.^2;
            mesh(X, Y, Z)  % XY是meshgrid得到的网格点，Z是网格顶点，c是用色矩阵可省略
        %三维表面图
            x = -5:0.1:5;  
            y = -4:0.1:4;
            [X, Y] = meshgrid(x, y);
            Z = X.^2+Y.^2;  % 以上部分同上
            surf(X, Y, Z)  % 与上一个类似

第4部分：多项式

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
%%
%多项式
    %创建
        p = [1, 2, 3, 4];  % 系数向量，按x降幂排列，最右边是常数
        f1 = poly2str(p, 'x');  % 生成好看的字符串 f1 = x^3 + 2 x^2 + 3 x + 4，不被认可的运算式
        f2 = poly2sym(p);  % 生成可用的符号函数 f2 = x^3 + 2*x^2 + 3*x + 4
    %求值
        x = 4;
        y1 = polyval(p, x);  % 代入求值；若x1为矩阵，则对每个值单独求值
    %求根
        r = roots(p); % p同上，由系数求根，结果为根植矩阵
        p0 = poly(r);  % 由根求系数，结果为系数矩阵
        
%%
%数据插值
    %一维插值
        %yi = interp1(X, Y, xi, 'method')
        X = [-3, -1, 0, 1, 3];
        Y = [9, 1, 0, 1, 9];  % XY为已知点横纵坐标向量
        y2 = interp1(X, Y, 2);  % 差值预估在x=2的y的值，x不能超过已知范围(此处x<3)
        y2m = interp1(X, Y, 2, 'spline');  % 用spline方法(三次样条)差值预估在x=2的y的值
    %二维插值
        %zi = interp1(X, Y, Z, xi, yi, 'method')
        
%%
X = [2, 3, 9, 15, 6, 7, 4];
A = [1, 7, 2; 9, 5, 3; 8, 4 ,6];
B = [1, 7, 3; 9, 5, 3; 8, 4 ,6];
%数据统计
    %矩阵最大最小值
        y = max(X);  % 求矩阵X的最大值，min最小值
        [y, k] = max(X);  % 求最大值，k为该值的角标
        [y, k] = max(A, [], 2);  % A是矩阵，'2'时返回y每一行最大元素构成的列向量，k元素所在列；'1'时与上述相同
    %均值和中值
        y = mean(X);  % 均值
        y = median(X);  % 中值
        y = mean(A, 2);  % '2'时返回y每一行均值构成的列向量；'1'时与上述相同
        y = median(A, 2);  % '2'时返回y每一行中值构成的列向量；'1'时与上述相同
    %排序
        Y = sort(A, 1, 'ascend');  % sort(矩阵, dim, 'method')dim为1按列排序，2按行排序；ascend升序，descend降序
        [Y, I] = sort(A, 1, 'ascend');  % I保留了元素之前在A的位置
    %求和求积累加累乘
        y = sum(X);  % 求和
        y = prod(X);  % 求积
        y = cumsum(X);  % 累加
        y = cumprod(X);  % 累乘
        
%%
%*数值计算
    %最(极)值
        %多元函数在给定初值附近找最小值点
        x = fminsearch(fun, x0);
    %函数零点
        x = fzero(fun, x0);  % 在给定初值x0附近找零点

第5部分：符号函数

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
%%
%符号对象创建
    %sym函数
        p = sin(pi/3);
        P = sym(p, 'r');  % 用数值p创建符号常量P；'d'浮点数'f'有理分式的浮点数'e'有理数和误差'r'有理数
    %syms函数
        syms x;  % 声明符号变量
        f = 7*x^2 + 2*x+9;  % 创建符号函数
    %符号运算
        % 加减乘除外
        % '转置 ； ==相等 ； ~=不等
        % sin, cos, tan; asin, acos, atan 三角反三角
        % sinh, cosh, tanh; asinh, acosh, atanh 双曲反双曲
        % conj复数共轭；real复数实部；imag复数虚部；abs复数模；angle复数幅角
        % diag矩阵对角；triu矩阵上三角；tril矩阵下三角；inv逆矩阵；det行列式；rank秩；poly特征多项式；
        % |----expm矩阵指数函数；eig矩阵特征值和特征向量；svd奇异值分解；
    %符号对象精度转换
        digits;  % 显示当前用于计算的精度
        digits(16);  % 将计算精度改为16位，降低精度有时可以加快程序运算速度或减少空间占用
        a16 = vpa(sqrt(2));  % vpa括起的运算使sqrt(2)运算按照规定的精度执行
        a8 = vpa(sqrt(2), 8);  % 在vpa内控制精度，离开这一步精度恢复
        
%%
%符号多项式函数运算
    %*符号表达形式与相互转化
        %多项式展开整理
            g = expand(f);  % 展开
            h = collect(g);  % 整理(默认按x整理)
            h1 = collect(f, x);  % 按x整理（降幂排列）
        %因式分解展开质因数
            fac = factor(h);  % 因式分解
            factor(12);  % 对12分解质因数
    %符号多项式向量形式与计算
        syms a b c;
        n = [a, b, c];
        roots(n);  % 求符号多项式ax^2+bx+c的根
        n = [1, 2, 3];
        roots(n);  % 求符号多项式带入a=1, b=2, c=3的根
    %*反函数
        fi = finverse(f, x);  % 对f中的变量x求反函数
        
%%
%符号微积分
    %函数的极限和级数运算
        % 常量a，b
        %极限
            limit(f, x, 4);  % 求f(x), x->4
            limit(f, 4);  % 默认变量->4
            limit(f);  % 默认变量->0
            limit(f, x, 4, 'right');  % 求f(x), x->4+, 'left' x->4-
        %*基本级数运算
            %求和
                symsum(s, x, 3, 5);  % 计算表达式s变量x从3到5的级数和，或symsum(s, x, [a b])或symsum(s, x, [a;b])
                symsum(s, 3, 5);  % 计算s默认变量从3到5的级数和
                symsum(s);  % 计算s默认变量从0到n-1的级数和
            %一维泰勒展开
                taylor(f, x, 4);  % f在x=4处展开为五阶泰勒级数
                taylor(f, x);  % f在x=0处展开为五阶泰勒级数
                taylor(f);  % f在默认变量=0处展开为五阶泰勒级数
    %符号微分
        %单变量求导（单变量偏导）
            n = 1;  % 常量n
            fn = diff(f, x, n);  % f对x的n阶导
            f1 = diff(f, x);  % f对x的1阶导
            diff(f, n);  % f对默认变量的n阶导
            diff(f);  % 默认变量1阶导
        %多元偏导
            fxy = diff(f, x, y);  % 先求x偏导，再求y偏导
            fxyz = diff(f, x, y, z);  % 先求x偏导，再求y偏导,再求z偏导
    %符号积分
        %积分命令
            int(f, x, 1, 2);  % 函数f变量x在1~2区间定积分
            int(f, 1, 2);  % 函数f默认变量在ab区间定积分
            int(f, x);  % 函数f变量x不定积分
            int(f);  % 函数f默认变量不定积分
            % 傅里叶，拉普拉斯，Z变换
            
%%
%*符号方程求解
    %符号代数方程
        %一元方程
            eqn1 = a*x==b;
            S = solve(eqn1);  % 返回eqn符号解
        %多元方程组
            eqn21 = x-y==a;
            eqn22 = 2*x+y==b;
            [Sx, Sy] = solve(eqn21, eqn22, x, y);  % [Svar1,...SvarN]=solve(eqn1,...eqnM, var1,...varN),MN不一定相等
            [Sxn, Syn] = solve(eqn21, eqn22, x, y, 'ReturnCondition', true);  % 加上参数ReturnCondition可返回通解及解的条件
            % 其他参数(参数加上true生效)
                % IgnoreProperty，忽略变量定义时一些假设
                % IgnoreAnalyticConstraints，忽略分析限制；
                % MaxDegree，大于3解显性解；
                % PrincipleValue，仅主值
                % Real，仅实数解
        %非线性fsolve
            X = fsolve(fun, X0, optimset(option));  % fun函数.m文件名；X0求根初值；option选项如('Display','off')不显示中间结果等

作者：glanny

原文地址

出处： bilibili

坚持原创技术分享，您的支持将鼓励我继续创作！