Graphics Programs Reference
In-Depth Information
for tau = -1.5*taup:.01:1.5*taup
ii = ii + 1;
val1 = 1. - abs(tau) / taup;
val2 = pi * taup * (1.0 - abs(tau) / taup);
val3 = (fd + mu * tau);
val = val2 * val3;
x(ii) = abs( val1 * (sin(val+eps)/(val+eps)));
end
figure(1)
plot(taux,x)
grid
xlabel ('Delay - seconds')
ylabel ('Uncertainty')
figure(2)
plot(taux,x.^2)
grid
xlabel ('Delay - seconds')
ylabel ('Ambiguity')
Listing 4.7. MATLAB Function Ðtrain_ambg.mÑ
function x = train_ambg (taup, n, pri)
if( taup > pri / 2.)
'ERROR. Pulsewidth must be less than the PRI/2.'
return
end
gap = pri - 2.*taup;
eps = 0.000001;
b = 1. / taup;
ii = 0.;
for q = -(n-1):1:n-1
tauo = q - taup ;
index = -1.;
for tau1 = tauo:0.0533:tauo+gap+2.*taup
index = index + 1;
tau = -taup + index*.0533;
ii = ii + 1;
j = 0.;
for fd = -b:.0533:b
j = j + 1;
if (abs(tau) <= taup)
val1 = 1. -abs(tau) / taup;
val2 = pi * taup * fd * (1.0 - abs(tau) / taup);
val3 = abs(val1 * sin(val2+eps) /(val2+eps));
Search WWH ::




Custom Search