1
2
code = randi([0,1],1,N);   %N为码元个数 1或者0

1
2
3
4
code = randi([0,1],1,N);   %N为码元个数 1或者0
I = find(code == 0);
code(I) = -1;

1
2
3
rng default
x = randn(36,1); 

1
2
signal2 = awgn(carrier2,-50,'measured');  %SNR = -50dB

1
2
3
4
5
6
7
8
9
10
11
12
13
14
%noisegen：信号加噪函数
%X：信号输入，SNR：信噪比
%Y:加了噪声的信号，NOISE:所加的噪声
function [Y,NOISE] = noisegen(X,SNR)  

NOISE=randn(size(X));                                  %生成均值为0方差为1的正态分布白噪声
NOISE=NOISE-mean(NOISE);                               %但在任意长度下x=randn(1,N)，x不一定是均值为0方差为1（有些小小的偏差），因此在这里做一点修正
signal_power = 1/length(X)*sum(X.*X);                  %计算信号的功率
noise_power_base = 1/length(NOISE)*sum(NOISE.*NOISE);  %当前噪声序列的功率
noise_power = signal_power / ( 10^(SNR/10) );          %实际需要的噪声功率
NOISE=sqrt(noise_power/noise_power_base)*NOISE;        %修正噪声序列的幅值，满足SNR的要求
Y=X+NOISE;
end

1
2
3
4
5
6
7
8
9
10
% 计算信噪比函数
% I :原始信号
% In:加噪信号（原始信号+噪声信号）
function snr=SNR_singlech(I,In)
snr=0;
Ps=sum((I.^2));                %信号功率
Pn=sum((In - I).^2);           %噪声功率
snr=10*log10(Ps/Pn);
end

1
2
参考之前GNSS中NOISE_GEN的处理方式，对噪声信号和噪声需要先进行滤波，去除带外能量