Matlab--微积分问题的计算机求解

news2024/11/19 15:22:51

目录

1.单变量函数的极限问题

1.1.公式例子

1.2.对应例题 1

2.多变量函数的极限问题

3.函数导数的解析解

4.多元函数的偏导数

5.Jacobian函数

6.Hessian矩阵

 7.隐函数的偏导

 8.不定积分问题的求解

9.定积分的求解问题

10. 多重积分的问题求解


1.单变量函数的极限问题

1.1.公式例子

%%
%3.1.1.单变量函数的极限
%直接趋近
syms x;
y=2*x+2;
L=limit(y,x,2)
%左右趋近
L1=limit(y,x,2,'left')
L2=limit(y,x,2,'right')

1.2.对应例题 1

%%
%例题1
syms x;
y=(exp(x^3)-1)/(1-cos(sqrt(x-sin(x))))
limit(y,x,0,'right')

1.3. 对应例题2

%%
%例题2
syms t;
f=tan(t);
L1=limit(f,t,pi/2,'left')
L2=limit(f,t,pi/2,'right')
%L1 =Inf L2 =-Inf 

2.多变量函数的极限问题

%多变量函数的极限
%需要嵌套使用limit()函数
%L=limit(limit(f,y,y0),x,x0)
syms x y a;
f=exp(-1/(y^2+x^2))*sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2)
L=limit(limit(f,x,1/sqrt(y)),y,inf)

3.函数导数的解析解

%%
%函数导数的解析解
%函数的导数和高阶导数
%y=diff(fun,x) && y=diff(fun,x,n)
syms x;
f=sin(x)/(x^2+4*x+3);
f1=diff(f)
x1=0:0.01:5;
y=subs(f,x,x1);
y1=subs(f1,x,x1);
plot(x1,y,x1,y1,':')
saveas(gcf,[ 'C:\Users\Zeng Zhong Yan\Desktop\MATLAB\matlab_code','matalb2', '.png']);

 

4.多元函数的偏导数

%%
%多元函数的偏导数
%f=diff(diff(f,x,m),y,n)
%f=diff(diff(f,y,n),x,m)
%例题3.8
syms x y;
z=(x^2-2*x)*exp(-x^2-y^2-x*y);
zx=simplify(diff(z,x)),
zy=diff(z,y)

[x0,y0]=meshgrid(-3:.2:3,-2:0.2:2);
z=(x0.^2-2*x0).*exp(-x0.^2-y0.^2-x0.*y0);
surf(x0,y0,z),
axis([-3 3 -2 2 -0.7 1.5])
saveas(gcf,[ 'C:\Users\Zeng Zhong Yan\Desktop\MATLAB\','偏导数', '.png']);

 

对应例题:

%%
%例题
syms x y z;
f=sin(x^2*y)*exp(-x^2*y-z^2);
f4=diff(diff(diff(f,x,2),y,1),z,1);
f4=simplify(f4)

5.Jacobian函数

%%
%Jacobian()函数
syms r theta phi
x=r*sin(theta)*cos(phi);
y=r*sin(theta)*sin(phi);
z=r*cos(theta);
J=jacobian([x;y;z],[r theta phi])

6.Hessian矩阵

%%
%Hessian偏导数矩阵
%Hessian矩阵就是两次的雅可比矩阵
H=jacobian(jacobian(f,x),x)

对应例题:

%%
%Hessian矩阵例题
syms x y;
f=(x^2-2*x)*exp(-x^2-y^2-x*y);
H=jacobian(jacobian(f,[x,y]),[x,y])

 7.隐函数的偏导

%%
%隐函数的偏导数
%已知f(x1,x2,x3....x4)=0求解偏导数
F=-diff(f,xj)/diff(f,xi)

 对应例题:

%%
%隐函数求导例题
syms x y;
f=(x^2-2*x)*exp(-x^2-y^2-x*y);
d=-diff(f,x)/diff(f,y);
d=simplify(d)

 8.不定积分问题的求解

%%
%不定积分的求解
F=int(fun,x);
%多重嵌套
%F=int(.......int(fun,x))

 对应例题:

%%
%不定积分例题
syms x;
y=sin(x)/(x^2+4*x+3);
y1=diff(y);
y0=int(y1);
y1,y0
y4=diff(y,4);
y0=int(int(int(int(y4))))

9.定积分的求解问题

%%
%定积分与无穷运算
I=int(fun,x,a,b);
I=int(fun,x,a,inf);

对应例题1:

%%
syms x;
y=exp(-x^2/2)
I1=int(y,x,0,1.5)
vpa(I1,70)
I1=int(y,x,0,inf)

对应例题2:

%%
%定积分例题
f=(-2*x^2+1)/(2*x^2-3*x+1)^2
I=simplify(int(f,x,cos(t),exp(-2*t)))

10. 多重积分的问题求解

%%
%多重积分的matlab问题求解
%积分公式差不多,只要注意积分的顺序
syms x y z
k=4*x*z*exp(-x^2*y-z^2);
int(int(int(k,x,0,2),y,0,pi),z,0,pi)
%ans =-(exp(-pi^2) - 1)*(eulergamma + log(4*pi) 

 

 

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1017799.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【Vue】快速入门案例与工作流程的讲解

🎉🎉欢迎来到我的CSDN主页!🎉🎉 🏅我是Java方文山,一个在CSDN分享笔记的博主。📚📚 🌟在这里,我要推荐给大家我的专栏《Vue快速入门》。&#x1f…

运维Shell牛刀小试(十一):for循环读取多个命令行参数|read重定向读取文件内容

运维Shell脚本小试牛刀(一) 运维Shell脚本小试牛刀(二) 运维Shell脚本小试牛刀(三)::$(cd $(dirname $0); pwd)命令详解 运维Shell脚本小试牛刀(四): 多层嵌套if...elif...elif....else fi_蜗牛杨哥的博客-CSDN博客 Cenos7安装小火车程序动画 运维Shell脚本小试…

【系统架构】系统架构设计基础知识

导读:本文整理关于系统架构设计基础知识来构建系统架构知识体系。完整和扎实的系统架构知识体系是作为架构设计的理论支撑,基于大量项目实践经验基础上,不断加深理论体系的理解,从而能够创造新解决系统相关问题。 目录 1、软件架…

高速信号处理板资料保存:383-基于kintex UltraScale XCKU060的双路QSFP+光纤PCIe 卡设计原理图

基于kintex UltraScale XCKU060的双路QSFP光纤PCIe 卡 一、板卡概述 本板卡系我司自主研发,基于Xilinx UltraScale Kintex系列FPGA XCKU060-FFVA1156-2-I架构,支持PCIE Gen3 x8模式的高速信号处理板卡,搭配两路40G QSFP接口&#xf…

PyQt5入门2——添加一个画布并且显示特定的图片

PyQt5入门2——添加一个画布并且显示特定的图片 学习前言使用到的PyQt5类实例使用1、窗口构建a、构建基础类b、读取已有的图片并且显示 2、主程序运行 全部代码 学习前言 搞搞可视化界面哈,虽然不一定有用,但是搞一下。 使用到的PyQt5类 创建画布需要…

MapReduce YARN 的部署

1、部署说明 Hadoop HDFS分布式文件系统,我们会启动: NameNode进程作为管理节点DataNode进程作为工作节点SecondaryNamenode作为辅助 同理,Hadoop YARN分布式资源调度,会启动:ResourceManager进程作为管理节点NodeM…

蓝桥杯2023年第十四届省赛真题-更小的数--题解

目录 蓝桥杯2023年第十四届省赛真题-更小的数 题目描述 输入格式 输出格式 样例输入 样例输出 提示 【思路解析】 【代码实现】 蓝桥杯2023年第十四届省赛真题-更小的数 时间限制: 3s 内存限制: 320MB 提交: 895 解决: 303 题目描述 小蓝有一个长度均为 n 且仅由数字…

【Qt-17】Qt调用matlab生成的dll库

matlab生成dll库 1、matlab示例代码 function BDCube(x,y)[x,y,z] cylinder(x,y);t1 hgtransform;s1 surf(3*x,3*y,4*z,Parent,t1);grid onview(3)shading interp end 2、matlab环境配置 首先检查自己的mcc编译器是否可用,输出以下命令: &#x…

无头单向非循环链表(详解)

无头单向非循环链表 链表的种类的简单介绍链表的定义typedef的运用 节点的创建遍历头插尾插头删尾删查找节点位置任意位置插入任意位置删除 链表的种类的简单介绍 链表的种类有很多 1.分有头链表和无头链表 2.分循环链表和非循环链表 3.分单向链表和双向链表 这里的话&#xff…

RecyclerView源码解析(四):RecyclerView对ViewHolder的回收

RecyclerView源码解析(四):RecyclerView对ViewHolder的回收 导言 前面几篇文章我们已经介绍了RecyclerView绘图的三大流程和其四层缓存机制,不过对于我来说还有一个疑问,那就是RecyclerView是如何将分离的ViewHolder进…

MATLAB:电机控制(Motor Control)

目录 前言1、Overview1.1、Specific objectives1.2、Resources1.3 、Industrial context1.4 、Architecture of the system1.5 、Approach/Steps 2、MODELLING THE MOTOR AND THE LOAD2.1Simulate a DC motor by its physical model2.2、Modelling the motor / load by its mat…

堆内存与栈内存

文章目录 1. 栈内存2. 堆内存3. 区别和联系参考资料 1. 栈内存 栈内存是为线程留出的临时空间 每个线程都有一个固定大小的栈空间,而且栈空间存储的数据只能由当前线程访问,所以它是线程安全的。栈空间的分配和回收是由系统来做的,我们不需…

《C++ primer plus》精炼(OOP部分)——对象和类(5)

“学习是照亮心灵的火炬,它永不熄灭,永不止息。” 文章目录 类的自动和强制类型转换原始类型转换为自定义类型将自定义类型转换为原始类型 类的自动和强制类型转换 原始类型转换为自定义类型 可以用一个参数的构造函数来实现,例如&#xff…

【Linux操作系统】信号的产生捕获

🔥🔥 欢迎来到小林的博客!!       🛰️博客主页:✈️林 子       🛰️博客专栏:✈️ Linux       🛰️社区 :✈️ 进步学堂       &#x1f6f0…

destoon关于archiver归档的性能优化

今天在处理一个项目时候发现archiver单个模块归档超过百万数据,打开速度就特慢,所以打开archiver下index.php文件进行分析,发现有句sql作怪,查询需要三四分钟,所以要修改这个。 $result $db->query("SELECT …

login (mobile email usercode)

手机、邮箱、用户编码登录 package login;/*** 常量类(改造拉到你自己项目里面)** author ZengWenFeng* date 2023.09.17* email 117791303qq.com* mobile 13805029595*/ public class Constant {// 登录界面的登录方式:手机、邮箱、用户名public static final Str…

YOLOv8学习

1 YOLOv8训练自己的数据集 1.1 部署和运行 1.1.1 第一步,下载代码 源码地址 1.1.2 第二步,创建环境 默认已经安装好conda conda create -n yolopy38 python3.8conda activate yolopy38安装所需要的包,先cd到代码目录下 pip install -r…

深眸科技自研工业AI视觉检测设备,检测精度99.9%加速智造进程

随着机器视觉技术的持续升级,国内制造行业不断发展,工艺水平持续优化,产品的数量和种类与日俱增。当前社会将产品质量标准提高,在满足正常的性能使用外,还需要具有良好的表面外观质量。但在工业制造过程中,…

【C语言】【数据存储】用%d打印char类型数据,猜结果是啥

题目代码如下&#xff1a; #include <stdio.h> int main() {char a -1;signed char b-1;unsigned char c-1;printf("a%d,b%d,c%d",a,b,c);return 0; }解题关键&#xff1a; 1.二进制存储&#xff1a;原码&#xff0c;反码&#xff0c;补码 互换 2.截断 3.整型…

密码学入门——Feistel网络

文章目录 参考书一、引入二、Feistel网络计算流程三、Feistel网络的解密四、Feistel网络的性质 参考书 图解密码技术&#xff0c;第三版 一、引入 DES的基本结构是由Horst Feistel设计的&#xff0c;因此也被成为Feistel网络&#xff08;Feistel结构&#xff0c;Feistel密码…