python中numpy常用函数(使用Python的SymPy库解决数学运算问题的方法)
python中numpy常用函数
使用Python的SymPy库解决数学运算问题的方法摘要:在学习与科研中,经常会遇到一些数学运算问题,使用计算机完成运算具有速度快和准确性高的优势。python的numpy包具有强大的科学运算功能,且具有其他许多主流科学计算语言不具备的免费、开源、轻量级和灵活的特点。本文使用python语言的numpy库,解决数学运算问题中的线性方程组问题、积分问题、微分问题及矩阵化简问题,结果准确快捷,具有一定的借鉴意义。
1.sympy库简介
sympy一个用于符号型数学计算(symbolic mathematics)的python库。它旨在成为一个功能齐全的计算机代数系统(computer algebra system,cas),同时保持代码简洁、易于理解和扩展。sympy完全是用python写的,并不需要外部的库。
本文选择sympy库的原因在于:
- 免费:该库基于bsd开源许可,免费且开源;
- 基于python:该库完全是用python写就,并以python作为该库操作语言;
- 轻量级:为了使sympy简单易用,该库仅基于mpmath库(一个纯python库,用于浮点运算);
- 灵活性:除了用作交互工具,还可插入其他应用或软件拓展功能中。
具体说来,如果x与y未曾赋值,那么下列语句就会报错
|
#测试语句 print (x + y) |
而符号运算则不同,符号运算多用于公式推导,不需要赋值,此时使用sympy进行符号运算便具有方便快捷的优势,如下述语句便不会报错。
|
#测试语句 x = symbol( 'x' ) y = symbol( 'y' ) print (x + y) |
2 sympy库解决数学运算问题实现
2.1 求解线性方程组
解方程的功能主要是使用sympy中solve函数实现。以式(1)为例,求解过程如下:
(1) 符号表示
sympy库中使用symbol函数定义符号变量,
|
from sympy import * x = symbol( 'x' ) y = symbol( 'y' ) #或者用如下语句 x,y = symbol( 'x y' ) #第二个用空格隔开 |
(2)方程表示
使用代码表示数学符号与手写体的数学运算符号存在一定的差异,下面列举常用的运算符:
- 加号
- 加号 +
- 减号 -
- 除号 /
- 乘号 *
- 指数 **
- 对数 log()
- e的指数次幂 exp()
对于长的表达式,如果不确定运算符的优先级,可以加入小括号提升其优先级。由于需要将表达式都转化成右端等于0,这里把常数3和7移到等式左边。题目中表达式可表示为:
|
2 * x - y - 3 = 0 3 * x + y - 7 = 0 |
(3)使用solve函数解方程
在使用solve函数解方程之前,我们先来看一下solve函数的定义。solve函数的第一个参数是要解的方程,要求右端等于0,第二个参数是未知数。
对于式(1)的求解,代码如下:
|
from sympy import * x = symbol( 'x' ) y = symbol( 'y' ) print (solve([ 2 * x - y - 3 , 3 * x + y - 7 ],[x,y])) |
2.2 求解微积分问题
2.2.1 求解极限问题
在2.1中通过一个简单的二元一次方程组求解熟悉了该库求解数学问题的基本过程,下面本文通过示例,讲解使用sympy库求解微积分的过程。
求解式(2)所示的极限问题,需要用到limit函数求极限。
(1)符号及方程表示
引入sympy库并定义n为符号变量与2.1中一致。
|
from sympy import * n = symbol( 'n' ) s = ((n + 3 ) / (n + 2 )) * * n |
(2)利用limit函数求极限
首先我们介绍limit函数的调用格式:limit(e, z, z0, dir='+'),e为任意表达式,表示求取e(z)在点z0处的极限,dir='+'表示取右极限,die='-'则表示取左极限。则上式的求解代码可表示如下:
|
from sympy import * n = symbol( 'n' ) s = ((n + 3 ) / (n + 2 )) * * n |
print(limit(s,n,oo)) #无穷的表示方法是两个小写的字母o
2.2.2 求解定积分
(1)符号表示
|
from sympy import * t = symbol(t) x = symbol(x) |
(2)方程表示
|
m = integrate(sin(t) / (pi - t),(t, 0 ,x)) n = integrate(m,(x, 0 ,pi)) |
完整代码如下:
|
from sympy import * t = symbol(t) x = symbol(x) m = integrate(sin(t) / (pi - t),(t, 0 ,x)) n = integrate(m,(x, 0 ,pi)) print (n) |
2.2.3 求解微分问题
如求取的通解
(1)符号表示
这里与之前不同的是增加了函数的表示(用f(x)表示y),即例题中的y还有微分表示
|
from sympy import * f = function( 'f' ) x = symbol( 'x' ) |
y'的表示方法由以下代码组成
|
diff(f(x),x) |
这里对diff函数稍作介绍:
上面是求一阶导的方法,求解高阶导的方法如下所示:
|
>>> diff(x * * 3 ,x) 3 * x * * 2 >>> diff(x * * 3 ,x, 1 ) 3 * x * * 2 >>> diff(x * * 3 ,x, 2 ) 6 * x >>> diff(x * * 3 ,x, 3 ) 6 >>> diff(x * * 3 ,x, 4 ) 0 |
即改变第三个参数即可。
下面继续我们的解题过程。
|
#左端 diff(f(x),x) #看一下 print (diff(f(x),x)) #result #d #--(f(x)) #dx #右端 2 * f(x) * x |
用dsolve函数解微分方程
dsolve函数是用来解决微分方程(differential equation)的函数。
函数的一个用法为:
|
dsolve(eq, f(x)) |
第一个参数为微分方程(要先将等式移项为右端为0的形式)。第二个参数为要解的函数(在微分方程中)
举个例子:
|
>>> from sympy import * >>> f = function( 'f' ) >>> x = symbol( 'x' ) >>> pprint( 2 * x - diff(f(x),x)) d 2 * x - - - (f(x)) dx >>> dsolve( 2 * x - diff(f(x),x), f(x)) #result #eq(f(x), c1 + x**2) |
这样,我们可以将我们要解的题目,用以下代码表示。
|
dsolve(diff(f(x),x) - 2 * f(x) * x, f(x)) |
结果为:
eq(f(x), c1*exp(x**2))
#即f(x) = c1*exp(x**2)
对比答案可以发现正确。
完整代码:
|
from sympy import * f = function( 'f' ) x = symbol( 'x' ) print (dsolve(diff(f(x),x) - 2 * f(x) * x,f(x)) |
2.2.4 矩阵化简
平时线性代数问题中我们会遇到化简问题,虽然不算难,但着实麻烦。而且,出一点错就会导致
结果出错。不过好运的是sympy提供了相关的支持。
例题:
符号表示与矩阵表示
|
from sympy import * x1,x2,x3 = symbols( 'x1 x2 x3' ) a11,a12,a13,a22,a23,a33 = symbols( 'a11 a12 a13 a22 a23 a33' ) m = matrix([[x1,x2,x3]]) n = matrix([[a11,a12,a13],[a12,a22,a23],[a13,a23,a33]]) v = matrix([[x1],[x2],[x3]]) |
注意m的表示,需要有两个中括号
化简实现
|
print (m * n * v) |
得到的是:
matrix([[x1*(a11*x1 + a12*x2 + a13*x3) + x2*(a12*x1 + a22*x2 + a23*x3) + x3*(a13*x1 + a23*x2 + a33*x3)]])
使用
|
f = m * n * v print f[ 0 ] |
可以进一步得到化简后的式子
也许你要问我要化简后在计算怎么办?下面我就举个例子。
如果上式中x1,x2,x3均等于1,则可这样代入。
|
from sympy import * x1,x2,x3 = symbols( 'x1 x2 x3' ) a11,a12,a13,a22,a23,a33 = symbols( 'a11 a12 a13 a22 a23 a33' ) m = matrix([[x1, x2, x3]]) n = matrix([[a11, a12, a13], [a12, a22, a23], [a13, a23, a33]]) v = matrix([[x1], [x2], [x3]]) f = m * n * v print f[ 0 ].subs({x1: 1 , x2: 1 , x3: 1 }) |
可得
a11 + 2*a12 + 2*a13 + a22 + 2*a23 + a33
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持开心学习网。
原文链接:https://blog.csdn.net/whcxytj/article/details/75270873
- python 验证码识别模块([机器视觉]使用python自动识别验证码详解)
- python函数参数讲解(Python高级特性与几种函数的讲解)
- python栈和队列(Python 实现数据结构中的栈队列)
- python判断对象是否是某一类型(Python判断对象是否相等及eq函数的讲解)
- python3语法规则(详解Python3注释知识点)
- python3循环使用教程(Python3.4学习笔记之 idle 清屏扩展插件用法分析)
- python基础教学之125 装饰器简介(python3 property装饰器实现原理与用法示例)
- 虚拟环境python详细教程(python创造虚拟环境方法总结)
- python 调钉钉接口(python3实现钉钉消息推送的方法示例)
- python线程安全队列(Python 限制线程的最大数量的方法Semaphore)
- python eval函数原理(浅谈Python中eval的强大与危害)
- python改变字体颜色指令(使用Python自动化破解自定义字体混淆信息的方法实例)
- python简单代码画图(python+pyqt5实现图片批量缩放工具)
- python入门练习网页(详解python项目实战:模拟登陆CSDN)
- python生成密码字典教程(Python利用字典破解WIFI密码的方法)
- python编程将一个三位数反序输出(python实现整数的二进制循环移位)
- 二胎家庭老大爱闹情绪,用这招很有效(二胎家庭老大爱闹情绪)
- 一个30岁男人外遇失败的全过程(一个30岁男人外遇失败的全过程)
- 《无敌破坏王2》 不聊彩蛋,聊聊我从动画里看到的现实那些事儿(无敌破坏王2不聊彩蛋)
- 《寄生虫》 三观不正 人类悲欢从来不相通,感同身受也并非本能(寄生虫三观不正)
- 这部动漫中的女孩子,可比101女孩更加励志(这部动漫中的女孩子)
- 《白狐的人生》热拍 贾征宇偶像包袱难自弃 图(白狐的人生热拍)
热门推荐
- php获取网页内容的几种方法(PHP实现的文件浏览器功能简单示例)
- C#插入排序
- python如何判断文件是否结束(python判断文件是否存在,不存在就创建一个的实例)
- 在linux系统中切换成其他用户(如何在Linux命令行与其他用户通信)
- mysql读写分离代码实现(MySQL Router实现MySQL的读写分离的方法)
- set statistics time on的理解
- css设置背景图片在角落(详解使用CSS固定页面背景图片位置的方法)
- css3动画效果图(css3实现波纹特效、H5实现动态波浪效果)
- pythonselenium查看当前页面元素(python+selenium 定位到元素,无法点击的解决方法)
- django参数判断(Django保护敏感信息的方法示例)
排行榜
- 1
- 2
- 3
- 4
- 5
- 6
- 7
- 8
- 9