Performance Comparison of Various Discrete Transforms

Total Page:16

File Type:pdf, Size:1020Kb

Performance Comparison of Various Discrete Transforms Appendix A Performance Comparison of Various Discrete Transforms This will not include the fast algorithms, separability, recursivity, orthogonality and fast algorithms (complexity of implementation). These topics are described else- where in detail (see Chapter 3). The focus is on their various properties. We will consider the random vector x is generated by I-order Markov process. When T x ¼ ðÞx0; x1; ...; xNÀ1 , ðx0; x1; ...; xNÀ1 are the N random variables) correlation matrix ½Rxx is generated by the I-order Markov process, hi jjÀkj ½Rxx ¼ r ; r ¼ adjacent correlation coefficient |ffl{zffl} |fflfflffl{zfflfflffl} ðÞNÂN ðÞNÂN j; k ¼ 0; 1; ...; N À 1 (A.1) The covariance matrix ½S in the data domain is mapped into the transform domain ½S~ as is Âà à S~ ¼ ½DOT ½S ½DOT T (A.2) |{z} |fflffl{zfflffl} |{z} |fflfflfflfflfflfflffl{zfflfflfflfflfflfflffl} ðÞ ðÞNÂN ðÞNÂN N N ðÞNÂN DOT stands for discrete orthogonal transform. Superscripts T and * denote trans- pose and complex conjugateÂà respectively. When DOT is KLT, S~ is a diagonal matrix as all the transform coefficients in the KLT domain are uncorrelated. For all the other DOTs, residual correlation (correlation left undone in the DOT domain) is defined as [G6, G10] ! ! XNÀ1 XNÀ1 XNÀ1 XNÀ1 1 2 ~ 2 1 2 ~ 2 r ¼ kkS À Snn ¼ jjSmn À Snn (A.3) N n¼0 N m¼0 n¼0 n¼0 XNÀ1 XNÀ1 2 2 1 2 where jjSjj is the HilbertÀSchmidt norm defined as kkS ¼ jjSmn . N m¼0 n¼0 Note that N is the size of discrete signal. Smn is the element of ½S in row m and column n (m, n ¼ 0, 1, ..., NÀ1). 317 318 Appendix A Performance Comparison of Various Discrete Transforms For a 2D-random signal such as an image assuming that row and column statistics are independent of each other, the variances of the ðÞN  N samples can be easily obtained. This concept is extended for computing the variances of the ðÞN  N transform coefficients. A.1 Transform Coding Gain ½S ¼ Correlation or covariance matrix in data domain (see Section 5.6) ðÞNÂN Âà à à À S~ ¼ ½A ½S ð½A TÞ (Note that ð½A T Þ ¼ ½A 1 for unitary transforms) ðÞNÂN ðÞNÂN ðÞNÂN ðÞNÂN ¼ Correlation or covariance matrix in transform domain The transform coding gain, GTC is defined as N -1 1 σ~ 2 å kk N k=0 Arithmetic Mean GTC = = æ N -1 ö1/ N Geometric Mean ç σ~ 2 ÷ çÕ kk ÷ è k=0 ø s~2 ð¼; ; ...; À Þ where kk is the variance of the kth transform coefficient k 0 1 N 1 .As the sum of the variances in any orthogonal transform domain is invariant (total energy is preserved), GTC can be maximized by minimizing the geometric mean [B23]. The lower bound on the gain is 1 (as seen in Fig. A.1), which is attained only if all the variances are equal. 16-point regular DCT 9 Proposed integer 8 DCTs for H.264 and 7 AVS in [LA14] TC G 6 Ma’s integer DCT 5 Wien-Kuo’s integer DCT 4 Coding gain 3 2 1 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 Correlation coefficient ρ Fig. A.1 Comparing coding gain of orthogonal transforms [LA14] # 2009 IEEE A.4 Rate versus Distortion (Rate-Distortion) [B6] 319 A.2 Variance Distribution in the Transform Domain It is desirable to have few transform coefficients with large variances (this implies the remaining coefficients will have small variances, as the sum of the variances is invariant). This variance distribution can be described graphically or in a tabular form for N ¼ 8, 16, 32, ...and r ¼ 0.9, 0.95, ..., etc. The compaction of the energy in few transform coefficients can be represented by the normalized basis restriction error [B6] defined as NPÀ1 s~2 kk ¼ J ðÞ¼r k m ; m ¼ 0; 1; ...; N À 1Eq: ð5.179) in½ B6 ð5:72Þ m NPÀ1 s~2 kk k¼0 s~2 where kk have been arranged in decreasing order. See Table 5.2 and Fig. 5.33 about variance distribution of transform coefficients (See also Table 5.2 and Fig. 5.18 in [B6]). A.3 Normalized MSE (see Figs. 5.35À5.39 and see also Figs. 5.21À5.23 in [B6]) PP 2 vk; l k; l2 stopband J ¼ s PPNÀ1 2 vk; l k; l¼0 Js ¼ Energy in stopband/total energy, vk; l ðk; l ¼ 0; 1; ...; N À 1Þ are the trans- form coefficients of an ðÞN  N image. A.4 Rate versus Distortion (Rate-Distortion) [B6] The rate distortion function RD is the minimum average rate (bits/sample) for coding a signal at a specified distortion D (mean square error) [B6]. Let x0; x1; ...; xNÀ1 be Gaussian random variables encoded independently and ^ x^0; x^1; ...; x^NÀ1 be their reproduced values. Xk and Xk k ¼ 0, 1, ..., N À 1 are the corresponding transform coefficients. Then the average mean square distortion is 320 Appendix A Performance Comparison of Various Discrete Transforms XNÀ1 hiXNÀ1 hiÀÁ 1 2 1 ^ 2 D ¼ ExðÞn À x^n ¼ EXk À Xk N n¼0 N k¼0 For a fixed average distortion D, the rate distortion function RD is XNÀ1 1 1 s~2 R ðyÞ ¼ max ; kk Eq: ð2:118Þ in½ B6 D 0 log2 y N k¼0 2 where threshold y is determined by solving XNÀ1 ÀÁ ðyÞ¼1 y; s~2 ; fs~2 gy fs~2 g : ð : Þ ½ D min kk min kk max kk Eq 2 119 in B6 kk kk N k¼0 Select a value for Û to get a point in the plot of versus . Develop RD versus D for various discrete transforms based on I-order Markov process given N and r adjacent correlation coefficient. Plot RD vs. D, for N ¼ 8, 16, 32, ...and r ¼ 0.9, 0.95, ... For I-order Markov process (Eq. 2.68 in [B6]): ½S ¼ s2 ¼ rjjÀkj ; ¼ ; ; ...; À jk jk j k 0 1 N 1 maximum achievable coding gain is ðÞ= ½S ÀÁÀðÞÀ = ðÞ¼r 1 N tr ¼ À r2 1 1 N GN = 1 ðÞdet½S 1 N where tr ¼ trace of the matrix, det ¼ determinant of the matrix (see Appendix C in [B9]). A.5 Residual Correlation [G1] While the KLT completely decorrelates a random vector [B6], other discrete transforms fall short of this. An indication of the extent of decorrelation can be gauged by the correlation left undone by the discrete transform. This can be measured by the absolute sum of the cross covariance in the transform domain i.e., A.5 Residual Correlation [G1] 321 Data domain ⇔ Transform domain ∼ [Σ] [Σ] Set off diagonal elements to zero ∼ [Σ]ˆ [Σkk] Âà Fig. A.2 The relation between ½S and S~ XNÀ1 XNÀ1 js~2 j ij (A.4) i¼0 j¼0 i6¼j for N ¼ 8, 16, 32, ...as a function of r (Fig. A.2). Âà à Given S~ ¼ ½A ½S ½A T (A.5) ðÞNÂN ðÞNÂN ðÞNÂN ðÞNÂN Âà à Âà ~ T ~ obtain S ¼ ½A Skk ½A (A.6) ðÞNÂN ðÞNÂN ðÞNÂN ðÞNÂN Âà S~ ÂÃwhere kk is a diagonal matrix whose diagonal elements are the same as those of S~ , i.e., Âà S~ ¼ s~2 ; s~2 ; ...; s~2 kk diag 00 11 ðNÀ1ÞðÞNÀ1 It should be recognized that the conjugate appears in (A.5) that is derived in (5.42a), whereas the 2-D discrete transform of ½S is ½A ½S ½A T defined in (5.6a) and has no conjugate. Thus (A.5) can be regarded as a separable two-dimensional unitary transform of ½S for purposes of computation. Plot residual correlation versus r for DCT, DFT, KLT and ST [B23]. Fractional correlation (correlation left undone by a transform – for KLT this is zero, as KLT diagonalizes a covariance or correlation matrix) is defined as kk½ÀS ½S 2 (A.7) kk½ÀS ½l 2 P P 2 ½ ðÞ kk½2 ¼ NÀ1 NÀ1 ½ where IN is an N N unit matrix and A j¼0 k¼0 A jk . Note that theÂà measuresÂà (A.3), (A.4) and (A.7) are zeros respectively for the KLT ~ ~ as S ¼ Skk . 322 Appendix A Performance Comparison of Various Discrete Transforms x + n X N xˆ + Inverse Corrupted Transform Filtered transform signal signal [G] Filter matrix n Additive noise Fig. A.3 Scalar Wiener filtering abc Identity Hadamard defUnitary DFT DST of type 1 DCT of type 2 Karhunen−Loéve Fig. A.4 Magnitude displays of Wiener filter matrices ½G for a vector length of 16 elements ðÞN ¼ 16 . Dark pixels represent zero values, light pixels represent one values, and gray pixels represent values in between. Signal-to-noise ratio is 0.3 and r ¼ 0:9. Dynamic ranges of the Wiener filter magnitudes in the figure are compressed via the logarithmic transformation defined in (5.26) A.6 Scalar Wiener Filtering Filter matrix ½G is optimized for a specific transform, such that the noise can be filtered (Fig. A.3)[G5]. Evaluate MSE ¼ Eðjjx À ^xjj2Þ for the discrete transforms (appearing in Fig. A.4, plus Haar and slant transforms defined in [B6]) for N ¼ 4, 8, 16, 32, and r ¼ 0.9 and 0.95. A.7 Geometrical Zonal Sampling (GZS) 323 Plot magnitude displays of various discrete transforms referring to Fig. A.4. Comparing the filter planes, the filter characteristic changes drastically for differ- ent unitary transforms. For the KLT the filtering operation is a scalar multiplica- tion, while for the identity transform most elements of the filter matrix are of relatively large magnitude.
Recommended publications
  • Arxiv:2105.00793V3 [Math.NA] 14 Jun 2021 Tubal Matrices
    Tubal Matrices Liqun Qi∗ and ZiyanLuo† June 15, 2021 Abstract It was shown recently that the f-diagonal tensor in the T-SVD factorization must satisfy some special properties. Such f-diagonal tensors are called s-diagonal tensors. In this paper, we show that such a discussion can be extended to any real invertible linear transformation. We show that two Eckart-Young like theo- rems hold for a third order real tensor, under any doubly real-preserving unitary transformation. The normalized Discrete Fourier Transformation (DFT) matrix, an arbitrary orthogonal matrix, the product of the normalized DFT matrix and an arbitrary orthogonal matrix are examples of doubly real-preserving unitary transformations. We use tubal matrices as a tool for our study. We feel that the tubal matrix language makes this approach more natural. Key words. Tubal matrix, tensor, T-SVD factorization, tubal rank, B-rank, Eckart-Young like theorems AMS subject classifications. 15A69, 15A18 1 Introduction arXiv:2105.00793v3 [math.NA] 14 Jun 2021 Tensor decompositions have wide applications in engineering and data science [11]. The most popular tensor decompositions include CP decomposition and Tucker decompo- sition as well as tensor train decomposition [11, 3, 17]. The tensor-tensor product (t-product) approach, developed by Kilmer, Martin, Bra- man and others [10, 1, 9, 8], is somewhat different. They defined T-product opera- tion such that a third order tensor can be regarded as a linear operator applied on ∗Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong, China; (liqun.qi@polyu.edu.hk). †Department of Mathematics, Beijing Jiaotong University, Beijing 100044, China.
    [Show full text]
  • Session Siptool: the 'Signal and Image Processing Tool'
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by DigitalCommons@CalPoly Session SIPTool: The ‘Signal and Image Processing Tool’ An Engaging Learning Environment Fred DePiero1 Abstract The ‘Signal and Image Processing Tool’ is a With the SIPTool, students create an integrated system multimedia software environment for demonstrating and that includes their processing routine along with developing Signal & Image Processing techniques. It has image/signal acquisition and display. This integrated system been used at CalPoly for three years. A key feature is is a very different result than the ‘haphazard’ line-by-line extensibility via C/C++ programming. The tool has a processing steps that students may or may not successfully minimal learning curve, making it amenable for weekly stumble through in a MatLab environment, as they follow a student projects. The software distribution includes given example. The SIPTool-based implementation is much multimedia demonstrations ready for classroom or more like a complete, commercial product. laboratory use. SIPTool programming assignments strengthen the skills needed for life-long learning by A variety of learning objectives can be readily requiring students to translate mathematical expressions addressed with the SIPTool including: time/frequency into a standard programming language, to create an relationships, 1-D and 2-D Fourier transforms, convolution, integrated processing system (as opposed to simply using correlation, filtering, difference equations, and pole/zero canned processing routines). relationships [1]. Learning objectives associated with image processing can also be presented, such as: gray scale Index Terms Image Processing, Signal Processing, resolution, pixel resolution, histogram equalization, median Software Development Environment, Multimedia Teaching filtering and frequency-domain filtering [2].
    [Show full text]
  • Towards a Fully Automated Extraction and Interpretation of Tabular Data Using Machine Learning
    UPTEC F 19050 Examensarbete 30 hp August 2019 Towards a fully automated extraction and interpretation of tabular data using machine learning Per Hedbrant Per Hedbrant Master Thesis in Engineering Physics Department of Engineering Sciences Uppsala University Sweden Abstract Towards a fully automated extraction and interpretation of tabular data using machine learning Per Hedbrant Teknisk- naturvetenskaplig fakultet UTH-enheten Motivation A challenge for researchers at CBCS is the ability to efficiently manage the Besöksadress: different data formats that frequently are changed. Significant amount of time is Ångströmlaboratoriet Lägerhyddsvägen 1 spent on manual pre-processing, converting from one format to another. There are Hus 4, Plan 0 currently no solutions that uses pattern recognition to locate and automatically recognise data structures in a spreadsheet. Postadress: Box 536 751 21 Uppsala Problem Definition The desired solution is to build a self-learning Software as-a-Service (SaaS) for Telefon: automated recognition and loading of data stored in arbitrary formats. The aim of 018 – 471 30 03 this study is three-folded: A) Investigate if unsupervised machine learning Telefax: methods can be used to label different types of cells in spreadsheets. B) 018 – 471 30 00 Investigate if a hypothesis-generating algorithm can be used to label different types of cells in spreadsheets. C) Advise on choices of architecture and Hemsida: technologies for the SaaS solution. http://www.teknat.uu.se/student Method A pre-processing framework is built that can read and pre-process any type of spreadsheet into a feature matrix. Different datasets are read and clustered. An investigation on the usefulness of reducing the dimensionality is also done.
    [Show full text]
  • Fourier Transform, Convolution Theorem, and Linear Dynamical Systems April 28, 2016
    Mathematical Tools for Neuroscience (NEU 314) Princeton University, Spring 2016 Jonathan Pillow Lecture 23: Fourier Transform, Convolution Theorem, and Linear Dynamical Systems April 28, 2016. Discrete Fourier Transform (DFT) We will focus on the discrete Fourier transform, which applies to discretely sampled signals (i.e., vectors). Linear algebra provides a simple way to think about the Fourier transform: it is simply a change of basis, specifically a mapping from the time domain to a representation in terms of a weighted combination of sinusoids of different frequencies. The discrete Fourier transform is therefore equiv- alent to multiplying by an orthogonal (or \unitary", which is the same concept when the entries are complex-valued) matrix1. For a vector of length N, the matrix that performs the DFT (i.e., that maps it to a basis of sinusoids) is an N × N matrix. The k'th row of this matrix is given by exp(−2πikt), for k 2 [0; :::; N − 1] (where we assume indexing starts at 0 instead of 1), and t is a row vector t=0:N-1;. Recall that exp(iθ) = cos(θ) + i sin(θ), so this gives us a compact way to represent the signal with a linear superposition of sines and cosines. The first row of the DFT matrix is all ones (since exp(0) = 1), and so the first element of the DFT corresponds to the sum of the elements of the signal. It is often known as the \DC component". The next row is a complex sinusoid that completes one cycle over the length of the signal, and each subsequent row has a frequency that is an integer multiple of this \fundamental" frequency.
    [Show full text]
  • 7 Image Processing: Discrete Images
    7 Image Processing: Discrete Images In the previous chapter we explored linear, shift-invariant systems in the continuous two-dimensional domain. In practice, we deal with images that are both limited in extent and sampled at discrete points. The results developed so far have to be specialized, extended, and modified to be useful in this domain. Also, a few new aspects appear that must be treated carefully. The sampling theorem tells us under what circumstances a discrete set of samples can accurately represent a continuous image. We also learn what happens when the conditions for the application of this result are not met. This has significant implications for the design of imaging systems. Methods requiring transformation to the frequency domain have be- come popular, in part because of algorithms that permit the rapid compu- tation of the discrete Fourier transform. Care has to be taken, however, since these methods assume that the signal is periodic. We discuss how this requirement can be met and what happens when the assumption does not apply. 7.1 Finite Image Size In practice, images are always of finite size. Consider a rectangular image 7.1 Finite Image Size 145 of width W and height H. Then the integrals in the Fourier transform no longer need to be taken to infinity: H/2 W/2 F (u, v)= f(x, y)e−i(ux+vy) dx dy. −H/2 −W/2 Curiously, we do not need to know F (u, v) for all frequencies in order to reconstruct f(x, y). Knowing that f(x, y)=0for|x| >W/2 and |y| >H/2 provides a strong constraint.
    [Show full text]
  • The Discrete Fourier Transform
    Tutorial 2 - Learning about the Discrete Fourier Transform This tutorial will be about the Discrete Fourier Transform basis, or the DFT basis in short. What is a basis? If we google define `basis', we get: \the underlying support or foundation for an idea, argument, or process". In mathematics, a basis is similar. It is an underlying structure of how we look at something. It is similar to a coordinate system, where we can choose to describe a sphere in either the Cartesian system, the cylindrical system, or the spherical system. They will all describe the same thing, but in different ways. And the reason why we have different systems is because doing things in specific basis is easier than in others. For exam- ple, calculating the volume of a sphere is very hard in the Cartesian system, but easy in the spherical system. When working with discrete signals, we can treat each consecutive element of the vector of values as a consecutive measurement. This is the most obvious basis to look at a signal. Where if we have the vector [1, 2, 3, 4], then at time 0 the value was 1, at the next sampling time the value was 2, and so on, giving us a ramp signal. This is called a time series vector. However, there are also other basis for representing discrete signals, and one of the most useful of these is to use the DFT of the original vector, and to express our data not by the individual values of the data, but by the summation of different frequencies of sinusoids, which make up the data.
    [Show full text]
  • Quantum Fourier Transform Revisited
    Quantum Fourier Transform Revisited Daan Camps1,∗, Roel Van Beeumen1, Chao Yang1, 1Computational Research Division, Lawrence Berkeley National Laboratory, CA, United States Abstract The fast Fourier transform (FFT) is one of the most successful numerical algorithms of the 20th century and has found numerous applications in many branches of computational science and engineering. The FFT algorithm can be derived from a particular matrix decomposition of the discrete Fourier transform (DFT) matrix. In this paper, we show that the quantum Fourier transform (QFT) can be derived by further decomposing the diagonal factors of the FFT matrix decomposition into products of matrices with Kronecker product structure. We analyze the implication of this Kronecker product structure on the discrete Fourier transform of rank-1 tensors on a classical computer. We also explain why such a structure can take advantage of an important quantum computer feature that enables the QFT algorithm to attain an exponential speedup on a quantum computer over the FFT algorithm on a classical computer. Further, the connection between the matrix decomposition of the DFT matrix and a quantum circuit is made. We also discuss a natural extension of a radix-2 QFT decomposition to a radix-d QFT decomposition. No prior knowledge of quantum computing is required to understand what is presented in this paper. Yet, we believe this paper may help readers to gain some rudimentary understanding of the nature of quantum computing from a matrix computation point of view. 1 Introduction The fast Fourier transform (FFT) [3] is a widely celebrated algorithmic innovation of the 20th century [19]. The algorithm allows us to perform a discrete Fourier transform (DFT) of a vector of size N in (N log N) O operations.
    [Show full text]
  • Circulant Matrix Constructed by the Elements of One of the Signals and a Vector Constructed by the Elements of the Other Signal
    Digital Image Processing Filtering in the Frequency Domain (Circulant Matrices and Convolution) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering 2 Toeplitz matrices • Elements with constant value along the main diagonal and sub-diagonals. • For a NxN matrix, its elements are determined by a (2N-1)-length sequence tn | (N 1) n N 1 T(,)m n t mn t0 t 1 t 2 t(N 1) t t t 1 0 1 T tt22 t1 t t t t N 1 2 1 0 NN C. Nikou – Digital Image Processing (E12) 3 Toeplitz matrices (cont.) • Each row (column) is generated by a shift of the previous row (column). − The last element disappears. − A new element appears. T(,)m n t mn t0 t 1 t 2 t(N 1) t t t 1 0 1 T tt22 t1 t t t t N 1 2 1 0 NN C. Nikou – Digital Image Processing (E12) 4 Circulant matrices • Elements with constant value along the main diagonal and sub-diagonals. • For a NxN matrix, its elements are determined by a N-length sequence cn | 01nN C(,)m n c(m n )mod N c0 cNN 1 c 2 c1 c c c 1 01N C c21 c c02 cN cN 1 c c c c NN1 21 0 NN C. Nikou – Digital Image Processing (E12) 5 Circulant matrices (cont.) • Special case of a Toeplitz matrix. • Each row (column) is generated by a circular shift (modulo N) of the previous row (column). C(,)m n c(m n )mod N c0 cNN 1 c 2 c1 c c c 1 01N C c21 c c02 cN cN 1 c c c c NN1 21 0 NN C.
    [Show full text]
  • Arxiv:1801.05832V2 [Cs.DS]
    Efficient Computation of the 8-point DCT via Summation by Parts D. F. G. Coelho∗ R. J. Cintra† V. S. Dimitrov‡ Abstract This paper introduces a new fast algorithm for the 8-point discrete cosine transform (DCT) based on the summation-by-parts formula. The proposed method converts the DCT matrix into an alternative transformation matrix that can be decomposed into sparse matrices of low multiplicative complexity. The method is capable of scaled and exact DCT computation and its associated fast algorithm achieves the theoretical minimal multiplicative complexity for the 8-point DCT. Depending on the nature of the input signal simplifications can be introduced and the overall complexity of the proposed algorithm can be further reduced. Several types of input signal are analyzed: arbitrary, null mean, accumulated, and null mean/accumulated signal. The proposed tool has potential application in harmonic detection, image enhancement, and feature extraction, where input signal DC level is discarded and/or the signal is required to be integrated. Keywords DCT, Fast Algorithms, Image Processing 1 Introduction Discrete transforms play a central role in signal processing. Noteworthy methods include trigonometric transforms— such as the discrete Fourier transform (DFT) [1], discrete Hartley transform (DHT) [1], discrete cosine transform (DCT) [2], and discrete sine transform (DST) [2]—as well as the Haar and Walsh-Hadamard transforms [3]. Among these methods, the DCT has been applied in several practical contexts: noise reduction [4], watermarking methods [5], image/video compression techniques [2], and harmonic detection [2], to cite a few. In fact, when processing signals modeled as a stationary Markov-1 type random process, the DCT behaves as the asymptotic case of the optimal Karhunen–Lo`eve transform in terms of data decorrelation [2].
    [Show full text]
  • + Wireless World
    ELECTRONICS Denmark DKr. 70.00 Germany DM 15.00 Greece Dra.760 Holland DFI. 14 Italy L. 7300 IR £3.30 Spain Pts. 780 WORLD Singapore SS 12.60 USA $6.70 + WIRELESS WORLD SOR DISTRIBUTION A REED BUSINESS PUBLICATION MARCH 1993 £1.95 APPLICATIONS Intelligent charging for f 4A 4 r NiCds B 7-1/8 RF DESIGN ATV receiver for 23cm SYSTEMS Radio circuitry for GPS systems SIGNAL PROCESSING Experimenting with DSP REVIEW Dadisp-the clumsy alternative ? IMF an Hickman's new hook READE na oguee for just £22.40* *Overseas f27.40..This month Normal price L 34.90 19.91.1) OFFER I 833004 THE WORLDS No.1 BEST SELLING UNIVERSAL PROGRAMMING AND TESTING SYSTEM. The PC82 Universal Programmer and Tester is a *More sold worldwide than any other of *Software supplied to write own test PC -based development tool designed to its type. vectors for custom ICs and ASICs etc. program and test more than 1500 ICs. The latest *UK users include BT, IBM, MOD, THORN *Protection circuitry to protect against version of the PC82 is based on the experience EMI, MOTOROLA, SANYO, RACAL wrong insertion of devices. gained after a 7 year production run of over *High quality Textool or Yamaichi zero 100,000 units. insertion force sockets. *Ground control circuitry using relay switching. *Rugged screened cabling. The PC82 is the US version of the Sunshine *High speed PC interface card designed *One model covers the widest range of Expro 60, and therefore can be offered at a very devices, at the lowest cost.
    [Show full text]
  • Comparison of FFT, DCT, DWT, WHT Compression Techniques on Electrocardiogram & Photoplethysmography Signals
    Special Issue of International Journal of Computer Applications (0975 – 8887) International Conference on Computing, Communication and Sensor Network (CCSN) 2012 Comparison of FFT, DCT, DWT, WHT Compression Techniques on Electrocardiogram & Photoplethysmography Signals Anamitra Bardhan Roy Debasmita Dey Devmalya Banerjee Programmer Analyst Trainee, B.Tech. Student, Dept. of CSE, JIS Asst. Professor, Dept. of EIE Cognizant Technology Solution, College of Engineering College of Engineering, WB, India Kolkata, WB, India Kolkata, WB, India, Bidisha Mohanty B.Tech. Student, Dept. of ECE, JIS College of Engineering, WB, India. ABSTRACT media requires high security and authentication [1]. This Compression technique plays an important role indiagnosis, concept of telemedicine basically and foremost needs image prognosis and analysis of ischemic heart diseases.It is also compression as a vital part of image processing techniques to preferable for its fast data sending capability in the field of manage with the data capacity to transfer within a short and telemedicine.Various techniques have been proposed over stipulated time interval. The therapeutic remedies are derived years for compression. Among those Discrete Cosine from the minutiae’s of the biomedical signals shared Transformation(DCT), Discrete Wavelet worldwide through telelinking. This communication of Transformation(DWT), Fast Fourier Transformation(FFT) exchanging medical data through wireless media requires high andWalsh Hadamard Transformation(WHT) are mostly level of compression techniques to cope up with the restricted used.In this paper a comparative study of FFT, DCT, DWT, data transfer capacity and rate. Compression of signals and and WHT is proposed using ECG and PPG signal, images can cause distortion in them. As the medical signals whichshows a certain relation between them as discussed in and images contains and relays information required for previous papers[1].
    [Show full text]
  • Dadisp / Filters Digital Filter Design Module
    DADiSP / Filters Digital Filter Design Module KEY FEATURES DADiSP/Filters is a menu-driven module for FIR and IIR digital filter design and analysis. From Simple User Interface any DADiSP Worksheet, you can quickly design, view and analyze both FIR (Finite Impulse Lowpass, Highpass, Bandpass, Bandstop and Response) and IIR (Infinite Impulse Response) Multiband Filters filters. Once you have designed a filter, you may filter the raw data, then view and analyze the Finite Impulse Response (FIR) Filter Design filtered signal. Through the easy-to-use dialog FIR Hilbert Transforms and Differentiators boxes or simple one line functions, you can tune the filter iteratively, re-filtering the data until you FIR Remez Exchange and Kaiser Window have separated the signal from the noise cleanly. Design Algorithms Infinte Impulse Response (IIR) Filter Design IIR Bessel, Butterworth, Chebychev I, Chebychev II and Elliptic Filters IIR Bilinear Transform and Matched Z Design Algorithms Magnitude, Phase, Group Delay and Impulse Response Output Coefficent Form Conversion Output Coefficent Quantization Time and Frequency Domain Filtering Functions Pole-Zero Plots www.dadisp.com DSP Development Corporation Toll Free: 1-800-424-3131 New Features DADiSP/Filters Version 5.0 includes a completely redesigned user interface to streamline the process of designing and applying digital filters. Straightforward dialog boxes with automatic option validation simplifies both the design and analysis of filters. Filter coefficients can be easily converted to various filter structures and quantization routines are included to help simulate DSP chipsets. IIR Bessel filters and the Matched Z Transform design method have been added. Linear phase FIR Kaiser filters have been expanded and enhanced.
    [Show full text]