I have multiple signal plots in one graph and i want to add their amplitudes to get one graph only

1 view (last 30 days)
This is the code:
n1=1.45;
n2=2.65;
c=299792458;
a=259*10^-9;
b=142*10^-9;
G=a+b;
Lambda0=1.5*10^-6;
Omega0=2*pi/Lambda0;
k1=2*pi*n1/Lambda0;
k2=2*pi*n2/Lambda0;
K=acos((-1/2).*((n1./n2)+(n2./n1)));
A=(exp(1i*k1*a))*(cos(k2*b)+((1i/2).*((n1./n2)+(n2./n1)).*sin(k2*b)));
B=(exp(-1i*k1*a))*((1i/2).*((n2./n1)-(n1./n2)).*sin(k2*b));
C=(exp(1i*k1*a))*((-1i/2).*((n2./n1)-(n1./n2)).*sin(k2*b));
D=(exp(-1i*k1*a))*(cos(k2*b)-((1i/2).*((n1./n2)+(n2./n1))).*sin(k2*b));
M=[D -B;-C A];
a0=B;
b0=(exp(1i.*K))-A;
N=transpose([a0 b0]);
z=0:0.5:10;
E=zeros(1,length(z));
for j=1:10
v=(M^j)*N;
w=inv([1 1;1i*k2 -1i*k2])*[exp(1i*k1*a) exp(-1i*k1*a);1i*k1*exp(1i*k1*a) -1i*k1*exp(1i*k1*a)]*v;
if z>(j*G-a) & z<j*G
E(j,:)=(v(1)*(exp(1i*k1*((j.*G)-z))))+(v(2)*(exp(-1i.*k1.*((j.*G)-z))));
hold on
else
E(j,:)=(w(1).*(exp(1i.*k2*((j.*G)-(a+z)))))+(w(2).*(exp(-1i*k2*((j.*G)-(a+z)))));
hold on
end
hold off
plot(z,abs(E));
end

Accepted Answer

Mathieu NOE
Mathieu NOE on 6 Apr 2021
hello
here you are :
n1=1.45;
n2=2.65;
c=299792458;
a=259*10^-9;
b=142*10^-9;
G=a+b;
Lambda0=1.5*10^-6;
Omega0=2*pi/Lambda0;
k1=2*pi*n1/Lambda0;
k2=2*pi*n2/Lambda0;
K=acos((-1/2).*((n1./n2)+(n2./n1)));
A=(exp(1i*k1*a))*(cos(k2*b)+((1i/2).*((n1./n2)+(n2./n1)).*sin(k2*b)));
B=(exp(-1i*k1*a))*((1i/2).*((n2./n1)-(n1./n2)).*sin(k2*b));
C=(exp(1i*k1*a))*((-1i/2).*((n2./n1)-(n1./n2)).*sin(k2*b));
D=(exp(-1i*k1*a))*(cos(k2*b)-((1i/2).*((n1./n2)+(n2./n1))).*sin(k2*b));
M=[D -B;-C A];
a0=B;
b0=(exp(1i.*K))-A;
N=transpose([a0 b0]);
z=0:0.5:10;
E=zeros(1,length(z));
for j=1:10
v=(M^j)*N;
w=inv([1 1;1i*k2 -1i*k2])*[exp(1i*k1*a) exp(-1i*k1*a);1i*k1*exp(1i*k1*a) -1i*k1*exp(1i*k1*a)]*v;
if z>(j*G-a) & z<j*G
E(j,:)=(v(1)*(exp(1i*k1*((j.*G)-z))))+(v(2)*(exp(-1i.*k1.*((j.*G)-z))));
else
E(j,:)=(w(1).*(exp(1i.*k2*((j.*G)-(a+z)))))+(w(2).*(exp(-1i*k2*((j.*G)-(a+z)))));
end
end
figure(1),plot(z,abs(E)); % original (invidual) plot
figure(2),plot(z,sum(abs(E),1)); % all curves summed plot

More Answers (0)

Categories

Find more on 2-D and 3-D Plots in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!