6716914
Newbie level 5
%% Simulate the DSM in time
Nfft = 2^13;
tone_bin = 31; what is tone_bin?
t = 0:Nfft-1;
u = 1*(nLev-1)*sin(2*pi*tone_bin/Nfft*t); why the input signal has to be this, what this means
v = simulateDSM(u,Ha,nLev);
figure()
n=1:1000;
plot(t, u, 'g');
hold on; grid on;
stairs(t, v, 'b');
title('DSM time-domain Simulation');
%% Plot the Spectrum
spec = fft(v.*ds_hann(Nfft))/(Nfft*(nLev-1)/4); what this means
snr = calculateSNR(spec(1:ceil(Nfft/(2*OSR))+1), tone_bin);
Neff = (snr-1.76)/6.02;
NBW = 1.5/Nfft; why 1.5
f = linspace(0, 0.5, Nfft/2+1);
Sqq = 4*(evalTF(Ha, exp(2i*pi*f))/(nLev-1)).^2/3;
figure()
plot(f, dbv(spec(1:Nfft/2+1)), 'b')
hold on; grid on;
plot(f, dbp(Sqq*NBW), 'm', 'Linewidth', 1);
title('DSM Output Spectrum');
Nfft = 2^13;
tone_bin = 31; what is tone_bin?
t = 0:Nfft-1;
u = 1*(nLev-1)*sin(2*pi*tone_bin/Nfft*t); why the input signal has to be this, what this means
v = simulateDSM(u,Ha,nLev);
figure()
n=1:1000;
plot(t, u, 'g');
hold on; grid on;
stairs(t, v, 'b');
title('DSM time-domain Simulation');
%% Plot the Spectrum
spec = fft(v.*ds_hann(Nfft))/(Nfft*(nLev-1)/4); what this means
snr = calculateSNR(spec(1:ceil(Nfft/(2*OSR))+1), tone_bin);
Neff = (snr-1.76)/6.02;
NBW = 1.5/Nfft; why 1.5
f = linspace(0, 0.5, Nfft/2+1);
Sqq = 4*(evalTF(Ha, exp(2i*pi*f))/(nLev-1)).^2/3;
figure()
plot(f, dbv(spec(1:Nfft/2+1)), 'b')
hold on; grid on;
plot(f, dbp(Sqq*NBW), 'm', 'Linewidth', 1);
title('DSM Output Spectrum');