复数的运算

Yasen.M

于 2024-10-31 10:33:53 发布

阅读量2.3k

点赞数 21

分类专栏：算法 python 文章标签：算法 python

本文链接：https://blog.csdn.net/qq_45449625/article/details/143381057

版权

python 同时被 2 个专栏收录

9 篇文章

订阅专栏

算法

2 篇文章

订阅专栏

1.介绍

复数是数学中一种扩展了实数概念的数字类型，通常表示为 ( z = a + bi )，其中：

( a ) 是实部（real part），
( b ) 是虚部（imaginary part），
( i ) 是虚数单位，满足 ( i^2 = -1 )。

复数的一些基本特性：

复平面：复数可以在二维平面上表示，称为复平面或阿根伯平面。实部沿水平轴（x轴），虚部沿垂直轴（y轴）。
模：复数的模（或绝对值）是其与原点的距离，计算公式为：

$|z| = \sqrt{a^2 + b^2}$

共轭：复数的共轭是将虚部的符号反转，表示为 $\overline{z} = a - bi$ 。
极坐标表示：复数可以用极坐标形式表示，形式为 $( r(\cos \theta + i \sin \theta) )$ 或用欧拉公式表示为 $( re^{i\theta} )$ ，其中 ( r ) 是模， $( \theta )$ 是幅角。

2.复数的应用

复数在多个领域中有广泛的应用，包括但不限于以下几个方面：

1.电工程

复数用于交流电路分析，特别是在描述电压和电流的相位差、阻抗等方面。通过使用复，可以将电路中的正弦波形简化为代数运算。

例子：交流电路中的阻抗计算

假设我们有一个交流电路，包含一个电阻 $( R = 4 , \Omega )$ 和一个电感 $( L = 2 , \text{H} )$ 。交流电源的频率为 $( f = 50 , \text{Hz} )$ 。

首先，我们需要计算电感的阻抗 $( Z_L )$ ：

$[ Z_L = j \omega L ]$
其中 ( $\omega = 2\pi f$ )，所以：
$[ \omega = 2 \pi \times 50 = 100\pi \approx 314.16 , \text{rad/s} ]$
因此，电感的阻抗为：
$[ Z_L = j(314.16)(2) \approx j628.32 , \Omega ]$

电路的总阻抗 ( Z ) 为：
$[ Z = R + Z_L = 4 + j628.32 , \Omega ]$

计算模和相位：

模：
$[ |Z| = \sqrt{R^2 + (Z_L)^2} = \sqrt{4^2 + 628.32^2} \approx \sqrt{16 + 394,889.41} \approx 628.36 , \Omega ]$
相位角：
$[ \phi = \tan^{-1}\left(\frac{\text{Im}(Z)}{\text{Re}(Z)}\right) = \tan^{-1}\left(\frac{628.32}{4}\right) \approx \tan^{-1}(157.08) \approx 89.63^\circ ]$

2.控制理论

在控制系统的分析与设计中，复数用于描述系统的动态特性，尤其是在根轨迹和频域分析中。系统的稳定性和响应特性可以通过复数平面上的极点和零点来分析。

3.信号处理：

复数广泛应用于数字信号处理（DSP），例如在傅里叶变换中。通过复数，可以更方便地处理和分析信号的频率成分。

例子：傅里叶变换

考虑一个简单的信号 $( x(t) = \cos(2 \pi f_0 t) )$ ，其中 $( f_0 = 1 , \text{Hz} )$ 。我们想要计算这个信号的傅里叶变换。

傅里叶变换定义为： $[ X(f) = \int_{-\infty}^{\infty} x(t) e^{-j 2 \pi f t} dt ]$

对于 $( x(t) = \cos(2 \pi f_0 t) )$ ，我们可以利用欧拉公式：
$[ \cos(2 \pi f_0 t) = \frac{e^{j 2 \pi f_0 t} + e^{-j 2 \pi f_0 t}}{2} ]$

所以傅里叶变换可以写成：

$[ X(f) = \int_{-\infty}^{\infty} \frac{e^{j 2 \pi f_0 t} + e^{-j 2 \pi f_0 t}}{2} e^{-j 2 \pi f t} dt ]$

这个积分可以分解为两个部分：

$[ X(f) = \frac{1}{2} \left( \int_{-\infty}^{\infty} e^{j 2 \pi (f_0 - f) t} dt + \int_{-\infty}^{\infty} e^{-j 2 \pi (f_0 + f) t} dt \right) ]$

根据傅里叶变换的性质，结果为：

$[ X(f) = \frac{1}{2} \left( \delta(f - f_0) + \delta(f + f_0) \right) ]$
这表明，信号 ( x(t) ) 在频域上具有两个冲击，分别位于 ( $f = 1 , \text{Hz}$ ) 和 ( $f = -1 , \text{Hz}$ )。

4.量子物理：

在量子力学中，波函数通常用复数表示，复数的模平方表示粒子存在的概率密度。这是量子力学的一个核心概念。

5.流体动力学：

复数在流体动力学中用于分析流场和潜流，尤其是在潜流的复数势流理论中，复数帮助简化问题的复杂性。

6.图像处理：

在图像处理和计算机视觉中，复数可用于图像的频域变换，例如使用傅里叶变换进行图像滤波和增强。

7.数学理论：

复数是数学分析、复变函数、代数和几何等多个领域的重要工具。复数理论在许多数学定理和公式的证明中发挥了关键作用，例如柯西-黎曼方程和留数定理。

例子：复变函数的积分

考虑复函数 ( $f(z) = z^2$ ) 在 ( z = 1 ) 和 ( z = i ) 之间的积分。

通过参数化路径，我们可以计算：

$[ \int_C f(z) , dz ]$
其中 ( C ) 是从 ( 1 ) 到 ( i ) 的直线。我们可以参数化为 $( z(t) = 1 + (i - 1)t )$ 其中 ( t ) 从 0 到 1。

计算 ( dz )：

$[ dz = (i - 1) dt ]$

因此，积分变为：

$[ \int_0^1 f(z(t)) \cdot dz = \int_0^1 (1 + (i - 1)t)^2 \cdot (i - 1) , dt ]$

展开 $( (1 + (i - 1)t)^2 )$ ：

$[ = (1 + (i-1)t)(1 + (i-1)t) = 1 + 2(i-1)t + (i-1)^2 t^2 = 1 + 2(i-1)t - t^2 \quad (\text{since } (i-1)^2 = -2i) ]$

因此，计算积分：

$[ \int_0^1 (1 + 2(i-1)t - t^2)(i-1) , dt ]$

这个过程比较复杂，但它展示了复变函数的积分计算过程。

3.复数的功能

复数在数学和工程领域具有多种功能和应用：

加法和减法：复数可以直接相加或相减，方法是分别加或减其实部和虚部。例如：
- ( $(a + bi) + (c + di) = (a + c) + (b + d)i$ )
乘法：复数的乘法遵循分配律和结合律。乘法的结果为：
- ( $(a + bi)(c + di) = (ac - bd) + (ad + bc)i$ )
除法：复数的除法通常通过乘以共轭（conjugate）来实现，以消去分母中的虚部。即：
- ( $\frac{a + bi}{c + di} = \frac{(a + bi)(c - di)}{c^2 + d^2}$ )
模：复数的模（或绝对值）是其与原点的距离，计算公式为：
- ( $|z| = \sqrt{a^2 + b^2}$ )
共轭：复数的共轭是将虚部的符号反转，表示为 ( $\overline{z} = a - bi$ )。共轭在复数运算中经常使用，尤其是在除法中。
极坐标形式：复数可以用极坐标表示，形式为 ( $r(\cos \theta + i \sin \theta)$ ) 或用欧拉公式表示为 ( $re^{i\theta}$ )，其中 ( r ) 是模，( $\theta$ ) 是幅角。极坐标形式在处理旋转和波动问题时非常有用。
复数的指数：复数可以通过欧拉公式与指数函数结合，表示为 ( $e^{i\theta} = \cos \theta + i \sin \theta$ )，这在傅里叶分析和信号处理中非常重要。

4.实现复数的运算类

import math

class Complex:
    def __init__(self, real=0.0, imag=0.0):
        self.real = real
        self.imag = imag

    def __add__(self, other):
        return Complex(self.real + other.real, self.imag + other.imag)

    def __sub__(self, other):
        return Complex(self.real - other.real, self.imag - other.imag)

    def __mul__(self, other):
        return Complex(self.real * other.real - self.imag * other.imag,
                       self.real * other.imag + self.imag * other.real)

    def __truediv__(self, other):
        denom = other.real**2 + other.imag**2
        if denom == 0:
            raise ValueError("Cannot divide by zero")
        return Complex((self.real * other.real + self.imag * other.imag) / denom,
                       (self.imag * other.real - self.real * other.imag) / denom)

    def modulus(self):
        return math.sqrt(self.real**2 + self.imag**2)

    def conjugate(self):
        return Complex(self.real, -self.imag)

    def __str__(self):
        return f"{self.real} + {self.imag}i"

    def __repr__(self):
        return f"Complex({self.real}, {self.imag})"

    def argument(self):
        return math.atan2(self.imag, self.real)

    def to_polar(self):
        r = self.modulus()
        theta = self.argument()
        return (r, theta)

    def from_polar(r, theta):
        return Complex(r * math.cos(theta), r * math.sin(theta))

    def __eq__(self, other):
        return self.real == other.real and self.imag == other.imag

# 示例
if __name__ == "__main__":
    z1 = Complex(3, 4)
    z2 = Complex(1, 2)

    print("z1:", z1)
    print("z2:", z2)
    print("z1 + z2:", z1 + z2)
    print("z1 - z2:", z1 - z2)
    print("z1 * z2:", z1 * z2)
    print("z1 / z2:", z1 / z2)
    print("z1 的模:", z1.modulus())
    print("z1 的共轭:", z1.conjugate())
    print("z1 的极坐标:", z1.to_polar())

功能说明

加法：__add__ 实现复数的加法。
减法：__sub__ 实现复数的减法。
乘法：__mul__ 实现复数的乘法。
除法：__truediv__ 实现复数的除法，包含了对分母为零的处理。
模：modulus 计算复数的模。
共轭：conjugate 返回复数的共轭。
字符串表示：__str__ 和 __repr__ 方法用于打印复数。
幅角：argument 计算复数的幅角（弧度）。
极坐标：to_polar 方法返回复数的极坐标形式。
极坐标转换：from_polar 类方法（静态方法）根据极坐标创建复数。
比较：__eq__ 比较两个复数是否相等。

可以根据需求扩展更多的功能，例如复数的指数函数、对数等。

当：

#当：z1和z2
z1 = Complex(8, 1)
z2 = Complex(1, 2)

运行结果：

z1: 8 + 1i
z2: 1 + 2i
z1 + z2: 9 + 3i
z1 - z2: 7 + -1i
z1 * z2: 6 + 17i
z1 / z2: 2.0 + -3.0i
z1 的模: 8.06225774829855
z1 的共轭: 8 + -1i
z1 的极坐标: (8.06225774829855, 0.12435499454676144)