Ideal Filter - Windowed - DTFT/Highpass

  • Thread starter Thread starter DSRadin
  • Start date Start date
  • Tags Tags
    Filter
Click For Summary
SUMMARY

The discussion focuses on the analysis of the discrete-time Fourier transform (DTFT) of a windowed function derived from the high-pass differentiator filter H_{dd}(e^{j\omega})=j\omega e^{\frac{-j\omega}{2}}. The user successfully synthesized h_{dd}(n) as -sin(π(n-1/2)) and derived the windowed DTFT result H_3(e^{j\omega}) using integration by parts. A critical error in MATLAB code was identified, where the loop index was incorrectly set, leading to constant output. The correct loop should iterate from 1 to length(n), resolving the issue and confirming the expected behavior of the filter as N approaches infinity.

PREREQUISITES
  • Understanding of discrete-time Fourier transform (DTFT)
  • Familiarity with MATLAB programming and debugging
  • Knowledge of window functions in signal processing
  • Basic concepts of high-pass filters and their characteristics
NEXT STEPS
  • Study the properties of window functions in signal processing
  • Learn about MATLAB for signal processing applications
  • Explore the derivation of high-pass differentiators in DTFT
  • Investigate the implications of increasing N in filter design
USEFUL FOR

Signal processing engineers, MATLAB users, and students studying digital filters and Fourier analysis will benefit from this discussion.

DSRadin
Messages
12
Reaction score
1

Homework Statement


Given:
H_{dd}\left(e^{j\omega}\right)=j\omega e^{\frac{-j\omega}{2}}, \left|\omega\right|\le\pi

Find: H_{3}\left(e^{j\omega}\right) where
H_{3}\left(e^{j\omega}\right) is the spectrum of h_{dd}\left(n\right)\left(W_N\left(n\right)\right) and W_N\left(n\right)=1 for \frac{-N}{2}\le n \le \frac{N}{2} , 0 else


Homework Equations


DTFT Synthesis: \frac{1}{2\pi}\int_{-\pi}^{\pi} H\left(e^{j\omega}\right)e^{j\omega n}d\omega
DTFT Analysis: \sum_{n=\frac{-N}{2}}^{\frac{N}{2}} h(n)e^{-j\omega n}

The Attempt at a Solution



Step 1: Synthesis h_{dd}\left(n\right). This is done through integration by parts and my result is:

h_{dd}(n)=-sin\left(\pi\left(n-\frac{1}{2}\right)\right)

Step 2: Window - ok. -\frac{N}{2}\le n \le \frac{N}{2} is the new range.

Step 3: DTFT windowed function result:
H_3\left(e^{j\omega}\right) = -\sum_{n=-N/2}^{N/2} \frac{sin\left(\pi n - \frac{\pi}{2}\right)}{\pi \left(n-1/2\right)^2}e^{-j\omega n}

Really cool - but when plotted versus frequency, I get a constant, regardless of the size of N.
The goal of this problem was to prove that as N--> big that the filter approaches the ideal high-pass differentiator. I must have made a mistake somewhere but I'm not sure where, if anyone could see if they receive a different result I would be much obliged. Thanks!

-DR
 
Physics news on Phys.org
Found my mistake - it turns out that the above is actually correct and corresponds to Sum( (-1)^(n-1/2)/(denom) * e^-jwn).

There was an error in my MATLAB code (ridiculous error) where my for loop looked like:

for i=length(n)

instead of

for i=1:length(n)

hence the reason I was only getting one constant value... You think it's some important mistake in your math and it turns out to be a typo.

Oh well, Go Bears.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
Replies
19
Views
3K
Replies
3
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K
Replies
24
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K