杭电3G移动通信实验(实验6)-伪随机序列生成编程(二)(共4页)

1、杭 州 电 子 科 技 大 学实 验 报 告学生姓名： xx 学 号： 13xxxxxx 指导教师：xxx实验地点：1教108 实验时间2016.5.19一、实验室名称： 1教108二、实验项目名称： 伪随机序列生成编程（二）三、实验学时：3四、实验内容：画出m序列的自相关函数和互相关函数并进行分析；五、实验结果：%±à³Ì²úÉúÁ½¸öÖÜÆÚ³¤¶ÈΪ31µ&

2、#196;mÐòÁУ»£¨±¾Ô­¶àÏîʽΪ£¨45£©8£¬£¨75£©8¡£%(45)8=100101,(75)8=111101 clear all; clc; r=5; N=2r-1;%ÖÜÆÚ T=30; Tmin=1; a

3、1=ones(1,r);%ÒÆλ¼Ä´æÆ÷ m1=zeros(1,2r-1);%Éú³ÉµÄmÐòÁÐ for i=1:(2r-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);%ÐòÁ&

4、#208;1 axis(0 ,2r, 0,2 ); title('±¾Ô­¶àÏîʽm1Ϊ£¨45£©8'); a2=ones(1,r);%ÒÆλ¼Ä´æÆ÷ m2=zeros(1,2r-1);%Éú³ÉµÄmÐòÁÐ

5、 for i=1:(2r-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 ,2r, 0,2 ); title('±¾Ô­¶àÏîʽm2Ϊ£¨75£©8

6、9;); 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(k>N) 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=

7、zeros(1,T+1); for i=0:T Rbb(i+1)=0; for j=1:N; k=i+j; while(k>N) 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(k>N) 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»¥Ï