通信原理第十次实验报告.docVIP

通信原理第十章实验报告 10.1 %m序列发生器及其自相关 mseq.m clear all; close all; g = 19; % G = 10011 state = 8;% state=1000 L = 1000; %m序列产生 N = 15; mq = mgen(g,state,L); %求序列自相关 ms = conv( 1-2*mq, 1-2*mq(15:-1:1) )/N; figure(1) subplot(222) stem(ms(15:end)); axis([0 63 -0.3 1.2]);title(m序列自相关序列) %m序列构成的信号(矩形脉冲） N_sample=8; Tc = 1; dt = Tc/N_sample; t = 0:dt:Tc*L-dt; gt = ones(1,N_sample); mt = sigexpand(1-2*mq,N_sample); mt = conv(mt,gt); figure(1) subplot(221); plot(t,mt(1:length(t))); axis([0 63 -1.2 1.2]);title(m序列矩形成形信号) st = sigexpand( 1-2*mq(1:15),N_sample ); s = conv(st,gt); st = s(1:length(st)); rt1 = conv(mt,st(end:-1:1))/(N*N_sample); subplot(223) plot(t,rt1(length(st):length(st)+length(t)-1) ); axis([0 63 -0.2 1.2]);title(m序列矩形成形信号的自相关);xlabel(t); %sinc脉冲 Tc = 1; dt = Tc/N_sample; t = -20:dt:20; gt = sinc(t/Tc); mt = sigexpand(1-2*mq,N_sample); mt = conv(mt,gt); st2 = sigexpand( 1-2*mq(1:15),N_sample ); s2 = conv(st2,gt); st2 = s2; rt2 = conv(mt,st2(end:-1:1))/(N*N_sample); subplot(224); t1 = -55+dt:dt:Tc*L-dt; %plot(t,mt(1:length(t)) ); plot(t1,rt2(1:length(t1))); axis([0 63 -0.5 1.2]);title(m序列sinc成形信号的自相关);xlabel(t); %m序列发生器, mgen.m function [out] = mgen(g,state,N) %输入 g: m序列生成多项式 (10进制输入） % state: 寄存器初始状态 （10进制输入） % N: 输出序列长度 % test g =11; state=3; N=15; gen = dec2bin(g) - 48; M = length(gen); curState = dec2bin(state,M-1) - 48; for k=1:N out(k) = curState(M-1); a = rem( sum( gen(2:end).*curState ),2 ); curState = [a curState(1:M-2)]; end 10.2 %码分多址复用示意 cdm.m clear all; close all; Ts = 1; N = 4; %用户数 %产生用户数据, 1000bits randn(state,sum(100*clock)); d1 = sign(randn(1,100000)); d2 = sign(randn(1,100000)); d3 = sign(randn(1,100000)); d4 = sign(randn(1,100000)); dd = [d1;d2;d3;d4]; %产生4阶walsh码 w = hadamard(4); %调用Matlab函数 w = w/2; %能量归一化 %用户数据复合 s = [d1;d2;d3;d4]*w; %复合 ss = reshape(s,1,4*100000); stairs(ss); %A点波形 axis([0 40 -2.2 2.2]); xlabel(t/Tc); yl