首页 > 其他分享 >微积分相关

微积分相关

时间:2024-01-20 13:44:41浏览次数:26  
标签:dots cos frac 微积分 theta 相关 underline lambda

拉格朗日乘数法

对于多元函数 \(f(x_1,x_2,\dots,x_n)\),有若 \(m\) 个约束条件形如:\(g_i(x_1,x_2,\dots,x_n)=0\)。

我们要求 \(f\) 在约束条件下的极值。

首先,对与一元情况,我们只要找到所有导数为 \(0\) 的点即可。

对于多元和约束,我们构造函数 \(L(x_1, x_2, \dots , x_n, \lambda)\),\(\lambda\) 是实数:

\[L(x_1, x_2, \dots , x_n, \lambda)=f(x_1,x_2,\dots,x_n)+\lambda(\sum_{i=1}^mg(x_1, x_2, \dots , x_n)) \]

对于这个函数,我们求出对 \(x_1,x_2,\dots,x_n\) 求出他们的偏微分,要求每个偏微分都等于 \(0\) 并且 \(\lambda\) 的偏微分也是 \(0\)(这样就相当于 \(g\) 为 \(0\))。

即满足方程组:

\[\begin{cases} \frac{\partial L}{\partial x_1}=0\\ \frac{\partial L}{\partial x_2}=0\\ \,\,\,\,\,\vdots\\ \frac{\partial L}{\partial x_n}=0\\ \frac{\partial L}{\partial \lambda}\,\,=0\\ \end{cases} \]

求出这个方程组的所有解就是可能的极值。

例1: 过河问题:有一条河,分成三段,每段水流速度分别为 \(v_1,v_2,v_3\),宽度分别为 \(w_1, w_2, w_3\),现在要从东岸到西岸,以恒定的速度 \(v\) 过河,要求上岸点和出发地点在同一垂直于河流的直线上,求最少时间。

思路:

显然这个问题中的变量是各段出发的角度,不妨设偏移的角度分别是 \(\theta_1,\theta_2,\theta_3\)。

通过对速度分解,表示过河时间的函数为:

\[f(\theta_1,\theta_2,\theta_3)= \frac{w_1}{v\cos\theta_1}+ \frac{w_2}{v\cos\theta_2}+ \frac{w_3}{v\cos\theta_3} \]

约束为:

\[g(\theta_1,\theta_2,\theta_3)= \frac{w_1}{v\cos\theta_1}(v_1-v\sin\theta_1)+ \frac{w_2}{v\cos\theta_2}(v_2-v\sin\theta_2)+ \frac{w_3}{v\cos\theta_3}(v_3-v\sin\theta_3) \]

于是我们可以用拉格朗日乘数法,引入实数 \(\lambda\),构造函数 \(L(\theta_1,\theta_2,\theta_3,\lambda)\)。

我们考虑对 \(\theta_1,\theta_2,\theta_3,\lambda\) 求偏微分。

首先,求导中有结论:

\[(\frac{u}{v})'=\frac{u'v-uv'}{v^2} \]

先对 \(\theta_1\) 求导,显然只有 \(\frac{w_1}{v\cos\theta_1}\) 和 \(\frac{w_1}{v\cos\theta_1}(v_1-v\sin\theta_1)\lambda\) 有影响。

第一个,\(\cos\) 的导数是 \(-\sin\),再根据结论:

\[(\frac{w_1}{v\cos\theta_1})'=\frac{w_1\sin \theta_1}{v\cos^2 \theta_1} \]

第二个,同理:

\[\begin{aligned} [\frac{w_1}{v\cos\theta_1}(v_1-v\sin\theta_1)\lambda]' &= \frac{w_1\sin \theta_1}{v\cos^2 \theta_1}\lambda v_1 -\frac{\lambda w_1}{\cos^2\theta_1}\\ \end{aligned} \]

所以它的偏微分就是:

\[\frac{w_1\sin \theta_1}{v\cos^2 \theta_1}(1+\lambda v_1) -\frac{\lambda w_1}{\cos^2\theta_1} \]

要求其等于 \(0\),化简就可以得到:

\[\sin \theta_1(1+\lambda v_1)=\lambda v \]

其他同理,我们就可以得到一个方程组,解开就行了。

怎么解还不会,咕咕咕咕。

例2: 平面上有若干个点,距离原点距离分别为 \(r_1,r_2,\dots,r_n\),要求形成凸包面积最大。

思路:

同理,也是以角度为变量,以和为 \(2\pi\) 为约束即可。

最终可以得出方程组,与 \(\lambda\) 有单调性,二分 \(\lambda\) 即可求出答案。

例3: 三位平面上 \(n\) 个点 \((a_1,b_1,c_1),\dots,(a_n,b_n,c_n)\),要求求一条直线 \(L\) 使得投影在这条直线上后所有值的方差最大。

思路:

学了线代再回来做。

离散微积分

当 \(f(x)\) 是定义在 \(\mathbb{Z}\) 上的函数时,我们定义差分(也就是导数):

\[g(x)=\Delta f(x)=f(x+1)-f(x) \]

离散积分就是:

\[\sum_{a}^{b}g(x)\delta x=\sum_{k=a}^{b-1}g(k)=f(b)-f(a) \]

与微积分基本定理类似。

接着,我们还要知道两个东西。

上阶乘幂:\(x^{\overline{m}}=x(x+1)(x+2)\dots(x+m-1)\)。

下阶乘幂:\(x^{\underline{m}}=x(x-1)(x-2)\dots(x-m+1)\)。

我们对下阶乘幂求离散导数就会得到:

\[\Delta x^{\underline{m}}=mx^{\underline{m-1}} \]

与幂函数求导类似。

运用这个我们还可以得到:

\[\sum_{a}^b x^{\underline{m}}\delta x=\frac{1}{m+1}x^{\underline{m+1}}|^{b}_{a} \]

例1: 求出 \(\sum_{0 \le k < n}k^2\) 的封闭形式。

思路:

可以用离散积分来做,有 \(x^2=x^{\underline{2}}+x^{\underline{1}}\),于是我们可以根据积分得出封闭形式。

这种方法可以解决任意次幂。

例2: 证明 \((a+b)^{\underline{n}}=\sum_{r=0}^n\binom{n}{r}a^{\underline{n-r}}b^{\underline{r}}\)

思路:

\(a+b\) 个球中选 \(n\) 个排成一行,double counting.

标签:dots,cos,frac,微积分,theta,相关,underline,lambda
From: https://www.cnblogs.com/rlc202204/p/17976381

相关文章

  • 数学和CNN里面的卷积和互相关
    卷积和互相关nndl上CNN这章的互相关讲的比较晦涩,简单辨析一下书上的互相关A.1数学意义上的卷积就是将卷积核进行翻转之后再进行我们熟悉CNN上的卷积运算同时互相关就是不将卷积核翻转直接CNN卷积运算说到这里就明白了,做如下总结\(数学上的互相关=CNN卷积\)$数学上的卷积......
  • 机器视觉 - yolo 相关工具
    模型可视化netron网站netron源码标签格式转换文章:https://zhuanlan.zhihu.com/p/461488682代码:https://github.com/KKKSQJ/DeepLearning/blob/master/others/label_convert/README.mdyolo标注文件可视化源码:https://github.com/KKKSQJ/DeepLearning/blob/master/oth......
  • 栈和相关算法
    栈栈是一种抽象数据结构(ADT),其主要特性是后进先出LIFO(LastinFirstout)实现方式可以用数组、链表实现,本质就是对一个列表进行后进先出的操作操作栈的操作主要有push入栈、pop出栈、isEmpty判空、getTop获取栈顶元素数组实现首先进行最基本的数据结构和操作定义://栈......
  • 教学相关工具
    从超星导出压缩包提取某些类型文件,如.docx、.doc等文件源码删除无用目录,如android源码中的build、.idea、.gradle等目录源码对文件及目录统一命名源码对实验报告统一打分源码教务系统期末成绩导入源码......
  • 多项式相关
    板子:1constintN=4e5+5;23structNTT{4constintmod=998244353,G=3,invG=332748118;5intrev[N],res[N],Inv[N],Res[N],eres[N],Eres[N];67voidpreinv(intn){8Inv[0]=Inv[1]=1;9......
  • Java 秘钥对相关操作
    生成JKS(JavaKeyStore)文件keytool-genkeypair-keystoremercury.jks-keyalgRSA-validity180-aliasmercury参数说明keytool:这是JavaKeytool工具,用于管理密钥和证书。-genkeypair:指示Keytool生成一个密钥对(公钥和私钥)。-aliasmercury:设置密钥对的别......
  • 常见的自动化测试相关框架
    Appium:一个开源、跨平台的自动化测试工具,用于测试原生和轻量移动应用,支持iOS、Android和FirefoxOS平台。Carina:一款Java自动测试框架,实现很完善、功能齐全,但文档较少,对于测试人员学习难度有要求。Galen:一个开放源码的测试网页布局和响应设计的开源工具。Gauge:一种相对较新的测试自......
  • Kubernetes之云原生相关
    什么是云原生可以简单看做就是K8S。将项目全部都通过K8S部署。实际上,云原生是一条最佳路径或者最佳实践。更详细的说,云原生为用户指定了一条低心智负担的、敏捷的、能够以可扩展、可复制的方式最大化地利用云的能力、发挥云的价值的最佳路径。因此,云原生其实是一套指导进行软件......
  • 最优订单执行算法相关Paper介绍
    更多精彩内容,欢迎关注公众号:数量技术宅,也可添加技术宅个人微信号:sljsz01,与我交流。随着量化交易、高频交易的竞争日益激烈,事实证明,交易执行显着影响量化策略的投资绩效。因此,许多从业者开始将交易执行视为一项独立的任务。本期文章分享,我们搜寻了一些用于最小化交易成本的最佳......
  • 【python网络编程相关】 ----操作系统相关了解
    title:【python网络编程相关】----操作系统相关了解date:2024-01-1615:54:06updated:2024-01-1616:20:00description:【python网络编程相关】----操作系统相关了解cover: https://www.cnblogs.com/YZL2333/p/10444200.htmlhttps://home.cnblogs.com/u/......