silas/imdd_silas
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

minor changes büro

Browse Source
This commit is contained in:
Silas Oettinghaus
2023-06-07 13:50:56 +02:00
parent b640c013cb
commit 2868887a15
5 changed files with 186 additions and 84 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
Classes/00_signals/Signal.m
View File

@@ -168,19 +168,22 @@ classdef Signal
fig = findall(groot, 'Type', 'figure', 'Name', 'power density'); fig = findall(groot, 'Type', 'figure', 'Name', 'power density');
if isvalid(fig) if isvalid(fig)
fig = get(fig); fig = get(fig);
ax = gca;
hold on hold on
else else
figure('name','power density'); figure('name','power density');
ax = gca;
end end
else else
fig = findall(groot, 'Type', 'figure', 'Name', options.figurename); fig = findall(groot, 'Type', 'figure', 'Name', options.figurename);
if isvalid(fig) if isvalid(fig)
fig = get(fig); ax = fig.CurrentAxes;
hold on hold on
else else
figure('name',options.figurename); figure('name',options.figurename);
ax = gca;
hold on
end end
end end
@@ -200,10 +203,12 @@ classdef Signal
psd = psd/length(Fsignal); psd = psd/length(Fsignal);
%smoothing %smoothing
psd = smooth(psd,100); psd = smooth(psd,50);
psd_plot = 20*log10(psd); psd_plot = 20*log10(psd);
psd_plot(psd_plot<-120) = -120;
testParseval = 1; testParseval = 1;
if testParseval == 1 if testParseval == 1
@@ -225,9 +230,9 @@ classdef Signal
if ~isempty(options.displayname) if ~isempty(options.displayname)
plot(freq_vec*1e-9,psd_plot,'Linewidth',0.5,'DisplayName',options.displayname); plot(freq_vec*1e-9,psd_plot,'Linewidth',0.5,'DisplayName',options.displayname,'Parent',ax);
else else
plot(freq_vec*1e-9,psd_plot,'Linewidth',0.5); plot(freq_vec*1e-9,psd_plot,'Linewidth',0.5,'Parent',ax);
end end
xlabel('Frequency [GHz]') xlabel('Frequency [GHz]')
@@ -235,9 +240,9 @@ classdef Signal
%Wavelength Axis %Wavelength Axis
freq_vec = physconst('LightSpeed')*linspace(-fsamp/2,fsamp/2,length(psd))./((physconst('LightSpeed')/1550e-9)^2); freq_vec = physconst('LightSpeed')*linspace(-fsamp/2,fsamp/2,length(psd))./((physconst('LightSpeed')/1550e-9)^2);
if ~isempty(options.displayname) if ~isempty(options.displayname)
plot(freq_vec*1e9,psd_plot,'Linewidth',0.5,'DisplayName',options.displayname); plot(freq_vec*1e9,psd_plot,'Linewidth',0.5,'DisplayName',options.displayname,'Parent',ax);
else else
plot(freq_vec*1e9,psd_plot,'Linewidth',0.5); plot(freq_vec*1e9,psd_plot,'Linewidth',0.5,'Parent',ax);
end end
xlabel('Wavelength [nm]') xlabel('Wavelength [nm]')

6
Classes/02_etc/Filter.m
View File

@@ -175,14 +175,14 @@ classdef Filter
% hold on % hold on
% plot(w/(2*pi)*obj.fsamp*1e-9,20*log10(abs(H)),'LineWidth',3); % plot(w/(2*pi)*obj.fsamp*1e-9,20*log10(abs(H)),'LineWidth',3);
% ax = gca; % ax = gca;
% ylim([-6 0]) % %ylim([-6 0])
% grid on % grid on
% xlabel('Freq in GHz') % xlabel('Freq in GHz')
% ylabel('Magnitude (dB)') % ylabel('Magnitude (dB)')
% freqz(B, A) % freqz(B, A, 2048, obj.fsamp);
%hfvt = fvtool(B,A); % hfvt = fvtool(B,A);
end end
end end

1
Classes/03_receive/Photodiode.m
View File

@@ -56,6 +56,7 @@ classdef Photodiode
% Magnitude squared detection % Magnitude squared detection
yout = sum( abs(xin) .^2*obj.responsivity, 2 ) ; yout = sum( abs(xin) .^2*obj.responsivity, 2 ) ;
% Shot Noise % Shot Noise
shot_noise = sqrt(k * obj.fsimu .* yout) .* randn(size(yout,1),1); shot_noise = sqrt(k * obj.fsimu .* yout) .* randn(size(yout,1),1);

70
filter_debug.m Normal file
View File

@@ -0,0 +1,70 @@
clear
%% Set Simulation Variables
O = 18; %order of prbs
N = 2^(O-1); %length of prbs
[~,seed] = prbs(O,1); %initialize first seed of prbs
% Modulation
M = 4; %PAM-M
bitpattern = zeros(N,log2(M));
% Symbol Rate
fsym = 112e9;
% DAC Rate
fdac = 120e9;
% Simulation oversampling rate "k";
kover = 16;
% ADC Rate
fadc = 256e9;
% Simulation frequency in "analog domain"
fsimu = kover * fdac ;
%% CONSTRUCT ALL CLASSES
digimod = PAMmapper(M,0);
pulsef = Pulseformer("pulseform","rrc","fdac",fdac,"fsym",fsym,"pulselength",32,"rrcalpha",0.05);
awg = AWG('preset','M8199B','fdac',fdac,'kover',kover,'lpf_active',1,'f_cutoff',56e9,'lpf_type',filtertypes.gaussian,'bit_resolution',5.5);
lp_laser = Filter('filtdegree',1,"f_cutoff",30e9,"fsamp",fdac,"filterType",filtertypes.bessel_inp);
%% PROCESS
% PRBS Generation
for i = 1:log2(M)
[bitpattern(:,i),seed] = prbs(O,N,seed);
end
% Build Inf. signal class
bits = Informationsignal(bitpattern);
% Digi Mod
mod_out = digimod.map(bits);
% merken für EQ training
reference = mod_out;
% shape shape
mod_out = pulsef.process(mod_out);
test = applyPulseShaping(reference.signal,fsym,fdac);
% AWG -> ELECTRICAL DOMAIN
X = lp_laser.process(mod_out);
%X.spectrum(fsimu,"displayname",'AWG out','figurename','after AWG');

162
setup_simulation.m
View File

@@ -1,8 +1,10 @@
clear
rng(2023);
%% Set Simulation Variables %% Set Simulation Variables
O = 18; %order of prbs O = 17; %order of prbs
N = 2^(O-1); %length of prbs N = 2^(O-1); %length of prbs
[~,seed] = prbs(O,1); %initialize first seed of prbs [~,seed] = prbs(O,1); %initialize first seed of prbs
@@ -26,7 +28,6 @@ fadc = 256e9;
fsimu = kover * fdac ; fsimu = kover * fdac ;
%% CONSTRUCT ALL CLASSES %% CONSTRUCT ALL CLASSES
digimod = PAMmapper(M,0); digimod = PAMmapper(M,0);
@@ -39,13 +40,13 @@ lp_laser = Filter('filtdegree',1,"f_cutoff",60e9,"fsamp",fdac,"filterType",fi
u_pi = 4.6; u_pi = 4.6;
vbias = (0.5*u_pi)-u_pi; vbias = (0.5*u_pi)-u_pi;
extmodlaser = EML("mode",eml_mode.im_cosinus,"power",5,"fsimu",fsimu,"lambda",1550,"bias",vbias,"u_pi",u_pi,"linewidth",10000); extmodlaser = EML("mode",eml_mode.im_cosinus,"power",15,"fsimu",fsimu,"lambda",1550,"bias",vbias,"u_pi",u_pi,"linewidth",10000);
amp = Amplifier("amp_mode","ideal_no_noise","amplification_db",0,"gain_mode","output_power"); amp = Amplifier("amp_mode","ideal_no_noise","amplification_db",0,"gain_mode","output_power");
fib = Fiber("fsimu",fdac*kover,"fiber_length",0,"alpha",0.2,"D",17,"lambda0",1550,"gamma",0); fib = Fiber("fsimu",fdac*kover,"fiber_length",1,"alpha",0.2,"D",17,"lambda0",1550,"gamma",0);
optatten = Amplifier("amp_mode","ideal_no_noise","gain_mode","output_power","amplification_db",5); optatten = Amplifier("amp_mode","ideal_no_noise","gain_mode","output_power","amplification_db",0);
edfaamp = Amplifier("amp_mode","edfa_increase_nase","gain_mode","output_power","amplification_db",0,"nase_mode","generate_ase","noifig",5); edfaamp = Amplifier("amp_mode","edfa_increase_nase","gain_mode","output_power","amplification_db",0,"nase_mode","generate_ase","noifig",5);
@@ -54,82 +55,107 @@ phdiode = Photodiode("fsimu",fdac*kover,"dark_current",2e-08,"responsivity",
lp_diode = Filter('filtdegree',1,"f_cutoff",80e9,"fsamp",fdac,"filterType",filtertypes.bessel_inp); lp_diode = Filter('filtdegree',1,"f_cutoff",80e9,"fsamp",fdac,"filterType",filtertypes.bessel_inp);
scp = Scope("fsimu",fdac*kover,"fadc",fadc,... scp = Scope("fsimu",fdac*kover,"fadc",fadc,...
"delay",0,"fixed_delay",0,"lpf_bw",113e9,"filtertype",filtertypes.butterworth,... "delay",0,"fixed_delay",0,"lpf_bw",113e9,"filtertype",filtertypes.butterworth,...
"samplingdelay",0,"rand_samplingdelay",0,"freq_offset",0,"samp_jitter",0,... "samplingdelay",0,"rand_samplingdelay",0,"freq_offset",0,"samp_jitter",0,...
"adcresolution",6,"quantbuffer",0.1); "adcresolution",6,"quantbuffer",0.1);
eq = EQ("K",2,"plottrain",0,"plotfinal",1,...
"training_length",4096,"training_loops",5,...
"Ne",[50,0,0],"Nb",[0,0,0],...
"DCmu",0.005,"DDmu",[0.0004 0.0004 0.0004 0.0004 ],"DFEmu",0.005,"FFEmu",0.000,...
"dd_loops",2,"epsilon",[10 100 1000 ],"M",4,...
"thres",[0.005 0.004 0.0005 ],"l1act",0,"delay",1,"rho",0.0005,"ideal_dfe",0,"DB_aim",0);
eq = EQ("K",2,"plottrain",0,"plotfinal",0,...
"training_length",4096,"training_loops",5,...
"Ne",[30,0,0],"Nb",[0,0,0],...
"DCmu",0.005,"DDmu",[0.0004 0.0004 0.0004 0.0004 ],"DFEmu",0.002,"FFEmu",0.00,...
"dd_loops",2,"epsilon",[10 100 1000 ],"M",4,...
"thres",[0.005 0.004 0.0005 ],"l1act",0,"delay",0,"rho",0.0005,"ideal_dfe",0,"DB_aim",0);
%% PROCESS %% PROCESS
% PRBS Generation output_pow = 0
for i = 1:log2(M) for lp1 = 1:length(output_pow)
[bitpattern(:,i),seed] = prbs(O,N,seed);
% change parameters
optatten.amplification_db = output_pow(lp1);
% PRBS Generation
for i = 1:log2(M)
[bitpattern(:,i),seed] = prbs(O,N,seed);
end
bitpattern = prms_out';
% Build Inf. signal class
bits = Informationsignal(bitpattern);
% Digi Mod
mod_out = digimod.map(bits);
% merken für EQ training
reference = mod_out;
% shape shape
mod_out = pulsef.process(mod_out);
test = applyPulseShaping(reference.signal,fsym,fdac);
% AWG -> ELECTRICAL DOMAIN
X = awg.process(mod_out);
X = lp_laser.process(X);
X = X.normalize;
% Laser; Modulation -> OPTICAL DOMAIN
X = extmodlaser.process(X);
%X.spectrum(fsimu,"displayname",'laser out','figurename','after laser');
%X = amp.process(X);
% Fiber Propagation
X = fib.process(X);
X = optatten.process(X);
X = edfaamp.process(X);
% Photo Diode -> ELECTRICAL DOMAIN
X = phdiode.process(X);
X = lp_diode.process(X);
%X.spectrum(fsimu,"displayname",'diode out','figurename','after diode');
% Oscilloscope (Sampling to f_adc; Quantization; Bandwidth Limitation)
X = scp.process(X);
% Resample to Symbol Rate
X = X.resample("fs_out",2*fsym,"fs_in",fadc);
% INFORMATION SIGNAL
X = X.normalize;
% Equalizer
eq_out = eq.process(X,reference);
rx_series = eq_out.signal;
% Digi Demod
demap_out = digimod.demap(eq_out);
% BER
[bits,errors,BER] = calc_ber(demap_out.signal(:,1:end-1),bitpattern(1:end,:)',0);
disp(['BER: ', sprintf('%2E',BER), ' ERRORS: ' ,num2str(sum(errors))]);
bercurve(lp1) = BER;
errorscurve(lp1) = sum(errors);
end end
% Build Inf. signal class
bits = Informationsignal(bitpattern);
% Digi Mod
mod_out = digimod.map(bits);
% merken für EQ training
reference = mod_out;
% shape shape
mod_out = pulsef.process(mod_out);
test = applyPulseShaping(reference.signal,fsym,fdac);
% AWG -> ELECTRICAL DOMAIN
X = awg.process(mod_out);
X = lp_laser.process(X);
X = X.normalize;
% Laser; Modulation -> OPTICAL DOMAIN
X = extmodlaser.process(X);
X.spectrum(fsimu,"displayname",'laser out','figurename','after laser');
%X = amp.process(X);
% Fiber Propagation
X = fib.process(X);
X = optatten.process(X);
X = edfaamp.process(X);
% Photo Diode -> ELECTRICAL DOMAIN
X = phdiode.process(X);
X = lp_diode.process(X);
X.spectrum(fsimu,"displayname",'diode out','figurename','after diode');
% Oscilloscope (Sampling to f_adc; Quantization; Bandwidth Limitation)
X = scp.process(X);
% Resample to Symbol Rate
X = X.resample("fs_out",2*fsym,"fs_in",fadc);
% INFORMATION SIGNAL
X = X.normalize;
% Equalizer
eq_out = eq.process(X,reference);
rx_series = eq_out.signal;
% Digi Demod
demap_out = digimod.demap(eq_out);
% BER
[bits,errors,BER] = calc_ber(demap_out.signal(:,1:end-1),bitpattern(1:end,:)',0);
disp(['BER: ', sprintf('%2E',BER)]);
% disp(demap_out.logbook);
@@ -139,13 +165,13 @@ if 0
figure(22) figure(22)
clf clf
% subplot(3,1,1) % subplot(3,1,1)
% hold on % hold on
% plot(reference.signal(10000:end-20),'DisplayName','Tx','Color',col(1,:),'LineWidth',3); % plot(reference.signal(10000:end-20),'DisplayName','Tx','Color',col(1,:),'LineWidth',3);
% plot(rx_series(10000:end-20),'DisplayName','Rx','Color',col(6,:)); % plot(rx_series(10000:end-20),'DisplayName','Rx','Color',col(6,:));
% title('Modulated Sequence Tx - Rx'); % title('Modulated Sequence Tx - Rx');
% legend % legend
% hold off % hold off
subplot(2,1,1) subplot(2,1,1)
hold on hold on
plot(reference.signal(4150:4175),'DisplayName','Tx','Color',col(1,:),'LineWidth',3); plot(reference.signal(4150:4175),'DisplayName','Tx','Color',col(1,:),'LineWidth',3);