correct my integral in line43
2 次查看(过去 30 天)
显示 更早的评论
please correct line 43 that is about integral function.
tnx a million.
clc;
clear;
close all;
% define sellmier equation to obtain refractive index
T=25; % tempreture in c
F = (T-24.5)*(T+570.5);
A1=4.5820; A2=0.09921; A3=0.21090; A4=0.021940; B1=5.2716E-8; B2=4.9143E-8; B3=2.2971E-7;
lambda_sh=0.775e-6; % SH wavelength
nlithium_sh = (sqrt(A1)+((A2+B1*F)/((lambda_sh^2)-(A3+B2*F)^2))+(B3*F)-(A4*(lambda_sh)^2));
lambda_ff=1.55e-6; % FF wavelength
nlithium_ff = (sqrt(A1)+((A2+B1*F)/((lambda_ff^2)-(A3+B2*F)^2))+(B3*F)-(A4*(lambda_ff)^2));
% parameters
k_sh=(2*pi/lambda_sh)*nlithium_sh;
k_ff=(2*pi/lambda_ff)*nlithium_ff;
capital_lambda = 20e-6;
Delta_kk = k_sh-2*k_ff;
l=pi/Delta_kk; % for first order QPM
Delta_k=(Delta_kk-pi)/l;
lc = pi/Delta_kk; % lc (coherence length);
L = 10e-3; % Length;
N=L/(capital_lambda/2); % number of domains
z=L;
% for d(z) as fourier serie
Km=2*pi/20e-6;
Gm=2/pi;
d33=27; % unit=pm/V
deff=(2/pi)*d33;
d0=(2/pi)*deff; % unit=pm/V
D=lc/capital_lambda;
% for SH field
syms z
GG = symsum(Gm * exp(1i * Km * z), 1, inf);
dz = GG * d0; % periodic modulation
gzz = (pi/2*L)*(dz)/d0;
syms z
f = gzz * exp(1i * Delta_kk * z);
int(f,z,0,L);
y=int(f,z,0,L);
figure(1)
plot(y)
xlabel('Delta_k');
ylabel('SHG efficiency');
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!