%使用1KHz正弦信号调制20KHz载波,绘制已调信号时域波形与频谱
%采用DSB-SC,另根据Nuquist低通采样定理,选定采样频率为150KHz
T_start=0;%开始时间
T_stop=1e-2;%结束时间
T=T_stop-T_start;%持续时间
fs1=150e3;%调制信号采样频率
T_sample=1/fs1;%调制信号采样间隔
N_sample=T/T_sample;%调制信号采样点数
n=0:N_sample-1;
fm=1e3;%调制信号频率
fn=sin(2pifmnT_sample);%调制信号
fc=20e3;%载波频率
cn=cos(2pifcnT_sample);%载波信号
sn=fn.cn;%DSB-SC调制
Nd_sample=N_sample;
nd=n;
fsd=fs1;
SN_O=abs(fft(sn))/N_sample;
SN(1:N_sample/2)=SN_O(N_sample/2+1:N_sample);
SN(N_sample/2+1:N_sample)=SN_O(1:N_sample/2);
subplot(2,1,1)
plot(nT_sample,sn,'b');
title('已调信号波形')
axis([0 1e-2 -2 2])
hold on
plot(nT_sample,fn,'y--','linewidth',3);%画出DSB信号的包络
subplot(2,1,2)
plot((n-N_sample/2)/N_samplefs1,SN);
title('已调信号频谱')
axis([-30e3 30e3 0 0.3])
以上为代码
下面三段的意思是什么
SN_O=abs(fft(sn))/N_sample;
SN(1:N_sample/2)=SN_O(N_sample/2+1:N_sample);
SN(N_sample/2+1:N_sample)=SN_O(1:N_sample/2);
第一句,用FFT求幅度谱后面两句正负频率频谱互换