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
| Po= 4; % [mW] initial peak power of signal source
C = 0; % Chirping Parameter
m = 1; % Super Gaussian parameter (m=1 ==> Gaussian)
t0 = 100e-12; %[sec] initial pulse width
vgroup= 2e8;
alpham=0.24;
attn=0.24;
taum = 10000;
dtau = 2*taum/2^12; %1.953
tunit= 1e-12; % make time unit in psec
tau = (-1000:dtau:(taum-dtau))*tunit; %
% fs = 1/(dtau*tunit);
at = sqrt(Po)*exp(-0.5*(1+i*C)*(tau/t0).^(2*m)); % Input field in the time domain
subplot(2,2,1);
plot(tau,at);
title('input wave');
for x=0:0.1:16;
N=length(x);
%bt = sqrt(Po)*exp(-0.5*((tau-((2*x)/vgroup))/t0).^2)*exp(3.84);
%bt = sqrt(Po)*exp(-0.5*((tau-((2*x)/vgroup))/t0).^2)*attn;
bt = sqrt(Po)*exp(-0.5*((tau*N)/t0).^2)*attn;
subplot(2,2,2)
plot(N,bt,'-rx');
end
l=trapez(N,bt);
subplot(2,2,3);
plot(tau,l,'-rx'); |
Partager