Clear Filters
Clear Filters

How to get to Power Spectral Density from Power spectrum as shown in figure

41 views (last 30 days)
fsine=10e6;
L=10;
PS=0.2;
A=7.85e-11;
c=3e8;
dt=2e-12;
V=1;
n=1.45; %Index of refraction
eps0=8.854e-12; % [F/m] Vacuum permittivity
T=10*2*L*n/c; %Total time
Nt=round(T/dt);
t = (-T/2/dt:1:T/2/dt)*dt; %time axis
nu=(-1/2/dt:1/T:1/2/dt); %frequency axis
I1_0=PS/A;
sine = V*sin(2*pi*fsine*t); %phase modulation condition
phi=sine;
ES_0t=sqrt(I1_0/2/n/c/eps0)*exp((1i*(pi)*phi)); % Original signal in time
Power=trapz(t,2*n*c*eps0*A*abs(ES_0t).^2)/T; % Area under the curve in time domain
FFt_EL0t=fftshift(abs(fft(ES_0t/Nt))); % fourier transform of the original signal
Power_FFt=T*trapz(nu,2*n*c*eps0*A*(FFt_EL0t).^2); % Area under the curve in frequency domain
figure;plot(nu,10*log10(2*n*c*eps0*A*(FFt_EL0t).^2)./1e-3);
xlim([-100e6 100e6]);
ylabel('Power[dBm]');
xlabel('frequency[Hz]');
The idea was to obtain PSD exactly like the second figure but I cannot understand what to change in the above code.I have converted Power units from W to dBm but evidently that is of no use.
Any suggestions will be greatly appreciated.
Thank you

Answers (1)

dpb
dpb on 19 Sep 2024 at 20:14
Edited: dpb on 19 Sep 2024 at 21:51
fsine=10e6;
L=10;
PS=0.2;
A=7.85e-11;
c=3e8;
dt=2e-12;
V=1;
n=1.45; %Index of refraction
eps0=8.854e-12; % [F/m] Vacuum permittivity
T=10*2*L*n/c; %Total time
Nt=round(T/dt);
t = (-T/2/dt:1:T/2/dt)*dt; %time axis
nu=(-1/2/dt:1/T:1/2/dt); %frequency axis
I1_0=PS/A;
sine = V*sin(2*pi*fsine*t); %phase modulation condition
phi=sine;
ES_0t=sqrt(I1_0/2/n/c/eps0)*exp((1i*(pi)*phi)); % Original signal in time
Power=trapz(t,2*n*c*eps0*A*abs(ES_0t).^2)/T; % Area under the curve in time domain
FFt_EL0t=fftshift(abs(fft(ES_0t/Nt))); % fourier transform of the original signal
Power_FFt=T*trapz(nu,2*n*c*eps0*A*(FFt_EL0t).^2); % Area under the curve in frequency domain
figure;
semilogy(nu/1E6,2*n*c*eps0*A*(FFt_EL0t).^2);
xlim([-1 1]*1E2)
ylabel('Amplitude');
xlabel('frequency [MHz]');
will scale to log y axis. Nothing can do about the shape of the spectrum; that's wholly dependent upon the characteristics of the time trace.

Categories

Find more on Fourier Analysis and Filtering 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!