3G实验六精要.docxVIP

课程实验报告课 程3G移动通信实验学院通信工程学院专业信息工程班生姓名樊嘉鑫实验六伪随机序列生成编程一、实验目的加深对伪随机序列性能的理解；能够使用Matlab（或者C语言）生成伪随机序列；比较这些伪随机序列的自相关函数和互相关函数特性。二、实验内容编程产生两个周期长度为31的m序列；（本原多项式为（45）8，（75）8。画出两个m序列各自的自相关函数，以及这两个序列的互相关函数。（本原多项式为（45）8，（75）8。）代码：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); enda1(1)=temp; m1(i)=a1(r); endsubplot(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);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); enda2(1)=temp; m2(i)=a2(r);endsubplot(2,1,2) stem(m2); axis([0 ,2^r, 0,2 ]); title(本原多项式m2为（75）8); m1=m1*2-1; m2=m2*2-1; % 求m1自相关函数figureRaa=zeros(1,T+1); for i=0:T Raa(i+1)=0; for j=1:N; k=i+j; while(kN) k=k-N; endRaa(i+1)=Raa(i+1)+m1(j)*m1(k); endRaa(i+1)=Raa(i+1)/N; endsubplot(3,1,1) stem(0:T,Raa); title(m1自相关函数); % 求m2自相关函数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; endRbb(i+1)=Rbb(i+1)+m2(j)*m2(k); endRbb(i+1)=Rbb(i+1)/N; endsubplot(3,1,2) stem(0:T,Rbb); title(m2自相关函数); % 求m1和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; endRab(i+1)=Rab(i+1)+m1(j)*m2(k); endRab(i+1)=Rab(i+1)/N; endsubplot(3,1,3) stem(0:T,Rab); title(m1,m2互相关函数); 实验结果PN码是一种具有类似白噪声性质的码。瞬态值服从正态分布，功率谱在很宽的频带内是均匀的，而且有良好的相关特性实验七直接序列扩频信号产生与分析一、实验目的加深对直接序列扩频信号产生的理解；能够使用Matlab语言（或者C语言）生成直接序列扩频信号；观察和比较发送端和接收端的星座图。二、实验内容使用使用Matlab语言（或者C语言）采用m序列作为扩频序列分别产生BPSK调制和QPSK调制时的扩频信号；clc; clear all; r=6; N=2^r-1; a=ones(1,r); m=zeros(1,N); for i=1:(2^r-1) temp= mod((a(5)+a(2)),2); for j=r:-1:2 a(j)=a(j-1); enda(1)=temp; m(i)=a(r); endm=m*2-1; %产生随机发送信息，并进行BPSK/QPSK调制n=100; j=sqrt(-1); s_real=randint(1,n); s_imag=randint(1,n); s=s_real+j*s_imag; s_BPSK=2*s_real-1; s_QPSK=s*2-(1+j); %产生扩频信号：for z=1:1:100 y_QPSK(63*(z-1)+(1