MATLAB 在线性代数和高等数学中的应用
线性代数
调用格式
目的 | 命令 |
---|---|
求方阵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 在线性代数和高等数学中的应用相关推荐
- 从零开始学AI(数学基础之线性代数和高等数学)
这周又开了一次线下课程,这次一点代码没讲,只是讲了一下数学基础,讲的内容很简单,听着一点也不吃力,都是大学的数学,内容也比较少,自我感觉吧,这点数学基础肯定是不够的,深度学习的东西虽然我还没入门,但是 ...
- MATLAB(四)在高等数学中的应用
MATLAB在高等数学中的应用 (一)求极限 matlab中求极限的命令为 limit(expr, x, a) limit(expr, a) limit(expr) limit(expr, x, a, ...
- 用Matlab求解高等数学中的问题(不定积分,定积分,多元微积分,微分方程,级数问题)
用Matlab求解高等数学中的问题 一.求解不定积分 利用matlab符号工具箱中的求积函数int, 可求函数的不定积分,int函数的调用格式如下: int(S) int(S,v) 说明: (1)in ...
- matlab选修结课作业,matlab在高等数学中的应用结课作业
matlab在高等数学中的应用结课作业 MatLab软件在作图中的应用 1.描点作图(MatLab及在电子信息课程中的应用131页) a=0:0.01:3; >> y=a: plot(y) ...
- matlab语法手册下载,MATLAB及其在理工课程中的应用指南 第4版.pdf
作 者 :陈怀琛 出版发行 : 西安:西安电子科技大学出版社 , 2018.02 ISBN号 :978-7-5606-4737-1 页 数 : 262 原书定价 : 39.00 开本 : 16开 主题 ...
- 计算机编程概率论,程序员的数学:线性代数和概率统计
程序员的数学:线性代数和概率统计 语音 编辑 锁定 讨论 上传视频 <程序员的数学:线性代数和概率统计>是2018年9月1日机械工业出版社出版的图书,作者是[美] 欧内斯特 戴维斯 .[1 ...
- 【数字信号处理】线性常系数差分方程 ( 使用 matlab 求解 “ 线性常系数差分方程 “ 示例 | A 向量分析 | B 向量分析 | 输入序列分析 | matlab 代码 )
文章目录 一.使用 matlab 求解 " 线性常系数差分方程 " 示例 1.B 向量元素 : x(n) 参数 2.A 向量元素 : y(n) 参数 3.输入序列 4.matlab ...
- 【数字信号处理】线性常系数差分方程 ( 卷积 与 “ 线性常系数差分方程 “ | 使用 matlab 求解 “ 线性常系数差分方程 “ )
文章目录 一.卷积 与 " 线性常系数差分方程 " 二.使用 matlab 求解 " 线性常系数差分方程 " 一.卷积 与 " 线性常系数差分方程 & ...
- Matlab实现 线性动态电路可视化分析
Matlab实现 线性动态电路可视化分析 版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明. 编程思路 这个编程总体采用面向过程的编程思想,将人在 ...
- matlab在电力电子的应用,MATLAB在电力电子仿真中的应用
MATLAB在电力电子仿真中的应用 MATLAB在电力电子仿真中的应用 田伟 于海生 (青岛大学自动化工程学院,青岛 266071) 摘 要:MATLAB 软件是当今控制系统的设计与仿真中重要的工具软 ...
最新文章
- 强烈不建议你用 a.equals(b) 判断对象相等!
- 用PyCharm Profile分析异步爬虫效率
- 这些产品基本功了解一下?
- 【小白学习C++ 教程】九、C++中字符型、字符串和转义字符
- 京东AI NLP高阶实训学习笔记
- 虚幻引擎发布4.8版,支持VR虚拟设备
- Linux目录的可写意味着,Linux权限分析 - osc_h5427nyq的个人空间 - OSCHINA - 中文开源技术交流社区...
- Windows 平台下的C++代理类(供C # 使用)
- python编程入门经典-总算理解python编程入门经典教程
- 电商设计素材:15000款淘宝促销水印元素模版,手慢的人领不到!...
- YUV420转RGB888
- 使用Windows命令行reg控制注册表键值
- CorelDRAW教程大全集
- 全国计算机考试等级考务管理系统,全国计算机等级考试考务管理系统:https://ncre-bm.neea.edu.cn/...
- “智多星”手机销售网
- POCO C++ 在IOS上的使用
- 阿里云oss文件服务器
- torch.diag() 取矩阵对角线元素,torch.diag_embed() 指定值变成对角矩阵
- 大数据征信助力金融科技腾飞
- 浪潮信息做5G服务器,看来浪潮和运营商在5G布局非常深,中标移动开了个好头...
热门文章
- 动物识别系统 c语言代码_C++程序员避不开虚函数的,就像C语言程序员避不开指针一样...
- C语言编程QQ管理系统,c语言制作学生管理系统srrpqq67.doc
- oracle永久表设置,Oracle 表空间简单管理永久表空间
- c语言linux系统宏,Linux下C语言中的预定义宏
- mysql分日汇总填充0_将记录按小时或按日分组,并在mysql中填充零或空的空白
- java泛型要用t,Java泛型中T和問號(通配符)的區別
- curlopt_ssl_verifypeer后https还是验证不过_2020CNAS能力验证小总结
- VS2015sql本地服务器为空,详解VS2015自带LocalDB数据库用法实例
- Angr安装与使用之使用篇(八)
- java程序编写需注意的问题