Python花式解方程

BakerJ 發布于2019-07-31 11:16 / 1964人閱讀

摘要：用來解方程的話有點復雜，需要用到矩陣的思維我矩陣沒學好再加上不能解非線性方程組，所以我也不會這玩意兒遜色于和，但解方程也是非常不錯的既能解線性方程組，又能解非線性方程組，堪稱解方程界的神器，但是表達式不支持位運算，比如與或非，取余以及

numpy

numpy 用來解方程的話有點復雜，需要用到矩陣的思維！我矩陣沒學好再加上 numpy 不能解非線性方程組，所以...我也不會這玩意兒！

sympy

遜色于 sage 和 z3，但解方程也是非常不錯的！

from sympy import *
x = symbols("x")
y = symbols("y")
res = solve([x+y-3,x-y-1],[x,y])[0]
print(res)

sage

sage 既能解線性方程組，又能解非線性方程組，堪稱解方程界的神器，但是表達式不支持位運算，比如：與或非，取余以及異或。出現位運算的方程就只能用 z3 創建約束求解！sage 的優點也很明顯：表達式簡單易寫，運算速度快！
在線sage求解

var("x y")
solve([x**3+y**2+666==142335262,x**2-y==269086,x+y==1834],[x,y])

z3 也叫約束求解器，用來解任何方程都沒有問題！但是 windows 不太好裝，所以我基本上是在linux上跑，python2 和 python3 都支持！使用的思路非常簡單：

先創建你所需類型的符號變量

再初始化一個約束器，

添加約束

最后判斷約束是否有解以及求解變量

下面列舉常用的函數，順便給個 z3-solver文檔

# 符號變量類型
Int("x")
Real("x")
Bool("x")
BitVec("x",N) # N bit的符號變量，用于位操作
BitVecVal(num,N) # N bit的數據 num
# 初始化約束器
solver = Solver()
# 添加約束
solver.add(x+y==10,x-y==0)
# 求解約束
solver.check()
ans = solver.mode()

# 初始化多個符號變量
x = [Int("x%d" % i) for i in range(n)]
# 取結果中某個變量的值
value = ans[x].as_long()

END

云服務器 GPU云服務器解方程用雅克比解方程c語言 python求解方程 python微分方程

文章版權歸作者所有，未經允許請勿轉載,若此文章存在違規行為，您可以聯系管理員刪除。

轉載請注明本文地址：http://specialneedsforspecialkids.com/yun/44894.html

貌離神合的RNN與ODE：花式RNN簡介

摘要：事實上，我記得確實有一些教程是直接通過微分方程來定義函數的。歐拉的解法來源很簡單，就是用來近似導數項。這樣一來，我們就知道的歐拉解法實際上就是的一個特例罷了。作者丨蘇劍林單位丨廣州火焰信息科技有限公司研究方向丨NLP，神經網絡個人主頁丨kexue.fm本來筆者已經決心不玩 RNN 了，但是在上個星期思考時忽然意識到 RNN 實際上對應了 ODE（常微分方程）的數值解法，這為我一直以來想做的...

darcrand 2019-04-25 18:28 評論0 收藏0
用Python學數學之Sympy代數符號運算

摘要：的符號運算如果之前是學數學相關專業了解計算機代數系統，就會對數學符號的運算比較熟悉，而如果之前是程序員，可能會有點不太明白，下面我們就來了解一下。在我們初、高中和大學近10年的學習時間里，數學一直占據著非常大的分量，但是回憶過去可以發現，我們把大量的時間都花在反復解題、不斷運算上，計算方法、運算技巧、筆算能力以及數學公式的記憶仿佛成了我們學習數學的全部。這些記憶和技巧沒幾年就忘掉了，...

Jackwoo 2019-07-31 10:13 評論0 收藏0
三對角線性方程組(tridiagonal systems of equations)的求解

摘要：三對角線性方程組三對角線性方程組對于熟悉數值分析的同學來說，并不陌生，它經常出現在微分方程的數值求解和三次樣條函數的插值問題中。三對角線性方程組(tridiagonal systems of equations) ??三對角線性方程組,對于熟悉數值分析的同學來說，并不陌生，它經常出現在微分方程的數值求解和三次樣條函數的插值問題中。三對角線性方程組可描述為以下方程組：$$a_{i}x_...

yimo 2019-07-30 16:28 評論0 收藏0