Returns the source normalization spectrum used to normalize data in the cwnorm state for standard fourier transform quantities. See the Units and normalization chapter for more information. If the source time signal of the jth source in the simulation is sj(t), and N is the number of active sources then
$$ s(\omega)=\operatorname{sourcenorm}(\omega)=\frac{1}{N} \sum_{s o u r c e s} \int \exp (i \omega t) s_{j}(t) d t $$
Syntax |
Description |
---|---|
out = sourcenorm(f); |
Returns the source normalization used to normalize data in the cwnorm state at the vector of frequency points f. (f is the frequency in Hz) |
out = sourcenorm(f, name); |
This function makes it possible to perform the normalization using the spectrum of one source, rather than the sum of all the sources. |
Example
This example shows how to reproduce the source spectrum figures shown in the Frequency/Wavelength tab of the source property window.
lambda1 = 0.4e-6; # start wavelength
lambda2 = 0.7e-6; # end wavelength
f=linspace(c/lambda2,c/lambda1,1000);
# get the source spectrum
spectrum=sourcenorm(f);
spectrum=abs(spectrum)^2;
spectrum=spectrum/max(spectrum);
# get the source time domain signal
time = getdata("source","time");
time_signal = getdata("source","time_signal");
plot(c/f*1e6,spectrum, "wavelength (um)","spectrum vs wavelength");
plot(f/1e12,spectrum, "frequency (THz)","spectrum vs frequency");
plot(time*1e15,time_signal,"time (fs)","amplitude","Source time signal");
See Also
sourcenorm2_avg, sourcenorm2_pavg, sourcepower, cwnorm, nonorm, Units and normalization