线性代数

调用格式

目的 命令
求方阵A的行列式 det(A)
求矩阵A的秩 rank(A)
求矩阵A的逆矩阵 inv(A)
求矩阵A的行最简行 rref(A)
求特征值和特征向量 [V,D]=eig(A)
求模最大的特征值和特征向量 [V,D]=eigs[A,1]

例题


MATLAB程序实现:

clc, clear, syms s
a=[3*s,-2;-6,s]; b=[4;1];
c=det(a) %计算系数矩阵的行列式
s0=solve(c) %求行列式等于零的点
x=inv(a)*b %当 a 可逆时,求方程组的唯一解
x=simplify(x) %对符号解进行化简

其输出为:

c =3*s^2 - 12s0 =-22x =(4*s)/(3*(s^2 - 4)) + 2/(3*(s^2 - 4))s/(s^2 - 4) + 8/(s^2 - 4)x =(2*(2*s + 1))/(3*(s^2 - 4))(s + 8)/(s^2 - 4)

符号微积分

极限

调用格式

命令 效果
limit(f,x,a) 计算当自变量 x 趋近于常数 a 时,符号函数 f(x)的极限值.
limit(f,x,a,‘left’) 计算当 x 从左侧趋近于 a 时,符号函数 f(x)的左极限值.
limit(f,x,a,‘right’) 计算当 x 从右侧趋近于 a 时,符号函数 f(x)的右极限值.

例题


MATLAB程序实现:

clc, clear, syms x
a=limit(sin(x)/x,x,0) %求 x 趋近于 0 时的极限
b=limit((1+1/(2*x))^x,x,inf)

输出结果为:

a =1b =exp(1/2)

微分

在 MATLAB 中,采用函数 diff 进行符号函数求导运算,gradient 求多元函数的梯度,利用函数 jacobian 求 Jacobian 矩阵等.

调用格式

命令 效果
diff(expr) 求符号表达式 expr 对第 1 个符号变量(symvar(expr,1)确定)的一阶导数.
diff(expr,n) 求符号表达式 expr 对第 1 个符号变量的 n 阶导数.
diff(expr,v,n) 求符号表达式 expr 对符号变量 v 的 n 阶导数.
gradient(f) 求多元函数f的梯度
hessian(f) 求f的Hessian 矩阵
jacobian(F,V) 计算向量函数F关于V的Jacobian矩阵,当F和V为标量时,jacobian(F,V)等价于 diff(F,V).

例题1


MATLAB程序实现:

clc, clear, syms x y
f(x,y)=x^3+x*y+y^2;
dx=diff(f) %求 f 关于 x 的一阶导数
d2x=diff(f,x,2) %求 f 关于 x 的二阶导数
d2y=diff(f,y,2) %求 f 关于 y 的二阶导数

输出结果:

dx(x, y) =3*x^2 + yd2x(x, y) =6*xd2y(x, y) =2

例题2


MATLAB实现:

clc, clear, syms y(x) dy
eq1=y^5+2*y-x-3*x^7;
eq2=diff(eq1,x) %求关于 x 的一阶导数
eq3=subs(eq2,diff(y(x), x),dy) %把 diff(y(x),x)替换为 dy,否则无法解代数方程
dy2=solve(eq3,dy) %解代数方程,求 y 关于 x 的导数
dy3=subs(dy2,{x,y(x)},{0,0}) %代入具体的数值

输出结果为:

eq2(x) =5*y(x)^4*diff(y(x), x) + 2*diff(y(x), x) - 21*x^6 - 1eq3(x) =2*dy + 5*dy*y(x)^4 - 21*x^6 - 1dy2 =(21*x^6 + 1)/(5*y(x)^4 + 2)dy3 =1/2

例题3


MATLAB实现:

clc, clear, syms x y
f(x,y)=exp(x)*sin(y)+x^2+x*cos(y);
gradf=gradient(f)
Hf=hessian(f)

输出结果为:

gradf(x, y) =2*x + cos(y) + exp(x)*sin(y)exp(x)*cos(y) - x*sin(y)Hf(x, y) =[     exp(x)*sin(y) + 2,     exp(x)*cos(y) - sin(y)]
[exp(x)*cos(y) - sin(y), - exp(x)*sin(y) - x*cos(y)]

例题4


MATLAB实现:

clc, clear, syms x y
f(x,y)=exp(x)*sin(y)+x^2+x*cos(y);
gradf=gradient(f)
Hf1=jacobian(gradf) %求梯度向量的 Jacobian 阵
Hf2=hessian(f) %求 f 的 Hessian 阵

输出结果为:

gradf(x, y) =2*x + cos(y) + exp(x)*sin(y)exp(x)*cos(y) - x*sin(y)Hf1(x, y) =[     exp(x)*sin(y) + 2,     exp(x)*cos(y) - sin(y)]
[exp(x)*cos(y) - sin(y), - exp(x)*sin(y) - x*cos(y)]Hf2(x, y) =[     exp(x)*sin(y) + 2,     exp(x)*cos(y) - sin(y)]
[exp(x)*cos(y) - sin(y), - exp(x)*sin(y) - x*cos(y)]

积分

在 MATLAB 中,提供了 int 函数计算符号表达式的不定积分和定积分,函数的调用格式为:

调用格式

命令 效果
int(expr,v) 求符号表达式 expr 关于符号变量 v 的不定积分.
int(expr,v,a,b) 求符号表达式 expr 关于 v 的定积分,积分区间为[a,b].

例题


MATLAB实现如下:

clc, clear, syms x
I1=int(x*exp(x)/(1+exp(x))^2) %求符号不定积分
pretty(I1) %书写习惯的显示方式
I2=int(x^2/(x^4+x^2+1),x,0,inf) %求符号定积分
pretty(I2) %书写习惯的显示方式
I2=double(I2) %把符号数据转换为 double 类型数据

输出结果如下所示:


I1 =(x*exp(x))/(exp(x) + 1) - log(exp(x) + 1)x exp(x)
---------- - log(exp(x) + 1)
exp(x) + 1I2 =(pi*3^(1/2))/6pi sqrt(3)
----------6I2 =0.9069

MATLAB 在线性代数和高等数学中的应用相关推荐

  1. 从零开始学AI(数学基础之线性代数和高等数学)

    这周又开了一次线下课程,这次一点代码没讲,只是讲了一下数学基础,讲的内容很简单,听着一点也不吃力,都是大学的数学,内容也比较少,自我感觉吧,这点数学基础肯定是不够的,深度学习的东西虽然我还没入门,但是 ...

  2. MATLAB(四)在高等数学中的应用

    MATLAB在高等数学中的应用 (一)求极限 matlab中求极限的命令为 limit(expr, x, a) limit(expr, a) limit(expr) limit(expr, x, a, ...

  3. 用Matlab求解高等数学中的问题(不定积分,定积分,多元微积分,微分方程,级数问题)

    用Matlab求解高等数学中的问题 一.求解不定积分 利用matlab符号工具箱中的求积函数int, 可求函数的不定积分,int函数的调用格式如下: int(S) int(S,v) 说明: (1)in ...

  4. matlab选修结课作业,matlab在高等数学中的应用结课作业

    matlab在高等数学中的应用结课作业 MatLab软件在作图中的应用 1.描点作图(MatLab及在电子信息课程中的应用131页) a=0:0.01:3; >> y=a: plot(y) ...

  5. matlab语法手册下载,MATLAB及其在理工课程中的应用指南 第4版.pdf

    作 者 :陈怀琛 出版发行 : 西安:西安电子科技大学出版社 , 2018.02 ISBN号 :978-7-5606-4737-1 页 数 : 262 原书定价 : 39.00 开本 : 16开 主题 ...

  6. 计算机编程概率论,程序员的数学:线性代数和概率统计

    程序员的数学:线性代数和概率统计 语音 编辑 锁定 讨论 上传视频 <程序员的数学:线性代数和概率统计>是2018年9月1日机械工业出版社出版的图书,作者是[美] 欧内斯特 戴维斯 .[1 ...

  7. 【数字信号处理】线性常系数差分方程 ( 使用 matlab 求解 “ 线性常系数差分方程 “ 示例 | A 向量分析 | B 向量分析 | 输入序列分析 | matlab 代码 )

    文章目录 一.使用 matlab 求解 " 线性常系数差分方程 " 示例 1.B 向量元素 : x(n) 参数 2.A 向量元素 : y(n) 参数 3.输入序列 4.matlab ...

  8. 【数字信号处理】线性常系数差分方程 ( 卷积 与 “ 线性常系数差分方程 “ | 使用 matlab 求解 “ 线性常系数差分方程 “ )

    文章目录 一.卷积 与 " 线性常系数差分方程 " 二.使用 matlab 求解 " 线性常系数差分方程 " 一.卷积 与 " 线性常系数差分方程 & ...

  9. Matlab实现 线性动态电路可视化分析

    Matlab实现 线性动态电路可视化分析 版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明. 编程思路 这个编程总体采用面向过程的编程思想,将人在 ...

  10. matlab在电力电子的应用,MATLAB在电力电子仿真中的应用

    MATLAB在电力电子仿真中的应用 MATLAB在电力电子仿真中的应用 田伟 于海生 (青岛大学自动化工程学院,青岛 266071) 摘 要:MATLAB 软件是当今控制系统的设计与仿真中重要的工具软 ...

最新文章

  1. 强烈不建议你用 a.equals(b) 判断对象相等!
  2. 用PyCharm Profile分析异步爬虫效率
  3. 这些产品基本功了解一下?
  4. 【小白学习C++ 教程】九、C++中字符型、字符串和转义字符
  5. 京东AI NLP高阶实训学习笔记
  6. 虚幻引擎发布4.8版,支持VR虚拟设备
  7. Linux目录的可写意味着,Linux权限分析 - osc_h5427nyq的个人空间 - OSCHINA - 中文开源技术交流社区...
  8. Windows 平台下的C++代理类(供C # 使用)
  9. python编程入门经典-总算理解python编程入门经典教程
  10. 电商设计素材:15000款淘宝促销水印元素模版,手慢的人领不到!...
  11. YUV420转RGB888
  12. 使用Windows命令行reg控制注册表键值
  13. CorelDRAW教程大全集
  14. 全国计算机考试等级考务管理系统,全国计算机等级考试考务管理系统:https://ncre-bm.neea.edu.cn/...
  15. “智多星”手机销售网
  16. POCO C++ 在IOS上的使用
  17. 阿里云oss文件服务器
  18. torch.diag() 取矩阵对角线元素,torch.diag_embed() 指定值变成对角矩阵
  19. 大数据征信助力金融科技腾飞
  20. 浪潮信息做5G服务器,看来浪潮和运营商在5G布局非常深,中标移动开了个好头...

热门文章

  1. 动物识别系统 c语言代码_C++程序员避不开虚函数的,就像C语言程序员避不开指针一样...
  2. C语言编程QQ管理系统,c语言制作学生管理系统srrpqq67.doc
  3. oracle永久表设置,Oracle 表空间简单管理永久表空间
  4. c语言linux系统宏,Linux下C语言中的预定义宏
  5. mysql分日汇总填充0_将记录按小时或按日分组,并在mysql中填充零或空的空白
  6. java泛型要用t,Java泛型中T和問號(通配符)的區別
  7. curlopt_ssl_verifypeer后https还是验证不过_2020CNAS能力验证小总结
  8. VS2015sql本地服务器为空,详解VS2015自带LocalDB数据库用法实例
  9. Angr安装与使用之使用篇(八)
  10. java程序编写需注意的问题