- #1

- 1

- 0

it=xt.*cos(2*pi*f0*t);

ht=9900*sinc((9900)*(t)).*exp(((20000*pi-200*pi)/2+200*pi)*j*t);

ift=conv(it, ht*ts);

qt=xt.*(cos(2*pi*f0*t)-sin(2*pi*f0*t));

qft=conv(qt, ht*ts);

cft(ift, ts,1)

The above is my lowpass filter. Then I used numerical convolution to calculate an approximate output of the filter output i(t) and q(t) when the input is i(t) and q(t) respectively:

for k=1:50

yt=ift+j*qft;

yk =yt(1:length(t)).* ((t>=.02 .*k - .02) - (t>= .02 .*k + .02)) ;

[Yw, w] = cft(yk, ts, 1 ); % calculate the Fourier transform, where

[maxY w_index] = max(abs(Yw)); % find the peak in ||xw||

omega_m(k) = w(w_index); % omega_m(k) is the Doppler shift frequency

end

but when I multiply y(t) by p(t-t0) I get 0 for my w0 and blood velocity is zero everywhere..What am I supposed to do?