移动通信课程设计实验报告-利用matlab进行m序列直接扩频仿真

目录

一、背景 .................................................................................................... 4 二、基本要求 ............................................................................................ 4 三、设计概述 .......................................................................................... 4 四、Matlab设计流程图 ........................................................................... 5 五、Matlab程序及仿真结果图 ............................................................... 6

1、生成m序列及m序列性质 ............................................................................ 6 2、生成50位随机待发送二进制比特序列，并进行扩频编码......................... 7 3、对扩频前后信号进行BPSK调制，观察其时域波形................................... 9 4、计算并观察扩频前后BPSK调制信号的频谱............................................. 10 5、仿真经awgn信道传输后，扩频前后信号时域及频域的变化 .................. 11 6、对比经信道前后两种信号的频谱变化......................................................... 12 7、接收机与本地恢复载波相乘，观察仿真时域波形..................................... 14 8、与恢复载波相乘后，观察其频谱变化......................................................... 15 9、仿真观察信号经凯萨尔窗低通滤波后的频谱............................................. 16 10、观察经过低通滤波器后无扩频与扩频系统的时域波形........................... 17 11、对扩频系统进行解扩，观察其时域频域 ................................................... 18 12、比较扩频系统解扩前后信号带宽............................................................... 19 13、比较解扩前后信号功率谱密度................................................................... 20 14、对解扩信号进行采样、判决....................................................................... 21 15、在信道中加入2040~2050Hz窄带强干扰并乘以恢复载波 ...................... 24 16、对加窄带干扰的信号进行低通滤波并解扩............................................... 25 17、比较解扩后信号与窄带强干扰的功率谱................................................... 27

六、误码率simulink仿真 ..................................................................... 28

1、直接扩频系统信道模型................................................................................. 28

2、加窄带干扰的直扩系统建模......................................................................... 29 3、用示波器观察发送码字及解扩后码字......................................................... 30 4、直接扩频系统与无扩频系统的误码率比较................................................. 31 5、不同扩频序列长度下的误码率比较............................................................. 32 6、扩频序列长度N=7时，不同强度窄带干扰下的误码率比较 ................... 33

七、利用Walsh码实现码分多址技术 ................................................. 34

1、产生改善的walsh码 ..................................................................................... 35

2、产生两路不同的信息序列............................................................................. 36 3、用两个沃尔什码分别调制两路信号............................................................. 38

- 1 -

4、两路信号相加，并进行BPSK调制............................................................. 39 5、观察调制信号频谱，并经awgn信道加高斯白噪和窄带强干扰 .............. 40 6、接收机信号乘以恢复载波，观察时域和频域............................................. 42 7、信号经凯萨尔窗低通滤波器......................................................................... 43 8、对滤波后信号分别用m1和m2进行解扩 .................................................. 44 9、对两路信号分别采样，判决......................................................................... 45

八、产生随机序列Gold码和正交Gold码 ......................................... 47

1、产生Gold码并仿真其自相关函数............................................................... 48

2、产生正交Gold码并仿真其互相关函数....................................................... 50

九、实验心得体会 .................................................................................. 51

- 2 -

直接序列扩频系统仿真

一、背景

直接序列扩频通信系统(DSSS)是目前应用最为广泛的系统。在发送端，直扩系统将发送序列用伪随机序列扩展到一个很宽的频带上去，在接受端又用相同的扩频序列进行解扩，回复出原有信息。由于干扰信息与伪随机序列不相关，扩频后能够使窄带干扰得到有效的抑制，提高输出信噪比。系统框图如下图所示：

二、基本要求：

1. 通过matlab建模，对直扩系统进行仿真，数据调制方式可以自由选择，可以

使用基带信号，但最好能使用频带信号，信道为高斯白噪信道。要仿真出扩频前的信号的频偏，扩频后的信号频谱，过信道之后的频谱以及解扩之后的频谱。

2. 研究并仿真产生m序列，写出生成m序列的算法。

3. 验证直扩系统对窄带干扰的抑制能力，在信道中加入一个窄带强干扰，仿真

出加了干扰后的频谱图和解扩后的频谱图，给出误码率等仿真图。

4. 在以上基础上仿真实现码分多址技术，使用Walsh码进行复用，实现多个信

号同时传输。(选做) 可选项：

1.在信道中加入多径，使用rake接收来抗多径效应。

2.产生除m序列之外的其他随机序列，如Gold码，正交Gold码等等。 3.对比无扩频的系统的误码率。

三、设计概述

本次课设完成基本要求，并选作了可选项码分多址，Gold码及误码率对比。通过matlab建模仿真了直扩系统BPSK调制的各点频偏及时域信号，并仿真了窄带强干扰对直扩系统的影响以及利用改善的WALSH码实现码分多址技术。另外，通过matlab的simulink工具盒bertool工具仿真对比了直扩系统和无扩频系统的误码率。

- 3 -

四、matlab设计流程图

基本扩频系统仿真流程图

100/7Hz 二进制比特信息 100Hz 7位双极性m序列 100Hz 扩频序列 2000Hz 载波cos4000πt BPSK调制信号 高斯白噪声 恢复载波cos4000πt 凯萨尔滤波器低通滤波 100Hz 7位双极性m序列 采样、判决 - 4 -

五、matlab程序及仿真结果图

1、生成m序列及m序列性质

实验产生7位m序列，频率100Hz，模拟线性反馈移位寄存器序列，原理图如下：

clear all; clc;

X1=0;X2=0;X3=1;

m=350; %重复50遍的7位单极性m序列 for i=1:m

Y3=X3; Y2=X2; Y1=X1; X3=Y2; X2=Y1; X1=xor(Y3,Y1); L(i)=Y1; end

for i=1:m

M(i)=1-2*L(i); %将单极性m序列变为双极性m序列 end

k=1:1:m; figure(1)

subplot(3,1,1) %做m序列图 stem(k-1,M); axis([0,7,-1,1]); xlabel('k');

ylabel('M序列');

title('移位寄存器产生的双极性7位M序列') ; subplot(3,1,2) ym=fft(M,4096); magm=abs(ym); %求双极性m序列频谱 fm=(1:2048)*200/2048;

plot(fm,magm(1:2048)*2/4096);

title('双极性7位M序列的频谱') axis([90,140,0,0.1]);

[a,b]=xcorr(M,'unbiased'); subplot(3,1,3) %求双极性m序列自相关函数 plot(b,a);

axis([-20,20,-0.5,1.2]);

title('双极性7位M序列的自相关函数');

由上图可以看出，7位m序列为1，-1，-1，-1，1，-1，1。另外，自相关函数的图形比较尖锐，最大值为1，最小值为-1/7，符合理论结果。

2、生成50位随机待发送二进制比特序列，并进行扩频编码

生成的信息码频率为100/7Hz，利用m序列编码后，频率变为100Hz。

N=50;a=0;

x_rand=rand(1,N); %产生50个0与1之间随机数 for i=1:N

if x_rand(i)>=0.5 %大于等于0.5的取1，小于0.5的取0 x(i)=1;a=a+1; else x(i)=0; end end

- 6 -

t=0:N-1; figure(2) %做信息码图 subplot(2,1,1) stem(t,x);

title('扩频前待发送二进制信息序列'); tt=0:349; subplot(2,1,2) l=1:7*N; y(l)=0; for i=1:N k=7*i-6;

y(k)=x(i);

k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i); end s(l)=0;

for i=1:350 %扩频后，码率变为100/7*7=100Hz s(i)=xor(L(i),y(i)); end

tt=0:7*N-1; stem(tt,s);

axis([0,350,0,1]);

title('扩频后的待发送序列码');

- 7 -

3、对扩频前后信号进行BPSK调制，观察其时域波形

BPSK调制采用2kHz信号cos(2*2000*t)作为载波

figure(3)

subplot(2,1,2) fs=2000;

ts=0:0.00001:3.5-0.00001;%为了使信号看起来更光滑，作图时采样频率为100kHz % ps=cos(2*pi*fs*ts);

s_b=rectpulse(s,1000); %将冲激信号补成矩形信号

s_bpsk=(1-2.*s_b).*cos(2*pi*fs*ts);%扩频后信号BPSK调制时域波形 plot(ts,s_bpsk); xlabel('s');

axis([0.055,0.085,-1.2,1.2])

title('扩频后bpsk信号时域波形'); subplot(2,1,1)

s_bb=rectpulse(x,7000);

s_bpskb=(1-2.*s_bb).*cos(2*pi*fs*ts);%无扩频信号BPSK调制时域波形 plot(ts,s_bpskb); xlabel('s');

axis([0.055,0.085,-1.2,1.2]);

title('扩频前bpsk信号时域波形')

可以看出，100/7Hz的无扩频信号每0.07s时由于序列极性变换产生相位变换，100Hz的扩频后调制信号每0.01s由于序列极性变换产生相位变换。

- 8 -

4、计算并观察扩频前后BPSK调制信号的频谱

对信号采用400000点fft计算，得到频谱 figure(4) N=400000;

ybb=fft(s_bpskb,N); %无扩频信号BPSK调制频谱 magb=abs(ybb);

fbb=(1:N/2)*100000/N; subplot(2,1,1)

plot(fbb,magb(1:N/2)*2/N); axis([1700,2300,0,0.8]);

title('扩频前调制信号频谱'); xlabel('Hz'); subplot(2,1,2)

yb=fft(s_bpsk,N); %扩频信号BPSK调制频谱 mag=abs(yb);

fb=(1:N/2)*100000/N; plot(fb,mag(1:N/2)*2/N); axis([1700,2300,0,0.8]);

title('扩频后调制信号频谱'); xlabel('Hz');

如图，扩频前信号主瓣宽度约为2*100/7=28Hz，扩频后，信号频谱展宽，主瓣1900~2100Hz约为200Hz，为无扩频信号频谱宽度的N=7倍，符合理论推算。

- 9 -

5、仿真经awgn信道传输后，扩频前后信号时域及频域的变化

awgn信道模拟了真实的信道，为传输信号增加了高斯白噪声。在本次仿真中，设定信道信噪比为3dB，即信噪比约为2。

figure(5)

subplot(2,2,1)

s_bpskba=awgn(s_bpskb,3,'measured');%经过信道加高斯白噪,信噪比为3dbw plot(ts,s_bpskb,ts,s_bpskba); axis([0,0.005,-1.2,1.2]); xlabel('t');

title('经过信道加噪后的信号与原信号时域波形对比');

subplot(2,2,3)

s_bpska=awgn(s_bpsk,3,'measured'); plot(ts,s_bpsk,ts,s_bpska);

title('扩频后经加噪过信道后的信号与原信号时域波形对比'); xlabel('t');

axis([0.0675,0.0725,-1.2,1.2]);

subplot(2,2,2)

ybba=fft(s_bpskba,N); %无扩频调制信号经信道后频谱分析 magba=abs(ybba);

plot(fbb,magba(1:N/2)*2/N);

title('扩频前经信道调制信号频谱'); axis([1700,2300,0,0.8]); xlabel('Hz');

subplot(2,2,4)

yba=fft(s_bpska,N); %扩频调制信号经信道后频谱分析 maga=abs(yba);

fb=(1:N/2)*100000/N; plot(fb,maga(1:N/2)*2/N); axis([1700,2300,0,0.8]); xlabel('Hz');

title('扩频后经信道调制信号频谱');

- 10 -

10、观察经过低通滤波器后无扩频与扩频系统的时域波形

figure(10) subplot(2,1,1)

yrebl=real(ifft(bs.*yreb,400000)); %对无扩频系统频谱做ifft变换 tm=(1:N)/N*4; plot(tm,yrebl); xlabel('t');

title('扩频前经过凯萨尔窗函数滤波后时域波形'); subplot(2,1,2)

yrel=real(ifft(bs.*yre,400000)); %对扩频系统频谱做ifft变换 plot(tm,yrel); xlabel('t');

title('扩频后经过凯萨尔窗函数滤波后时域波形');

如图，经过低通滤波器后，高频分量基本消失，剩下的信号已经能够进行采样判决，时域波形与原信息基本吻合。

- 16 -

11、对扩频系统进行解扩，观察其时域频域

figure(11) subplot(2,1,1)

jj=rectpulse(M,1000); %扩频信号乘以解扩序列 yrej=jj.*yrel(1:350000); plot(ts(1:350000),yrej); xlabel('t');

axis([0,4,-0.5,0.5]);

title('解扩后信号波形'); subplot(2,1,2) yj=fft(yrej,N); magj=abs(yj);

plot(freb,magj(1:N/2)*2/N); axis([0,500,0,0.2]);

title('解扩后信号频谱'); xlabel('Hz');

由于扩频信号与m序列具有良好的相关性，故乘以m序列以后，能基本还原出原信号波形。同时可以看出，频谱已经由扩展带宽再次缩短，还原出原信号频谱。

- 17 -

12、比较扩频系统解扩前后信号带宽

figure(12)

title('解扩前后信号频偏对比'); subplot(2,1,1)

plot(freb,magrel(1:N/2)*2/N); axis([0,200,0,0.4]);

title('解扩前信号频偏'); subplot(2,1,2)

plot(freb,magj(1:N/2)*2/N); axis([0,200,0,0.4]);

title('解扩后信号频偏');

可以清楚看出，解扩前信号主瓣约为100Hz，解扩后恢复为100/7Hz，与发送信息吻合。

- 18 -

13、比较解扩前后信号功率谱密度

figure(13) subplot(2,1,1) yjb=fft(yrel,N);

prelb=yjb.*conj(yjb)/N; plot(freb,prelb(1:N/2)*2/N); axis([0,200,0,0.01]);

title('解扩前信号功率谱');xlabel('Hz'); subplot(2,1,2) yj=fft(yrej,N);

prel=yj.*conj(yj)/N;

plot(freb,prel(1:N/2)*2/N); axis([0,200,0,0.01]);

title('解扩后信号功率谱'); xlabel('Hz');

如图，解扩后信号的频谱被压缩，功率幅度增加，符合理论分析结果。

- 19 -

14、对解扩信号进行采样、判决

figure(14) subplot(2,1,1) for i=1:1:350

ij=i*1000-500; ss(i)=yrej(ij); end

stem(ss);

title('解扩信号采样'); subplot(2,1,2) for i=1:1:350 %判决信号算法 if ss(i)>0.2 ss(i)=1; elseif ss(i)<-0.2 ss(i)=-1; else ss(i)=0; end end

for i=1:1:50 ij=7*i-6; if ss(ij)==0

ss(ij)=ss(ij+4); end end

for i=1:1:348 if ss(i)==0

ss(i)=ss(i+2); end end

for i=1:1:50

S(i)=ss(i*7-3); if S(i)==0 S(i)=S(i)+1; end

S(i)=(1-S(i))/2; end stem(S);

title('判决后的最终信号');

17、比较解扩后信号与窄带强干扰的功率谱

figure(17)

prelnz=ynzj.*conj(ynzj)/N; prelz=yzj.*conj(yzj)/N;

plot(freb,prelnz(1:N/2)*2/N,freb,prelz(1:N/2)*2/N); axis([0,100,0,0.007]); xlabel('Hz');

title('信号与窄带干扰经解扩后的功率谱');

比较功率谱。蓝色部分为信号功率谱，绿色部分为窄带强干扰功率谱，可以发现窄带强干扰已经完全淹没在信号功率里。由理论计算可知，信噪比与N=1无扩频相比能够提高N2倍，抗干扰能力明显增强。

六、误码率simulink仿真

误码率仿真采用了MATLAB里面的SIMULINK及BERTOOL工具。SIMULINK是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。SIMULINK具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点SIMULINK已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于SIMULINK。

1、直接扩频系统信道模型

2、加窄带干扰的直扩系统建模

仿真过程中，仍然采用了100Hz的扩频序列，100/7Hz的数据码流，2kHz的BPSK调制并加入了与中心频点频偏20Hz的窄带强干扰。误码率仿真采用了SIMULIMK自带的Error Rate Calculation模块，来对比发送与接收的码流。建模的过程中，考虑了信号的频率，极性，窄带干扰和awgn信道的高斯白噪声，完全符合一般的通信系统的组成结构。仿真时间为100s。

3、用示波器观察发送码字及解扩后码字

上图为仿真过程中截取的部分发送与接收码字，上半部分为发送码字序列，下半部分为最终解扩后的码字序列。可以看出，大多数时间信号吻合，在第25.86s左后出现两个误码。

4、直接扩频系统与无扩频系统的误码率比较

采用SIMULINK下的BERTOOL工具，可以轻松地仿真出信号的误码率。在误码率计算中，我分别仿真了不同m序列长度和不同窄带干扰强度下，误码率与awgn高斯信道信噪比的关系图。

下图为无窄带干扰时，无扩频系统与N=7的m序列直接扩频BPSK系统的误码率比较，横坐标为信道信噪比。

可以看出，扩频序列误码率与无扩频系统的理论误码率基本相同，说明在没有窄带干扰的情况下扩频与否对于误码率影响不大。这说明，在干扰为高斯白噪声的情况下，扩频系统与无扩频系统的抗干扰能力相同。

- 30 -

5、不同扩频序列长度下的误码率比较

如图，扩频系统的误码率与扩频序列的长度有关，可以看出，在信噪比很小的情况下，不同扩频序列长度的误码率基本相同。当信噪比逐渐增大时，系统的误码情况有所好转。扩频码越长，误码率越低，抗干扰能力越强。

6、扩频序列长度N=7时，不同强度窄带干扰下的误码率比较

如图，当窄带强干扰的振幅为信号幅度的10倍时，扩频系统的误码率不变，说明扩频系统对窄带强干扰有很强的抑制能力，符合之前的matlab仿真结果。当窄带强干扰的振幅增至信号幅度的250倍和500倍时，系统误码率随干扰振幅增大而增大。窄带干扰幅度为信号幅度250倍时，误码率小于0.01，可以接受。当窄带干扰幅度为信号幅度的500倍时，误码率较大，不能忍受，但这种情况基本不会出现。

- 32 -

七、利用Walsh码实现码分多址技术

沃尔什码的自相关特性和互相关特性都不理想，这意味着CDM信号经过多径信道传输时，每个用户的不同径之间会产生严重干扰，同时不同用户不同径之间也会产生严重干扰。为此，可以用相关性较好的伪随机序列，即m序列与沃尔什码模二加，得到改善的沃尔什码。这样，既保留了沃尔什码的正交特性，又大大改善了其相关特性。

12.5Hz二进制比特信息1 12.5Hz二进制比特信息2 100Hz8位改善walsh码m1 100Hz8位改善walsh码m2 100Hz 扩频序列 2000Hz 载波cos4000πt BPSK调制信号 高斯白噪声 2040~2050Hz 窄带强干扰 恢复载波cos4000πt 凯萨尔滤波器低通滤波 100Hz8位改善walsh码1 采样、判决信息1 100Hz8位改善walsh码m2 采样、判决信息2 1、产生改善的walsh码

在本实验中，我采用了前面的7位m序列改善的沃尔什码。产生方法为在双极性m序列后加“1”，并与正交沃尔什码相乘，并用改善后的沃尔什码对两路信息序列分别进行扩频，相加后共同传输，实现码分多址技术。

由于m序列为7位，补位后0和1出现概率相等，成为了真正的随机序列。两路沃尔什函数采用了八阶沃尔什函数的的第一个和第五个，分别为1,1,1,1,1,1,1,1和1,1,1,1,0,0,0,0，为正交码。

clear all; clc;

X1=0;X2=0;X3=1; m=7; for i=1:m

Y3=X3; Y2=X2; Y1=X1; X3=Y2; X2=Y1; X1=xor(Y3,Y1); L(i)=Y1; end

L(8)=0; %单极性7位m序列末尾补0 m=m+1; for i=1:m

M(i)=1-2*L(i); end

k=1:1:m; figure(1)

subplot(3,1,1) stem(k-1,M); axis([0,9,-1,1]); xlabel('k');

title('移位寄存器产生的双极性7位M序列,末尾加一个1') ; subplot(3,1,2)

for i=1:m %生成改善的沃尔什码m1 m1(i)=M(i); end

stem(k-1,m1); axis([0,8,-1,1]); xlabel('k');

title('改善的沃尔什码m1') ; subplot(3,1,3)

for i=1:4 %生成改善的沃尔什码m2 m2(i)=M(i); end

for i=5:8

m2(i)=-M(i); end

stem(k-1,m2); axis([0,8,-1,1]); xlabel('k');

title('改善的沃尔什码m2') ;

如图，改善的沃尔什码m1与m2正交。下面将用这两路沃尔什码分别对两路信息进行扩频调制。

2、产生两路不同的信息序列

figure(2) N=50;a=0;

x_rand=rand(1,N); for i=1:N

if x_rand(i)>=0.5 x1(i)=1;a=a+1;

- 35 -

else x1(i)=0; end end

t=0:N-1;

subplot(2,1,1) stem(t*0.08,x1); xlabel('t/s');

title('扩频前待发送二进制信息序列1'); x_rand=rand(1,N); for i=1:N

if x_rand(i)>=0.5 x2(i)=1;a=a+1; else x2(i)=0; end end

subplot(2,1,2) stem(t*0.08,x2);

title('扩频前待发送二进制信息序列2'); xlabel('t/s');

- 36 -

3、用两个沃尔什码分别调制两路信号

figure(3)

y1=rectpulse(x1,8);y2=rectpulse(x2,8); for i=1:1:50

M1(8*i-7)=m1(1);M2(8*i-7)=m2(1); M1(8*i-6)=m1(2);M2(8*i-6)=m2(2); M1(8*i-5)=m1(3);M2(8*i-5)=m2(3); M1(8*i-4)=m1(4);M2(8*i-4)=m2(4); M1(8*i-3)=m1(5);M2(8*i-3)=m2(5); M1(8*i-2)=m1(6);M2(8*i-2)=m2(6); M1(8*i-1)=m1(7);M2(8*i-1)=m2(7); M1(8*i)=m1(8);M2(8*i)=m2(8); end

subplot(2,1,1) tt=0:8*N-1; for i=1:400;

x11(i)=1-2*y1(i);x12(i)=1-2*y2(i);s1(i)=M1(i)*x11(i);s2(i)=M2(i)*x12(i); end

subplot(2,1,1) stem(tt/100,s1);

title('扩频后的待发送序列码1'); xlabel('t/s'); subplot(2,1,2) stem(tt/100,s2);

title('扩频后的待发送序列码2'); xlabel('t/s');

4、两路信号相加，并进行BPSK调制 figure(4)

subplot(2,1,1) s=s1+s2;

stem(tt/100,s);

title('两路信息相加,即进行walsh码复用'); xlabel('t/s'); subplot(2,1,2) fs=2000;

ts=0:0.00001:4-0.00001; s_b=rectpulse(s,1000);

s_bpsk=s_b.*cos(2*pi*fs*ts); plot(ts,s_bpsk); xlabel('t/s');

axis([0.065,0.095,-2.2,2.2])

title('walsh码分复用做BPSK变换');

对比walsh码进行复用后的信号与无码分多址系统的信号，可以看出，无码分多址系统的信息码只有两个取值-1和1，但是复用后存在三个值-1,0,1，所以BPSK调制信号存在为零的时刻。

5、观察调制信号频谱，并经awgn信道加高斯白噪和窄带强干扰

figure(5)

subplot(3,1,1) N=400000;

yb=fft(s_bpsk,N); mag=abs(yb);

fb=(1:N/2)*100000/N; plot(fb,mag(1:N/2)*2/N); axis([1000,3000,0,0.25]);

title('码分复用扩频后调制信号频谱'); xlabel('f/Hz'); subplot(3,1,2)

s_bpska=awgn(s_bpsk,3,'measured');%经过awgn信号，信噪比3dB plot(ts,s_bpsk,ts,s_bpska);

title('码分复用扩频后经加噪过信道后的信号与原信号时域波形对比'); xlabel('t/s');

axis([0.0775,0.0825,-2.4,2.4]); subplot(3,1,3)

yba=fft(s_bpska,N); maga=abs(yba);

fb=(1:N/2)*100000/N; fd=200000;

Wp1=2*2040/fd;%%%%%%%%%%带通滤波器 Wp2=2*2050/fd; Wc1=2*2030/fd; Wc2=2*2060/fd; Ap=1; As=100;

W1=(Wp1+Wc1)/2; W2=(Wp2+Wc2)/2;

wdth=min((Wp1-Wc1),(Wc2-Wp2)); Nd=ceil(11*pi/wdth)+1; bd=fir1(Nd,[W1 W2]); zd(1)=1;%单位冲激 for i=2:1:350000 zd(i)=0; end

ds=abs(freqz(bd,1,400000,fd))'; ybz=fft(zd,N)*40000; magz=abs(ybz); dz=ds.*magz;

dsz=maga+dz;%将窄带干扰叠加到信号上

plot(fb,dsz(1:N/2)*2/N); axis([1000,3000,0,0.4]); xlabel('f/Hz');

title('码分复用扩频后经信道调制加窄带强干扰信号频谱');

- 40 -

本文来源：https://www.bwwdw.com/article/j0wr.html