现代数字信号处理05.pptVIP

现代数字信号处理 作业 构造两个m码序列， C1 = 453，C2 = 703 用上面两个m序列5个构造戈登码序列 求它们的相关特性 程序 corr_xy.m % cross-correlation x y % QZD/11.03/Nokia Tampere function z = corr_xy(x, y) N = length(x); b = [y, y]; for k = 1:N z(k) = x * b(k:k+N-1); end 程序 m_1023.m function c = m_1023(k) % 1023 m-code, QZD/05.94/Nokia Oulu % k is for different polynomial R = ones(1,10); % initial for n = 1:1023 switch k case 1 B = xor(R(10),R(7)); % 2011 case 2 B1 = xor(R(9), R(10)); B2 = xor(R(6), R(7)); B = xor(B1, B2); % 2033 otherwise % some more, add later B = xor(R(10),R(7)); % 2011 end c(n)= R(10); if(c(n)1), c(n) = -1; end for r=10:-1:2, R(r)=R(r-1); end, R(1)=B; end % 两个m码的相关特性 clear all c1 = m_1023(1); c2 = m_1023(2); r11 = corr_xy(c1, c1); r22 = corr_xy(c2, c2); r12 = corr_xy(c1, c2); subplot(311), plot(r11/1023); axis([1 1023 -0.1 1]) subplot(312), plot(r22/1023); axis([1 1023 -0.1 1]) subplot(313), plot(r12/1023); axis([1 1023 -0.1 1]) * 上次的作业 考虑一个10个码片的周期序列（ 1001101010 ），采样率为100 Hz(每个码片一个采样点）。a, 对该序列进行谱分析; b, 对该序列的自相关序列进行谱分析。画出谱分布。 2. 对上面的周期序列用一个2.3Hz的Doppler信号进行调制，观察Dopper对自相关峰（t = 0) 的影响。 exp(j2pft) = exp(j2p2.3[0:9]/100) 3. 调整Doppler频率，观察其变化。 4. 预习编码过程。 M序列（最长序列）码产生器 反馈式移位寄存器的输出是一个周期序列； 周期由反馈系数和寄存器初始状态决定； 寄存器初始状态为全0，输出为全0，与反馈系数无关； 对于莫种反馈系数，寄存器初始状态为非全0时，产生 最长周期 N = 2m-1, 称为m序列，也称为最长周期序列。 产生m序列的充要条件： 满足这四个条件的多项式，为本原多项式。 m码的基本性质 1. 由于全部状态都经历（除全0）周期N=2n-1, 在一个周期内，0出现的次数与1出现的次数仅仅差1; 在一个周期内，共有2n-1个游程，其长度为k的游程数的比例为2-k (1 = k = m-2), +一个长度为n的1的游程和长度为n-1的0的游程； 具有双值自相关特性，归一化自相关函数为 主瓣中集中了约90%的功率，定义信号带宽为1/Tc m码的优缺点 伪随机特征好，均衡性一般 良好的自相关 简单，实现方便，硬件少 不同特征值的不同码之间非相关性能好 有小的直流分量 必威体育官网网址性能差 码长受限制（2n-1) 序列的数量受限制 m构造码 戈德码（Golden code) 时频空多维码 软件码（Software codes) *