PLRandom.m

for ii=1:100
   a0=0;
a1=8e-6;
k=0.5;
j=sqrt(-1);
%di=[150];
di=150*rand(1,10);
%gi=[0.4];
gi=0.2*(rand(1,10)-0.5);
f=[1:1000]/1000*30e6;
for i=1:1000
   Hf(i)=sum(gi.*exp(-(a0+a1*f(i)^k).*di).*exp(-2e-6/150*pi*j*f(i).*di));
end
Hfcj=real(Hf(1000:-1:1))-j*imag(Hf(1000:-1:1));
HHf=[Hf Hfcj];
imp=real(ifft(HHf));
figure(1)
plot(f,20*log10(abs(Hf)))
figure(2)
plot(imp)
pause
end