david753
Full Member level 1
In system, we know the input signal X(t) convolution H(t) equals Y(t).
Mathmatic notation is X(t)*H(t)=Y(t).
In my lab, I do an interesting experience.
I input a unit step signal, ch1 into a system(a band pass filter), and measuring its output signal, ch2.
Then, I use deconvolution formula to retrieve system's Impulse response.
In matlab, it is presented as [q,r]=deconv(y,x).
Then, FFT(q), and got the frequency response.
Am I right?
But, unfortunately, the result is wrong.
Why???
Please refer to as following matlab code.
ch1--> **broken link removed**
ch2--> **broken link removed**
===============================
load ch1.txt;
plot(ch1);
hold on;
load ch2.txt
plot(ch2,'r');
hold on;
[q,r]=deconv(ch2,ch1);
plot(q,'g');
q_fft=abs(fft(q(1:250))).^2;
figure;
plot(q_fft);
=================================
P.S. scope sampling frequency is 1GHZ.
Mathmatic notation is X(t)*H(t)=Y(t).
In my lab, I do an interesting experience.
I input a unit step signal, ch1 into a system(a band pass filter), and measuring its output signal, ch2.
Then, I use deconvolution formula to retrieve system's Impulse response.
In matlab, it is presented as [q,r]=deconv(y,x).
Then, FFT(q), and got the frequency response.
Am I right?
But, unfortunately, the result is wrong.
Why???
Please refer to as following matlab code.
ch1--> **broken link removed**
ch2--> **broken link removed**
===============================
load ch1.txt;
plot(ch1);
hold on;
load ch2.txt
plot(ch2,'r');
hold on;
[q,r]=deconv(ch2,ch1);
plot(q,'g');
q_fft=abs(fft(q(1:250))).^2;
figure;
plot(q_fft);
=================================
P.S. scope sampling frequency is 1GHZ.