[傅里葉變換及其應用學習筆記] 二十四. 級聯,脈沖響應

?

我們上節課學習了

  1. 在離散有限維空間中,任何線性系統都是通過矩陣間的相乘得到的
  2. 在連續無限維空間中,任何線性系統都是通過對核函數的積分得到的

?

?

脈沖響應(impulse response)

?

級聯線性系統(Cascading linear system)

?

如果$L$與$M$都是線性的,有

$w=MLv$

在連續無限維空間中

$\begin{align*}
MLv
&=M\left( \int_{-\infty}^{\infty}k(x,y)v(y)dy \right )\\
&\approx M\left( \sum_{i=-\infty}^{\infty}k(x,y_i)v(y_i)\Delta y_i \right )\\
&=\sum_{-\infty}^{\infty}M\left( k(x,y_i)v(y_i)\Delta y_i \right )\\
&=\sum_{-\infty}^{\infty}M_x \left( k(x,y_i)v(y_i)\Delta y_i \right ) \qquad M\ deal\ with\ the\ function\ that\ Lv\ ouput\ based\ on\ x \\
&\approx \int_{-\infty}^{\infty}M_x(k(x,y))v(y)dy
\end{align*}$

?

?

脈沖響應的定義

上述關于級聯的討論是為了引出下面的這個結論:

任何線性系統都由對核(函數)的積分得到,核(函數)就是該線性系統對脈沖函數的響應。(Any linear system is given by integration against a kernel (impulse response).)

?

推導過程如下:

$\begin{align*}
v(x)
&=(\delta * v)(x)\\
&=\int_{-\infty}^{\infty}\delta(x-y)v(y)dy \qquad(\delta\ shift\ property)
\end{align*}$

那么線性系統有如下表示

$\begin{align*}
Lv(x)
&=L\left( \int_{-\infty}^{\infty}\delta(x-y)v(y)dy \right )\\
&=\int_{-\infty}^{\infty}L_x\delta(x-y)v(y)dy
\end{align*}$

令$h(x,y) = L_x\delta(x-y)$,則有,

$\displaystyle{ Lv(x) = \int_{-\infty}^{\infty}h(x,y)v(y)dy }$

?

其中$h(x,y)$是該線性系統的核函數,它由$L_x\delta(x-y)$得到,同時他也是該線性系統的脈沖響應。

脈沖響應的定義如下

  • $\delta(x-y)$是位置在$y$上的脈沖,$L_x\delta(x-y)$表示了把脈沖輸入到該線性系統,此時系統會做出響應,并輸出脈沖響應$h(x,y)$。

?

?

?

Schwartz核函數定理

如果$L$是廣義函數(分布)的一個線性算符,即$L$在符合疊加性原則的基礎上將一個廣義函數變換為另一個廣義函數,那么就會存在唯一的核$k$,使得$Lv = <k,v>$

?

?

傅里葉變換的脈沖響應

當輸入脈沖函數$\delta(x-y)$,傅里葉變換會輸出

$h(x,y) = \mathcal{F}(\delta(x-y)) = e^{-2\pi ixy}$

另外,傅里葉變換的公式如下

$\mathcal{F}f(x) = \displaystyle{ \int_{-\infty}^{\infty}e^{-2\pi ixy}f(y)dy }$

它的核函數為

$k(x,y) = e^{-2\pi ixy}$

我們注意到,核函數與脈沖響應式一樣的,其中有如下關系:

  • 如果一個線性算符能表示成對于$k(x,y)$與輸入函數乘積的積分形式,那么$k(x,y)$就是脈沖響應了。
  • 反過來說,如果我們能得到線性系統的脈沖響應,就能通過對脈沖響應和輸入函數的乘積進行積分來表達該線性系統的線性算符。

?

?

離散有限維線性系統的脈沖響應

在連續無限維的線性系統中,脈沖響應是線性系統對輸入脈沖$\delta(x-y)$的響應。在離散有限維線性系統也同樣是對輸入脈沖序列的響應,用矩陣乘法的表達如下:

$A\cdot \left[ \underline{\delta}_0,\underline{\delta}_1,\underline{\delta}_2,...,\underline{\delta}_{n-1} \right]
=A\cdot\begin{bmatrix}
1 &0? &0? &...? &0 \\
0 &1? &0? &...? &0 \\
0 &0? &1? &...? &0 \\
\vdots &\vdots? &\vdots? &...? &\vdots \\
0 &0? &0? &...? &1
\end{bmatrix}=A$

?

?

脈沖響應的例子

開關

$Lv = \Pi v$

它的脈沖響應為

$h(x,y) = L\delta(x-y) = \Pi(x)\delta(x-y) = \Pi(y)\delta(x-y) \qquad(\delta\ sampling\ property)$

對脈沖響應與輸入函數乘積的積分會得到開關的線性算符

$\begin{align*}
\int_{-\infty}^{\infty}h(x,y)v(y)dy
&=\int_{-\infty}^{\infty}\Pi(y)\delta(x-y)v(y)dy\\
&=\int_{-\infty}^{\infty}\delta(x-y)\left( \Pi(y)v(y) \right )dy\\
&=\left(\delta * (\Pi v) \right )(x) \qquad (\delta\ shift\ property)\\
&=\Pi(x)v(x)
\end{align*}$

這個結果證明了前面的結論是正確的

?

?

?

卷積,連續無限維線性時不變系統

引入時移符號$\tau$

$\tau_a v(x) = v(x-a)$。

?

假設有用卷積表達的線性系統如下

$Lv = h*v$

如果我們對輸入$v$進行延時$a$

$\begin{align*}
L\tau_a v(x)
&= (h*\tau_av)(x)\\
&= \int_{-\infty}^{\infty}h(x-y)v(y-a)dy \\
&= \int_{-\infty}^{\infty}h(x-z-a)v(z)dz \qquad (letting\ z=y-a) \\
&= \int_{-\infty}^{\infty}\left(h(x-z)*\delta(x-a) \right )v(z)dz \qquad (\delta\ shift\ property)\\
&= \left(\int_{-\infty}^{\infty}h(x-z)v(z)dz \right )*\delta(x-a)\\
&= (h*v)(x-a)\\
&= \tau_a(h*v)(x)\\
&= \tau_aLv(x)
\end{align*}$

結果顯示該線性系統輸出的延時與輸入的延時同為$a$,這被稱為線性時不變系統(Linear Time Invariant System)。

結論是:

  • 如果一個線性系統是由卷積給定的,那么他就是時不變的。

?

反過來也是成立的:

  • 如果一個線性系統是時不變的,那么它一定是由卷積給定的。

證明過程如下:

任何連續無限維線性系統都有如下表示

$Lv = \displaystyle{ \int_{-\infty}^{\infty}L_x\delta(x-y)v(y)dy }$

我們令$h(x) = L_x\delta(x)$,就是線性系統對$\delta_0$進行脈沖響應。則有

$L_x\delta(x-y) = L_x(\tau_y\delta(x))$

如果$L$是時不變的,則輸出與輸入會有同一延時

$L_x\delta(x-y) = L_x(\tau_y\delta(x)) = \tau_y(L_x\delta(x)) = \tau_yh(x) = h(x-y)$

$Lv = \displaystyle{ \int_{-\infty}^{\infty}h(x-y)v(y)dy }$

轉載于:https://www.cnblogs.com/TaigaCon/p/5138144.html

本文來自互聯網用戶投稿,該文觀點僅代表作者本人,不代表本站立場。本站僅提供信息存儲空間服務,不擁有所有權,不承擔相關法律責任。
如若轉載,請注明出處:http://www.pswp.cn/news/375157.shtml
繁體地址,請注明出處:http://hk.pswp.cn/news/375157.shtml
英文地址,請注明出處:http://en.pswp.cn/news/375157.shtml

如若內容造成侵權/違法違規/事實不符,請聯系多彩編程網進行投訴反饋email:809451989@qq.com,一經查實,立即刪除!

相關文章

WPF如何實現TreeView節點重命名

我們經常看到一些軟件比如酷狗音樂&#xff0c;在對列表右鍵進行重命名的時候&#xff0c;當前列表會泛白并且進入可編輯狀態&#xff0c;當我們更改完成后就會并進入非編輯狀態&#xff0c;這些具體是怎么實現的呢&#xff1f;下面的方法也許會提供一些思路&#xff0c;下面的…

8086地址傳送指令LEA

LEA REG,OPRD ;操作數OPRD必須是一個存儲器操作數 LEA AX,IDATA ;把IDATA的偏移地址傳送到AX寄存器中DATA SEGMENTIDATA DW 1,2,3,4 DATA ENDS CODE SEGEMNT BEG:MOV AX,OFFSET IDATA ;AXIDATA的偏移地址LEA AX,IDATA ;AXIDATA的偏移地址LEA AX,DS:[IDATA] ;把ds:[IDA…

Shell --- 批量修改文件后綴腳本

for f in *.$1; dofilenamebasename $fmv $f "${filename%.*}".$2; done; Usage:&#xff1a; rename suffix rename_suffix eg: rename dat txt > ls > a.dat > rename dat txt > ls > a.txt 轉載于:https://www.cnblogs.com/RookieCoder/p/5140265.…

8086標志操作指令

標號傳送指令LAHF 把FLAG低八位送入AH。不影響FLAG的任何位LAHF ;把CF,PF,AF,ZF,SF送入AH的相應位即0,2,4,6,7位SAHF 把AH送入FLAG低八位。根據AH中的內容FLAG的低八位受到影響,高位不受影響MOV AH,11111111B SAHF ;(執行后CF,PF,AF,ZF,SF等于1)PUSHF 把FLAG壓入棧中。不影…

最大子序列和

一、遍歷所有可能情況 (1234...N)所有可能子序列如下&#xff1a; 1&#xff1b;12&#xff1b;123...... 2&#xff1b;23&#xff1b;234...... ...... N 共N趟&#xff0c;沒趟可能的情況由N&#xff0c;N-1...&#xff0c;1依次遞減。 時間復雜度O(N3)的算法&#xff1a; i…

!--注釋--和%--注釋--%有什么區別

轉載&#xff1a;版權所有&#xff1a;基礎軟件。作者郵箱&#xff1a;s.j.l.studiohotmail.com&#xff0c;sun.j.l.studiogmail.com。本文首發于 http://www.cnblogs.com/FoundationSoft。文章轉載請保持此版權信息并注明出處。 使用ASP.NET編輯頁面時&#xff0c;有兩種注釋…

團隊開發——用戶需求調研報告

用戶需求調研報告 項目名稱&#xff1a; 躲避小球 項目編號&#xff1a;001 調研主題&#xff1a; 用戶需求 訪談時間&#xff1a;2015.4.10 調研地點&#xff1a; 石家莊鐵道大學圖書館 訪談部門&#xff1a; 三個人行 參與人員&#xff1a; 林彥汝 1. 訪談目的 1、讓用…

8086地址傳送指令LDS,LES

LDS REG,OPRE;REG是一個16位寄存器,OPRD是一個存儲器操作數,為雙字項,高16位送入DS,低16位送入REG LDS BX,[SI](ds:si取得存儲器中32位(雙字),高16位作為段值送入DS中,低十六位送入指定的BX寄存器中)data segmentidata dd 12345678H data ends code segment assume cs:code,ds:…

設計模式(十五):解釋器模式

一、定義 在設定環境中&#xff0c;定義一種規則或者語法&#xff0c;通過解釋器來解釋規則或者語法的含義. 二、實例&#xff1a;將 二十一 —> 21 2.1 設定我們的環境 Context public class Context{public string Input { get; set; }public int Output { get; se…

8086地址傳送指令LEA,LDS,LES

LEA BX,[SI] ;BX<--DS:[SI]內存單元的偏移地址,與該單元中的數據無關 LDS BX,[SI] ;BX<--DS:[si]內存單元中的數據的低16位,DS<--內存單元中數據的高16位 LES BX,[SI] ;BX<--DS:[SI] ,ES<--DS:[SI2] 顯然LDS和LES的操作數是32位的(兩個字 四個字節 32位) 寫個…

MySQL 5.7.10 免安裝配置

# 配置環境&#xff1a;windows 64bit # 安裝版本&#xff1a;mysql-5.7.10-win32&#xff08;zip archive版本&#xff09; 1. ZIP Archive版是免安裝的&#xff0c;只需把mysql-5.7.10-win32.zip解壓到安裝目錄即可。 2. 在D:\Program Files\mysql-5.7.10-win32文件夾下新建配…

Cacti安裝教程

CentOS 6.0架設流量監控及集中日志系統第一章、cacti的安裝1、 系統的基本設置2、 設置主機名3、 [rootlocalhost ~]# vi /etc/sysconfig/network4、 [rootlocalhost ~]# exit5、 重新登陸6、 [rootGZ-Cacti ~]#7、 安裝screen后臺運行程序yum -y install screen8、 安裝優秀的…

8086條件轉移指令JE,JZ

JZ/JE 指令格式 轉移條件 說明 JZ/JE 標號 ZF1 等于0或相等轉移code segment assume cs:codemov ax,2mov bx,2cmp ax,bxje ok ;如果axbx 進行轉移(其實是檢測ZF的值進行轉移操作);根據ZF1 得出axbx ,所以才有axbx 相等則轉移 iout: mov ax,4c00Hint 21H ok: add a…

單片機C語言下LCD多級菜單的一種實現方法

摘要: 介紹了在C 語言環境下,在LCD 液晶顯示屏上實現多級嵌套菜單的一種簡便方法,提出了一個結構緊湊、實用的程序模型。 關鍵詞: 液晶顯示屏; 多級菜單; 單片機; C 語言; LCD  中圖分類號:TP311. 1    文獻標識碼:B 收稿日期:2005 - 11 - 210  引言 由于C 語…

fortran語法筆記

1&#xff0c;數據類型&#xff0c;fortran支持整形&#xff0c;real型&#xff0c;logical型&#xff0c;char型&#xff0c;復數型。整形分為為長整形和短整形定義長整形的方法 同時聲明多個變量的話可以用逗號隔開。 加兩個冒號的話可以直接在聲明的時候賦值。 fortran是唯一…

8086條件轉移指令JNE,JNZ

JNZ/JNE 指令格式 轉移條件 說明 JNZ/JNE 標號 ZF0 不等于0或不相等轉移code segment assume cs:codemov ax,5mov bx,2cmp ax,bx ;兩數比較 不相等則轉移 jne ok iout: mov ax,4c00Hint 21H ok: add ax,bxjmp iout code ends endAX0005 BX0002 CX0013 DX0000 SP00…

8086條件轉移指令JE/JZ,JNE/JNZ

JE/JZ,JNE,JNZ檢測的只是ZF的值,只要ZF滿足則進行轉移 JE的邏輯含義是相等則轉移 但是JE進行操作的是ZF1才執行所有相等則轉移是通過CMP指令來配合使用體現出的CMP AX,BX ;ax-bx 結果是否為0 {true ,ZF1} {false ,ZF0} JE OK (ZF1說明兩數相等) (ZF…

git 使用手冊

git 作為開源項目首選版本管理軟件&#xff0c;必然其優勢不容小覷&#xff0c;下面總結一下初次解除用到的一些技能 1. 基本命令介紹 git help branch 查看branch幫助文檔 git branch -a 列出當前所有分支 git branch good_proj 創建分支good_proj git checkout good_p…

WZY社區

WZY社區是我自己做的一個網站&#xff0c;后面會詳細更新&#xff0c;敬請關注&#xff01;轉載于:https://www.cnblogs.com/jxldjsn/p/5146682.html

條件轉移指令和無條件轉移指令練習

無條件轉移指令 JMP條件轉移指令 JE JNE JB 低于轉移(小于) JNB 不低于轉移(大于等于) JA 高于轉移(大于) JNA不高于轉移(小于等于)查找出數據中0&#xff5e;9的無符號數寫到itable單元中 0&#xff5e;9 大于等于0 小于等于9 用JNB和JNA最合適不過了 貼代碼~assume cs:code,d…