基于Matlab的CDMA通信系统仿真(3)

- 11 -

通信系统综合设计与实践

致谢

在本次仿真设计中，感谢指导老师在论文选题、设计任务书中对我的无私帮助，加深了我对专业基础理论的理解，拓宽了我的专业知识面，实在是受益匪浅。

感谢我组同学的团结协作，在设计过程中，共同商讨，共同去图书馆查阅相关参考资料，使我组能够抓紧宝贵的时间，完成论文的写作。

在此，谨对所有在论文写作中帮助过我的老师、同学表示衷心的感谢和由衷的敬意！谢谢大家！

- 12 -

通信系统综合设计与实践

参考文献

[1] 刘学勇,编著.通信系统建模与仿真(电子工业出版社) [2] 邓薇，编著.MATALB函数速查手册（人民邮电出版社）

[3] 【美】William C.Y.Lee,著.移动通信工程理论和应用（第二版）（人民邮电出版社）

[4] 王华奎，李艳萍等编著.移动通信原理与技术(清华大学出版社)

- 13 -

通信系统综合设计与实践

附录（各部分编程）：

1、M序列

function [mout] = mseq(n, taps, inidata, num)

% **************************************************************** % n : m序列的阶数n % taps : 反馈寄存器的连接位置 % inidata : 寄存器的初始值序列

% num : 输出的m序列的个数

% mout : 输出的m序列，如果num>1,则每一行为一个m序列

% ****************************************************************

if nargin < 4 num = 1; end

mout = zeros(num,2^n-1); fpos = zeros(n,1);

fpos(taps) = 1;

for ii=1:2^n-1

mout(1,ii) = inidata(n); % 寄存器的输出值 temp = mod(inidata*fpos,2); % 计算反馈数据

inidata(2:n) = inidata(1:n-1); % 寄存器移位一次 inidata(1) = temp; % 更新第1个寄存器的值 end

if num > 1 %如果要输出多个m序列，生成其他m序列

for ii=2:num

mout(ii,:) = shift(mout(ii-1,:),1);

- 14 -

通信系统综合设计与实践

end end

关于shift函数的代码：

function [outregi] = shift(inregi,shiftr)

% **************************************************************** % inrege : 输入序列

% shiftr : 循环右移的位数

% outregi : 输出序列

% ****************************************************************

v = length(inregi); outregi = inregi;

shiftr = rem(shiftr,v);

if shiftr > 0

outregi(:,1:shiftr) = inregi(:,v-shiftr+1:v); %循环移位 outregi(:,1+shiftr:v) = inregi(:,1:v-shiftr); elseif shiftr < 0

outregi(:,1:v+shiftr) = inregi(:,1-shiftr:v); outregi(:,v+shiftr+1:v) = inregi(:,1:-shiftr); end

%********************* end of file ********************************

- 15 -

通信系统综合设计与实践

2、Gold序列：

function [gout] = goldseq(m1, m2, num)

% **************************************************************** % m1 : m序列1 % m2 : m序列2

% num : 生成的Gold序列个数 % gout : 生成的Gold序列输出

% ****************************************************************

if nargin < 3 %如果没有指定生成的Gold序列个数，默认为1 num = 1; end

gout = zeros(num,length(m1));

for ii=1:num %根据Gold序列生成方法生成Gold序列 gout(ii,:) = xor(m1,m2); m2 = shift(m2,1); end

%********************** end of file ********************************

其shift函数代码同M序列