1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66
| **********Doppler function***************************
function [fm,am,iflaw]=doppler(N,Fs,f0,d,v,t0,c);
// Parameters
// N : number of points.
// FS : sampling frequency (in Hertz).
// F0 : target frequency (in Hertz).
// D : distance from the line to the observer (in meters).
// V : target velocity (in m/s)
// T0 : time center (default : N/2).
// C : wave velocity (in m/s) (default : 340).
// FM : Output frequency modulation.
// AM : Output amplitude modulation.
// IFLAW : Output instantaneous frequency law.
// Description
// Returns the frequency modulation (FM), the amplitude
// modulation (AM) and the instantaneous frequency law (IFLAW)
// of the signal received by a fixed observer from a moving target
// emitting a pure frequency f0.
else
tmt0=((1:N)'-t0)/Fs;
dist=sqrt(d^2+(v*tmt0).^2);
fm = exp(j*2.0*pi*f0*(tmt0-dist/c));
if (nargout>=2),
if abs(f0)<eps,
am=0;
else
am= 1.0 ./ sqrt(dist);
end
end;
if (nargout==3), iflaw=(1-v^2*tmt0./dist/c)*f0/Fs; end;
end ;
********************Doppler function************************* |
Partager