t=0:255;%信号时间范围
x1=sin(8*pi*t/100)+1;%信号源
sigPower = sum(abs(x1).^2)/length(x1) ; %求出信号功率
k=sigPower;
subplot(221);
plot(t,x1);
title('原始输入波形')
x2=sin(40*pi*t/100);
x3=x1.*x2; %AM调制
subplot(222);
plot(t,x3);
title('与载波相乘后波形')
SNR=exp(25*log(10)/10);
x3=awgn(x3,SNR*k*k*k/27,'measured','linear');%加入AWGN噪声
subplot(223);
plot(t,x3);
title('加入噪声后波形')
x4=envelop(x3,1);%包络检波解调
subplot(224);
plot(t,x4);
title('包络检波后输出波形')
- 1
- 2
前往页