MATLAB基本操作
1. 对象定义 使用sym定义单个对象、使用syms定义多个对象
2. 使用limit
求极限
$$ \lim_{v \rightarrow a} f(x) $$
limit(f,v,a) % 使用limit(f,v,a,'left')可求左极限
3. 导数 使用diff(f,v,n)
对$ f(v)=v^{t-1} $求 $ n $ 阶导 $ \frac{d^nf}{d^nv} $,n缺省时,默认为1,diff(f)默认求一阶导数。
4. 定积分和不定积分 使用int(f,v)
求f对变量v的不定积分,使用int(f,v,a,b)
求f对变量v的定积分,a、b为积分上下标。$ \int{f(v)dv} $、$ \int^{a}_{b}{f(v)dv} $。
5. matlab函数文件定义形式
function [输出形参列表] = 函数名(输入形参列表)
函数体
function spir_len = spirallength(d, n, lcolor)
% SPIRALLENGTH plot a circle of radius as r in the provided color and calculate its area
% 输入参数:
% d: 螺旋的旋距
% n: 螺旋的圈数
% lcolor:画图线的颜色
% 输出参数:
% spir_len:螺旋的周长
% 调用说明:
% spirallength(d,n):以参数d,n画螺旋线,螺旋线默认为蓝色
% spirallength(d,n,lcolor):以参数d,n,lcolor画螺旋线
% spir_len = spirallength(d,n):计算螺旋线的周长,并以蓝色填充螺旋线
% spir_len = spirallength(d,n,lcolor):计算螺旋线的周长,并以lcolor颜色填充螺旋线
% 版本号V1.0,编写于1999年9月9号,修改于1999年9月10号,作者:亚索
if nargin > 3
error('输入变量过多!');
elseif nargin == 2
lcolor = 'b'; % 默认情况下为蓝色
end
j = sqrt(-1);
phi = 0 : pi/1000 : n*2*pi;
amp = 0 : d/2000 : n*d;
spir = amp .* exp(j*phi);
if nargout == 1
spir_len = sum(abs(diff(spir)));
fill(real(spir), imag(spir), lcolor);
elseif nargout == 0
plot(spir, lcolor);
else
error('输出变量过多!');
end
axis('square');
6. matlab程序设计语句
% for循环
for 循环变量=初值:步长:终值
循环体
end
% while循环
while 条件
循环体
end
% if语句
if 条件
语句组1
elseif
语句组2
else
语句组3
end
% switch语句
switch 表达式
case 表达式1
语句组1
case 表达式2
语句组2
... ...
case 表达式m
语句组m
otherwise
语句组
end
% try语句
try
语句组1 %语句组1若正确则跳出该结构
catch
语句组2
end
7. 矩阵操作
操作 | 作用 |
---|---|
size(A) | 求矩阵A的行数和列数 |
length(x) | 返回向量x的长度 |
A' | A的转置 |
A(:,n) | 取矩阵A第n列数,A(n,:)取第n行 |
det(A) | 求矩阵A的行列式 |
inv(A) | 求A的逆 |
rank(A) | 求A的秩 |
trace(A) | 求A的迹 |
max(A)、min(A) | 求A的各列最大、最小元素 |
mean(A) | 求A各列的平均值 |
sum(A) | 求A各列元素之和 |
8. matlab简单绘图
plot函数是MATLAB中最核心的二维绘图函数,有诸多语法格式,可实现多种功能。常用格式有:
- plot(x):缺省自变量的绘图格式,x可为向量或矩阵。
- plot(x, y):基本格式,x和y可为向量或矩阵。
- plot(x1, y1, x2, y2,…):多条曲线绘图格式,在同一坐标系中绘制多个图形。
- plot(x, y,‘s’):开关格式,开关量字符串s设定了图形曲线的颜色、线型及标示符号(见下表)。
无约束优化问题求解
fminbnd、fminunc函数输出变量解释
变量 | 描述 |
---|---|
x | 由优化函数求得的值. 若exitflag>0,则x为解; 否则,x不是最终解, 它只是迭代制止时优化过程的值 |
fval | 解 x 处的目标函数值 |
exitflag | 描述退出条件:exitflag>0,表目标函数收敛于解x处;exitflag=0,表已达到函数评价或迭代的最大次数;exitflag<0,表目标函数不收敛 |
output | 包含优化结果信息的输出结构。Iterations:迭代次数;Algorithm:所采用的算法;FuncCount:函数评价次数 |
一元函数无约束优化问题-fminbnd
常用格式
$$ min f(x), x_1<x<x_2 $$
(1)x= fminbnd (fun, x1, x2) (2)x= fminbnd (fun, x1, x2 , options) (3)[x , fval]= fminbnd(...) (4)[x , fval , exitflag]= fminbnd(...) (5)[x , fval , exitflag , output]= fminbnd(...) 函数fminbnd的算法基于黄金分割法和二次插值法,它要求目标函数必须是连续函数,并可能只给出局部最优解
例子
求函数 $ f(x)=2e^{-x}sin(x) $ 在 $ 0<x<8 $ 时的最小值
% 如果求最大需要对f取反
f = @(x) (2*exp(-x)*sin(x));
[x,fval] = fminbnd(f,0,8);
x
fval
多元函数无约束优化问题-fminunc
常用格式
$$ min f(X),这里X为n维变量 $$ fminunc常用格式为: (1)x= fminunc(fun, X0); (2)x= fminunc(fun, X0,options); (3)[x,fval]= fminunc(...); (4)[x,fval,exitflag]= fminunc(...); (5)[x,fval,exitflag,output]= fminunc(...) 其中 X0为初始值
例子
求函数$ f(x_1,x_2)=(4x_1^2+2x_2^2+4x_1x_2+2x_2^2+1)e^x $的最小值,$ X_0=[-1,1] $
f = @(x) (4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1)*exp(x(1));
x0 = [-1,1];
[x,fval] = fminunc(f, x0);
x
fval
线性规划问题求解
使用linprog求解一般线性规划问题
常见问题(linprog默认求最小值) $$ minz=cX $$
$$ s.t. \begin{cases} AX\leq{b}\ Aeq\cdot{X}=beq\ VLB\leq{X}\leq{VUB} \end{cases}$$
求解命令
[x,fval] = linprog(c,A,b,Aeq,beq,VLB,VUB)
例子
$$ min z=13x_1+9x_2+10x_3+11x_4+12x_5+8x_6 $$
$$ s.t.\left{ \begin{aligned} & x_1+x_2=400\ & x_2+x_5=600\ & x_3+x_6=500\ & 0.4x_1+1.1x_2+x_3\leq{800}\ & 0.5x_4+1.2x_5+1.3x_6\leq{900}\ & x_i\geq0,i=1,2,...,6 \end{aligned} \right. $$
f = [13 9 10 11 12 8];
A = [0.4 1.1 1 0 0 0
0 0 0 0.5 1.2 1.3];
b = [800; 900];
Aeq=[1 0 0 1 0 0
0 1 0 0 1 0
0 0 1 0 0 1];
beq=[400 600 500];
vlb = zeros(6,1);
vub=[];
[x,fval] = linprog(f,A,b,Aeq,beq,vlb,vub)
使用bintprog求解0-1规划问题
matlab2014以上版本使用intlinprog
求解0-1规划问题
$$ minz=cX $$
$$ s.t. \begin{cases} AX\leq{b}\ Aeq\cdot{X}=beq\ X为0-1变量 \end{cases}$$
% 命令
[x,fval] = bintprog(c,A,b,Aeq,beq)
例子
$$ min z=3x_1+7x_2-x_3+x_4 $$ $$ s.t. \begin{cases} 2x_1-x_2+x_3-x_4\geq{1}\ x_1-x_2+6x_3+4x_4\geq{8}\ 5x_1+3x_2+x_4\geq{5}\ x_i=0或1(i=1,2,3,4) \end{cases} $$
z = [3;7;-1;1];
A = [-2 1 -1 1;
-1 1 -6 -4;
-5 -3 0 -1];
b = [-1;-8;-5];
Aeq = [];
beq = [];
[x,fval] = bintprog(z,A,b,Aeq,beq)
数据插值与拟合
数据插值,使用interpl进行一维插值
matlab命令
yi = interpl(X,Y,xi,method)
该命令用指定的算法找出一个一元函数,然后以该函数给出xi处的值。其中x=[x1,x2,…,xn]’和 y=[y1,y2,…,yn]’两个向量分别为给定的一组自变量和函数值,用来表示已知样本点数据;xi为待求插值点处横坐标,可以是一个标量,也可以是一个向量,是向量时,必须单调;yi得到返回的对应纵坐标。
- method可以选取以下方法之一:
- ‘nearest’:最近邻点插值,直接完成计算;
- ‘spline’:三次样条函数插值;
- ‘linear’:线性插值(缺省方式),直接完成计算;
- ‘cubic’:三次函数插值;
例子
作函数$ y=(x^2-3x+7)e^{-4x}sin(2x) $在[0,1]取间隔为0.1的点图,用插值进行实验
x=0:0.1:1;
y=(x.^2-3*x+7).*exp(-4*x).*sin(2*x); %产生原始数据
subplot(1,2,1);
plot(x,y,x,y,'ro') %作图
xx=0:0.02:1; %待求插值点
yy=interp1(x,y,xx,'spline'); %此处可用nearest,cubic,spline分别试验
subplot(1,2,2)
plot(x,y,'ro',xx,yy,'b') %作图
曲线拟合
拟合函数polyfit
p=polyfit(x,y,n)
[p,s]= polyfit(x,y,n)
说明:x,y为数据点,n为多项式阶数,返回p为幂次从高到低的多项式系数向量p。p是n+1维参数向量p(1),p(2)….那么拟合后对应的多项式即为: $$ p(1)x^n+p(2)x^{n-1}+\cdot\cdot\cdot+p(n)x+p(n+1) $$
x必须是单调的。矩阵s用于生成预测值的误差估计
多项式求值函数polyval
y=polyval(p,x)
[y,DELTA]=polyval(p,x,s)
说明:y=polyval(p,x)为返回对应自变量x在给定系数p的多项式的值; [y,DELTA]=polyval(p,x,s) 使用polyfit函数的选项输出s得出误差估计DELTA。它假设polyfit函数数据输入的误差是独立正态的,并且方差为常数。则DELTA将至少包含50%的预测值。
例子
求如下给定数据的拟合曲线 x=[0.5,1.0,1.5,2.0,2.5,3.0],y=[1.75,2.45,3.81,4.80,7.00,8.60]文章来源:https://www.toymoban.com/news/detail-838619.html
x=[0.5,1.0,1.5,2.0,2.5,3.0];
y=[1.75,2.45,3.81,4.80,7.00,8.60];
plot(x,y,‘*r’) %先观察数据点的大致形态
p=polyfit(x,y,2) %用二次多项式拟合
x1=0.5:0.05:3.0; % 步长0.05
y1=polyval(p,x1);
plot(x,y,'*r',x1,y1,'-b')
本文由博客群发一文多发等运营工具平台 OpenWrite 发布文章来源地址https://www.toymoban.com/news/detail-838619.html
到了这里,关于matlab 基操~的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!