UWB帯域のガウスモノパルスのパワースペクトル密度(PSD)を推定します。私はfftを使用して2つのコードを持っています。しかし、その寸法がわからないため、y軸に問題があります(dBm / MHzのPSDが必要です)。また、y軸に1つの値しか表示されないため、最初のコードでは間違いであるはずです。
Code1
fs=1e11;
g=0.1e-9;
tmax=1e-9;
fftl=2048;
t=(-tmax:1./fs:tmax)';
s=t./(g.^3.*sqrt(2.*pi)).*exp(-t.^2./(2.*g.^2))./2.5e19;
figure(1)
plot(t,s);
xlabel('Time, s');
ylabel('Amplitude, V');
ffts=abs(fft(s,fftl));
ffts=2.*ffts./fftl;
fftp=abs(ffts.*conj(ffts))./2;
fftps=(fftp-30)./1e-6;
f=0:fs./fftl:fs./2-fs./fftl;
figure(2)
plot(f,fftps(1:length(f))),grid;
xlabel('Frequency, Hz');
Code2
fs=1e11;
g=0.1e-9;
tmax=1e-9;
t=(-tmax:1./fs:tmax)';
s=t./(g.^3.*sqrt(2.*pi)).*exp(-t.^2./(2.*g.^2))./2.5e19;
figure(1)
plot(t,s);
xlabel('Time, s');
ylabel('Amplitude, V');
S=fft(s,8192);
f=fs.*(0:4095)./8192;
Pss=S.*conj(S)./8192;
figure(2)
plot(f,Pss(1:4096));
助けてくれてありがとう!