clear all;
format long;
N = 4; %设置沃尔什-哈达玛序列阶数，可供2^n个用户使用。
x_num = 50; %设置传输信息长度
UserCode = hadamard(N); %生成沃尔什-哈达玛序列
x1 = bingen(x_num); %随机生成用户1要发送的信号（子程序见附录）
c1 = UserCode(4,:); %取出哈达玛矩阵第4行作为用户1的识别码
y1 = ds_mod(c1(:),x1); %对用户1的信号进行ds扩频（子程序见附录）
x2 = bingen(x_num); %随机生成用户2要发送的信号
c2 = UserCode(2,:); %取出哈达玛矩阵第2行作为用户2的识别码
y2 = ds_mod(c2(:),x2); %对用户2的信号进行ds扩频
x3 = bingen(x_num); %随机生成用户3要发送的信号
c3 = UserCode(3,:); %取出哈达玛矩阵第3行作为用户3的识别码
y3 = ds_mod(c3(:),x3); %对用户2的信号进行ds扩频
y = y1+y2+y3; %把三路信号加在一起形成基带信号
figure(1),
subplot(3,1,1);plot(x1);title('用户1');
subplot(3,1,2);plot(x2);title('用户2');
subplot(3,1,3);plot(x3);title('用户3');
figure(2),
subplot(3,1,1);plot(y1);title('用户1扩频后');
subplot(3,1,2);plot(y2);title('用户2扩频后');
subplot(3,1,3);plot(y3);title('用户3扩频后');
z = awgn(y,1); %加入高斯白噪声
figure(3),
subplot(2,1,1);plot(y);title('基带信号');
subplot(2,1,2);plot(z);title('加入高斯白噪声后的基带信号');
x1_de = ds_demod(c1(:),z); %使用用户1的识别码对加入噪声的基带信号进行解扩
x1_de(find(x1_de < 0)) = -1; %(解扩子程序见附录)
x1_de(find(x1_de >=0)) = 1; %对解扩后的信号进行判决，从而恢复出用户1的信号
x2_de = ds_demod(c2(:),z); %以下同上
x2_de(find(x2_de < 0)) = -1;
x2_de(find(x2_de >=0)) = 1;
x3_de = ds_demod(c3(:),z);
x3_de(find(x3_de < 0)) = -1;
x3_de(find(x3_de >=0)) = 1;
figure(4),
subplot(3,1,1);plot(x1_de);title('解扩后用户1');
subplot(3,1,2);plot(x2_de);title('解扩后用户2');
subplot(3,1,3);plot(x3_de);title('解扩后用户3');
评论0