channelpass.m

来自「就是自适应的一写程序」· M 代码 · 共 20 行

function receive=channelpass(sig_in,channel,snr_db)
%信道：多径加噪声
avg_sig_power = sum(abs(sig_in).^2)/length(sig_in);
%求信号的平均功率
snr = power(10,snr_db/10);
%换算信噪比
avg_noise_power = avg_sig_power/(2*snr);    
% 实部或虚部的噪声功率
real_noise = randn(length(sig_in)+length(channel)-1,1)*sqrt(avg_noise_power);
image_noise = randn(length(sig_in)+length(channel)-1,1)*sqrt(avg_noise_power);
noise_out = real_noise+j*image_noise;       
%加高斯白噪声
% channel1=[0.26,
%           0.93,
%           0.26];                             %针对p496页的信道情况
receive=conv(sig_in,channel)+noise_out;  
%将映射后的信号和信道做卷积，再加高斯白噪
% scatterplot(receive)                           %画星座图