problem with IFFT and Convolution

8 views (last 30 days)
Maxilyn Tan
Maxilyn Tan on 19 Aug 2015
Commented: Maxilyn Tan on 28 Aug 2015
I’m having difficulty with executing IFFT and convolution
I want to see the output of a system with an input sine wave x(t) and transfer function H(f). In order to get the output y(t), I used IFFT with the transfer function to convert it from the frequency domain to the time domain. After which I convoluted both h(t) and x(t) to obtain y(t).
However when I plot y(t) I get a distorted signal. Could someone please instruct me as to what I did wrong? Your help would be highly appreciated
%Transfer function (multipathmodel)
f = 0:1*10^6:30*10^6; %frequency of multipathmodel
a = 0+(7.8*10^-10)*f.^1;
l = [200 225]; %length of wire
g = [0.55 0.45];
Vp = 150*10^6;
h = 0;
for i = 1:2
temp = (g(i)*exp(-a.*l(i))).*exp((-f.*j*2*pi)*l(i)/Vp);
h = h + temp;
%h is in the frequency domain
d1 = 10*log(abs(h)); %H(f)/dB
d2 = ifft(h); %Impulse response
%sine wave
fo = 3*10^6; %frequency of the sinewave
Fs = 2^nextpow2((fo*10*2)+fo); %samplingrate
Ts = 1/Fs; %sampling time interval
E = 1/fo; %duration of sinewave (t)
t = 0:Ts:E; %xaxis
x = sin(2*pi*fo*t); %sinewave equation
%convolution of sine wave and transfer function
y = conv2(x,d2);
%plot transfer function
title('Transfer function of Twopathmodel')
grid on
%plot the sine wave
title('Sine wave')
xlabel('Time secs')
grid on
%plot the impulse response of the transfer function
t2 = (0:E/(length(d2)-1):E);
%axis([0 3.5*10^-6 -1 1])
title('Impulse response')
%plot the convolution
t1 = (0:E/(length(y)-1):E);
%axis([0 3.5*10^-6 -1 1])
Maxilyn Tan
Maxilyn Tan on 28 Aug 2015
This is what i'm expecting to get
Here is my input, x(t)
my transfer function H(f)
I converted the transfer function to the time domain by using ifft, this was the result
after which I convoluted both x(t) and the impulse response to get the output, then I used fft. The graph is supposed to be the same as the first graph above but i'm getting something different

Sign in to comment.

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by