DSPRelated.com

David Shaw (@dgshaw6)

I am a Senior level DSP software and algorithms engineer. I have many, many years of DSP product development experience applied to world class products, including modems, speakerphones, network and acoustic echo cancellation for end user products and carrier grade infrastructure. I have also been involved in wireless communications techniques including some SDR analysis of unknown signals. My current work is in the analysis of LTE radio signals to derive the time delays to allow for triangulation to determine location.

There is a fairly complete C source library from ITU that has lots of basic functions that deal with saturation and overflow protection.G.191

Re: fft process gain revisited

Reply posted 1 year ago (03/10/2023)
There are many interesting answers in the thread.I have believed for a long time that the effective quantization noise in individual bins after an FFT is > the...

Re: Comparing Utterances

Reply posted 2 years ago (08/29/2022)
Is there any implicit assumption that the utterances came from the same person?If so, then some of what you have suggested might work.  You will probably need some...

Re: #OFDM material

Reply posted 2 years ago (07/25/2022)
One of the best books I have found recently, is Digital Signal Processing in Modern Communications Systems, by Andreas Schwarzinger.It is self published and comes...

Re: #OFDM material

Reply posted 2 years ago (07/25/2022)
Have you tried looking at the code for GNU Radio?

Re: #OFDM sample rate

Reply posted 2 years ago (07/18/2022)
The earlier DJI drones used proprietary modulation, but the newer ones use an LTE chipset.There is a company that you may want to look at, that is doing what you...

Re: 'near IQ' sampling

Reply posted 2 years ago (07/18/2022)
The term I have heard in the past is "quadrature sampling".

Measuring RSRP in 5G

New thread started 2 years ago
Hi all,In the 3GPP standard 38.215 Section 5.1.19, there is a reference to the definition of the RSRP for SRS.It appears to simply be the average of the power in...

Re: How to separate these orthogonal inputs

Reply posted 2 years ago (06/08/2022)
I think is it is as simple as performing the IFFT on the combined received signal, selecting only the impulse response for the UE of interest, by zeroing out the...

Re: How to separate these orthogonal inputs

Reply posted 2 years ago (06/08/2022)
Hi Kaz,The cyclic shift (as defined actually) amounts to a time delay.If you go to the time domain, and then window the specific "impulse" of interest, you can separate...

Re: correlation of two complex signals

Reply posted 2 years ago (04/18/2022)
Hi Kaz,I've been working a lot with the Zadoff Chu sequences for 5G lately.For my case, I am able to retrieve the frequency domain bins from the receiver where the...
Very cool, and I bet he is writing backwards on a large transparent panel!The reason why, is that he is probably right handed, and I watched him erase some things.Also,...
Thanks to all for the kind contributions.It turns out that I found an extremely promising scheme (since the realistic number of modes is low) that is referred to...
Thanks Chuck.We are using Matlab for the "offline" processing right now, and it does the job we need.In fact, we are able to run in "real-time" for a significantly...
Thanks for your idea.I will have to see if there is some promise there.
Thanks for your insight.For our Matlab implementation, we have used svds() which allow us to specify that the processing should find only the first N singular values.For...
I have a matrix that is complex but truly symmetric (not conjugate symmetric).I'm trying to find a reference that will help me build an extremely efficient implementation...

Re: Complex IF

Reply posted 2 years ago (09/17/2021)
...

Re: Complex IF

Reply posted 3 years ago (09/17/2021)
Many excellent answers here, but 1 item I think has not been mentioned.The word reflection is used in several of them, and that is the essence of my comment.An exchange...
If you only care about signals up to 1kHz, you can filter and decimate the sample data down to a lower rate, and then run 2048 point FFTs on the lower sampling rate...

Re: DTFT optimization question

Reply posted 3 years ago (06/03/2021)
Again, I agree Tim.I view Goertzel as a "sympathetic" driven oscillator.The poles are on the unit circle (zero damping as you said) as they would be if you chose...

Re: DTFT optimization question

Reply posted 3 years ago (06/03/2021)
I just looked up the details of the computing platform described (Jetson Xavier AGX) in the original problem statement, and it appears to have several floating...

Re: DTFT optimization question

Reply posted 3 years ago (06/03/2021)
Thanks Tim,I totally agree, and you still have to run the correct number of samples through it to get to the right result.Just thought it might be an option for...

Re: DTFT optimization question

Reply posted 3 years ago (06/03/2021)
A Goertzel filter is effectively a 1 bin FT.Have you looked into it for what you need?It is effectively implemented using a biquad section for each desired frequency.All...

Re: DSP principles question

Reply posted 3 years ago (04/27/2021)
Hello Rich,The concept in the book is related to requirements for sampling relative to the content of the signal.If you are working in audio, then the hardware will...

Re: Impulse response in OFDM system using IFFT

Reply posted 3 years ago (04/20/2021)
I confess that you have confused me a little.You started with some recorded data that had a certain sample rate.You performed an FFT of that data at some specific...

Re: Impulse response in OFDM system using IFFT

Reply posted 3 years ago (04/20/2021)
The plot you show, of the time domain response should also have an fftshift applied to it to bring the impulse into the middle where the new 0 delay will be.The...

Re: Loop optimization in C6416T Processor

Reply posted 3 years ago (04/13/2021)
We don't know the format of your data, but it seems that the name "knownsym" implies data communications, so I would guess that fixed point math should be the most...
In the work that I was doing, the problem showed up in very narrow bandpass filters, which are almost always either a single biquad or combination of biquads.I have...
I believe that the size of the internal variables inside the fixed point IIR filters is actually the harder problem to deal with.Most IIR filters are > order...

Re: Autocorrelation of OFDM Signal ?..

Reply posted 3 years ago (12/20/2020)
One source of other correlation components is pilot signals.Almost every OFDM system has pilot symbols to assist with certain timing acquisition.Depending on your...

Re: LTE delay measurement resolution?

Reply posted 3 years ago (12/16/2020)
OK.  I stand corrected.  Thanks.

Re: LTE delay measurement resolution?

Reply posted 3 years ago (12/16/2020)
Hi Kaz,I think you may have your IFFT length wrong.1 Ts is 1/30.72e6, and the 20 MHz band is assumed to have a symbol length of 2048 samples, and so is the length...

Re: short time subband coding

Reply posted 3 years ago (10/21/2020)
The implementation details of G.722 may give you some insight in terms of split band adpcm coding.

Re: short time subband coding

Reply posted 3 years ago (10/20/2020)
It is not easy to see from your plots, but there are clearly some repeating stripes in the lower spectral plot.16e3 * 0.003 gives 48 samples per frame, but you...

Re: IFFT Synthesis Amplitude Estimation

Reply posted 4 years ago (08/13/2020)
This problem is studied extensively in the context of OFDM modulation in data communications, where the power amplifiers are non-linear as you push closer to full...

Re: IFFT Synthesis Amplitude Estimation

Reply posted 4 years ago (08/13/2020)
Which is a train of impulses.

Re: understanding no perfect Zero-IF

Reply posted 4 years ago (07/25/2020)
I'm going to try to understand what you "meant", and comment based on those assumptions.First, in my mind, zero-IF implies that correct downshifting would give...
I think the "driven" oscillator approach will work well.Do you have access to a copy of the power-line noise independently?That could be used as a reference for...

Re: Adaptive filter and noise reduction

Reply posted 4 years ago (07/23/2020)
I saw you additional question in another post, so I understand you problem a little better now.I incorrectly assume that you would be working on audio, sampled at...
I think that your best approach is in fact an adaptive filter, which I now see you posted a question about in a different posting.It is not that hard to remove a...

Re: Adaptive filter and noise reduction

Reply posted 4 years ago (07/23/2020)
Hi sh_ee,I'm not sure that the information that you have provided is enough to answer the question.Lets assume that you are talking audio, and that you sample rate...
Hi napierm,I went through almost the exact same idea several years ago for  measuring signal power in an audio application.The only thing I would mention that relates...

Re: Question regarding CIC filter

Reply posted 4 years ago (05/28/2020)
As you suggest Rick, it certainly looks like a 4 AM symbol sequence.I'm guessing that the sample rate of the captured data is 53 MHz, but that the symbol rate is...
I'm confused by your results too.10^-4 at baseband (presumably without fading) and 10^-5 at passband (with fading) is the opposite of what we would expect.Did I...

Re: Computing the delay from a room response FIR filter

Reply posted 4 years ago (05/06/2020)
Aha! Nice work.My bet is that the 500 or so samples of delay are dominated by the buffering.  If we guess 480 samples, that is 10 msec at 48k sampling.My first...

Re: Filtering a sine wave with ALE

Reply posted 4 years ago (05/06/2020)
...

Re: Computing the delay from a room response FIR filter

Reply posted 4 years ago (05/05/2020)
Hi planb_buzz,If your measured delay is only 6 samples, and you are looking at a room impulse response of, probably, hundreds of samples, then I don't understand...

Re: Filtering a sine wave with ALE

Reply posted 4 years ago (05/04/2020)
Hi Artur,If the interference is a pure sine wave, even fluctuating slightly, then consider a driven oscillator, in which you have a second order IIR filter that...
Should not be any different if you are starting with time domain samples.  If you up-convert the baseband to 800 MHz (interpolated by 20) for the real and imaginary...
Hello Amit,Sorry to be too basic if this is such.Just a quick question as to whether you are carefully choosing a sampling rate that can handle the real only signal...

Re: Setting 25KHz playback sampling rate in sound card

Reply posted 4 years ago (03/27/2020)
Hi Rick,Farrow's original paper.All the plots in the paper were generated by Bill using HTBasic, which was his favorite platform.00015483.pdfOne of the references...

Re: Setting 25KHz playback sampling rate in sound card

Reply posted 4 years ago (03/26/2020)
Hi Rick,I started my career at Bell Labs in 1981 and worked for several AT&T elements and derivatives including AT&T Information Systems, AT&T Microelectronic, Lucent,...

Re: Setting 25KHz playback sampling rate in sound card

Reply posted 4 years ago (03/25/2020)
You can use a Farrow interpolator to resample from 25 kHz to whatever rate is supported by your sound card. (22.050, 24, 44.1, 48 ksps etc) However, if you are...
Very good overview Rick,@amitjonak Very often, the design of the transmit and received filter pair are designed as root-raised-cosine pulse-shaping filters.The TX...

Re: How to obtain impulse response of a room acoustics

Reply posted 4 years ago (01/22/2020)
Many, many years back I worked with Bill Farrow of Farrow interpolator fame.  He worked out an incredibly clever scheme to get an LMS filter to converge extremely...

Re: How to obtain impulse response of a room acoustics

Reply posted 4 years ago (01/22/2020)
For humor's sake.In the EE department at Columbia University, there was a great black and white picture, in the hallway, of a faculty member and some students on...

Re: ADC sampling rate matlab/simulink M-PSK

Reply posted 4 years ago (12/24/2019)
Is there a compelling reason why you want to prove the theorem with M-PSK signal?It is far easier to prove the theorem using a simple sine wave that is close to...

Re: Practical digital communications book

Reply posted 4 years ago (11/13/2019)
One of the best books I have found recently, that is very practical, but not down to the FPGA level, is Digital Signal Processing in Modern Communications Systems,...

Re: phase change

Reply posted 4 years ago (11/09/2019)
Hi Sharan,Several good suggestions here so far from drmike, ChuckMcM and Neil.However, your question is vague, so it is not clear what you are trying to achieve.Some...

Re: concept of bandwidth

Reply posted 4 years ago (10/23/2019)
Thanks Rick.I agree completely.  The bandwidth of the signal is not changed based on the number of bits per sample.  Only the signal to quantization noise ratio...

Re: concept of bandwidth

Reply posted 4 years ago (10/22/2019)
Is your question related to the comments I left about bandwidth in the OFDM discussion the other day?I may be wrong in my understanding of your question here, but...

Re: OFDM questions

Reply posted 4 years ago (10/18/2019)
As suggested by SlartibartFast, the usage of the subcarriers and symbols within a frame are a shared resource between the user equipment (UE) and the base station...

Re: OFDM questions

Reply posted 4 years ago (10/18/2019)
The parameters you defined sound exactly like the 10 MHz LTE channel commonly used for cellular telephony.1) It does matter which subcarriers are used.  So, for...
20 ppm is even better then, and I re-enforce my previous statement.  If you think about it, then it takes 50,000 samples to have a one sample shift in timing. I...
A questions related to your estimation process.What is the worst case difference between the TX and the RX clock in PPM?Remember, that you have an uncertainty problem...
First, I assume in my technique that the stretch/shrink factor is consistent across the entire capture.  If this is true, then no windowing is necessary, and the...
The most typical application for the sample rate conversion technique has been in modem signal analysis/reception. This is often referred to as timing recovery and...
Thanks for your kind support!
I'm not sure that you care necessarily to make the signals the same length.  I think that the technique I propose does not care about the relative length of the...
I have had to do this several times.  I have used a Farrow interpolator to do fractional rate sample rate conversion on only one of the signals (presumably the...
I have to confess that I am a little confused by the question.You mention random hopping over a range of F0 + [500:300:5000]The you mention another type of hopping...

Re: Frequency estimation in between the bins ...

Reply posted 5 years ago (07/19/2019)
Can you give us a little more information?What is the size if the FFT?  How much zero padding did you have to do, since the % of zero padding defines the width...

Re: Clock drift and compensation

Reply posted 5 years ago (06/28/2019)
Fun problem.Many, many years ago, we had this same problem when synchronizing the internals of the first Sirius Satellite Radio receivers to the audio play-out.We...

Re: Timing Recovery Loop Filter Rate

Reply posted 5 years ago (05/11/2019)
Hi Rick,Thanks for your confidence in assuming I know something.The simple answer is yes.  The signal sees an analog low pass filter. The process is a little more...

Re: Timing Recovery Loop Filter Rate

Reply posted 5 years ago (05/10/2019)
Many years ago, I worked on the development of the modem for the AT&T VideoPhone 2500.We had an oversampling rate that allowed for the equivalent of 3 sampling...

Re: Todo List: Improvements to the Related Sites

Reply posted 5 years ago (05/02/2019)
I agree with lamabrew.However, I do find that even on my phone, the current format is readable if I turn my phone to landscape mode.  Unfortunately in this mode,...

Re: SNR of a pulsed signal in noise

Reply posted 5 years ago (04/27/2019)
I have read the responses from Dr. Mike and Slartibartfast, and their ideas are completely valid, and should work.Here are some other ideas.How much do you know...

Re: Time delay

Reply posted 5 years ago (04/04/2019)
I agree with the general principles espoused in the other comments here, but I think that things can be made less complicated.However, there is an important question we...
I agree with the single polyphase structure you suggests, but I would implement it so that you get 1 new output for every 3 or 4 inputs.So you have a delay line...

Re: fixing sample rate error/mismatch

Reply posted 5 years ago (01/15/2019)
fred harris' response is spot on, as stated.The control signal he mentions can be run into a PLL to smooth the process significantly, and the phase difference derived...

Re: Folded FIR filter in Matlab or Octave

Reply posted 5 years ago (01/15/2019)
In addition to the responses already given, I believe that MatLab is likely to have a very optimized compiled routine that is called for filter() or conv(), that...

Re: Matlab-implementation of dynamic range compressor

Reply posted 5 years ago (01/14/2019)
This idea, as a general concept, is the essence of mu-law or a-law coding.https://www.itu.int/rec/T-REC-G.711-198811-I/en
One more piece of advice, that no-one has mentioned here.You will have much more success with this whole thing if you go to sub-band processing using a polyphase...
I totally agree with achesir (hello by the way).However, I would also caution, that the recording into the mic will probably be in a space that has some significant...

Re: OFDM Cyclic Prefix Understanding

Reply posted 5 years ago (12/06/2018)
I understand exactly what you are saying, and you make a completely valid point.Thanks.

Re: OFDM Cyclic Prefix Understanding

Reply posted 5 years ago (12/06/2018)
I disagree slightly here.The dispersion of the channel is alway causal, so that the ISI from one symbol to the next (if there is any) will always creep into the...

Re: OFDM Cyclic Prefix Understanding

Reply posted 5 years ago (12/05/2018)
That makes sense now.However, the code looks wrong now.  What you really need is the last part of the FFT segment rather than the earlier part of the cyclic prefix.You...

Re: OFDM Cyclic Prefix Understanding

Reply posted 5 years ago (12/05/2018)
Not sure I quite understand the situation.If symOffset is set to 0.75, how can it be used as an index into an array in the line: postCPRemoval = x([CPLen+1:FFTLen+symOffset,...
Thanks for the correction to an obvious error in my first response. I have attempted to make the correction above.I took a trip down a deep rabbit hole as I was...
Fairly simple answer for waveforms with relatively clean tones.Simply demodulate successive samples against each other:y = x(n)*conj(x(n-1));If you normalize the...

Re: Magnitude of "wav" files

Reply posted 5 years ago (11/27/2018)
The way in which MatLab views wave files is simply to scale them so that they fall in the range +- 1.0.If you had access to the original samples, they would fall...
You are using a 1st order farrow, which is nothing more than a linear interpolation.From your diagram, I see a link from a scaled (by the left interpolator coefficient)...
The 3gpp standard for wide-band speech coders, has an example C code reference that includes the full sub-band separation for the noise reduction algorithm.  This...
Unfortunately, am not familiar with MCLT, so can't comment on their use here.Thanks for the reference though.
For cancellation, and noise reduction, the most common frequency domain techniques use FFTs, but usually as part of a sub-band processing scheme using filter-banks. ...

Re: Singularity issue due to Zero-frequency component

Reply posted 5 years ago (10/29/2018)
Sorry to say I don't have an answer, but rather a question to consider.From a practical point of view, can your forcing function, applied to your non-linear elastic...
There are so many places where this synchronization problem occurs, that it is in fact a very common problem.  Some applications places identified below: Modems,...
In my experience with telecom, where we often use similar filters for telephone bandwidth speech, we typically use a set of elliptic IIR filters.  One HP and one...

Re: MCU signal processing PDM to PCM

Reply posted 5 years ago (10/03/2018)
I assume that you are decimating after the filtering operation.  In order to bring the signal back up to a reasonable level, you probably need to scale up by your...

Re: MCU signal processing PDM to PCM

Reply posted 5 years ago (10/02/2018)
Replace the data with signed values of -1 and +1 instead of 0 and 1, by converting your uint8_t to int8_t, multiply by 2 and subtract 1.No need for high pass.

Re: OFDM Constellation.

Reply posted 5 years ago (10/01/2018)
Plot your constellations as a 3D plot.plot3(symbol1:symbolN,real(symbol1:symbolN),imag(symbol1:symbolN))I even go so far as to plot two to form a stereoscopy image.---------figure(1)subplot(1,2,1)plot3(symbol1:symbolN,real(symbol1:symbolN),imag(symbol1:symbolN))camproj('perspective')view(50,10)grid...

Re: Machine learning and detection in noise

Reply posted 5 years ago (09/23/2018)
Before you throw out the conventional techniques, it seems appropriate to ask a few questions.1) Are you trying to just determine the existence of certain known...

Re: How to cancel REALLY heavy echo?

Reply posted 5 years ago (09/20/2018)
The sub-band separation requires polyphase filtering and a DFT.In our case, we were processing 8ksps data for 3.4 kHz audio, and we used a very odd number selection. ...

Re: How to cancel REALLY heavy echo?

Reply posted 5 years ago (09/20/2018)
In the field of acoustics, there is a concept of T60 time.  This is how long an impulse takes to decay by 60 dB as a result of the reverberation within an acoustic...

Re: Hot DSP Topics 2018

Reply posted 5 years ago (09/19/2018)
Fun to hear about your experience with noise cancellation.  I worked for more than 30 years on adaptive filtering techniques for audio and data communications cancellation...
My understanding is that the design of the IIR filter you want will define the delay.  A parallel or cascade implementation in the form of biquads are just different...
An FIR filter with minimum phase response can - at best - match the delay of a set of biquads.I believe that you won't be able to do any better WRT delay.A typical...
Looks like I may have missed the essence of your issue.I re-read you comment and see that you are trying to align your RMS estimate with the signal itself.If this...
For decimation, you will have to have some anti-aliasing filtering.Decimation will add to delay, but it is un-avoidable.You have expressed a desire to measure RMS...
Can you decimate a couple of times first, and then do the Biquad?Consider a cascade of half-band filters to get the sampling rate down first.Also, if you are measuring...
I'm not sure exactly where to jump in here.I have read the discussion with @Slartibartfast, and I'm not as pessimistic as he is about what you can achieve.If you...
Hehe.  Been there and done that!!No worries :-)
You make it sound simple, but a complex multiply by what?How much has to be done to figure out the values of the complex coefficients for the multiplies?  Table...
The Goertzel is very efficient.It is effectively the left side of a biquad with only one coefficient requiring a multiply and a couple of added, and two data element...
Your are most welcome.  I have built these detectors more times than I care to mention.Goertzel detectors require some experimentation to get them right.  Depends...
The detection of the "tring tring" is a standard call progress detection problem.There are three tones to look for: 350, 440, 480 and 620 Hz.Each call progress is...

Re: IIR filters

Reply posted 6 years ago (11/29/2017)
In the design of anti-aliasing filters, we have very often used IIR filters.  They are almost alway elliptic designs.  They are also most often higher order than...

Re: Quick question on Deconvolution

Reply posted 6 years ago (11/28/2017)
Hi Steve and nelsona,I have responded to some of nelsona's previous postings - about his thesis - a long time back.One thing that has not been mentioned in any of...

Re: DSBSC Demodulator

Reply posted 6 years ago (11/02/2017)
Just a little less than half of the 45 msec?That is pretty good.  I suspect that atan2 will bring you a lot of savings too.

Re: DSBSC Demodulator

Reply posted 6 years ago (11/02/2017)
Hi,What percentage of your total real time is 6000 cycles?I don't have a tutorial for you, I'm afraid.  Most of the hard learn comes from experience.Overflows are...

Re: DSBSC Demodulator

Reply posted 6 years ago (10/31/2017)
Hi again.Again angle takes on values between 0x0000 and 0x7fff.0x7fff represents 2*pi - epsilon.Remember that angle "value" = "value" & 0x7f80 + "value" &...

Re: DSBSC Demodulator

Reply posted 6 years ago (10/31/2017)
First, of course you are right that 256/4 is 64.  My mistake.Next.  The table entries are not calculated from the index, but rather from a fraction of a cycle.table[0x00]...

Re: DSBSC Demodulator

Reply posted 6 years ago (10/31/2017)
My advice, is not to try to convert from float to fixed.  Rather, start in the fixed point domain.For the sin and cos, start with the identity sin(a+b) = sin(a)*cos(b)+cos(a)*sin(b)....

Re: DSBSC Demodulator

Reply posted 6 years ago (10/31/2017)
Wow! Your code is ripe for some clever stuff.Don't be scared of fixed point math.Sin and cos in 4 multiplies and 2 adds etc.Keep angle as binary between 0 and 0x7FFF...

Re: DSBSC Demodulator

Reply posted 6 years ago (10/30/2017)
Have you run a profile of your code to see where the major real time is consumed?Can you "compromise" the algorithms to reduce their real time and possibly their...

Re: Reverberation rejection filter

Reply posted 6 years ago (10/02/2017)
I have to ask the question.  How can the first 10 chunks be reverberation only?  Almost by definition, reverberation is delayed compared to the original signal.Am...
Y. Ephraim and D. Malah, "Speech Enhancement Using a Minimum Mean-Square Error Short-Time Spectral Amplitude Estimator", IEEE Trans. on Acoustics, Speech, and Signal...
In order to create a signal that is shifted by 90 degrees, the Hilbert filter is an interpolator. This virtually requires some sort of low pass filtering.  That...

Re:

Reply posted 7 years ago (09/13/2017)
Hi NatiyalTD,Best way to show you is to dump my code on you.Code to get the plots I used.

Re:

Reply posted 7 years ago (09/12/2017)
@NatiyalTD,I have some specific results.I went about the process in the following order. Zoomed in on one of the longest bursts of tone. Used just the cleanest part...

Re:

Reply posted 7 years ago (09/11/2017)
Wav is great.

Re:

Reply posted 7 years ago (09/11/2017)
I'm not sure we can help further.If you post a signal segment, maybe we can look at it.David

Re:

Reply posted 7 years ago (09/09/2017)
You just updated your need by saying "My goal with this exercise is to clean my signal from noise and possible interference."If you are trying to remove some noise...

Re:

Reply posted 7 years ago (09/08/2017)
As a result of your post, I tried the goetzel in Matlab, and I cannot imagine how they could have made it more confusing.  It is not well documented.  Why would...

Re:

Reply posted 7 years ago (09/07/2017)
Thank you for setting up the community and making it run.

Re:

Reply posted 7 years ago (09/07/2017)
I tried to type a long response, and it was all lost, because the system thought I was not logged in, which I in fact was.I will reiterate only a very small portion...

Re: audio artifacts sleuthing

Reply posted 7 years ago (09/04/2017)
Which segments do you consider damaged?The earliest utterances are definitely lacking in higher frequency components (muffled).The four counting sequences in the...

Re: Removing a known noise component from a signal

Reply posted 7 years ago (09/04/2017)
I simpathize with Tim.  I don't think you can recover from the terrible non-linearity of the sample data.

Re:

Reply posted 7 years ago (09/04/2017)
Hi,No one has yet suggested a bank of closely space Goertzel detectors.  If you space then fairly close together around the expected frequency, then you should...

Re: Estimating SNR without windowing

Reply posted 7 years ago (06/17/2017)
I have used you fundamental approach several times.I have even used it for data communications, where I synthesized a transmit signal, and manually tinkered with...
I'm not a control expert, but I assume the speed requirements are most likely related to physical world activities.  For many of these, the time constants are...
As a first experiment, run the filtering in the frequency domain.  Overlap add, or overlap save etc..Do 4096 point FFTs and kill the first 8 bins to have an effective...

Re: DFT/FFT

Reply posted 7 years ago (06/05/2017)
I believe that you need the complex data input to do the shift correctly.  The Hilbert Filters can be very simple FIR implementation.  The delayed real input is...

Re: DFT/FFT

Reply posted 7 years ago (06/05/2017)
You might want to create a complex version of you signal by using a Hilbert filter which will then give you a real and imaginary part in the time domain, and then...

Re: complex values in Short-time Fourier Transform

Reply posted 7 years ago (06/01/2017)
Hi,I think that you simply need atan2(y,x) for the phase.  I may be wrong.Remember that atan has a range of only -pi/2 -> pi/2 if you don't have the full values...

Re: GFSK parameters

Reply posted 7 years ago (05/22/2017)
Thanks for your message.  I have seen several articles that mention the same issue of the non-linear relationship that you have mentioned.The problem I am trying...

Re: GFSK parameters

Reply posted 7 years ago (05/22/2017)
Thanks for your response.The frequency deviation I was defining is the deviation from the carrier for the mark and space frequencies, rather than the carrier offset...

GFSK parameters

New thread started 7 years ago
I'm trying to understand the exact relationship between the frequency deviation and the time-bandwidth product for GFSK modulation.The basic set of parameters for...

Re: Average FFT Magnitude in bins

Reply posted 7 years ago (05/17/2017)
Hi,Thanks for the correction to my understanding of the pwelch function.  It has caused me to go back and reconsider some of the work I have done using spectrogram...

Re: Average FFT Magnitude in bins

Reply posted 7 years ago (05/16/2017)
1) Yes each bin is 10 Hz DC, 10Hz, ... Fs/2.  DC and Fs/2 will have real content only.2) Yes. The pwelch function will pull blocks of data in that are 4410 samples...

Re: Average FFT Magnitude in bins

Reply posted 7 years ago (05/15/2017)
Hi,Sorry to take so much time to get back to you.By the time of this writing, you have many of your questions answered, and it sounds like you have the tools needed...

Re: DFT

Reply posted 7 years ago (05/15/2017)
Hi,In answer to the question about conj(f0) vs. 3*f0.  Remember, that in the digital domain, there is a limit of Fs/2 for the highest representable frequency.  Since...

Re: Average FFT Magnitude in bins

Reply posted 7 years ago (05/15/2017)
Just a couple of observations.1) It is clear from the spectra shown, that the recordings were actually captured at 32 kHz sampling before being converted to 44.1...

Re: DFT

Reply posted 7 years ago (05/15/2017)
Hi,1) You are correct.  If you have 4 time domain points, then you will have 4 frequency domain points. The frequency domain points will be complex even if the...

Re: FFT Speed, FIR Output

Reply posted 7 years ago (05/13/2017)
Hi Roger,The Goertzel detector is the most effective solution that I know of with respect to noise.The items that you really need to look at are essentially specs....

Re: FFT Speed, FIR Output

Reply posted 7 years ago (05/12/2017)
Hi Roger,Do you really need to reconstruct the signal or is it really a detection problem?My best suggestion for your detection problem is to use a Goertzel detector...

Re: OFDM sampling

Reply posted 7 years ago (05/12/2017)
I believe that you have all the information (and more) that you need.If the sampling rate that the author was using is 192kHz, and the sub-carrier frequency is 5...

Re: OFDM sampling

Reply posted 7 years ago (05/11/2017)
Sometimes the confusion is related to the fact that some/many of the frequency domain bins may not be occupied for the inverse FFT.In WiFi for example, only 52 bins...

Re: Semi-blind OFDM channel estimation.

Reply posted 7 years ago (04/26/2017)
Hi,I reviewed the document, but I'm not sure what advice you are looking for.  I don't think I can offer more unless I can look at the data myself.David

Re: Semi-blind OFDM channel estimation.

Reply posted 7 years ago (04/24/2017)
Hi SenSen,Your first diagram is exactly as I had suggested.  PR can be used for carrier offset and some phasing information, but not for EQ.  TR1 and TR2 should...

Re: Filtering of amplitude spectrum

Reply posted 7 years ago (04/21/2017)
What about the spectrum do you not like that makes you want to smooth it?Remember, that the spectrum of speech is always very dynamic and rough when comparing across...

Re: Semi-blind OFDM channel estimation.

Reply posted 7 years ago (04/21/2017)
Hi SenSen,Re: 3. above:  The second repetition of the very first sequence has the same spectral problems as the very first sequence. These two repetitions are...

Re: Semi-blind OFDM channel estimation.

Reply posted 7 years ago (04/20/2017)
One more thing.The frequency response of this early sequence has intentional ripples in the transmit signal.  If you equalize the rest of the burst with this early...

Re: Semi-blind OFDM channel estimation.

Reply posted 7 years ago (04/19/2017)
Hi SenSen,A couple of observations.You might be better off doubling your sampling rate and doubling your FFT length.  The reason is that the phase of the sample...
How about a simple sum(abs(coefficients)*32768)?This will tell you the largest result that you could ever get at the output, because it assumes that the samples...

Re: Deriving FIR coefficients for a higher sample rate

Reply posted 7 years ago (04/08/2017)
Hello Rick,Thanks for your kind words.  Coming from someone who has contributed so many extremely useful ideas, I take your comment as a great compliment.I am...

Re: How to calculate time delay estimation?

Reply posted 7 years ago (04/07/2017)
I was going to propose an adaptive filter so that you can run the data over and over to get a good averaged correlation.However, I was not even able to run you code,...

Re: Deriving FIR coefficients for a higher sample rate

Reply posted 7 years ago (04/06/2017)
That was quick.The principle is simple.Normally to up-sample, we would add zeros between existing samples, and run through the newly designed h1.Instead we keep...

Re: Deriving FIR coefficients for a higher sample rate

Reply posted 7 years ago (04/06/2017)
Sorry, I seem to have caused confusion.Use the two sets of coefficients to interpolate your existing FIR filter.1) Your empirical FIR filter, that needs to be resampled, we...

Re: Deriving FIR coefficients for a higher sample rate

Reply posted 7 years ago (04/06/2017)
Hi Probbie,How about a simple interpolation of the coefficients?Design a new FIR filter that is low pass with cutoff at quarter of the sampling rate.  Then divide...

Re: OFDM Burst Preamble.

Reply posted 7 years ago (03/25/2017)
Hi SenSen,The sequence is so short, that I cannot make much of it at all.I performed the standard test for FSK looking for differential phase (frequency):1) ann...

Re: OFDM Burst Preamble.

Reply posted 7 years ago (03/25/2017)
Hi SenSen.  David again.Is this the same data set that we discussed before?I don't believe that there is any FSK sequence for the preamble.I am posting a plot from...

Re: Frequency resolution

Reply posted 7 years ago (03/24/2017)
I like Fred's response.I have an additional thought to keep in mind.If there are an integer number of cycles of the sine/s in the length of you data sample, then...

Re: Bit width in a digital delta sigma

Reply posted 7 years ago (03/14/2017)
Hi Kaz,I can see your concern, but every accumulator is preceded by a subtraction from which the output magnitude is expected to be less than either of its inputs....

Re: Bit width in a digital delta sigma

Reply posted 7 years ago (03/13/2017)
Hi Deba,Looking at the block diagram that you have shared with us, I conclude that all the bit precision at each stage is set up to preserve the required number...

Re: Signal power at the receiver

Reply posted 7 years ago (03/10/2017)
The only difference for the baseband will be the signal to quantization ratio SQR.For most modulation schemes though, there is an optimal received signal level.For...

Re: Signal power at the receiver

Reply posted 7 years ago (03/10/2017)
My immediate response it in relation to the baseband signal to which you refer.The level of the analog signal will also define the signal to quantization ratio....

Re: Signal power at the receiver

Reply posted 7 years ago (03/09/2017)
At the RF front end of the receiver, there is a Low Noise Amplifier (LNA).Sometimes, there is both an analog and a digital AGC that adjusts the signal to bring into...

Re: Signal power at the receiver

Reply posted 7 years ago (03/09/2017)
dudelsound is exactly right about the issues of SNR.Everything about the receiver should be designed to maximize the SNR and in the "N" he has correctly added the...

Re: Who clones DSP56362?

Reply posted 7 years ago (03/08/2017)
I agree with Joe West on the programming issues.Have you considered a huge heat sink on top of the DSP to keep it cool.  Lots of wings.Maybe it requires some clever...

Re: FHSS Signal Analysis.

Reply posted 7 years ago (03/08/2017)
Hi SenSen,I apologize for taking so long to get back to you.Did you get everything figured out?David
Nelson,My appologies for the delay.I have just sent an email from my private account.David
Thanks for your kind response.I guess that you can tell that I have a passion for this stuff.Best wishes for your research.David
Nelson,What you require, is certainly the Nth level of exactness.However, if you just want to cancel enough so that the residual is pretty small and basically inperceptable....
Hi Nelson,Wow!  You have just set up a "contest" to see if we can do it.I love being a sleuth.  :-)1) To do what you said means that you can never send any of...
Hi Nelson,What you are after is related to a fairly standard system identification problem.  Very often this identification is made using an adaptive filter.  Echo...
Thanks for this reminder.When I use the RIR software for simulating a room, I assign a positive sign to one and and a negative sign to the other of two opposing...
Hi Tim,Part of my thinking relates to the way in which a recording studio is designed to avoid "modes".  The walls are never parallel, the glass window - used by...
Sorry Joe, but I'm now going to nit-pick.I believe Tim is correct. Maybe our definition of reflect and scatter are different, but reflect to me implies bounce off...
Hi Nelson,First answer is that you assertions are not correct.Only the amplitude of the various frequency components is changed by the reflections off the surfaces.Reflections...

Re: Interpolation and filter gain

Reply posted 7 years ago (02/21/2017)
I really like Tim's answer.I would like to add a couple of ideas related to the filter design.A typical implementation would use a polyphase decomposition of the...

Re: FHSS Signal Analysis.

Reply posted 7 years ago (01/31/2017)
Hi SenSen,I am working on a more formal reply, but for now, please note above that i was using fhss_2 for my experiments.  Your document indicates that you were...
Hi Tim,I realize that you used the term gut feeling, but I would like to know why your gut believes that 20 dB is the best you could get?David

Re: FHSS Signal Analysis.

Reply posted 7 years ago (01/23/2017)
Hi SenSen,I have decoded the first L2 correlation sequence and found 54 subcarriers each side of DC. So there are 108 occupied bins.  Each is carrying a BPSK constellation.I...

Re: Spurious when Fs/F is not an integer

Reply posted 7 years ago (01/23/2017)
Hi Loganathan N, Please be careful not to confuse the MatLab plot display with some real spurious modulation. The plot of the waveform shows an envelope due to...

Re: Upsampled input to an Adaptive filter?

Reply posted 7 years ago (01/21/2017)
Hello khurram6050,You said that you up-sampled and then interpolated.Then you said that you also tried after you low pass filtered the inputs.This lowpass filtering...

Re: FHSS Signal Analysis.

Reply posted 7 years ago (01/21/2017)
Hi SenSen,I have modulated the fhss_2 source by exp(2i*pi*32.99666e6) and resampled it at 5.6e6 Hz and find some very satisfactory results using a 512 point FFT.Attached...

Re: FHSS Signal Analysis.

Reply posted 7 years ago (01/19/2017)
Hi SenSen,You are correct that there appear to be three copies of the short sequence.Careful review will show that the third is in fact a reflection of the back...

Re: FHSS Signal Analysis.

Reply posted 7 years ago (01/16/2017)
Hi SenSen,A little more digging verifies that the correlation sequence pair that is usually allocated early - for fine frequency estimation and fine timing alignment...

Re: FHSS Signal Analysis.

Reply posted 7 years ago (01/16/2017)
Hi SenSen,Your modulation looks like OFDM.I simply looked at the fhss_1 data.  It is sitting at about 29e6 based on the 40e6 sampling rate that you mention.I...

Re: How to Convert IQ Signal to Real

Reply posted 7 years ago (12/15/2016)
I know that his tread is now really old, but I'm thinking that the steps needed are in fact a lot simpler than the rest of the thread.How about using the dsp.DigitalDownConverter...

Use this form to contact dgshaw6

Before you can contact a member of the *Related Sites:

  • You must be logged in (register here)
  • You must confirm you email address