FDMA频分多址通信课程设计

1、二九二一学年第 一 学期信息科学与工程学院课程设计报告书课程名称： 数字信号处理课程设计 班 级： 电子信息工程2007级 6 班 学 号： 姓 名： 张琦 指导教师： 李娟二一年 01 月一、课程设计目的综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB作为编程工具进行计算机实现，从而复习巩固了课堂所学的理论知识，提高了对所学知识的综合应用能力，并从实践上初步实现了对数字信号的处理。二、课程设计内容及要求设计内容：在Matlab 环境中，利用编程方法对FDMA通信模型进行仿真研究。载波信号1载波信号n载波信号n频分多址（FDMA）通信系统模型设计

2、要求：（1）Matlab支持麦克风，可直接进行声音的录制，要求录制3路不同人的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图。（2）将各路语音信号分别与各自的不同的高频载波信号相乘，由于各高频载波信号将各语音信号频谱移到不同频段，复用信号频谱为各信号频谱的叠加，因此，只需传输该复用信号便可在同一信道上实现各路语音信号的同时传输。画出复用信号的频谱图。（3）传输完成后，通过选择合适的带通滤波器，即可获得各个已调信号；给定滤波器的性能指标，设计合适数字滤波器，并画出带通滤波器的频率响应。（4）再进行解调，即将各个已调信号分别乘以各自的高频载波信号，这样，原始低频信号被移到低

3、频段。画出解调后3路信号各自的频谱图。（5）最后通过选择合适的低通滤波器恢复出各原始语音信号，从而实现FDMA通信传输。画出低通滤波器的频率响应，恢复信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化。回放语音信号。三、详细设计过程(1)录音模块：首先，Matlab中支持麦克风，因此可直接录制三段语音信号，并保存和播放。由时域采样定理知采样频率要大于信号高频率的两倍，而且调制时要将音频调到不同的频段，所以选取采样频率时，可选高一点，便于调制。可选fs为44100Hz, 为CD音质。(2)调制与解调模块：首先设计载波频率，各个载频的间隔既要大于2倍的声音的最高截止频率,而且最高的载

4、频与采样频率也必须大于2倍的声音的最高截止频率。因为fs为44100Hz，所以选择载波频率时，最好不超过fs/2，否则+ 容易与+fs-相混淆（为载频，为声音最高截止频率，f为采样频率）。综合考虑三个载频分别取为4000Hz,11000Hz,18000Hz。然后将三路调制信号相加得到复用信号，并对其进行频谱分析。接着经过带通滤波出三路信号后，再解调并分析各路信号的频谱。(3)滤波器模块：为达到较好效果，采用切比雪夫2型滤波器，并且取Rp=0.5;Rs=40; 三个带通滤波器分别要滤出三路信号，其通频带要依据载频而定，可以滤出上边频，也可以滤出下边频。综合最后的试听效果，分别同频带取430074

5、00Hz，760011000Hz，1460018000Hz。由于三个信号都是人的声音信号，有着相同的频率范围，所以可以共用一个低通滤波器，且低通滤波器的通带范围为03400Hz。四、代码及图形%(1)获取录音文件fs=44100; %声音的采样频率为44.1khzduration=3; %录音的时间fprintf('按任意键开始录音1:n');pausefprintf('录音中n');sd1=wavrecord(duration*fs,fs);fprintf('放音中n');wavplay(sd1,fs);fprintf('录音1播放完毕

6、。n');wavwrite(sd1,fs,'sound1.wav');fprintf('按任意键开始录音2:n');pausefprintf('录音中n');sd2=wavrecord(duration*fs,fs);fprintf('放音中n');wavplay(sd2,fs);fprintf('录音2播放完毕。n'); %将录音文件保存为wav格式的声音文件，下同 %duration*fs每次获得总的采样数为132300,保存声音 %文件名为s1,以下类同.wavwrite(sd2,fs,'so

7、und2.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=wavr

8、ead('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,

9、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');ylabe

10、l('幅度');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),&

11、#39;.');xlabel('单位:Hz');ylabel('幅度');title('复用信号的频谱分析');%(4)带通滤波器设计Rp=0.5;Rs=40; %用切比雪夫2型设计带通滤波器1； Wp1=4000 8000/22050; %fs/2=22050Ws1=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);W

12、p2=9000 13000/22050; %用切比雪夫2型设计带通滤波器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; %用切比雪夫2型设计带通滤波器3； Ws3=14000 19000/22050;n3,Wn3=cheb2ord(Wp3,Ws3,Rp,Rs);b3,a3=cheby2(n3,Rs,Wn3);h3,

13、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(ej

14、w)|');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);%(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)subplot(311

15、) %各个已调信号分别乘以各自 %的高频载波信号 %图五为解调后3路信号各自的频谱图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('单位

16、:Hz');ylabel('幅度');%(6)低通滤波Rp=0.5;Rs=40;Wp1=3400/(22050);Ws1=4000/(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);figure(6) %图六为低通滤波器的频率响应 %低通滤波器参数选择 %采用切比雪夫2型(cheby2)带通滤波器 plot(w1/pi,db1);axis(0 1 -50 20);xl

17、abel('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'

18、;);ylabel('幅度');subplot(313)plot(t,yy3);xlabel('单位：s');ylabel('幅度');title('恢复信号的时域波形');figure(8) %图八为恢复信号的频谱分析 subplot(311)stem(t,abs(fft(yy1);xlabel('单位:Hz');ylabel('幅度');subplot(312)stem(t,abs(fft(yy2);xlabel('单位:Hz');ylabel('幅度');sub

19、plot(313)stem(t,abs(fft(yy3);xlabel('单位:Hz');ylabel('幅度');title('恢复信号的频谱分析');wavplay(yy1,fs);wavplay(yy2,fs);wavplay(yy3,fs); %恢复声音信号的再现五、结果分析与体会 结果分析：经过不断的修改调试，在MATLAB上仿真频分多址通信技术取得了较好的效果。录音的声音在经过调试和解调后的信号与原来相比较为接近。我觉得仿真的成功关键在于载波频率的选择以及带通和低通滤波器的参数设置。（详细设计中已有介绍，不在赘述）另外在低通滤波阶段，得到的恢复信号与原始信号基本一致，但在t=0附近有所失真，这是由于