correct my integral in line43

2 次查看（过去 30 天）

显示更早的评论

asma 2024-4-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2108746-correct-my-integral-in-line43

编辑： Torsten 2024-4-18

在 MATLAB Online 中打开

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');