clear variables;
close all;
clc;
% 产生Walsh 64序列
walsh = hadamard(64);
% 二进制信息码
infoCode = rand(1,64);
infoCode(infoCode > 0.5) = 1;
infoCode(infoCode <= 0.5) = -1;
subplot(4,2,1);
stairs(infoCode);axis([1 64 -1.5 1.5]);title('信息码波形图');
infoCode(infoCode==-1) = 0; % 便于xor运算
% index = randperm(64); % 随选择扩频码
% 生成扩频码
% ssCode1 = walsh(index(1,1),:);
ssCode1 = walsh(8,:);
subplot(4,2,2);
stairs(ssCode1);axis([1 64 -1.5 1.5]);title('扩频码1');
% ssCode2 = walsh(index(1,2),:);
ssCode2 = walsh(16,:);
subplot(4,2,3);
stairs(ssCode2);axis([1 64 -1.5 1.5]);title('扩频码2')
% 观察扩频码之间的自相关性和互相关性
[a,b]=xcorr(ssCode1,'coeff');
subplot(4,2,4);
plot(b,a);title('自相关性');
[a,b]=xcorr(ssCode1,ssCode2,'coeff');
subplot(4,2,5);
plot(b,a);ylim([-1,1]);title('互相关性');
ssCode1(ssCode1 == -1) = 0; % 便于xor运算
ssCode2(ssCode2 == -1) = 0;
% 扩频
ssedCode = xor(infoCode,ssCode1); % 异或
ssedCode_disp = im2double(ssedCode);
ssedCode_disp(ssedCode_disp==0) = -1;
subplot(4,2,6);
stairs(ssedCode_disp);axis([1 64 -1.5 1.5]);title('扩频码1扩频');
% 解扩
dssCode = xor(ssedCode,ssCode1);
dssCode = im2double(dssCode);
dssCode(dssCode==0) = -1;
subplot(4,2,7);
stairs(dssCode);axis([1 64 -1.5 1.5]);title('扩频码1解扩');
% 解扩2
dssCode = xor(ssedCode,ssCode2);
dssCode = im2double(dssCode);
dssCode(dssCode==0) = -1;
subplot(4,2,8);
stairs(dssCode);axis([1 64 -1.5 1.5]);title('扩频码2解扩');
实现效果