Question

matlab 画频谱图

dt=1/100000;
t=0:dt:1;
ac=4;%载波幅度

fc=1700;%载波频率km=1000;%调幅指数
mt=cos(2*pi*20*t);%调制信号

ct=ac*cos(2*pi*fc*t);%载波信号
st=ac*cos(2*pi*1700*t+2*pi*km*cumsum(mt*dt));%已调信号
snt=awgn(st,20);%模拟加性高斯白噪声信道传输

求画snt的频谱！请把程序写在下面！

Answer 1

1、假设该信号的域被借出并且向量t是N维向量，则该信号离散化的采样周期是Ts = 1 / fs = lent /（N-1），其中fs是采样频率。

2、从上面的离散傅立叶公式，我们可以知道在使用fft函数之后，我们仍然得到一个N维向量。

3、频域的长度定义为lenf，其满足lent =（N-1）* Ts，solid lenf =（（N-1）*（N-1）/ N）/ lent = fs *（N-1） ）/ N.当N很大时，lenf = fs。

4、到目前为止在频域中获得的矢量是f =（0：N-1）/（N-1）* lenf =（0：N-1）* fs / N.根据奈奎斯特采样定理，采样频率高于信号中最高频率的两倍，因此不会失真。

5、因此，频谱中的（fs / 2，fs）是无用的频谱，因此域矢量和频域信号的范围矢量应该减半。即：f = f（1，N / 2）; y = abs（fft（x））; y = y（1. N / 2）;好吧，我们可以绘制信号的时域和频域图像。

Answer 2

fft就可以画频谱，工程上画图习惯还要取对数，完整如下：
fsnt = fftshift(fft(snt));
pfsnt = log10(abs(fsnt));
plot(pfsnt);

追问

横坐标的单位是啥？画出来非常复杂呢。感觉不太对呢

追答

这是很典型的调幅信号加上噪声的频谱，您觉得那里不对呢？
横坐标的单位，我只提醒一下，数字频率范围是[-pi, pi]; 数字频率2*pi对应于模拟频率中的的采样率。既然您是做通信的应该很清楚这些了。