通信原理Matlab仿真：模拟信号的幅度调制
【摘要】
用Matlab产生一个频率为1Hz、功率为1的余弦信源，设载波频率为10Hz，试画出：
DSB调制信号及其功率谱密度；
clc,clear,close all;dt=0.001; t=-10:dt:10; % 仿真时间fm=1;%基带信号频率fc=10;%载波频率mt=sqrt(2)*cos(2*pi*fm*t);%基带...
用Matlab产生一个频率为1Hz、功率为1的余弦信源，设载波频率为10Hz，试画出：
DSB调制信号及其功率谱密度；
-
clc,clear,close all;
-
dt=0.001;
-
t=-10:dt:10; % 仿真时间
-
fm=1;%基带信号频率
-
fc=10;%载波频率
-
mt=sqrt(2)*cos(2*pi*fm*t);%基带信号
-
c=cos(2*pi*fc*t);%载波信号
-
st=mt.*c;%DSB
-
%求功率谱
-
T=20;
-
df=1/T;% 抽样精度
-
N=length(st);%确定抽样个数
-
f=-N/2*df:df:N/2*df-df;%功率谱范围
-
sf=fft(st);%傅里叶变化
-
sf=T/N*fftshift(sf);%将零频点移到频谱的中间
-
psf=(abs(sf).^2)/T;%求功率
-
%绘制基带信号
-
subplot(311)
-
plot(t,mt)
-
xlabel('t');
-
ylabel('mt');
-
%plot(t,mt,'r--')
-
title('基带信号波形');
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%绘制DSB波形
-
subplot(312)
-
plot(t,st)
-
xlabel('t');
-
ylabel('st');
-
hold on;
-
plot(t,mt,'r--');
-
title('DSB调制波形');
-
legend('DSB','基带信号')
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%绘制DSB功率谱
-
subplot(313)
-
plot(f,psf);
-
xlabel('f');
-
ylabel('psf');
-
%plot(f,abs(st1));
-
axis([-2*fc 2*fc 0 max(psf)]);
-
title('DSB信号功率谱');
-
grid on;
-
gtext('图2-1 DSB调制信号及其功率谱密度波形')
- 将已调信号解调，在时域内将解调后的波形与原信号进行对比；
-
clc,clear,close all;
-
dt=0.001;
-
t=-10:dt:10;% 仿真时间
-
fm=1;%基带信号频率
-
fc=10;%载波频率
-
w_lpf = 1; % 低通滤波器截止频率（理想）
-
mt=sqrt(2)*cos(2*pi*fm*t);%基带信号
-
c=cos(2*pi*fc*t);%载波信号
-
st=mt.*c;%DSB
-
%解调
-
st1=st.*c;
-
sample = length(t);
-
st2 = fft(st1, sample)/sample*2; % 这里要/抽样次数*2 才是最后结果
-
f = 1/dt/sample*(0:1:sample-1);
-
st3 = st2.*(f<w_lpf); % 进过低通滤波
-
st4 = real(ifft(st3)); %解调后波形
-
%绘制基带信号
-
subplot(311)
-
plot(t,mt)
-
xlabel('t');
-
ylabel('mt');
-
%plot(t,mt,'r--')
-
title('基带信号波形');
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%绘制DSB波形
-
subplot(312)
-
plot(t,st)
-
xlabel('t');
-
ylabel('st');
-
hold on;
-
plot(t,mt,'r--');
-
title('DSB调制波形');
-
legend('DSB信号','基带信号')
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%绘制解调后信号与原信号波形对比
-
subplot(313)
-
plot(t,st4*sample);
-
hold on;
-
plot(t,mt)
-
xlabel('t');
-
ylabel('st1');
-
%plot(f,abs(st1));
-
axis([-5,5,-1.5,1.5]);
-
hold on;
-
%plot(t,mt,'r--');
-
title('解调后波形');
-
legend('解调后信号','基带信号')
-
grid on;
-
gtext('图2-2 时域内DSB解调后的波形与原信号进行对比')
A=2的AM调制信号及其功率谱密度；
-
clc,clear,close all;
-
dt=0.001;
-
t=-10:dt:10;% 仿真时间
-
A=2;
-
fm=1;%基带信号频率
-
fc=10;%载波频率
-
w_lpf = 1; % 低通滤波器截止频率（理想）
-
w_ipf = 1; % 带通滤波器截止频率（理想）
-
fm = 0.3; % 带通滤波器通带
-
mt=sqrt(2)*cos(2*pi*fm*t);%基带信号
-
c=cos(2*pi*fc*t);%载波信号
-
st=(mt+A).*c;%AM
-
%解调
-
st1=st.*c;
-
sample = length(t);
-
st2 = fft(st1, sample)/sample*2; % 这里要/抽样次数*2 才是最后结果
-
f = 1/dt/sample*(0:1:sample-1);
-
st3 = st2.*(f<w_lpf); % 进过低通滤波
-
st4 = real(ifft(st3)); %解调后波形
-
%求功率谱
-
T=20;
-
df=1/T;%抽样精度
-
N=length(st);%确定抽样个数
-
f=-N/2*df:df:N/2*df-df;%功率谱范围
-
sf=fft(st);%傅里叶变化
-
sf=T/N*fftshift(sf);%将零频点移到频谱的中间
-
psf=(abs(sf).^2)/T;%求功率
-
%绘制基带信号
-
subplot(411)
-
plot(t,mt)
-
xlabel('t');
-
ylabel('mt');
-
%plot(t,mt,'r--')
-
title('基带信号波形');
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%绘制AM波形
-
subplot(412)
-
plot(t,st)
-
xlabel('t');
-
ylabel('st');
-
hold on;
-
plot(t,mt+2,'r--');
-
title('AM调制波形');
-
legend('AM','基带信号')
-
axis([-5,5,-5,5]);
-
grid on;
-
%绘制AM功率谱
-
subplot(413)
-
plot(f,psf);
-
xlabel('f');
-
ylabel('psf');
-
%plot(f,abs(st1));
-
axis([-2*fc 2*fc 0 max(psf)]);
-
title('AM信号功率谱');
-
grid on;
-
%绘制解调后信号与原信号波形对比
-
subplot(414)
-
plot(t,st4*sample-2)
-
xlabel('t');
-
ylabel('st4');
-
%plot(f,abs(st1));
-
axis([-10,10,-3,3]);
-
hold on;
-
plot(t,mt,'r--');
-
title('解调后波形');
-
legend('解调后信号','基带信号')
-
grid on;
-
gtext('图2-3 A=2的AM调制信号及其功率谱密度波形')
SSB调制信号及其功率谱密度；
-
clc,clear,close all;
-
dt=0.001;
-
t=-10:dt:10;% 仿真时间
-
fm=1;%基带信号频率
-
fc=10;%载波频率
-
w_lpf = 1; % 低通滤波器截止频率（理想）
-
w_ipf = 1; % 带通滤波器截止频率（理想）
-
fm = 0.3; % 带通滤波器通带
-
mt=sqrt(2)*cos(2*pi*fm*t);%基带信号
-
c=cos(2*pi*fc*t);%载波信号
-
st=real(hilbert(mt).*exp(j*2*pi*fc*t));%SSB
-
%解调
-
st1=st.*c;
-
sample = length(t);
-
st2 = fft(st1, sample)/sample*2; % 这里要/抽样次数*2
-
f = 1/dt/sample*(0:1:sample-1);
-
st3 = st2.*(f<w_lpf); % 进过低通滤波
-
st4 = real(ifft(st3)); %解调后波形
-
%求功率谱
-
T=20;
-
df=1/T;%抽样精度
-
N=length(st);%确定抽样个数
-
f=-N/2*df:df:N/2*df-df;%功率谱范围
-
sf=fft(st);%傅里叶变化
-
sf=T/N*fftshift(sf);%将零频点移到频谱的中间
-
psf=(abs(sf).^2)/T;%求功率
-
%绘制基带信号
-
subplot(411)
-
plot(t,mt)
-
xlabel('t');
-
ylabel('mt');
-
%plot(t,mt,'r--')
-
title('基带信号波形');
-
axis([-5,5,-2,2]);
-
grid on;
-
%绘制SSB波形
-
subplot(412)
-
plot(t,st)
-
xlabel('t');
-
ylabel('st');
-
hold on;
-
plot(t,mt,'r--');
-
title('SSB调制波形');
-
legend('SSB','基带信号')
-
axis([-5,5,-2,2]);
-
grid on;
-
%绘制SSB功率谱
-
subplot(413)
-
plot(f,psf);
-
xlabel('f');
-
ylabel('psf');
-
%plot(f,abs(st1));
-
axis([-2*fc 2*fc 0 max(psf)]);
-
title('SSB信号功率谱');
-
grid on;
-
%绘制解调后信号与原信号波形对比
-
subplot(414)
-
plot(t,st4*sample)
-
xlabel('t');
-
ylabel('st4');
-
%plot(f,abs(st1));
-
axis([-10,10,-2,2]);
-
hold on;
-
plot(t,mt,'r--');
-
title('解调后波形');
-
legend('解调后信号','基带信号')
-
grid on;
-
gtext('图2-4 SSB调制信号及其功率谱密度波形')
在信道中各自加入经过带通滤波器后的窄带高斯白噪声，功率为0.1，解调各个信号，并画出解调后的波形。
-
clc,clear,close all;
-
A=2;
-
dt=0.001;
-
t=-10:dt:10;% 仿真时间
-
fm=1;%基带信号频率
-
fc=10;%载波频率
-
w_lpf = 1;% 低通滤波器截止频率（理想）
-
w_ipf=1; % 带通滤波器截止频率（理想）
-
mt=sqrt(2)*cos(2*pi*fm*t);%基带信号
-
c=cos(2*pi*fc*t);%载波信号
-
%获得各种信号
-
st1=mt.*c;%DSB信号
-
st2=(mt+A).*c;%AM信号
-
st3= 0.5*mt.* c + 0.5*hilbert(mt).*hilbert(c); %SSB信号
-
%给各种信号加高斯白噪声
-
st11=awgn(st1,-10);%DSB信号加高斯白噪声
-
st22=awgn(st2,-10);%AM信号加高斯白噪声
-
st33=awgn(st3,-10);%SSB信号加高斯白噪声
-
%解调
-
sample = length(t);
-
f = 1/dt/sample*(0:1:sample-1);
-
st11=st11.*c;%DSB信号解调
-
st111 = fft(st11, sample)/sample*2; % 这里要/抽样次数*2 才是最后结果
-
st111 = st111.*(f<w_lpf); % 进过低通滤波
-
st1111 = real(ifft(st111)); %解调后波形
-
st22=st22.*c;%AM信号解调
-
st222 = fft(st22, sample)/sample*2; %这里要/抽样次数*2 才是最后结果
-
st222 = st222.*((f>w_ipf-fm) - (f>w_ipf+fm) );%进过低通滤波
-
st2222 = real(ifft(st222)); %解调后波形
-
st33=st33.*c;%SSB信号解调
-
st333 = fft(st33, sample)/sample*2; % 这里要/抽样次数*2 才是最后结果
-
st333 = st333.*((f>w_ipf-fm) - (f>w_ipf+fm) ); % 进过低通滤波
-
st3333 = real(ifft(st333)); %解调后波形
-
%DSB经过高斯白噪声后解调波形与原波形比较
-
subplot(311);
-
plot(t,st1111*sample);
-
xlabel('t');
-
ylabel('mt');
-
hold on;
-
plot(t,mt,'r--');
-
title('经过高斯白噪声后波形比较');
-
legend('DSB','基带信号');
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%AM经过高斯白噪声后解调波形与原波形比较
-
subplot(312)
-
plot(t,st2222*sample)
-
xlabel('t');
-
ylabel('mt');
-
hold on;
-
plot(t,mt,'r--')
-
title('经过高斯白噪声后波形比较');
-
legend('AM','基带信号');
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
%SSB经过高斯白噪声后解调波形与原波形比较
-
subplot(313)
-
plot(t,st3333*sample)
-
xlabel('t');
-
ylabel('mt');
-
hold on;
-
plot(t,mt,'r--');
-
title('经过高斯白噪声后波形比较');
-
legend('SSB','基带信号')
-
axis([-5,5,-1.5,1.5]);
-
grid on;
-
gtext('图2-5 基带信号与经过高斯白噪声后DSB、SSB、AM解调后波形比较')
文章来源: blog.csdn.net，作者：渣渣ye，版权归原作者所有，如需转载，请联系作者。
原文链接：blog.csdn.net/yyfloveqcw/article/details/123811018
【版权声明】本文为华为云社区用户转载文章，如果您发现本社区中有涉嫌抄袭的内容，欢迎发送邮件进行举报，并提供相关证据，一经查实，本社区将立刻删除涉嫌侵权内容，举报邮箱：
cloudbbs@huaweicloud.com
- 点赞
- 收藏
- 关注作者
评论（0）