 2-H Relations among Filter Length N, Transition Band, and Accuracy

Slides:



Advertisements
Similar presentations
做中国梦 走特色路 —— 宁波电大业余党校时政课 林志标 四川雅安地震 2013 年 4 月 20 日 8 时 02 分四川省雅安市芦山县(北纬 30.3, 东 经 )发生 7.0 级地震。震源深度 13 公里。震中距成都约 100 公里。成都、重庆及陕西的宝鸡、汉中、安康等地均有较.
Advertisements

海南省疾病预防控制中心. (一)基本情况  工作用房面积: ㎡,其中实验室使用面积为 6500 ㎡  中心定编 213 人,其中全额预算编制 193 人,自筹编制 20 人  现有在职职工 320 名,其中专业技术人员占 84.3% 。 人性化的办公场所实验室区域 一、海南省疾病预防控制中心概况.
663 Chapter 14 Integral Transform Method Integral transform 可以表示成如下的積分式的 transform  kernel Laplace transform is one of the integral transform 本章討論的 integral.
仪 容. 一、化妆的技巧 眼部的化妆 唇部化妆 眉部化妆 鼻部化妆 根据脸型化妆 根据脸型选发型.
Final Review Chapter 1 Discrete-time signal and system 1. 模拟信号数字化过程的原理框图 使用 ADC 变换器对连续信号进行采样的过程 使用 ADC 变换器对连续信号进行采样的过程 x(t) Analog.
H7N9 禽流感. H7N9 流感确诊病例主要表现 1 、起病急; 2 、病程早期均有高热 (38 ℃以上 ) ,伴咳嗽等呼 吸道感染症状,起病 5-7 天出现呼吸困难; 3 、典 型的病毒性肺炎,重症肺炎并进行性加重,部分 病例可迅速发展为急性呼吸窘迫综合症并死亡。
第七章 获利能力分析. 第一节 获利能力分析概述 获利能力的内涵 获利能力(盈利能力)是指企业获取利润的能力。 评价方法: ①利润与销售收入之间的比率 ②利润与资产之间的比率.
粉尘作业和使用场所防范 粉尘爆炸安全大检查注意事项
人感染H7N9禽流感医院感染 预防与控制技术指南
传染病预检分诊工作要求 发热门诊管理要求.
五專醫護類科介紹 樹人醫專 職業教育組 李天豪 組長.
情景应变 V1 段超.
针灸治疗学讲稿 山东中医药大学 高树中.
做好学校甲型H1N1流感防控工作 确保师生身体健康
H7N9禽流感相关知识
甘肃4班面试专项练习4 应急应变 主讲: 凌宇 时间:6月3日.
对本书、视频等任何MATLAB问题,作者做到有问必答!
只要大家共同努力,禽流感是可以預防的疾病。
菏泽市初中历史水平考试备考研讨与交流 菏泽市教研室 张红霞.
第七章 认识大洲 第一讲 亚洲及欧洲.
挖掘市场预期分布 建立有效投资策略 权证市场2006年中期投资策略
国王赏麦的故事.
歡迎蒞臨 三年八班大家族 導師:陳冠諠老師 16個帥氣寶貝 16個漂亮寶貝.
学校突发公共卫生事件 的应急处置 主讲人:汕头市中小学卫生保健所 姚广荣.
情緒行為障礙之教學與輔導 新竹縣情緒障礙巡迴教師 陳弘念.
人力資源管理委員會 主席:魏麗香部長 執秘:董家檥督導 委員:林姿伶HN、黃士豪HN、潘秋華HN 林素琴專師組長、卓惠瑄、張維恩、王孟萱、
第五組 幼兒安全與衛生教育 組員: 譚郁馨 張喻晴 沈恩華
我班最喜愛的零食 黃行杰.
10.2 分子动理论的初步知识 蒙城县乐土中学 袁亮.
我的過動人生 圖.文: 吳沁婕.
指導教授:林君憶老師 報告者:江秀英 日期: (三)
《中华人民共和国传染病防治法》部分知识 河西区卫生局.

八桥初中九年级思想品德课复习导学案之五---
建議題.
第二章 控制系统的数 学模型 烟台大学光电学院.
XI. Hilbert Huang Transform (HHT)
III. Gabor Transform III-A Definition Standard Definition:
IV. Implementation IV-A Method 1: Direct Implementation 以 STFT 為例
Rate and Distortion Optimization for Reversible Data Hiding Using Multiple Histogram Shifting Source: IEEE Transactions On Cybernetics, Vol. 47, No. 2,February.
聲音檔和 Video 檔的讀與寫 (by Matlab)
V. Homomorphic Signal Processing
數 學 基 礎 2 ※ 本章主要目的 1. 介紹拉氏轉換的基本理論。 2. 舉例說明應用拉氏轉換求解線性常微分方程式的方法。
第三章 影像顯示 3.1 前言 看看空間解析度和量化會如何影響影像的顯示,尤其是影像品質的部分,各種影像的屬性如何影響到顯示的品質。
滤波器设计matlab相关函数.
Chapter 11 Orthogonal Functions and
On Some Fuzzy Optimization Problems
第五章 数字滤波器设计 Filtering Beijing Institute of Technology 数字信号处理.
6.5 数字高通、带通和带阻滤波器的设计.
II. Short-time Fourier Transform
冠词的特殊位置(续).
何清波 博士 副教授 中国科学技术大学 精密机械与精密仪器系 安徽合肥 电话:
Digtlal Signal Processing —— Using MATLAB
Advanced Digital Signal Processing 高等數位訊號處理
XIV. Orthogonal Transform and Multiplexing
VII. Data Compression (A)
期末考.
聖方濟各英文小學 升中派位結果(2002/2004) 入讀英文中學:95.9% 第一組別(Band 1)學生:80.2%
Chapter 15 檔案存取 LabVIEW中的檔案存取函數也可將程式中的資料儲存成Excel或Word檔。只要將欲存取的檔案路徑位址透過LabVIEW中的路徑元件告訴檔案存取函數後,LabVIEW便可將資料存成Excel或Word檔;當然也可以將Excel或Word檔的資料讀入LabVIEW的程式中。
第三章 預測.
IV. Implementation IV-A Method 1: Direct Implementation 以 STFT 為例
Introduction to Matlab
本講義為使用「訊號與系統,王小川編寫,全華圖書公司出版」之輔助教材
認識H1N1 盧亞人醫院 感控護士 劉秀屏.
本講義為使用「訊號與系統,王小川編寫,全華圖書公司出版」之輔助教材
新高中通識教育科課堂的 教學規劃和應試訓練
II. Short-time Fourier Transform
Principle and application of optical information technology
認識﹋禽流感*.
Section 1 Basic concepts of web page
Presentation transcript:

 2-H Relations among Filter Length N, Transition Band, and Accuracy  Suppose that we want:  passband ripple  1,  stopband ripple  2,  width of transition band  F (expressed by normalized frequency) F = (f1  f2)/fs = (f1  f2)T (fs: sampling frequency, T: sampling interval) Then, the estimated length N of the digital filter is: . When there are two transition bands, F = min(F1, F2) 犧牲 transition band 的 frequency response, 換取較高的 passband and stopband accuracies

[Ref] F. Mintzer and L. Bede, “Practical design rules for optimum FIR bandpass digital filter”, IEEE Trans. ASSP, vol. 27, no. 2, pp. 204-206, Apr. 1979. 問題:假設 1 = 2 =  ,N 為固定, 當 F 變為 A 倍時,  變為多少?

N = 17, k+2=10 transition band: 0.24 ~ 0.26 N = 17, boundary of transition band N = 17, transition band: 0.2 ~ 0.3

 2-I Relations between Weighting Functions and Accuracy If we treat the passband more important than the stop band W(f) = 1 in the passband, 0 < W(f) < 1 in the stopband If we treat the stop band more important than the pass band 0 < W(f) < 1 in the passband, W(f) = 1 in the stopband

Larger error near the transition band Hd(F) H(F)

 2-J FIR Filter in MSE Sense with Weighting Functions 可對照 pages 47~49 F = f / fs n = 0 ~ k when n   問題:

τ = 0 ~ k, n = 0 ~ k 可以表示成 (k+1)  (k+1) matrix operation B S = C S = B−1C

思考: Is it possible to apply the transition band to the FIR filter in the MSE sense?

 2-K Four Types of FIR Filter h[n] = 0 for n < 0 and n  N 點數為 N  Type 1 之前的方法只討論到Type 1 h[n1] = h[n2 − n] and N is odd. (even symmetric) k = (N1)/2

 Type 1: h[n] = h[N−1− n] (even symmetric) and N is odd.  Type 2: h[n] = h[N−1− n] (even symmetric) and N is even.  Type 3: h[n] = −h[N−1− n] (odd symmetric) and N is odd.  Type 4: h[n] = −h[N−1− n] (odd symmetric) and N is even. k = (N−1)/2

 Type 2: When h[n] = h[N−1− n] and N is even: (even symmetric) 令 r[n] = h[n + k], where k = (N1)/2 (注意此時 k 不為整數) (比較 page 41) h[0] h[1] h[k-3/2] h[k-1/2] h[k+1/2] h[k+3/2] h[N-2] h[N-1] r[-k] r[-k+1] r[-3/2] r[-1/2] r[1/2] r[3/2] r[k-1] r[k] s[1]/2 s[2]/2 s[k-1/2]/2 s[k+1/2]/2

當 n = 1, 2, …, k+1/2 設計出 s[n] 之後 r[n] = s[n+1/2]/2, h[n] = r[n − k],

Design Method for Type 2 由於 n 和 n+1 兩項相加可得 所以可以「判斷」 R(F) 能被改寫成 求 s1[n] 和 s[n] 之間的關係

(令 k1 + 1/2 = k) 比較係數可得 for n = 2, 3, …., k−1/2

只需將 pages 52-55 的方法當中,Hd(F) 換成 W(F) 換成 k 換成 k − 1/2 = N/2 − 1 注意s1[n] 和 s[n] 之間的關係即可

Design Method for Type 3 由於 n−1 和 n+1 兩項相減可得 所以「判斷」可將 R(F) 改寫為

令 k1 = k − 1, 比較係數可得 for n = 2, 3, …., k−2

將 page 52-55 的方法當中,Hd(F) 換成 W(F) 換成 k 換成 k − 1 注意s1[n] 和 s[n] 之間的關係即可 思考:Design Method for Type 4

附錄三:Matlab 寫程式需注意的地方 一、各種程式語言寫程式共通的原則 二、Matlab 寫程式特有的技巧 (1) 能夠不在迴圈內做的運算,則移到迴圈外 (2) 寫一部分即測試,不要全部寫完再測試 (縮小範圍比較容易 debug) (3) 先測試簡單的例子,成功後再測試複雜的例子 二、Matlab 寫程式特有的技巧 (1) 迴圈能避免就儘量避免 (2) 儘可能使用 Matrix 及 Vector operation Example: 由 1 加 到100,用 Matlab 一行就可以了 sum([1:100]) 完全不需迴圈

三、一些重要的 Matlab 指令 (1) function: 放在第一行,可以將整個程式函式化 (2) tic, toc: 計算時間 (3) find: 找尋一個 vector 當中不等於 0 的entry 的位置 範例:find([1 0 0 1]) = [1, 4] find(abs([-5:5])<=2) = [4, 5, 6, 7, 8] (因為 abs([-5:5])<=2 = [0 0 0 1 1 1 1 1 0 0 0]) (4)  : Hermitian (transpose + conjugation), .  : transpose (5) imread: 讀圖, image, imshow, imagesc: 將圖顯示出來, (註: 較老的 Matlab 版本 imread 要和 double 並用 A=double(imread(‘Lena.bmp’)); (6) imwrite: 製做圖檔

(7) xlsread: 由 Excel 檔讀取資料 (8) xlswrite: 將資料寫成 Excel 檔 (9) aviread: 讀取 video 檔 (10) dlmread: 讀取 *.txt 或其他類型檔案的資料 (11) dlmwrite: 將資料寫成 *.txt 或其他類型檔案