FDMA通信系统设计.doc

FDMA通信系统的设计一、实验目的 1、加深对调制与解调原理及过程的理解，能将调制与解调原理应用到FDMA通信系统中；2、加深对滤波器滤波特性的理解。 3、掌握FDMA通信系统的原理，并利用Matlab进行实现。 二、实验原理 1、频分复用的原理 复用是指将若干个彼此独立的信号合并成可在同一信道上传输的复合信号的方法，常见的信号复用采用按频率区分与按时间区分的方式，前者称为频分复用，后者称为时分复用。 通常在通信系统中，信道所提供的带宽往往比传输一路信号所需要的带宽宽得多，这样就可以将信道的带宽分割成不同的频段，每频段传输一路信号，这就是频分复用（frequency division multiple access）(FDMA)。为此，在发送端首先要对各路信号进行调制将其频谱函数搬移到相应的频段内，使之互不重叠。再送入信道一并传输。在接收端则采用不同通带的带通滤波器将各路信号分隔，然后再分别解调，恢复各路信号。调制的方式可以任意选择，但常用的是单边带调制。因为每一路信号占据的频段小，最节省频带，在同一信道中传送的路数可以增加。图1 频分复用系统的示意图图1给出了频分复用系统的示意图。如图所示，其中f1(t),f2(t),fn(t)为n路低频信号，通过调制器形成各路处于不同频段上的边带信号。频分复用的理论基础仍然是调制和解调。通常为防止邻路信号的相互干扰，相邻两路间还要留有防护频带，因此各路载频之间的间隔应为每路信号的频带与保护频带之和。以语音信号为例，其频谱一般在0.33.4kHz范围内，防护频带标准为900Hz，则每路信号占据频带为4.3kHz，以此来选择相应的各路载频频率，在接收端则用带通滤波器将各路信号分离再经同步检波即可恢复各路信号，为减少载波频率的类型，有时也用二次调制。频分复用技术除传统意义上的频分复用（FDMA）外，还有一种是正交频分复用（OFDM）。（1）传统的频分复用传统的频分复用典型的应用莫过于广电HFC网络电视信号的传输了，不管是模拟电视信号还是数字电视信号都是如此，因为对于数字电视信号而言，尽管在每一个频道（8 MHz）以内是时分复用传输的，但各个频道之间仍然是以频分复用的方式传输的。（2）正交频分复用OFDM（Orthogonal Frequency Division Multiplexing）实际是一种多载波数字调制技术。OFDM全部载波频率有相等的频率间隔，它们是一个基本振荡频率的整数倍，正交指各个载波的信号频谱是正交的。OFDM系统比FDMA系统要求的带宽要小得多。由于OFDM使用无干扰正交载波技术，单个载波间无需保护频带，这样使得可用频谱的使用效率更高。另外，OFDM技术可动态分配在子信道中的数据，为获得最大的数据吞吐量，多载波调制器可以智能地分配更多的数据到噪声小的子信道上。目前OFDM技术已被广泛应用于广播式的音频和视频领域以及民用通信系统中，主要的应用包括：非对称的数字用户环线（ADSL）、数字视频广播（DVB）、高清晰度电视（HDTV）、无线局域网（WLAN）和第4代（4G）移动通信系统等。频分复用系统最大的优点是信道复用率高，允许的复用路数较多，同时分路也很方便，是模拟通信中主要的一种复用方式，在有线和微波通信中应用十分广泛。频分复用的缺点是设备生产较为复杂，同时因滤波性能不够理想，及信道内存在的非线性容易产生路间干扰。2、FDMA通信系统的原理 FDMA通信系统模型如图8-2所示。同学们还可了解一下WDMA。WDMA和FDMA基本上都基于相同原理，所不同的是，WDMA应用于光纤信道上的数字化光波传输过程，而FDMA应用于模拟传输，诸如双绞线话路传输、电缆接入、峰窝、无线电以及 TV 通信等。一直以来， TDMA 、CDMA 也是结合 FDMA 共同作用的。图8-2 FDMA通信系统模型三、实验涉及的Matlab函数 1. fft():对信号进行谱分析。 2. ifft():对信号进行傅里叶反变换。 3. fir1():窗函数法设计FIR数字滤波器。 4. fir2():频率采样法设计FIR数字滤波器5. buttord():设计巴特沃思型的IIR数字滤波器。还有cheb1ord()、cheb2ord()、ellipord()设计切比雪夫和椭圆型的IIR数字滤波。 6.filter()：IIR数字滤波器实现滤波。7. fftfilt():FIR数字滤波器实现滤波。8.y,fs,nbits=wavread(file)：采样值放在向量y中，fs表示采样频率（Hz），nbits表示采样位数。y=wavread(file,N)：读取前N点的采样值放在向量y中。9. sound(x,fs,bits)：将x的数据通过声卡转化为声音。四、实验内容 在Matlab 环境中，利用编程方法对FDMA通信模型进行仿真研究。1、设计要求： （1）Matlab支持麦克风，可直接进行声音的录制，要求至少获取3路语音信号。（2）对各路语音信号进行频谱分析。（3）将各路语音信号分别与各自的高频载波信号相乘，由于各高频载波信号将各语音信号频谱移到不同频段，复用信号频谱为各信号频谱的叠加，因此，只需传输该复用信号便可在同一信道上实现各路语音信号的同时传输。（4）传输完成后，通过选择合适的带通滤波器，即可获得各个已调信号。（5）再进行解调，即将各个已调信号分别乘以各自的高频载波信号，这样，原始低频信号被移到低频段。 （6）最后通过选择合适的低通滤波器恢复出各原始语音信号，从而实现FDMA通信传输。2、实验代码1、获取录音文件fs=44100; %声音的采样频率为44.1khzduration=3; %录音的时间fprintf(按任意键开始录音1:n);pausefprintf(录音中n);sd1=wavrecord(duration*fs,fs); %duration*fs每次获得总的采样数为132300,保存声音fprintf(放音中n); %文件名为s1,以下类同.wavplay(sd1,fs);fprintf(录音1播放完毕。n);wavwrite(sd1,fs,sound1.wav); %将录音文件保存为wav格式的声音文件，下同fprintf(按任意键开始录音2:n);pausefprintf(录音中n);sd2=wavrecord(duration*fs,fs);fprintf(放音中n);wavplay(sd2,fs);fprintf(录音2播放完毕。n);wavwrite(sd2,fs,sound2.wav);fprintf(按任意键开始录音3:n);pausefprintf(录音中n);sd3=wavrecord(duration*fs,fs);fprintf(放音中n);wavplay(sd3,fs);fprintf(录音3播放完毕。n);wavwrite(sd3,fs,sound3.wav); 三个声音的时域波形2、声音样本的时域和频域分析fs=44100; %声音的采样频率为44.1khzduration=3; t=0:duration*fs-1; %总的采样数s1,fs=wavread(sound1.wav); %打开保存的录音文件s2,fs=wavread(sound2.wav);s3,fs=wavread(sound3.wav);figure(1) %图一为三个声音样本的时域波形subplot(311)plot(t,s1);xlabel(单位：s);ylabel(幅度);title(三个声音样本的时域波形);subplot(312)plot(t,s2);xlabel(单位：s);ylabel(幅度);subplot(313)plot(t,s3);xlabel(单位：s);ylabel(幅度);figure(2) %图二为三个声音样本的频谱分析subplot(311)stem(t,abs(fft(s1),.);xlabel(单位:Hz);ylabel(幅度); %fft对声音信号进行快速傅里叶变换title(三个声音样本的频谱分析);subplot(312)stem(t,abs(fft(s2),.);xlabel(单位:Hz);ylabel(幅度);subplot(313)stem(t,abs(fft(s3),.);xlabel(单位:Hz);ylabel(幅度);三个声音的频谱分析3、调制,将三个声音信号用高频载波x1=4*s1.*cos(2*pi*4000*t/fs);x2=4*s2.*cos(2*pi*11000*t/fs);x3=4*s3.*cos(2*pi*18000*t/fs);s=x1+x2+x3;%复用信号频谱为各信号频谱的叠加figure(3)%图三为复用信号的频谱分析stem(t,abs(fft(s),.);xlabel(单位:Hz);ylabel(幅度);title(复用信号的频谱分析);复用信号的频谱分析4、带通滤波器设计Rp=0.5;Rs=40; %用切比雪夫设计带通滤波器1； Wp1=4000 8000/22050; %fs/2=22050 Ws1=3800 8500/22050;n1,Wn1=cheb2ord(Wp1,Ws1,Rp,Rs);b1,a1=cheby2(n1,Rs,Wn1);h1,w1=freqz(b1,a1);mag1=abs(h1);db1=20*log10(mag1+eps)/max(mag1);Wp2=9000 13000/22050; %用切比雪夫设计带通滤波器2； Ws2=8000 14000/22050;n2,Wn2=cheb2ord(Wp2,Ws2,Rp,Rs);b2,a2=cheby2(n2,Rs,Wn2);h2,w2=freqz(b2,a2);mag2=abs(h2);db2=20*log10(mag2+eps)/max(mag2);Wp3=14500 18500/22050; %用切比雪夫设计带通滤波器3； Ws3=14000 19000/22050;n3,Wn3=cheb2ord(Wp3,Ws3,Rp,Rs);b3,a3=cheby2(n3,Rs,Wn3);h3,w3=freqz(b3,a3);mag3=abs(h3);db3=20*log10(mag3+eps)/max(mag3);figure(4);subplot(3,1,1);plot(w1/pi,db1);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|);title(用切比雪夫2型设计三个带通滤波器);subplot(3,1,2);plot(w2/pi,db2);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|);subplot(3,1,3);plot(w3/pi,db3);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|);y1=filter(b1,a1,s); %滤出三路未解调信号y2=filter(b1,a1,s);y3=filter(b1,a1,s); 利用切比雪夫设计的3个带通滤波器5、解调fs=44100y01=y1.*cos(2*pi*4000*t/fs); %各个已调信号分别乘以各自y02=y2.*cos(2*pi*11000*t/fs);%的高频载波信号y03=y3.*cos(2*pi*18000*t/fs);figure(5) %图五为解调后3路信号各自的频谱图subplot(311)stem(t,abs(fft(y01),.);xlabel(单位:Hz);ylabel(幅度);title(解调后3路信号各自的频谱图);subplot(312)stem(t,abs(fft(y02),.);xlabel(单位:Hz);ylabel(幅度);subplot(313)stem(t,abs(fft(y03),.);xlabel(单位:Hz);ylabel(幅度); 解调后的信号频谱分析6、低通滤波Rp=0.5; %低通滤波器参数选择Rs=40;Wp1=3400/(22050);Ws1=4000/(22050);n1,Wn1=cheb2ord(Wp1,Ws1,Rp,Rs); %采用切比雪夫2型(cheby2)带通滤波器b1,a1=cheby2(n1,Rs,Wn1);h1,w1=freqz(b1,a1);mag1=abs(h1);db1=20*log10(mag1+eps)/max(mag1);figure(6) %图六为低通滤波器的频率响应plot(w1/pi,db1);axis(0 1 -50 20);xlabel(w/pi);ylabel(20lg|H(ejw)|);title(低通滤波器的频率响应); 低通滤波器的频率响应7、回复信号的时域波形和频谱分析yy1=filter(b1,a1,y01);yy2=filter(b1,a1,y02);yy3=filter(b1,a1,y03);figure(7) %图七为恢复信号的时域波形subplot(311)plot(t,yy1);xlabel(单位：s);ylabel(幅度);subplot(312)plot(t,yy2);xlabel(单位：s);ylabel(幅度);subplot(313)plot(t,yy3);xlabel(单位：s);ylabel(幅度); title(恢复信号的时域波形);figure(8) %图八为恢复信号的频谱分析subplot(311)stem(t,abs(fft(yy1);xlabel(单位:Hz);ylabel(幅度);subpl