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

杭 州 电 子 科 技 大 学 实 验 报 告 学生姓名： xx 学 号： 13xxxxxx 指导教师：xxx 实验地点： 1教108 实验时间2016.5.19 一、实验室名称： 1教108 二、实验项目名称： 伪随机序列生成编程（二） 三、实验学时：3 四、实验内容：画出m序列的自相关函数和互相关函数并进行分析； 五、实验结果： %±à3ì2úéúá????ü?ú3¤?è?a31μ?mDòáD￡?￡¨±??-?à??ê??a￡¨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);%éú3éμ?mDòáD 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);%DòáD1 axis([0 ,2^r, 0,2 ]); title(±??-?à??ê?m1?a￡¨45￡?8); a2=ones(1,r);%ò?????′??÷ m2=zeros(1,2^r-1);%éú3éμ?mDòáD 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);%DòáD2 end subplot(2,1,2) stem(m2); axis([0 ,2^r, 0,2 ]); title(±??-?à??ê?m2?a￡¨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×??à1?oˉêy); %************************ 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×??à1?oˉêy); %********************************** 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?￥?à1?oˉêy); 六、总结及心得体会： 加深对伪随机序列性能的理解；能够使用Matlab（或者C语言）生成伪随机序列； 比较这些伪随机序列的自相关函数和互相关函数特性。