fdm 等Word下载.docx
《fdm 等Word下载.docx》由会员分享,可在线阅读,更多相关《fdm 等Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
画出复用信号的频谱图。
(3)传输完成后,通过选择合适的带通滤波器,即可获得各个已调信号;
给定滤波器的性能指标,设计合适数字滤波器,并画出带通滤波器的频率响应。
(4)再进行解调,即将各个已调信号分别乘以各自的高频载波信号,这样,原始低频信号被移到低频段。
画出解调后3路信号各自的频谱图。
(5)最后通过选择合适的低通滤波器恢复出各原始语音信号,从而实现FDMA通信传输。
画出低通滤波器的频率响应,恢复信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化。
回放语音信号。
三、详细设计过程
(1)录音模块:
首先,Matlab中支持麦克风,因此可直接录制三段语音信号,并保存和播放。
由时域采样定理知采样频率要大于信号高频率的两倍,而且调制时要将音频调到不同的频段,所以选取采样频率时,可选高一点,便于调制。
可选fs为44100Hz,为CD音质。
(2)调制与解调模块:
首先设计载波频率,各个载频的间隔既要大于2倍的声音的最高截止频率,而且最高的载频与采样频率也必须大于2倍的声音的最高截止频率。
因为fs为44100Hz,所以选择载波频率时,最好不超过fs/2,否则Ω+ω容易与ω+fs-Ω相混淆(ω为载频,Ω为声音最高截止频率,f为采样频率)。
综合考虑三个载频分别取为4000Hz,11000Hz,18000Hz。
然后将三路调制信号相加得到复用信号,并对其进行频谱分析。
接着经过带通滤波出三路信号后,再解调并分析各路信号的频谱。
(3)滤波器模块:
为达到较好效果,采用切比雪夫2型滤波器,并且取Rp=0.5;
Rs=40;
三个带通滤波器分别要滤出三路信号,其通频带要依据载频而定,可以滤出上边频,也可以滤出下边频。
综合最后的试听效果,分别同频带取4300~7400Hz,7600~11000Hz,14600~18000Hz。
由于三个信号都是人的声音信号,有着相同的频率范围,所以可以共用一个低通滤波器,且低通滤波器的通带范围为0~3400Hz。
四、代码及图形
%
(1)获取录音文件
fs=44100;
%声音的采样频率为44.1khz
duration=3;
%录音的时间
fprintf('
按任意键开始录音1:
\n'
);
pause
录音中……\n'
sd1=wavrecord(duration*fs,fs);
%duration*fs每次获得总的采样数为132300,保存声音
放音中……\n'
%文件名为s1,以下类同.
wavplay(sd1,fs);
录音1播放完毕。
wavwrite(sd1,fs,'
sound1.wav'
%将录音文件保存为wav格式的声音文件,下同
按任意键开始录音2:
sd2=wavrecord(duration*fs,fs);
wavplay(sd2,fs);
录音2播放完毕。
wavwrite(sd2,fs,'
sound2.wav'
按任意键开始录音3:
sd3=wavrecord(duration*fs,fs);
wavplay(sd3,fs);
录音3播放完毕。
wavwrite(sd3,fs,'
sound3.wav'
%
(2)声音样本的时域和频域分析
t=0:
duration*fs-1;
%总的采样数
[s1,fs]=wavread('
%打开保存的录音文件
[s2,fs]=wavread('
[s3,fs]=wavread('
figure
(1)%图一为三个声音样本的时域波形
subplot(3,1,1)
plot(t,s1);
xlabel('
单位:
s'
ylabel('
幅度'
title('
三个声音样本的时域波形'
subplot(3,1,2)
plot(t,s2);
subplot(3,1,3)
plot(t,s3);
figure
(2)%图二为三个声音样本的频谱分析
stem(t,abs(fft(s1)),'
.'
单位:
Hz'
%fft对声音信号进行快速傅里叶变换
三个声音样本的频谱分析'
stem(t,abs(fft(s2)),'
stem(t,abs(fft(s3)),'
%(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)),'
复用信号的频谱分析'
%(4)带通滤波器设计
Rp=0.5;
%用切比雪夫2型设计带通滤波器1;
Wp1=[40008000]/22050;
%fs/2=22050
Ws1=[38008500]/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=[900013000]/22050;
%用切比雪夫2型设计带通滤波器2;
Ws2=[800014000]/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=[1450018500]/22050;
%用切比雪夫2型设计带通滤波器3;
Ws3=[1400019000]/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([01-5020]);
w/pi'
20lg|H(ejw)|'
用切比雪夫2型设计三个带通滤波器'
subplot(3,1,2);
plot(w2/pi,db2);
subplot(3,1,3);
plot(w3/pi,db3);
y1=filter(b1,a1,s);
%滤出三路未解调信号
y2=filter(b2,a2,s);
y3=filter(b3,a3,s);
%(5)解调
fs=44100
y01=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路信号各自的频谱图
stem(t,abs(fft(y01)),'
解调后3路信号各自的频谱图'
stem(t,abs(fft(y02)),'
stem(t,abs(fft(y03)),'
%(6)低通滤波
%低通滤波器参数选择
Wp1=3400/(22050);
Ws1=4000/(22050);
%采用切比雪夫2型(cheby2)带通滤波器
figure(6)%图六为低通滤波器的频率响应
低通滤波器的频率响应'
%(7)回复信号的时域波形和频谱分析
yy1=filter(b1,a1,y01);
yy2=filter(b1,a1,y02);
yy3=filter(b1,a1,y03);
figure(7)%图七为恢复信号的时域波形
subplot(3,1,1)
plot(t,yy1);
plot(t,yy2);
plot(t,yy3);
恢复信号的时域波形'
figure(8)%图八为恢复信号的频谱分析
stem(t,abs(fft(yy1)));
stem(t,abs(fft(yy2)));
stem(t,abs(fft(yy3)));
恢复信号的频谱分析'
wavplay(yy1,fs);
%恢复声音信号的再现
wavplay(yy2,fs);
wavplay(yy3,fs);
五、结果分析与体会
结果分析:
经过不断的修改调试,在MATLAB上仿真频分多址通信技术取得了较好的效果。
录音的声音在经过调试和解调后的信号与原来相比较为接近。
我觉得仿真的成功关键在于载波频率的选择以及带通和低通滤波器的参数设置。
(详细设计中已有介绍,不在赘述)
另外在低通滤波阶段,得到的恢复信号与原始信号基本一致,但在t=0附近有所失真,这是由于频谱混叠所致,各信号频谱混叠部分均为高频部分,恢复信号在附近的波峰变化最快。
即为频率最高的区域,引起高频部分失真,这是因为录音期间引入频率高于语音信号的噪声,所以如果在完全无噪音的环境中进行录音,可得无失真的恢复信号。
最后仿真结果分析表明,信号在频分复用时还存在着频间干扰的问题。
对此,采用了适当加大采样频率的方法,在较大程度打夯使该问题得以解决。
至于完全消除频谱间的干扰,还有待进一步研究的研究与完善。
体会:
这次课程设计中,我不仅复习巩固了课堂所学的理论知识,提高了对所学知识的综合应用能力,并从实践上初步实现了对数字信号的处理。
而且通过这次较综合的课程设计,对以前在课本中问题没有弄懂的问题,我都有了更深入的理解。
比如时域采样定理、滤波器参数设计等。
总的来说,通过数字信号处理的这次课程设计,让我对《信号与系统》这门课程有了深刻的认识,收获很大。
上学期我们专注于信号的学习,涉及到的复杂的数学计算分析一度让我失去学习《信号与系统》的信心和兴趣。
下学期的系统方面或者说数字信号处理,以及通信相关知识的学习,我重新认识了信号处理这门课的地位,更重要的是对于将来学习通信技术我也增强了很多信心!
六、附录或参考资料
1)频分多址接入模型设计及MATLAB仿真计算
作者:
陈慧慧,郑宾中北大学
2)现代通信系统分析与仿真MATLAB通信工具箱
李建新,刘乃安,刘继平编西安电子科技大学出版社2000
3)MATLAB通信仿真及应用实例详解
邓华等编著人民邮电出版社2003