画出频谱、功率谱密度图。
dt=0.001;
fmax=1;
fc=10;
T=5;
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fmax*t);
A=0;
s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t));
[f,Xf]=FFT_SHIFT(t,s_ssb);
PSD=(abs(Xf).^2)/T;
figure(1)
subplot(211);
plot(t,s_ssb);hold on ;
title('SSB 调制信号');
subplot(212);
plot(f,PSD);
axis([-2*fc 2*fc 0 1.5*max(PSD)]);
title('SSB 信号功率谱');
xlabel('f');
xlabel('f');
00.20.40.60.81 1.2 1.4 1.6 1.82
-2-1
1
2
SSB 调制信号
-20-15-10-50
51015200
1
2
3
SSB 信号功率谱
f
画出频谱、功率谱密度图。
dt=0.001; %时间采样频谱
fmax=1; %信源最高频谱
fc=10; %载波中心频率
T=5; %信号时长
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fmax*t); %信源
A=0;
s_dsb=mt.*cos(2*pi*fc*t);
[f,Xf]=FFT_SHIFT(t,s_dsb); %调制信号频谱
PSD=(abs(Xf).^2)/T; %调制信号功率谱密度
figure(1)
subplot(211);
plot(t,s_dsb);hold on; %画出DSB 信号波形
plot(t,A+mt,'r--'); %表示DSB 包络
plot(t,-A-mt,'r--');
title('DSB 调制信号及其包络');
xlabel('t');
subplot(212); %画出功率谱图形
plot(f,PSD);
axis([-2*fc 2*fc 0 1.5*max(PSD)]);
title('DSB 信号功率谱');
xlabel('f');
xlabel('f');
00.51 1.52 2.53
3.54
4.55
-2-1
1
2
DSB 调制信号
及其包络
t
-20-15-10-50
51015200
0.2
0.4
0.6
0.8
DSB 信号功率谱f
画出频谱、功率谱密度图。
程序：dt=0.001
T=5
N=T/dt
T=[0:N-1]*dt
mt=sqrt(2)*cos(2*pi*fmax*t)
s_am=10-cos(10*t(2*pi)*(5/0.3)*sin(0.3*t)) [f,Xf]=fft_shift(t,s_am)
PSD=(abs(Xf).^2)/T
figure(1)
Subplot(211)
Plot(t,s_am)
hold on
Plot(t,mt,’r--’)
title(‘FM调制信号及包络’)
xlabel(‘t’)
subplot(212)
plot(f,PSD)
axis([-2*fc 2*fc 0 1.5*max(PSD)])
tltle(‘FM信号功率谱’)
xlabel(‘f’)
