[讨论] 高斯讯号乘上cos讯号做傅立叶转换

楼主: heizesun (RedSun)   2019-09-05 17:05:17
各位高手好,
小妹目前还是程式新手,
有些问题想问版上的高手们。
我有一个cos讯号要和高斯讯号相乘后做傅立叶转换
(cos2pi*fc*t)*((1/根号2)*exp(-t^2/2a^2))
fc=10^6Hz
a为常数
1.请问该如何设定dt的值?
2.如何确认a的值?
3.我做出来的图形正确吗?一直觉得哪里怪怪的.
请大家帮帮忙 谢谢!
以下是我的程式码
dt = 0.0000001;
fs = 1/dt;
fn = fs/2;
n = 100;
t = dt*(-n/2:n/2);
a = 0.0001; aa = a^2;
ft = 1/(sqrt(2*pi))*(exp(-t.^2/(2*aa)));
fc = 1000000;
xt = cos(2*pi*t*fc) .* ft;
subplot(2,1,1); plot(t,xt,'b');
title('Gaussian Pulse');
xlabel('Time(s)'); ylabel('Amplitude');
xf = fftshift(fft(xt));
f = fs*(-n/2:n/2)/n; %Frequency Vector
subplot(2,1,2); plot(f,abs(xf),'r');
title('Magnitude of FFT');
xlabel('Frequency (Hz)'); ylabel('Magnitude |X(f)|');
楼主: heizesun (RedSun)   2019-09-05 17:14:00
更正ft = 1/(sqrt(2))*(exp(-t.^2/(2*aa)));

Links booklink

Contact Us: admin [ a t ] ucptt.com