fs=1e006;
tf4=1.5e-006;a4=1/tf4;
tr4=0.9e-006;b4=1/tr4;
t=0:1/fs:200/fs;
e3=rand(size(t));
e4=zeros(size(t));
N=length(t);
for k=1:N
if t(k)>=50/fs & t(k)<=150/fs
e4(k)=exp(-(t(k)-(50)/fs)*a4)-exp(-(t(k)-(50)/fs)*b4);
end
end
s=e4+0.008*e3;
s=round(s*32767);
figure(2);
plot(t*fs,s);title('signal');xlabel('time(us)');