AWGN信道模型下的模拟调制解调matlab实现.zip_AWGN信道，解调_AWGN模型_matlab实现AWGN_matla

% MATLAB script for Illustrative Problem 3.5. % Demonstration script for DSB-AM demodulation. The message signal % is +1 for 0 < t < t0/3, -2 for t0/3 < t < 2t0/3, and zero otherwise. clear echo on f0=1e3; % 基带信号频率 fc=4e4; % 载波信号频率 fs=2e5; % 采样频率 ts=1/fs; % 采样周期 t0=1/f0; %信号周期 t=[0:ts:t0]; % 时间向量 df=100; % 所需的频率分辨率 m=cos(2*pi*f0*t);%消息信号 m_n=m/max(abs(m)); c=cos(2*pi*fc.*t); % 载波信号 %u1=m.*c; %已调信号 A2=24.5;%已调信号振幅 f=fc; %载波信号频率 u1=(1+0.85*m_n).*c; u=awgn(u1,1);%加入噪声，awgn（）函数的第二个参数是信噪比 y=u.*c;%同步相乘检波 [M,m,df1]=fftseq(m,ts,df); % 对基带信号进行FFT变换 M=M/fs; % 缩放 [U,u,df1]=fftseq(u,ts,df); % 对DSB已调信号进行FFT变换 U=U/fs; % 缩放 [Y,y,df1]=fftseq(y,ts,df); % 对同步相乘信号进行FFT变换 Y=Y/fs; % 缩放 f_cutoff=3e4; % 低通滤波器截止频率 n_cutoff=floor(f_cutoff/df1); % Design the filter.设置LPF参数 f=[0:df1:df1*(length(y)-1)]-fs/2; H=zeros(size(f)); H(1:n_cutoff)=2*ones(1,n_cutoff); H(length(f)-n_cutoff+1:length(f))=2*ones(1,n_cutoff); DEM=H.*Y; % 滤波器频谱输出，spectrum of the filter output dem=real(ifft(DEM))*fs; % 滤波器输出，filter output pause % 按任意键显示混频结果，Press a key to see the effect of mixing. clf subplot(3,1,1) %三行一列figure图中第一个图 plot(f,fftshift(abs(M))) title('基带信号的频谱') xlabel('Frequency') subplot(3,1,2) plot(f,fftshift(abs(U))) title('已调信号的频谱') xlabel('Frequency') subplot(3,1,3) plot(f,fftshift(abs(Y))) title('同步检波输出频谱') xlabel('Frequency') pause % Press a key to see the output. clf subplot(3,1,1) plot(f,fftshift(abs(Y))) title('同步检波输出') xlabel('Frequency') subplot(3,1,2) plot(f,fftshift(abs(H))) title('LPF低通特性') xlabel('Frequency') subplot(3,1,3) plot(f,fftshift(abs(DEM))) title('解调器输出频谱') xlabel('Frequency') pause % Press a key to compare the spectra of the message and the received signal. clf subplot(2,1,1) plot(f,fftshift(abs(M))) title('消息信号频谱') xlabel('Frequency') subplot(2,1,2) plot(f,fftshift(abs(DEM))) title('解调输出频谱') xlabel('Frequency') pause % Press a key to see the message and the demodulator output signals. subplot(2,1,1) plot(t,m(1:length(t))) title('信息信号') xlabel('Time') subplot(2,1,2) plot(t,unwrap(dem(1:length(t))-1)) axis([0,1e-3,-1,1]) title('在AWGN信道下信噪比为1dB的解调输出') xlabel('Time')