main title on subplot using for loop

1 view (last 30 days)
son
son on 28 Aug 2014
Commented: Image Analyst on 29 Aug 2014
hi, i created a whole programme but don't know how to create one sub title using for loop:
if type =1 => main title: G.652 fibre
if type =2 => main title: G.653 fibre
if type =3 => main title: G.655 fibre
clear all
close all
%%Insert value in Matlab
n=[10 40 80];
p=1*10^-3;
delta_lambda = 0.4*10^-9;
distance=1000;
type=1:3;
alphadbset=[0.20 0.22 0.23];
dispersionset=[16 0 3];
slopeset=[0.080 0.075 0.050];
Aset=[80 50 72];
j1 = 1:n(1);
S1 = zeros(1,length(j1));
if mod(n(1),2)==1
for k = 1:length(j1)
if mod(j1(k),2) == 0
S1(k) = (n(1)^2-1)/4 + n(1)*j1(k)/2 - j1(k)^2/2 - n(1) + j1(k)/2;
else
S1(k) = (n(1)^2+1)/4 + n(1)*j1(k)/2 - j1(k)^2/2 - n(1) + j1(k)/2;
end
end
else
for k = 1:length(j1)
S1(k) = n(1)^2/4 + n(1)*j1(k)/2 - j1(k)^2/2 - n(1) + j1(k)/2;
end
end
j2 = 1:n(2);
S2 = zeros(1,length(j2));
if mod(n(2),2)==1
for k = 1:length(j2)
if mod(j2(k),2) == 0
S2(k) = (n(2)^2-1)/4 + n(2)*j2(k)/2 - j2(k)^2/2 - n(2) + j2(k)/2;
else
S2(k) = (n(2)^2+1)/4 + n(2)*j2(k)/2 - j2(k)^2/2 - n(2) + j2(k)/2;
end
end
else
for k = 1:length(j2)
S2(k) = n(2)^2/4 + n(2)*j2(k)/2 - j2(k)^2/2 - n(2) + j2(k)/2;
end
end
j3 = 1:n(3);
S3 = zeros(1,length(j3));
if mod(n(3),2)==1
for k = 1:length(j3)
if mod(j3(k),2) == 0
S3(k) = (n(3)^2-1)/4 + n(3)*j3(k)/2 - j3(k)^2/2 - n(3) + j3(k)/2;
else
S3(k) = (n(3)^2+1)/4 + n(3)*j3(k)/2 - j3(k)^2/2 - n(3) + j3(k)/2;
end
end
else
for k = 1:length(j3)
S3(k) = n(3)^2/4 + n(3)*j3(k)/2 - j3(k)^2/2 - n(3) + j3(k)/2;
end
end
for type=1:3
alphadb=alphadbset(type);
dispersion=dispersionset(type);
slope=slopeset(type);
A=Aset(type);
alpha = alphadb/(10*log10(exp(1)));
s = slope*10^3;
d = 6;
lambda = 1.550*10^-6;
Aeff = A*10^-12;
c = 3*10^8;
n2 = 2.68*10^-20;
nsp=1.3;
h=6.626*10^-34;
B=12.5*10^9;
D = dispersion*10^-6;
delta_belta = 2 * pi * c * (delta_lambda^2) / (lambda^2)* (D+s*delta_lambda)*10^3;
L=30:1:100;
m=ceil(distance./L)-1;
Leff=(1-exp(-alpha.*L))/alpha;
eta = (alpha^2) / ((alpha^2) + (delta_belta^2))*(1+((4*exp(-alpha.*L).*(sin(delta_belta.*L/2)).^2)./((1-exp(-alpha.*L)).^2)));
gamma=(2*pi*n2)/(lambda*Aeff);
Pfwm= eta .*(d^2/9).*(gamma^2).*((Leff*10^3).^2).*(p^3).*exp(-alpha.*L);
Pfwmchannelmax10=max(S1).*(Pfwm.*((sin(m.*delta_belta.*L./2)).^2)./((sin(delta_belta.*L/2)).^2)+(eta .*(d^2/9).*(gamma^2).*((Leff*10^3).^2).*(p^3).*exp(-alpha.*(distance-m.*(distance-m.*L)))));
Pfwmchannelmax40=max(S2).*(Pfwm.*((sin(m.*delta_belta.*L./2)).^2)./((sin(delta_belta.*L/2)).^2)+(eta .*(d^2/9).*(gamma^2).*((Leff*10^3).^2).*(p^3).*exp(-alpha.*(distance-m.*(distance-m.*L)))));
Pfwmchannelmax80=max(S3).*(Pfwm.*((sin(m.*delta_belta.*L./2)).^2)./((sin(delta_belta.*L/2)).^2)+(eta .*(d^2/9).*(gamma^2).*((Leff*10^3).^2).*(p^3).*exp(-alpha.*(distance-m.*(distance-m.*L)))));
%%Calculation ASE noise:
G=exp(alpha.*L);
Pase=2*nsp*h*c/lambda*B .*(G-1).*m;
Pasedbm=30+(10*log(Pase));
%%Optical signal-to-noise ratio per channel:
Pnoisemax10=Pfwmchannelmax10+Pase;
Pnoisemax40=Pfwmchannelmax40+Pase;
Pnoisemax80=Pfwmchannelmax80+Pase;
Pout=p*exp(-alpha.*(distance-m.*L));
OSNR10=10*log(Pout./Pnoisemax10);
OSNR40=10*log(Pout./Pnoisemax40);
OSNR80=10*log(Pout./Pnoisemax80);
figure(type)
FigHandle = figure(type);
set(FigHandle, 'Position', [400, 350, 1400, 750]);
subplot(3,1,1);
plot(L,OSNR10,'m')
title('Optical signal-to-noise ratio for 10 channels','Color','b');
xlabel('Amplifier spacing (km)','Color','b');
ylabel('OSNR (dB)','Color','b');
grid
subplot(3,1,2);
plot(L,OSNR40,'r')
title('Optical signal-to-noise ratio for 40 channels','Color','b');
xlabel('Amplifier spacing (km)','Color','b');
ylabel('OSNR (dB)','Color','b');
grid
subplot(3,1,3);
plot(L,OSNR80,'k')
title('Optical signal-to-noise ratio for 80 channels','Color','b');
xlabel('Amplifier spacing (km)','Color','b');
ylabel('OSNR (dB)','Color','b');
grid
end
  1 Comment
Geoff Hayes
Geoff Hayes on 28 Aug 2014
son - I ran your above code, and a title was created for each subplot. What exactly do you mean by don't know how to create one sub title using for loop? Do you want a title for each figure (this is missing), as well as one title for each subplot (this works)?

Sign in to comment.

Sign in to answer this question.

Answers (1)

Image Analyst
Image Analyst on 28 Aug 2014
  2 Comments
son
son on 29 Aug 2014
How to use this file?
Image Analyst
Image Analyst on 29 Aug 2014
Download the file to a folder on your search path, and then call it like it probably tells you in the instructions. I haven't run it in a log time, if ever.

Sign in to comment.

Categories

Find more on Subplots in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!