初识代数方程一

最新推荐文章于 2021-03-24 20:16:08 发布

playwrighter

最新推荐文章于 2021-03-24 20:16:08 发布

阅读量2.6k

点赞数 1

分类专栏：人工智能 Python数据分析文章标签：代数方程

原文链接：https://github.com/shaiic

版权

Python数据分析同时被 2 个专栏收录

77 篇文章 8 订阅

订阅专栏

人工智能

59 篇文章 2 订阅

订阅专栏

解方程

简单的一元一次方程

$\begin{array}{l}x+16=-25 \\ x+16-16=-25-16 \\ x=-25-16 \\ x=-41\end{array}$

x = -41                # 验证方程的解
x + 16 == -25

True

带系数的方程

$\begin{array}{l}3 x-2=10 \\ 3 x-2+2=10+2 \\ 3 x=12 \\ x=4\end{array}$

x = 4                  # 代入 x = 4
3 * x - 2 == 10

True

系数是分数的方程

$\begin{array}{l}\frac{x}{3}+1=16 \\ \frac{x}{3}=15 \\ \frac{3}{1} \cdot \frac{x}{3}=15 \cdot 3 \\ x=45\end{array}$

x = 45
x/3 + 1 == 16

True

需要合并同类项的例子

$\begin{array}{l}3 x+2=5 x-1 \\ 3 x+3=5 x \\ 3=2 x \\ \frac{3}{2}=x \\ x=\frac{3}{2} \\ x=1 \frac{1}{2}\end{array}$

x = 1.5
3*x + 2 == 5*x -1

True

一元方程练习

$\begin{array}{l}4(x+2)+3(x-2)=16 \\ 4 x+8+3 x-6=16 \\ 7 x+2=16 \\ 7 x=14 \\ \frac{7 x}{7}=\frac{14}{7} \\ x=2\end{array}$

x = 2
4 * (x + 2) + 3 * (x - 2) == 16

True

线性方程

例子

$\begin{array}{l}2 y+3=3 x-1 \\ 2 y+4=3 x \\ 2 y=3 x-4 \\ y=\frac{3 x-4}{2}\end{array}$

import numpy as np
from tabulate import tabulate

x = np.array(range(-10, 11))  # 从 -10 到 10 的21个数据点
y = (3 * x - 4) / 2           # 对应的函数值

print(tabulate(np.column_stack((x,y)), headers=['x', 'y']))

from IPython.core.interactiveshell import InteractiveShell
InteractiveShell.ast_node_interactivity = "last_expr"
# %matplotlib inline

from matplotlib import pyplot as plt

plt.plot(x, y, color="grey", marker = "o")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-SfXUtcwM-1609746991704)(output_13_0.png)]

截距

plt.plot(x, y, color="grey")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.axhline()                # 画出坐标轴
plt.axvline()
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-8jXbqMiS-1609746991705)(output_15_0.png)]

x轴上的截距

$\begin{array}{l}0=\frac{3 x-4}{2} \\ \frac{3 x-4}{2}=0 \\ 3 x-4=0 \\ 3 x=4 \\ x=\frac{4}{3} \\ x=1 \frac{1}{3}\end{array}$

y轴上的截距

$\begin{array}{l}y=\frac{3 \cdot 0-4}{2} \\ y=\frac{-4}{2} \\ y=-2\end{array}$

plt.plot(x, y, color="grey")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.axhline()
plt.axvline()
plt.annotate('x',(1.333, 0)) # 标出截距点
plt.annotate('y',(0,-2))
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-GXOaeJrQ-1609746991706)(output_17_0.png)]

截距的作用是明显的：两点成一线。连接截距就可以画出函数曲线。

斜率

$\begin{array}{l}\text {slope }=\frac{\Delta y}{\Delta x} \\ m=\frac{y_{2}-y_{1}}{x_{2}-x_{1}} \\ m=\frac{7--2}{6-0} \\ m=\frac{7+2}{6-0} \\ m=1.5\end{array}$

plt.plot(x, y, color="grey")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.axhline()
plt.axvline()
m = 1.5
xInt = 4 / 3
yInt = -2
mx = [0, xInt]
my = [yInt, yInt + m * xInt]
plt.plot(mx, my, color='red', lw=5)  # 用红色标出
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-RiGsDfI7-1609746991707)(output_20_0.png)]

plt.grid()                             # 放大图
plt.axhline()
plt.axvline()
m = 1.5
xInt = 4 / 3
yInt = -2
mx = [0, xInt]
my = [yInt, yInt + m * xInt]
plt.plot(mx, my, color='red', lw=5) 
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-QO5795N9-1609746991708)(output_21_0.png)]

直线的斜率、截距形式

$\begin{array}{l}y=m x+b \\ y=\frac{3 x-4}{2} \\ y=1 \frac{1}{2} x+-2\end{array}$

m = 1.5
yInt = -2

x = np.array(range(-10, 11))
y2 = m * x + yInt                       # 斜率截距形式

plt.plot(x, y2, color="grey")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.axhline()
plt.axvline()

plt.annotate('y', (0, yInt))
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-GhrM3znT-1609746991710)(output_23_0.png)]

线形方程组

解的含义：直线的交点

\begin{equation}x + y = 16 \end{equation}
\begin{equation}10x + 25y = 250 \end{equation}

l1p1 = [16, 0]  # 线1点1
l1p2 = [0, 16]  # 线1点2


l2p1 = [25,0]   # 线2点1
l2p2 = [0,10]   # 线2点2


plt.plot(l1p1,l1p2, color='blue')
plt.plot(l2p1, l2p2, color="orange")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()

plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-zyyfpOPU-1609746991711)(output_25_0.png)]

解线形方程组（消去法）

$\begin{array}{l}x+y=16 \\ 10 x+25 y=250 \\ -10(x+y)=-10(16) \\ 10 x+25 y=250 \\ -10 x+-10 y=-160 \\ 10 x+25 y=250 \\ 15 y=90 \\ y=\frac{90}{15} \\ y=6 \\ x+6=16 \\ x=10\end{array}$

x = 10
y = 6
print ((x + y == 16) & ((10 * x) + (25 * y) == 250))

True

指数、根和对数

指数

$\begin{array}{l}2^{2}=2 \cdot 2=4 \\ 2^{3}=2 \cdot 2 \cdot 2=8\end{array}$

x = 5**3
print(x)

根

$\begin{array}{l}?^{2}=9 \\ \sqrt{9}=3 \\ \sqrt[3]{64}=4\end{array}$

import math

x = math.sqrt(9)            # 平方根
print (x)

cr = round(64 ** (1. / 3))   # 立方根
print(cr)

3.0
4

根是分数形式的指数

$\begin{array}{l}8^{\frac{1}{3}}=\sqrt[3]{8}=2 \\ 9^{\frac{1}{2}}=\sqrt{9}=3\end{array}$

print (9**0.5)
print (math.sqrt(9))

3.0
3.0

对数

对数是指数的逆运算

$\begin{array}{l}4^{?}=16 \\ \log _{4}(16)=2\end{array}$

x = math.log(16, 4)
print(x)

2.0

以10为底的对数

$\log (64)=1.8061$

自然对数

$log _{e}(64)=\ln (64)=4.1589$

print(math.log10(64))

print (math.log(64))

1.806179973983887
4.1588830833596715

幂运算（合并同类项）

$\begin{array}{l}2 y=2 x^{4}\left(\frac{x^{2}+2 x^{2}}{x^{3}}\right) \\ 2 y=2 x^{4}\left(\frac{3 x^{2}}{x^{3}}\right) \\ 2 y=2 x^{4}\left(3 x^{-1}\right) \\ 2 y=6 x^{3} \\ y=3 x^{3}\end{array}$

x = np.array(range(-10, 11))
y3 = 3 * x ** 3

print(tabulate(np.column_stack((x, y3)), headers=['x', 'y']))

plt.plot(x, y3, color="magenta")  # y3是曲线
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.axhline()
plt.axvline()
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-weu8HHRc-1609746991712)(output_41_0.png)]

看一个成指数增长的例子

$y=2^{x}$

y4 = 2.0**x
plt.plot(x, y4, color="magenta")
plt.xlabel('x')
plt.ylabel('y')
plt.grid()
plt.axhline()
plt.axvline()
plt.show()

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-YM6uMwDm-1609746991714)(output_43_0.png)]

复利计算

存款100，年利5%，20年后余额是多少（记复利）

$\begin{array}{l}y 1=100+(100 \cdot 0.05) \\ y 1=100 \cdot 1.05 \\ y 2=100 \cdot 1.05 \cdot 1.05 \\ y 2=100 \cdot 1.05^{2} \\ y 20=100 \cdot 1.05^{20}\end{array}$

year = np.array(range(1, 21))                      # 年份
balance = 100 * (1.05 ** year)                     # 余额
plt.plot(year, balance, color="green")
plt.xlabel('Year')
plt.ylabel('Balance')
plt.show()