语音信号调制解调实验报告.doc

大连理工大学实验报告学院（系）： 专业： 班级： 姓 名： 学号： 组： _ 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验五 语音信号的调制解调一、实验目的和要求（略去） 二、实验原理和内容（略去）三、实验程序和结果按照讲义中的实验要求，按步骤实现信号的调制解调。信息信号用3Hz的正弦信号，载波信号用50Hz的余弦信号代码如下：clccleardt=0.001;fs=1000;t=-1:dt:1;f1=sin(6*pi*t);figure;subplot(3,2,1)fy1 = fft(f1);t1=(0:length(t)-1)*fs/length(t)-fs/2;plot(t1,fftshift(abs(fy1);title(信息信号的频谱);subplot(3,2,2);plot(t,f1);title(信息信号的时域波形);f2=cos(100*pi*t);subplot(3,2,3)fy2 = fft(f2);t2=(0:length(t)-1)*fs/length(t)-fs/2;plot(t2,fftshift(abs(fy2);title(载波信号的频谱);subplot(3,2,4);plot(t,f2);title(载波信号的时域波形);f3=(2+f1).*f2;subplot(3,2,5)fy3 = fft(f3);t3=(0:length(t)-1)*fs/length(t)-fs/2;plot(t3,fftshift(abs(fy3);title(已调信号的频谱);subplot(3,2,6);plot(t,f3);title(已调信号的时域波形);y0=rand(1,2001);y1=y0+f3;figure;subplot(2,1,1);yy1=fft(y1);yt1=(0:length(t)-1)*fs/length(t)-fs/2;plot(yt1,fftshift(abs(yy1);title(添加噪声信号的已调波频谱);subplot(2,1,2);plot(t,y1);title(添加噪声信号的已调波时域波形);a=45,55;b=40,60;Wp=a/(fs/2);Ws=b/(fs/2);Rp=3; Rs=15;N,Wn= Buttord(Wp,Ws,Rp,Rs);B,A=Butter(N,Wn,bandpass);q=filtfilt(B,A,y1);figure;subplot(2,1,1);qq=fft(q);qt=(0:length(t)-1)*fs/length(t)-fs/2;plot(qt,fftshift(abs(qq);title(已调波带通滤波后的频谱);subplot(2,1,2);plot(t,q);title(已调波带通滤波后的时域波形);y2=q.*cos(100*pi*t);figure;subplot(2,1,1);yy2=fft(y2);yt2=(0:length(t)-1)*fs/length(t)-fs/2;plot(yt2,fftshift(abs(yy2);title(经过乘法器后的频谱);Wp=5/(fs/2);Ws=10/(fs/2);Rp=3; Rs=20;N,Wn= Buttord(Wp,Ws,Rp,Rs) ;B,A=Butter(N,Wn,low);m0=filtfilt(B,A,y2);subplot(2,1,2);mm=fft(m0);mt=(0:length(t)-1)*fs/length(t)-fs/2;plot(mt,fftshift(abs(mm);title(低通滤波后的频谱);figure;plot(t,2*m0-2);axis(-1,1,-1,1);title(最后解调出来的信号);出现的matlab图形如下：按顺序依次为：接着添加噪声信号，分别画出添加噪声信号后的频谱以及其时域波形：然后画出经过带通滤波后的频谱以及时域波形：对信号进行解调，先经过乘法器，再低通滤波，画出两个步骤后的频谱：低通滤波后，查看解调后的信号，观察其时域波形如下：可以看到解调出来的信号有微小的误差，但是大致上能反映出原有的信息信号，可知经过调制解调后的信号是成功的。四、实验体会本次实验用到的东西比较多，综合了前面的所学以及新的知识的参杂来实现信号的调制解调。做完本次实验有如下的问题：滤波器通带和阻带的设置不是很明确，感觉通带和阻带的设置如果不