成人女人看片免费视频放人_亚洲色精品三区二区一区_欧美亚洲国产精品久久_成人无遮挡裸免费视频在线观看_97SE亚洲国产综合在线_精品久久久久久777米琪桃花_天天躁日日躁很很很躁_色噜噜狠狠一区二区三区果冻_国产免费久久精品国产传媒_67194成是人免费无码

切換到寬版
  • 廣告投放
  • 稿件投遞
  • 繁體中文
    • 8542閱讀
    • 1回復(fù)

    [分享]求解光孤子或超短脈沖耦合方程的Matlab程序 [復(fù)制鏈接]

    上一主題 下一主題
    離線tianmen
     
    發(fā)帖
    58
    光幣
    15
    光券
    0
    只看樓主 倒序閱讀 樓主  發(fā)表于: 2011-06-12
    計算脈沖在非線性耦合器中演化的Matlab 程序 ~K96y$ DTE  
    E?mp6R]}%  
    %  This Matlab script file solves the coupled nonlinear Schrodinger equations of B|=maz:_  
    %  soliton in 2 cores coupler. The output pulse evolution plot is shown in Fig.1 of 5r<(Z0  
    %  Youfa Wang and Wenfeng Wang, “A simple and effective numerical method for nonlinear +I n"OR%  
    %   pulse propagation in N-core optical couplers”, IEEE Photonics Technology lett. Vol.16, No.4, pp1077-1079, 2004 zrO|L|F&P  
    ug,|'<G+  
    %fid=fopen('e21.dat','w'); V^z;^mdd  
    N = 128;                       % Number of Fourier modes (Time domain sampling points) D,j5k3< #  
    M1 =3000;              % Total number of space steps wKlCx  
    J =100;                % Steps between output of space  {Yk20Zn  
    T =10;                  % length of time windows:T*T0 } XU:DE  
    T0=0.1;                 % input pulse width --YUiNhh  
    MN1=0;                 % initial value for the space output location S1`0d9ds#  
    dt = T/N;                      % time step Iq]6]  
    n = [-N/2:1:N/2-1]';           % Index ~ p.W*skD  
    t = n.*dt;   "T%'Rp`j|  
    u10=1.*sech(1*t);              % input to waveguide1 amplitude: power=u10*u10 -!>ZATL<B  
    u20=u10.*0.0;                  % input to waveguide 2 %QgAilj,  
    u1=u10; u2=u20;                 |a$w;s>\  
    U1 = u1;   eKz~viM'  
    U2 = u2;                       % Compute initial condition; save it in U SWsv,  
    ww = 4*n.*n*pi*pi/T/T;         % Square of frequency. Note i^2=-1. *.c9$`s  
    w=2*pi*n./T; Q|KD/s??  
    g=-i*ww./2;                    % w=2*pi*f*n./N, f=1/dt=N/T,so w=2*pi*n./T XHM"agrhSQ  
    L=4;                           % length of evoluation to compare with S. Trillo's paper &[Zap6]  
    dz=L/M1;                       % space step, make sure nonlinear<0.05 zZc@;S#  
    for m1 = 1:1:M1                                    % Start space evolution T;r];Y(b*  
       u1 = exp(dz*i*(abs(u1).*abs(u1))).*u1;          % 1st sSolve nonlinear part of NLS 64' 涿州市| 富民县| 武汉市| 龙江县| 宁蒗| 吴桥县| 成武县| 平谷区| 从江县| 北宁市| 大埔县| 广南县| 七台河市| 孝昌县| 潞西市| 贡嘎县| 宁化县| 喀喇| 玉龙| 德惠市| 丹东市| 武穴市| 乐昌市| 农安县| 民乐县| 江华| 兴城市| 荔浦县| 吴川市| 荆门市| 新巴尔虎左旗| 高尔夫| 湘潭县| 托里县| 郁南县| 玉龙| 麻栗坡县| 分宜县| 南京市| 西贡区| 温宿县|