![ba575128180f65e930dea8b9ce0ef1d0.png](https://i-blog.csdnimg.cn/blog_migrate/926aba64172d47954b3ac1f5e4108d4b.jpeg)
实矩阵也可能碰到复特征值,因此无可避免地在矩阵运算中碰到复数矩阵当然也有可能包含复数,最重要的复矩阵是傅立叶矩阵,它用于傅立叶变换。一种特殊的傅立叶变换是快速傅立叶变换(fast Fourier transform),简称FFT,在计算机中很常用,特别是涉及到大数据时,FFT将把傅立叶变换中的n阶方正阵乘法的运算次数从
本文相关前置知识
复数和复平面、复平面上的旋转
傅立叶矩阵中w的形式是什么
正交向量和正交矩阵 (格拉姆-施密特正交化)
对称矩阵的特征值
1 复向量
先给出一个复向量,即向量的分量中至少有一个是复数:
![a508eec6a2f7a10713ba0b6b365b07c3.png](https://i-blog.csdnimg.cn/blog_migrate/d76f2d73b44e2ad4fedbdc98883b1710.png)
虽然这个向量在表达上和普通的实向量没什么区别,但这个向量不再属于实空间
1.1 模长
关于复向量的第一个问题是模长怎么计算?
由于向量中有复数分量,再用过
![46b4f512881d1b34031c9a27f4388070.png](https://i-blog.csdnimg.cn/blog_migrate/4c59c96050967b8e15766a693c276b92.jpeg)
但很明显,
![0e967f2d17d460269f3b69f17057d9e6.png](https://i-blog.csdnimg.cn/blog_migrate/bfb6fa4ffaa4795d2d59972a2534d2ac.png)
我们知道一个复数的模长的平方等于这个复数与它的共轭复数的乘积,因此可以通过下面的方式计算复向量的模长:
![e1faadf837a0968199999b76eaa2da40.png](https://i-blog.csdnimg.cn/blog_migrate/ceed0ae8f3e28800f2d93e609d5c5c9c.png)
通常用
![216005f52129af7ab84cc1db0d0af927.png](https://i-blog.csdnimg.cn/blog_migrate/6a51910b60e15fdc8feac197c3a1136a.png)
1.2 点积
与模长类似,如果有两个复向量p和q,它们的点积也不能简单地定义成
![ca83cca8252be7abd38ec7a4c557955d.png](https://i-blog.csdnimg.cn/blog_migrate/2be662e2dbdc74c03633a76e30da2ed2.png)
复平面上有两个向量
![e62e8c79cf8f721914dd2d19e46a4c04.png](https://i-blog.csdnimg.cn/blog_migrate/b2824adabc5d49b7c9435881404dccee.png)
二者的点积是0,因此可以判断两个复向量互相垂直:
![568086e39fb59edfe57d139e4c5252c4.png](https://i-blog.csdnimg.cn/blog_migrate/afc4c23baf75b8c8fc19672b25314b16.png)
2 复矩阵
我们曾讲过,对于一个矩阵A来说,如果
2.1 厄米特矩阵
如果一个复矩阵是对称矩阵,那么:
![f404e9435f23d141d9f3ca8074824b70.png](https://i-blog.csdnimg.cn/blog_migrate/3a86a6b08f91b4726fda1b39a9763d2a.png)
通常写作另一种方式:
![458279d5877f1b59077a589afa73d5ef.png](https://i-blog.csdnimg.cn/blog_migrate/d0865dff78f01e7d33a37ec105534160.jpeg)
这种对称复矩阵称为厄米特矩阵(或埃米特矩阵,Hermitian matrix),比如下面这个:
![1fe70476d72072d2d3f07652f9c2c291.png](https://i-blog.csdnimg.cn/blog_migrate/4d2888420e5b4b9b3810837030928912.png)
2.2 酉矩阵
上一章讲到,一个复对称矩阵的特征值仍然是实数,且可以找到互相垂直的特征向量,其对角元素都是实数。 假设有一个由n个标准正交向量组成的复矩阵
![be0c09c42c2d962c9bfa1d6602d1584a.png](https://i-blog.csdnimg.cn/blog_migrate/b68b0a7316731a58bfd62b089bf09675.jpeg)
这个复空间的正交矩阵Q称为酉矩阵(unitary matrix)。
3 傅立叶矩阵和傅立叶快速变化
傅立叶变换是一种分析信号的方法,它可分析信号的成分,也可用这些成分合成信号。许多波形可作为信号的成分,比如正弦波、方波、锯齿波等,傅立叶变换用正弦波作为信号的成分。
在电子工程或计算机中,n×n矩阵的行和列都是从0开始的,到n-1结束,由于傅立叶变换经常用在计算机上,所以我们在讨论傅立叶矩阵的时候遵从这种下标规则。
3.1 傅立叶矩阵
形如
![685f519c65cb5cad222236c93c2fbd23.png](https://i-blog.csdnimg.cn/blog_migrate/51dc580414e045539eb58d7dbd8fabd6.png)
矩阵中的每个元素都不为0,是个全矩阵。w是个特殊的值:
![89da3113f02bd6e04da6ebc420958cf4.png](https://i-blog.csdnimg.cn/blog_migrate/f62bc3554c7df1238a40ddf9d29e081f.png)
在计算w的乘方的时候,我们需要考虑用极坐标表示复平面。在极坐标下,w表示模长为1的向量从(1,0)开始,绕原点逆时针旋转了
![ef81debd1ccf0739dac734b9b0973c1b.png](https://i-blog.csdnimg.cn/blog_migrate/f02de1a9f38b7a9f42af271926ecae5d.png)
同理,n=4时,
![499f84ce0bdee7fa671d816e7b32c91f.png](https://i-blog.csdnimg.cn/blog_migrate/d266f9aee21c57a24d6368e75fbf148c.png)
傅立叶矩阵可以得到一个四点(离散的)傅立叶变换,它的逆矩阵可以得到傅立叶逆变换。此外,傅立叶矩阵的列向量是正交的,所以很容易求得逆矩阵。实际上傅立叶矩阵可以分解成一系列稀疏矩阵,这些矩阵有大量的0元素,所以相应的乘法和求逆都很简单。
![3d669b855891411933079c03923785d1.png](https://i-blog.csdnimg.cn/blog_migrate/02b9ea3c79565118640927dd8abfbd6c.png)
上一章介绍过,复向量的点积不是这么算的,正确算法应该是取共轭的转置:
![4f4210a7550961570a99aeeb338fd65b.png](https://i-blog.csdnimg.cn/blog_migrate/b482179dc494f7aad83d68c14f207dfa.png)
![248926801438892876b5953b98e36b64.png](https://i-blog.csdnimg.cn/blog_migrate/fbbe43d11ff03a79796e392866c60797.jpeg)
对于标准正交的实矩阵来说,矩阵的逆等于矩阵的转置,傅立叶矩阵可化简为标准正交的复矩阵,具有同样的性质,
![f63bfcd410211c16c9def201f818d831.png](https://i-blog.csdnimg.cn/blog_migrate/69f58b2352c259bc6f7f2a48ea4306b1.jpeg)
由于
3.2 快速傅立叶变换
什么是快速傅立叶变换呢?举个例子,
以
![52547d739738fb240d7ecf37b9af002d.png](https://i-blog.csdnimg.cn/blog_migrate/5bc56cd97c85538d8fff1436fec7dd37.png)
![014a73e86a17ca6f7a75c12916107b17.png](https://i-blog.csdnimg.cn/blog_migrate/e0e1765489c224aa3e48fdedabe52ae0.png)
既然如此,
![f86c6f13cc31ed4e03ec745c17f0106a.png](https://i-blog.csdnimg.cn/blog_migrate/ad54d173329eceab0c5b50bc0b8c459b.jpeg)
这种分解称为快速傅立叶变换。其中P是一个
![d370410025ace4ebfb899d405aa4070d.png](https://i-blog.csdnimg.cn/blog_migrate/03e906cce74e2eb061b294073d505a51.png)
P的效果是使得所乘行向量x中序号为奇数的分量
![da2f9ffb4c9655f2cdcaf355d4a35157.png](https://i-blog.csdnimg.cn/blog_migrate/8ac60d52a80dfb4b16fdc47c3e78f11a.png)
可以看到,快速傅立叶变换实际上使用的是分治算法。计算64阶傅立叶变换的计算量是