This calculates diffraction from a thick phase grating, as for ex. in an
acousto-optic modulator, in each order as functions of propagation distance (in
units of Talbot length), modulation index (in units of Pendell\"osung
length), and incidence angle (in units of Bragg angle) (see help above).
The output is expressed as a matrix giving the diffracted order and any of the
other 3 inputs.
// i/pi*dS_j/dz = j*(th-j)*S_j+a*(S_(j+1)+S_(j-1))
// sinusoidal grating diffraction: power in j-th mode (input in 0-th)
// rn(j=order, s=z/Ltalbot, a=Ltalbot/Lpendell., th=angle(1=bragg))
// (Ltalb=period^2*n/wvlgth0, Lpend=2*wlgth*sqrtepsilon/deltaepsilon,
// braggangle=wvlgth0/(2*period) ) (j may be array)
// input: row vectors, output: column vectors or matrix(j,s)
// (accuracy of the order of 0.001)
Upload date : 2011-06-15 11:40:30
MD5 : 0da461ae974d0ae887781caf3fc07f70
SHA1 : ba948d94c4b1989b3ee2692fd9995578fc7a4d81
Downloads : 708
I can not get script running
Send me instructions who to exec it
I have tried to call the functions as:
> nothing happens
just remove the ending ";"
script is not working
Nice function applied to physics.
But please remove the heading "clear"!
By the way, including the heading bloc of comments in the function, just after the
"function rana=..." line without any uncommented blank line in-between, would
using head_comments() to view these comments as help lines in a handy way, without editing
the script: http://help.scilab.org/docs/5.5.2/en_US/head_comments.html