日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程语言 > python >内容正文

python

第4章 Python 数字图像处理(DIP) - 频率域滤波4 - 单变量的离散傅里叶变换DFT

發布時間:2023/12/10 python 26 豆豆
生活随笔 收集整理的這篇文章主要介紹了 第4章 Python 数字图像处理(DIP) - 频率域滤波4 - 单变量的离散傅里叶变换DFT 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

目錄標題

    • 單變量的離散傅里葉變換
      • 由取樣后的函數的連續變換得到DFT
      • 取樣和頻率間隔的關系

單變量的離散傅里葉變換

由取樣后的函數的連續變換得到DFT

對原函數的變換取樣后的業的發展的變換F~(μ)\tilde F(\mu)F~(μ),但未給出取樣后的函數f~(t)\tilde f(t)f~?(t)的變換F~(μ)\tilde F(\mu)F~(μ)的表達式。
F~(μ)=∫?∞∞f~(t)e?j2πμtdt(4.39)\tilde F(\mu) = \int_{-\infty}^{\infty} \tilde f(t) e^{-j2\pi\mu t} dt\tag{4.39}F~(μ)=??f~?(t)e?j2πμtdt(4.39)

F~(μ)=∫?∞∞f~(t)e?j2πμtdt=∫?∞∞∑n=?∞∞f(t)δ(t?nΔT)e?j2πμtdt=∑n=?∞∞∫?∞∞f(t)δ(t?nΔT)e?j2πμtdt=∑n=?∞∞fne?j2πμnΔT(4.40)\begin{aligned} \tilde F(\mu) & = \int_{-\infty}^{\infty} \tilde f(t) e^{-j2\pi\mu t} dt = \int_{-\infty}^{\infty} \sum_{n = -\infty}^{\infty} f(t) \delta(t - n\Delta T) e^{-j2\pi\mu t} dt\\ & = \sum_{n = -\infty}^{\infty} \int_{-\infty}^{\infty} f(t) \delta(t - n\Delta T) e^{-j2\pi\mu t} dt \\ & = \sum_{n = -\infty}^{\infty} f_n e^{-j2\pi\mu n \Delta T} \end{aligned} \tag{4.40}F~(μ)?=??f~?(t)e?j2πμtdt=??n=??f(t)δ(t?nΔT)e?j2πμtdt=n=????f(t)δ(t?nΔT)e?j2πμtdt=n=??fn?e?j2πμnΔT?(4.40)

μ=mMΔT,m=0,1,2,?,M?1(4.41)\mu = \frac{m}{M\Delta T}, \quad m = 0, 1, 2, \cdots, M-1 \tag{4.41}μ=MΔTm?,m=0,1,2,?,M?1(4.41)

下面表達式就是我們所求的離散傅里葉變換
Fm=∑n=0M?1fne?j2πμnm/M,m=0,1,2,?,M?1(4.42)F_m = \sum_{n = 0}^{M - 1} f_n e^{-j2\pi\mu n m/M}, \quad m = 0, 1, 2, \cdots, M-1 \tag{4.42}Fm?=n=0M?1?fn?e?j2πμnm/M,m=0,1,2,?,M?1(4.42)

離散傅里葉反變換
fn=1M∑m=0M?1Fmej2πμnm/M,n=0,1,2,?,M?1(4.43)f_n = \frac{1}{M}\sum_{m = 0}^{M - 1} F_m e^{j2\pi\mu n m/M}, \quad n = 0, 1, 2, \cdots, M-1 \tag{4.43}fn?=M1?m=0M?1?Fm?ej2πμnm/M,n=0,1,2,?,M?1(4.43)

一般二維情況下,使用xxxyyy表示圖像坐標變量并使用uuuvvv表示頻率變量更為直觀。離散傅里葉變換對可以改寫為
F(u)=∑x=0M?1f(x)e?j2πux/M,u=0,1,2,?,M?1(4.44)F(u) = \sum_{x = 0}^{M - 1} f(x) e^{-j2\pi u x/M}, \quad u = 0, 1, 2, \cdots, M-1 \tag{4.44}F(u)=x=0M?1?f(x)e?j2πux/M,u=0,1,2,?,M?1(4.44)

離散傅里葉反變換
f(x)=1M∑u=0M?1F(u)ej2πux/M,x=0,1,2,?,M?1(4.45)f(x) = \frac{1}{M}\sum_{u = 0}^{M - 1} F(u) e^{j2\pi u x/M}, \quad x = 0, 1, 2, \cdots, M-1 \tag{4.45}f(x)=M1?u=0M?1?F(u)ej2πux/M,x=0,1,2,?,M?1(4.45)

取樣和頻率間隔的關系

# 例4.4 計算DFT # x = 0, 1, 2, 3 # f(x) = 1, 2, 4, 4 x = np.arange(4) y = np.array([1, 2, 4, 4]) fft = np.fft.fft(y) print('DFT') print(fft)# IDFT反變換 ifft = np.fft.ifft(fft) print('IDFT') print(ifft) DFT [11.+0.j -3.+2.j -1.+0.j -3.-2.j] IDFT [1.+0.j 2.+0.j 4.+0.j 4.+0.j]

總結

以上是生活随笔為你收集整理的第4章 Python 数字图像处理(DIP) - 频率域滤波4 - 单变量的离散傅里叶变换DFT的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。