matleb的使用

最新推荐文章于 2024-05-13 21:32:22 发布

HKTK-CWW

最新推荐文章于 2024-05-13 21:32:22 发布

阅读量649

点赞数

分类专栏： linux驱动嵌入式文章标签：线性代数矩阵算法

本文链接：https://blog.csdn.net/m0_54657729/article/details/129354122

版权

嵌入式同时被 2 个专栏收录

80 篇文章 0 订阅

订阅专栏

linux驱动

69 篇文章 3 订阅

订阅专栏

一、矩阵的运算与方程的求解

rem等同于mod
rem(5,3)=mod(5,3)=2

%重复单引号来输入含有单引号的字符串

randn函数：产生均值为0、方差为1的标准正态分布随
机数x

o 对角阵：只有对角线上有非零元素的矩阵
o 数量矩阵：对角线上的元素相等的对角矩阵
o 单位矩阵：对角线上的元素都为1的对角矩阵

diag既可以产生列向量也能根据列向量产生对角阵

diag(A,k)：提取第k条对角线的元素，产生一个列向量

diag(V)：产生一个m×m对角矩阵，其主对角线元素即为向量
V的元素。
diag(V,k)：以向量V为第k条对角线元素，产生对角矩阵

6.三角阵 (矩阵的对角线以上/下的元素全为零的矩阵)

triu(A)/tril(A)
提取矩阵A的主对角线及以上/下的元素

7.矩阵的转置

A.'转置
A'共轭转置
若矩阵的元素是实数，
则转置和共轭转置的结
果是一样的

8.矩阵旋转

rot90(A,k)：将矩阵A按逆时针方向旋转90º
的k倍，当k为1时可省略。

9.矩阵翻转

fliplr(A)：对矩阵A实施左右翻转
flipud(A)：对矩阵A实施上下翻转

10.求逆解线性方程组

在这里插入图片描述

A=[2 -1 3;3 1 -5;4 -1 1]
B=[5,5,9]’
X=inv(A)*B
或者X=A\B

11.求行列式值(要求是方阵)!

det（A）:求方阵A所对应的行列式值

12.求矩阵特征值与向量

[X,D] = eig(A)
如果只接收一个结果那么就是特征值

在这里插入图片描述

13.求解多项式方程 (3x5-7x4+5x2+2x-18=0)

1.求其伴随矩阵的特征值
p=[3,-7,0,5,2,-18];
A=compan(p); %A的伴随矩阵
x1=eig(A) %求A的特征值
2.直接求根
X2 = roots(p);

二、程序设计

1. 在这里插入图片描述

2. 在这里插入图片描述

在这里插入图片描述

4. 在这里插入图片描述

5. 在这里插入图片描述

6.函数

1.匿名函数

f = @(x,y) x^2+y^2
f(3,4)

2.全局变量

global 变量

四、二维绘图

1.二维绘图(plot)

绘制多条二维曲线
1.
x=0:pi/100:2*pi;
y1=sin(x);
y2=cos(x);
plot(x,y1,x,y2)
2.
x=0:pi/100:2*pi;
y=[sin(x);cos(x)];
plot(x,y)

在这里插入图片描述

2.自适应绘图（fplot）

绘制f(x)=cos(tan(πx))
fplot(@cos(tan(pi*x)),[ 0,1])

在这里插入图片描述

3.隐函数绘图

对于参数方程x = x(t)和y = y(t)，ezplot函数的调用
格式为：
ezplot(x,y)：在默认区间0<t<2π绘制x=x(t)和y=y(t)的
图形。
ezplot(x,y, [tmin,tmax])：在区间tmin < t < tmax绘制
x=x(t)和y=y(t)的图形。

3.图形标注

在这里插入图片描述

grid（设置网格命令）
legend(string1,string2,string3…)
gtext(‘string’)

4.控制坐标

axis([xmin xmax ymin ymax zmin zmax])
axis equal 纵、横坐标轴采用等长刻度
axis squar 产生正方形坐标系(缺省为矩形)

5.分割窗口

subplot(m,n,p)

6.新建窗口

xx = figure;

7.其它坐标下绘图

1.半对数与对数
semilogx(x1,y1,选项1,x2,y2,选项2,…)
semilogy(x1,y1,选项1,x2,y2,选项2,…)
loglog(x1,y1,选项1,x2,y2,选项2,…)
2.极坐标
polar(theta,rho,选项)

8.统计图

bar(x,y,选项) 条形图
stairs(x,y,选项) 阶梯图
stem(x,y,选项) 杆图
fill(x1,y1,选项1,x2,y2,选项2,…) 填充图色 :选项为基本颜色向量[r g b]表示

五、三维绘图

1.三维曲线

1.plot3

plot3(x1,y1,z1,x2,y2,z2,…,xn,yn,zn)

2.fplot3

fplot3(funx,funy,funz,tlims)

在这里插入图片描述

2.三维曲面

步骤：

1.先用meshgrid创建网格

2.再用mesh或surf

t = -2:0.2:2;
[X, Y] = meshgrid(t);
Z = X .* exp(-X.^2 - Y.^2);
mesh(X,Y,Z);

标准三维曲面

1.球面(R=1)

1.获取x,y,z
[x,y,z]=sphere(n) n为光滑程度
2.调用mesh或surf

2.柱面

[x,y,z]= cylinder(R,n)

三维曲面隐函数

funx = @(u,v) u.*sin(v);
funy = @(u,v) -u.*cos(v);
funz = @(u,v) v;
fsurf(funx,funy,funz,[-5 5 -5 -2])
hold on
fmesh(funx,funy,funz,[-5 5 -2 2])
hold of

三维统计图

在这里插入图片描述

六、数值计算

1.max/min

1.向量
[y,I]=max(X) y为最大值，i为最大值第一次出现的下标
2.矩阵
	1.求矩阵最大值
	max(A)返回一个行向量，为每列最大的元素
	max((max(A))) 能求出矩阵最大值
	2.两个矩阵最大值合并
	U=max(A,B)
	U的每个元素等于A,B对应元素的较大者

2.取平均值与中值

① mean(X)：返回向量X的算术平均值。
② median(X)：返回向量X的中值。
③ mean(A)：返回一个行向量，其第i个元素是A的第i列的
算术平均值。
④ median(A)：返回一个行向量，其第i个元素是A的第i列
的中值。

3.累加与累乘向量

cumsum(x)
cumprod(x)

4.求标准方差

std(x)
std(A)

5.求相关系数

corrcoef(A)
corrcoef(X,Y)=corrcoef([X,Y])

6.排序

[Y,I]=sort(A,dim,mode)
若dim=1，则按列排；若dim=2，则按行排。
mode指明按升序还是降 ‘descend’表示降序
默认:按列升序

7.数据插值

Y1=interp1(X,Y,X1,'method')
method是插值方法，允许的取值有‘linear’（线
性）、
‘nearest’（最近点）、
‘pchip’（3次多项
式）、
‘spline’（3次样条）

不同方法比较：
在这里插入图片描述

二维插值同理：

Z1=interp2(X,Y,Z,X1,Y1,'method')

8.曲线拟合

1.取得要以多少次最小二乘拟合的系数
[P,S]=polyfit(X,Y,m) m为次数 ,P为多项式系数
2.按照所得的系数进行逼近
Y2=polyval(P,X1)

案例：

X=linspace(0,2*pi,50);
Y=sin(X);
P=polyfit(X,Y,3)
X1=linspace(0,2*pi,20);
Y1=sin(X1);
Y2=polyval(P,X1)
plot(X1,Y1,’:o’,X1,Y2,’-*’)

9.多项式计算

1.乘除

[Q,r]=deconv(P1,P2)用于对多项式P1和P2作除法运算
Q为商，r为余数
deconv是conv的逆函数，即有P1=conv(P2,Q)+r。

2.多项式求导

p=polyder(P)：求多项式P的导函数
p=polyder(P,Q)：求P·Q的导函数
[p,q]=polyder(P,Q)：求P/Q的导函数，导函数的分子存入p，分母存入q。

10.多项式求值

1.以每个元素待人
polyval(P,x) P为多项式系数
2.以矩阵整体待人
polyvalm(P,x) 要求x为方阵

11.多项式求根与用根求多项式

P=[3,0,4,-5,-7.2,5];
X=roots(P)
G=poly(X)

12.求单变量非线性方程的根

z=fzero('fname',x0,tol,trace)
x0为搜索的起点。一个函数可能有多个根，但fzero函数只给出离x0最近的那个根。

案例：

求f(x)=x-10x+2=0在x0=0.5附近的根。
(1) 建立函数文件funx.m。
function fx=funx(x) 　　
	fx=x-10.^x+2;
(2) 调用fzero函数求根。
z=fzero('funx',0.5)

13。求解非线性方程组

fsolve

七、符号计算

创建符号变量

1.创建单个变量name = sym('name')
2.创建多个变量syms name1 name2 .

基本符号运算：

factor(S) —— 对S分解因式
expand(S) ——对S进行展开
collect(S) ——对S合并同类项
collect(S,v) —— 对S按变量v合并同类项
simplify(S) ——应用函数规则对S进行化简

符号求极限：

在这里插入图片描述

limit(f,x,a,'right')与limit(f,x,a,'left')

符号求导：

diff(f,x,n) 求函数f对变量x的n阶导数

符号求积分：

不定积分：
	int(f,x) 以x为自变量，对被积函数或符号表达式f求不定积分
定积分：
	int(f,x,a,b) 定积分的下限和上限

符号级数求和：

symsum(s,v,n0,nn) s为通项公式 v为变量 n0-nn

符号方程求解:

代数方程:
	solve(eq,v)求解符号表达式表示的代数方程eq，求解变量为v。
常微分方程：
	Dy y’
	D2y y’’
	Dy(5) y’(5)=0
	y’’’+y’’+y’-x+5=0
	D3y+D2y+Dy-x+5=0
	dsolve('eqn1','condition','var')

例题：

在这里插入图片描述