Я пытаюсь создать импульс Гаусса, который соответствует стандартам F CC в Matlab, это формула, которую я сейчас использую. Я не могу получить меньшее значение сигмы, как это должно быть в течение наносекунд, может кто-нибудь объяснить, что я делаю не так с генерацией гауссовского импульса?
Вот гауссов импульс:
dt=.0001;
n=1000;
t=dt*(-n/2:n/2); %time base
Sigma=0.001;
alpha=4*pi*Sigma;
A = 1*10^-20;
Gaussian = -(exp(-t.^2/(2*Sigma^2)))/(Sigma*sqrt(2*pi));
He5 = -(t.^5 - (10.*t.^3)*Sigma.^2 + 15.*t.*Sigma.^4)./Sigma.^10;
x5 = A*He5.*Gaussian;
subplot(2,1,1);
plot(t,x5,'DisplayName','5th Derivative');
legend;
title(['Gaussian Pulse \sigma=', num2str(Sigma),'s','\alpha=',num2str(alpha), 's']);
xlabel('Time(s)'); ylabel('Amplitude');
А вот PSD по сравнению с правилами F CC
subplot(2,1,2);
hold on;
[Pxx,g5] = periodogram(x5);
plot(g5,10*log10(Pxx),'DisplayName','PSD');
xplot = [-41.3 -41.3 -41.3 -75.3 -75.3 -53.3 -53.3 -51.3 -51.3 -41.3 -41.3 -41.3 -51.3 -51.3 -51.3];
yplot = [0 0 0.96 0.96 1.61 1.61 1.99 1.99 3.1 3.1 10.6 10.6 10.6 15 15];
plot(yplot,xplot,'DisplayName','FCC');
grid on;
legend
title('PSD');
xlabel('Frequency(GHz)'); ylabel('UWB EIRP(dBm/MHz)');
Большое спасибо.