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 : 1694
I can not get script running Send me instructions who to exec it I have tried to call the functions as: rana(1,1,1,1); nothing happens Thanks Luis
> rana(1,1,1,1); > > nothing happens just remove the ending ";" ;)
script is not working crap
Hello, 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 allow 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 Thanks Samuel
