- #1
renxiaoxie
- 1
- 0
Fiber Grating Using Finite Difference Method
I have to solve the Fiber grating coupled-equations using FDM. These equations are coupled each other and strange inital conditions. That means ode45 from the MATLAB library function doesn't work here!
I am totally fresh in numerical analysis and read some books, but unlukcy, i still confuse how to do this code:(
Anyone has the program that solve this problem, pls share with me! Thank you
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The probelm is:
Parameters described:
AF(z,t) is the amplitude of forward-propagation-waves;
AB(z,t) corresponds to the backward-propagation-waves;
z and t are space and time coordinates;
AFz = dAF/dz AFt = dAF/dt;
ABz = dAF/dz ABt = dAF/dt;
i=sqrt(-1);
vg, delta, kapa, gama are constant terms.
The coupled-equations is:
i*AFz + i*AFt/vg + delta*AF + kapa*AB + gama*(|AF|^2+2*|AB|^2)*AF = 0
-i*ABz + i*ABt/vg + delta*AB + kapa*AF + gama*(|AB|^2+2*|AF|^2)*AB = 0
and the inital conditons:
AF(0,t)=A(t);
AF(z,0)=0;
AB(L,t)=0;
AB(z,0)=0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I have to solve the Fiber grating coupled-equations using FDM. These equations are coupled each other and strange inital conditions. That means ode45 from the MATLAB library function doesn't work here!
I am totally fresh in numerical analysis and read some books, but unlukcy, i still confuse how to do this code:(
Anyone has the program that solve this problem, pls share with me! Thank you
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The probelm is:
Parameters described:
AF(z,t) is the amplitude of forward-propagation-waves;
AB(z,t) corresponds to the backward-propagation-waves;
z and t are space and time coordinates;
AFz = dAF/dz AFt = dAF/dt;
ABz = dAF/dz ABt = dAF/dt;
i=sqrt(-1);
vg, delta, kapa, gama are constant terms.
The coupled-equations is:
i*AFz + i*AFt/vg + delta*AF + kapa*AB + gama*(|AF|^2+2*|AB|^2)*AF = 0
-i*ABz + i*ABt/vg + delta*AB + kapa*AF + gama*(|AB|^2+2*|AF|^2)*AB = 0
and the inital conditons:
AF(0,t)=A(t);
AF(z,0)=0;
AB(L,t)=0;
AB(z,0)=0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%