fushe.rar_Image Encryption_encryption matlab_图像加密matlab_图像恢复matl

######################################################################## #This is an example of the input file for F2P calculations. #It demonstrates the negative refraction in a PC with a triangular # lattice of air holes. #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #------------------------------------------------Prepared by Dr. Min Qiu ######################################################################## ######################################################################## #Calculate the TM or TE modes?? # ==0: TE modes # ==1: TM modes 1 ######################################################################## ######################################################################## #Define the materials used in the calculations. #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #The number of materials 2 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #The dielectric constant, and conductivity of each material #Index---dielectric constant---conductivity 1 1.0000 0.00E+0 #air 2 10.5000 0.00E+0 #InP/GaInAsP (n=3.24) #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #End of material definition. ######################################################################## ######################################################################## #define the size of the computational domain #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #--- lattice constant (For normalized purpose) 1.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #The size of the computational domain. #---Sx-------Sy-- 60.0 30.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Background materials index 2 ######################################################################## ######################################################################## #Absorbing boundary condition (ABC) #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #ABC type: =0 All surrounded by PML; =1 Only the Y direction surrounded #by PML, while periodic BC in the X direction. #---ABC type---------------------------------------------------------- 0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #---number of PML layers in X----in Y---------The exponent index for PML 12 12 2.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ######################################################################## ######################################################################## #Define parameters for FDTD #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #----Dx---------Dy------ 0.05 0.05 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Dt coefficient. Dt=Dt_coe*1/sqrt(1/(Dx*Dx)+1/(Dy*Dy))/vc #Should be less than 1.0. #For example, 0.95 or 0.866 are all good choices. #---Dt_Coe--- 0.95 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Number of total time steps 10000 ######################################################################## ######################################################################## #Output the z-field component in some time steps. #For snapshot purpose. #If interval is less than 1, no output. # The current z-field distribution is always stored in "Fieldz_Bin.dat". #----Start Time steps---------End Time steps------Interval-------------- 20 10000 -20 ######################################################################## ######################################################################## #define the inclusions, could be waveguides, cylinders, ... #----------------------------------------------------------------------- #Inclusion type: # == 1: ellipse # == 2: quadrilateral # == 3: triangular #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #The total number of inclusions 2 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Inclusion No. 1, circles, layer down 1 #Inclusion Type ----Material index---- 1 1 #-----Rx---------Ry---------Xc---------Yc--------Theta-- 0.36 0.36 16.00 0.5 0.0 #-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat) 6 1.73205 30 1.00 0.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Inclusion No. 2, circles, layer down 2 #Inclusion Type ----Material index---- 1 1 #-----Rx---------Ry---------Xc---------Yc--------Theta-- 0.36 0.36 16.86603 1.0 0.0 #-----Mx--------Lx--------My--------Ly------Alpha-----(For Repeat) 6 1.73205 29 1.0 0.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Examples of Triangular #Inclusion Type ----Material index---- # 3 1 #---X1------Y1------X2------Y2------X3-------Y3--- # 10.5 0.0 24.5 3.964 16.0 9.964 #---Mx-----Lx-----My-----Ly----Alpha---(For Repeat) # 1 0.00 1 0.00 0.0 ######################################################################## ######################################################################## #Source part #----------------------------------------------------------------------- # For source type: ==1 point source # ==2 line source, direction, from P1 to P2, in P3 # ==3 waveguide mode source #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #----Number of Source------------ 1 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Point Source Examples #---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)- # 1 0.305 0.007 #----XS-----YS-------- (Point source position) # 4.0 6.06218 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Gaussian Line Source sample #For line source, direction == +1(X); -1(-X); +2(+Y); -2(-Y); #---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)- # 2 0.25 0.1 #----direction---XS-----YS-----Beam Width---Material Type----Theta------- # 1 9.00 6.06218 2.0 2 0.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Waveguide Source sample #For Waveguide source, direction == +1(X); -1(-X); +2(+Y); -2(-Y); #---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)- # 3 0.295 0.007 #----direction---XS-----YS-----Waveguide Width---Mode Number----------- # 1 1.00 6.06218 2.60 1 #---Material Up (Left)----Material Core----Material Down (Right)------- # 1 2 1 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: #Gaussian Line Source sample #For line source, direction == +1(X); -1(-X); +2(+Y); -2(-Y); #---Source Type------Source frequency (a/lambda)---Pulse Width (Delta_w)- 2 0.325 -0.1 #----direction---XS-----YS-----Beam Width---Material Type----Theta------- 1 1.00 15.0 2.0 2 0.0 #::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ######################################################################## ##########################################################