SLoscillator.zip_ SLoscillator_SLoscillator_agrawal_激光_激光器 模型

%边发激光器动态仿真 clear; close all; eel_parameter; global Anr B C L c beta Gama am aint a vg Ntr V e; landa=1.3e-6;%激光波长 Rf=0.3; Rb=0.3; Nth=(am+aint)/Gama/a+Ntr; Ith=(Anr*Nth+B*Nth*Nth+C*Nth*Nth*Nth)*e*V; T0=2;% ns delta_t=T0/100; t=0:delta_t:8*T0; options=optimset('Display','off'); I=10;%电流单位(mA) % x0=1.271709;%初始值 % x0=((am+aint)/Gama/a+ntr)/1e24; x0=2.127; x=fsolve('eel_I_N',x0,options,I,1.3e-6,Rf,Rb); N0=x; % N0=x*1e24; %------------算出光子密度--------- g=a*(N0*1e24-Ntr);%增益系数 Gs=exp((Gama*g-aint)*L); S0=1/vg*(beta*Gama*B/1e-15*x*x)/(am+aint-Gama*g); options = odeset('RelTol',1e-4,'AbsTol',[1e-10 1e-15]); UU=[0.6;0]; % UU=[1000;1000]; [t,y]=ode45(@eel_time_nidy,t,UU,options,I,Rf,Rb,landa); plot(t,y(:,1),t,x);figure plot(t,y(:,2)); plot(t,x-y(:,1),'r') % fid1=fopen('Carrier.dat','w'); % fid2=fopen('Photon.dat','w'); % for ii=1:1:length(t) % data_1=[t(ii);y(ii,1)]; % data_2=[t(ii);y(ii,2)*1e3]; % fprintf(fid1,'%10.8f %10.8f\n',data_1); % fprintf(fid2,'%10.8f %10.8f\n',data_2); % end % fclose(fid1); % fclose(fid2);