3g移动通信实验7.docVIP

实验报告 课程名称： 3G移动通信实验 实验名称：伪随机序列生成编程（2） 学生姓名： 章仕波 学生学号： 实验日期： 5月15日 实验六 伪随机序列生成编程（二） 【实验目的】 加深对伪随机序列性能的理解； 能够使用Matlab（或者C语言）生成伪随机序列； 比较这些伪随机序列的自相关函数和互相关函数特性。 【实验内容】 画出m序列的自相关函数和互相关函数并进行分析； 【实验步骤】 画出两个m序列各自的自相关函数，以及这两个序列的互相关函数。（本原多项式为（45）8，（75）8。） 【实验报告】 %编程产生两个周期长度为31的m序列；（本原多项式为（45）8，（75）8。 %(45)8=100101,(75)8=111101 clear all; clc; r=5; N=2^r-1;%周期 T=30; Tmin=1; a1=ones(1,r);%移位寄存器 m1=zeros(1,2^r-1);%生成的m序列 for i=1:(2^r-1) temp= mod((a1(5)+a1(2)),2); for j=r:-1:2 a1(j)=a1(j-1); end a1(1)=temp; m1(i)=a1(r); end subplot(2,1,1) stem(m1);%序列1 axis([0 ,2^r, 0,2 ]); title(本原多项式m1为（45）8); a2=ones(1,r);%移位寄存器 m2=zeros(1,2^r-1);%生成的m序列 for i=1:(2^r-1) temp= mod((a2(5)+a2(4)+a2(3)+a2(2)),2); for j=r:-1:2 a2(j)=a2(j-1); end a2(1)=temp; m2(i)=a2(r);%序列2 end subplot(2,1,2) stem(m2); axis([0 ,2^r, 0,2 ]); title(本原多项式m2为（75）8); m1=m1*2-1; m2=m2*2-1; %******************** % figure Raa=zeros(1,T+1); for i=0:T Raa(i+1)=0; for j=1:N; k=i+j; while(kN) k=k-N; end Raa(i+1)=Raa(i+1)+m1(j)*m1(k); end Raa(i+1)=Raa(i+1)/N; end subplot(3,1,1) stem(0:T,Raa); title(m1自相关函数); %************************ Rbb=zeros(1,T+1); for i=0:T Rbb(i+1)=0; for j=1:N; k=i+j; while(kN) k=k-N; end Rbb(i+1)=Rbb(i+1)+m2(j)*m2(k); end Rbb(i+1)=Rbb(i+1)/N; end subplot(3,1,2) stem(0:T,Rbb); title(m2自相关函数); %********************************** Rab=zeros(1,T+1); for i=0:T Rab(i+1)=0; for j=1:N; k=i+j; while(kN) k=k-N; end Rab(i+1)=Rab(i+1)+m1(j)*m2(k); end Rab(i+1)=Rab(i+1)/N; end subplot(3,1,3) stem(0:T,Rab); title(m1,m2互相关函数); 【背景知识】 【伪代码】 a = zeros(1,r); seq_m = zeros(1, 2^r-1); for j = 1: (2^r-1) temp = fun_feedback(a)； For i = r:-1:2 a(i) = a(i-1); end a(1) = temp; seq_m(j) = a(r); end