第二单元 用python学习微积分(十四)无穷小量和不定积分

本文主要是介绍第二单元 用python学习微积分(十四)无穷小量和不定积分,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

本文内容来自于学习麻省理工学院公开课:单变量微积分-无穷小量和不定积分-网易公开课

开发环境准备:CSDN

目录

一、无穷小量

二、不定积分

1、 ​,这个式子的含义是求谁的导数是sinx

2、 ​

3、 ​

4、 ​

5、 ​

6、 ​

三、例子

1、 ​

2、 ​

3、 ​

4、 ​

5、 ​

6、 ​


一、无穷小量

有函数y = f(x)

y的微分写作: dy = f'(x)dx ( \frac{dy}{dx} = f'(x)莱布尼茨的写法,把导数记作为两个无穷小量的比,按老师的说法这种比值的写法更流行,更接近直觉,更易于理解和使用)

dx 取代了 \Delta x

dy 取代了 \Delta y

添加图片注释,不超过 140 字(可选)

例子:

(64.1)^{\frac{1}{3}}\approx ?

y = x^{\frac{1}{3}}, dy = \frac{1}{3}x^{-\frac{2}{3}}dx

A+x = 64, y = 64^{\frac{1}{3}} = 4

dy = \frac{1}{3} (64)^{-\frac{2}{3}}dx = \frac{1}{3} \frac{1}{16}dx = \frac{1}{48} dx

x = 64, x+dx = 64.1\Rightarrow dx = \frac{1}{10}

(64.1)^{\frac{1}{3}}\approx y + dy\approx4 + \frac{1}{48}dx = 4 + \frac{1}{48*10}(值得注意的是这里的dy实际是 \Delta y ,而并不是无穷小量)

感觉这里老师是用另一种方式阐释了线性近似,因为dy/dx其实取得是x=64时的斜率,而后面dx用的其实是 \Delta x ( 64.1-64)

用线性近似 f(x) \approx f(x_0) + f'(x_0)(x-x_0) 公式算一遍

(64.1)^{\frac{1}{3}}\approx f(64) + f'(64)(x-64) = 4 + \frac{1}{3} (64)^{-\frac{2}{3}}(\frac{1}{10}) = 4 + \frac{1}{480}

二、不定积分

G(x) = \int_{}{}g(x)dx , G(x) 称做g的反导数或g的不定积分

1、 \int sinx\times dx,这个式子的含义是求谁的导数是sinx

from sympy import *
import numpy as np 
x = symbols('x')
y = integrate(1/(x*ln(x)))
y

\int sin(x) dx = -cos(x)

G(x) = -cos(x)

G'(x) = sin(x)

这里也可以写成 \int sin(x) dx = -cos(x) +c(c是常数)

G(x) = -cos(x) + c

所以叫做不定积分

2、 \int x^adx = \frac{1}{a+1}x^{a+1} + c......(a\ne-1)

因为... d(x^{a+1}) = (a+1) x^adx

3、 \int \frac{dx}{x} = ln(\left| x \right|) + c......(x\ne-0)

当x>0时,结果正确

当x<0时, ln(\left| x \right|) + c = ln(-x)+c

( ln(-x)+c)' = \frac{1}{-x} (-x)' = \frac{1}{x}

from sympy import *
import numpy as np import matplotlib.pyplot as plt fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.spines['left'].set_position('zero')
ax.spines['bottom'].set_position('zero')
ax.spines['right'].set_color('none')
ax.spines['top'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.yaxis.set_ticks_position('left')
ax.set_aspect(1 ) def DrawXY(xFrom,xTo,steps,expr,color,label,plt):yarr = []xarr = np.linspace(xFrom ,xTo, steps) for xval in xarr:yval = expr.subs(x,xval)yarr.append(yval)y_nparr = np.array(yarr) plt.plot(xarr, y_nparr, c=color, label=label)    def TangentLine(exprY,x0Val,xVal):diffExpr = diff(exprY)x1,y1,xo,yo = symbols('x1 y1 xo yo')expr = (y1-yo)/(x1-xo) - diffExpr.subs(x,x0Val)eq = expr.subs(xo,x0Val).subs(x1,xVal).subs(yo,exprY.subs(x,x0Val))eq1 = Eq(eq,0)solveY = solve(eq1)return xVal,solveYdef DrawTangentLine(exprY, x0Val,xVal1, xVal2, clr, txt):x1,y1 = TangentLine(exprY, x0Val, xVal1)x2,y2 = TangentLine(exprY, x0Val, xVal2)plt.plot([x1,x2],[y1,y2], color = clr, label=txt)def Newton(expr, x0):ret = x0 - expr.subs(x, x0)/ expr.diff().subs(x,x0)return retx = symbols('x')
y = ln(x)
DrawXY(0.1,4,100,y,'blue','y = ln(x) x>0',plt)
y = ln(0-x)
DrawXY(-0.1,-4,100,y,'blue','y = ln(-x) x<0',plt)plt.legend(loc='lower right')
plt.show()

添加图片注释,不超过 140 字(可选)

4、 \int sec^2(x) dx = tan(x) + c

5、 \int \frac{dx}{(1-x^2)^{\frac{1}{2}}} = sin^{-1}(x) +c

6、 \int \frac{dx}{1+x^2} = tan^{-1}(x) + c

积分的唯一性要看它的常数

F' = G', F(x) = G(x) + c

证明:如果 F' = G'

(F-G)' = F' - G' = 0

所以 F(x) -G(x) = c => F(x) = G(x) + c

三、例子

1、 \int x^3(x^4+2)^5dx = \frac{1}{24}(x^4+2)^6 +c

这里做积分有一个技巧 ,换元法:

u = x^4 + 2, du = 4x^3dx

把这个带入上面的积分式子

\int (x^4 + 2)^5 x^3dx = \int(u^5)\frac{du}{4} = \frac{1}{6}u^6\frac{du}{4} +c= \frac{1}{24}(x^4 +2)^6 +c

2、 \int \frac{xdx}{(1+x^2)^{\frac{1}{2}}} =

u = 1+x^2, du = 2x

\int \frac{du}{2}(u)^{\frac{-1}{2}} = 2u^{\frac{1}{2}} \frac{du}{2} +c= (1+x^2)^{\frac{1}{2}} +c

推荐的方法——提前猜测:

假设你一眼看出这个会最终导向 (1+x^2)^{\frac{1}{2}}就试一试

((1+x^2)^{\frac{1}{2}})' = \frac{1}{2} (1+x^2)^{-\frac{1}{2}}(2x) = x(1+x^2)^{-\frac{1}{2}}

所以 \int \frac{xdx}{(1+x^2)^{\frac{1}{2}}} = (1+x^2)^{\frac{1}{2}} +c

3、 \int e^{6x}dx

猜测法 (e^{6x})' = 6 e^{6x}

\int e^{6x}dx = \frac{1}{6} e^{6x} +c

4、 \int xe^{x^2}dx

猜测法 (e^{-x^2} )' = (e^{-2x})' = -2xe^{-x^2}

\int xe^{x^2}dx = -\frac{-1}{2}e^{-x^2} +c

5、 \int sin(x)cos(x)dx

换元法:u = sin(x) du= cos(x)

\int sin(x)cos(x)dx = \int udu = \frac{1}{2}u^2 +c= \frac{1}{2}sin^2(x)+c

猜测法 (sin^2(x))' = 2sin(x)cos(x)

\int sin(x)cos(x)dx = \frac{1}{2}sin^2(x) +c_1

(cos^2(x))' = -2sin(x)cos(x)

\int sin(x)cos(x)dx = \frac{-1}{2}cos^2(x) +c_2

\frac{1}{2}sin^2(x) +c1 = \frac{-1}{2}cos^2(x) +c2

c1 -c2 = \frac{-1}{2}cos^2(x) - \frac{1}{2}sin^2(x) = \frac{-1}{2}(cos^2(x) + sin^2(x)) = -\frac{1}{2}

6、 \int \frac{dx}{xln(x)}

u = ln(x) , du = \frac{1}{u}

\int \frac{dx}{xln(x)} = \int \frac{1}{ln(x)} \frac{dx}{x} = \int \frac{du}{u} =ln\left| u \right| +c = ln\left| ln(x) \right| +c

from sympy import *
import numpy as np 
x = symbols('x')
y = integrate(1/(x*ln(x)))
y

这篇关于第二单元 用python学习微积分(十四)无穷小量和不定积分的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/280778

相关文章

Python脚本实现自动删除C盘临时文件夹

《Python脚本实现自动删除C盘临时文件夹》在日常使用电脑的过程中,临时文件夹往往会积累大量的无用数据,占用宝贵的磁盘空间,下面我们就来看看Python如何通过脚本实现自动删除C盘临时文件夹吧... 目录一、准备工作二、python脚本编写三、脚本解析四、运行脚本五、案例演示六、注意事项七、总结在日常使用

Python将大量遥感数据的值缩放指定倍数的方法(推荐)

《Python将大量遥感数据的值缩放指定倍数的方法(推荐)》本文介绍基于Python中的gdal模块,批量读取大量多波段遥感影像文件,分别对各波段数据加以数值处理,并将所得处理后数据保存为新的遥感影像... 本文介绍基于python中的gdal模块,批量读取大量多波段遥感影像文件,分别对各波段数据加以数值处

python管理工具之conda安装部署及使用详解

《python管理工具之conda安装部署及使用详解》这篇文章详细介绍了如何安装和使用conda来管理Python环境,它涵盖了从安装部署、镜像源配置到具体的conda使用方法,包括创建、激活、安装包... 目录pytpshheraerUhon管理工具:conda部署+使用一、安装部署1、 下载2、 安装3

Python进阶之Excel基本操作介绍

《Python进阶之Excel基本操作介绍》在现实中,很多工作都需要与数据打交道,Excel作为常用的数据处理工具,一直备受人们的青睐,本文主要为大家介绍了一些Python中Excel的基本操作,希望... 目录概述写入使用 xlwt使用 XlsxWriter读取修改概述在现实中,很多工作都需要与数据打交

使用Python实现在Word中添加或删除超链接

《使用Python实现在Word中添加或删除超链接》在Word文档中,超链接是一种将文本或图像连接到其他文档、网页或同一文档中不同部分的功能,本文将为大家介绍一下Python如何实现在Word中添加或... 在Word文档中,超链接是一种将文本或图像连接到其他文档、网页或同一文档中不同部分的功能。通过添加超

Python MySQL如何通过Binlog获取变更记录恢复数据

《PythonMySQL如何通过Binlog获取变更记录恢复数据》本文介绍了如何使用Python和pymysqlreplication库通过MySQL的二进制日志(Binlog)获取数据库的变更记录... 目录python mysql通过Binlog获取变更记录恢复数据1.安装pymysqlreplicat

利用Python编写一个简单的聊天机器人

《利用Python编写一个简单的聊天机器人》这篇文章主要为大家详细介绍了如何利用Python编写一个简单的聊天机器人,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 使用 python 编写一个简单的聊天机器人可以从最基础的逻辑开始,然后逐步加入更复杂的功能。这里我们将先实现一个简单的

基于Python开发电脑定时关机工具

《基于Python开发电脑定时关机工具》这篇文章主要为大家详细介绍了如何基于Python开发一个电脑定时关机工具,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录1. 简介2. 运行效果3. 相关源码1. 简介这个程序就像一个“忠实的管家”,帮你按时关掉电脑,而且全程不需要你多做

Python实现高效地读写大型文件

《Python实现高效地读写大型文件》Python如何读写的是大型文件,有没有什么方法来提高效率呢,这篇文章就来和大家聊聊如何在Python中高效地读写大型文件,需要的可以了解下... 目录一、逐行读取大型文件二、分块读取大型文件三、使用 mmap 模块进行内存映射文件操作(适用于大文件)四、使用 pand

python实现pdf转word和excel的示例代码

《python实现pdf转word和excel的示例代码》本文主要介绍了python实现pdf转word和excel的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价... 目录一、引言二、python编程1,PDF转Word2,PDF转Excel三、前端页面效果展示总结一