MATLAB代码参考

1小时Matlab速成代码资料

学习 2020-10-15 21:213.7万阅读 · 1307喜欢 · 65评论

glanny

粉丝：7264文章：4

关注

由于需要Matlab代码的同学太多，以致于资源群日常满员。同时，因为部分同学表示github使用不方便，因此将代码直接放在专栏投稿，以便大家查阅。转载请注明出处。

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

%%
% 建议有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部分：程序结构

%%
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 = pirr;
p = 2pir;
disp(s)
disp(p)
%**********************************************
%以下是调用
%**********************************************
exp
%**********************************************
%函数文件：以function开头，有输入输出，变量为局部变量不返回工作空间，需要调用
%以下是函数文件
%**********************************************
function [s, p] = circ(r) % 文件命名应与函数名一致，系统找文件名circ.m
%CIRC 计算圆面积和周长 % 简单说明
%参数：输入参数r:圆半径；输出参数s:圆面积，p:周长 % 详细说明
s = pirr;
p = 2pir;
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部分：图像绘制

%%
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([-1pi, 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(2x)’;
ezplot(f, [0, 2pi]); % 绘制f并规定横坐标范围，也有[xmin, xmax, ymin, ymax]
x = ‘2cos(t)’;
y = ‘4sin(t)’;
ezplot(x, y); % 绘制x(t),y(t)在[0, 2pi]图像, 也可以在最后用[tmin, tmax]规定t的范围

%%
%三维曲线曲面绘制
%三维曲线
x = 0:0.1:2pi;
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部分：多项式

%%
%多项式
%创建
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 + 2x^2 + 3x + 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部分：符号函数

%%
%符号对象创建
%sym函数
p = sin(pi/3);
P = sym(p, ‘r’); % 用数值p创建符号常量P；’d’浮点数’f’有理分式的浮点数’e’有理数和误差’r’有理数
%syms函数
syms x; % 声明符号变量
f = 7x^2 + 2x+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 = ax==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’)不显示中间结果等

转载请注明出处

本文为我原创本文禁止转载或摘编

速成 Matlab