实验报告二高.docVIP

广州大学学生实验报告 开课学院及实验室：电子信息楼4月7日 学院 机械与电气学院 年级、专业、班 姓名 学号 实验课程名称 多媒体通信技术 成绩 实验项目名称 实验二： 综合所学理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用MATLAB作为变成工具进行计算机实现，从而加深对所学知识的理解。 二、实验要求 熟悉离散信号和系统的时域特性。 熟悉语音信号的特点。 掌握序列快速傅立叶变换方法。 学会MATLAB的使用，掌握MATLAB的程序设计方法。 掌握MATLAB设计各种数字滤波器的方法和对信号进行滤波的方法。 三、 多重回声滤波器，滤波器的系统函数为 无限个回声滤波器，滤波器的系统函数为 全通结构的混响器 4 用滤波器对信号进行滤波 用自己设计的滤波器对采集的信号进行滤波。在MATLAB中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。 5 分析得到的频谱 对得到的信号进行fft快速傅立叶变换，得到信号的频谱特性，与原始信号进行比较，并画出信号的时域波形和频谱。 6回放语音信号 利用sound函数可以对语音进行回放，听效果看看有没有达到加入延时和混响的目的 编写代码如下: n=200; y=wavread(1.wav); yy=fft(y,4096); figure subplot(121); plot(y);title(原信号时域波形’) subplot(122);plot(abs(yy));title(波形) y1=[zeros(n,2);y]; y2=[y;zeros(n,2)]; y3=y1+y2; y1=fft(y3); y1f=fft(y1,4096); figure subplot(121);plot(y3);title(时域波形)%wav of time delay subplot(122);plot(abs(y1));title(的频域波形)%wav of frequency domain after delay %[y,fs,bits]=wavread(1.wav); y=y(:,1); m=10000;%用来调整回声的间隔a=0.6;%用来调整回声的强度b=[1;zeros(m,1);a]; a=[1]; y2=filter(b,a,y); %wav of time domain through single echo filter figure y3=fft(y2); subplot(121);plot(y2) title(单回声滤波后的时域波形) subplot(122);plot(abs(y3)) title(的域波形)%wave of frequency domain through single echo filter sound (y2,fs,bits); %滤波器clear all [x,fs,bits]=wavread(1.wav); x=x(:,1); p=200;%用来调整回声的间隔a=0.6%用来调整回声的音量b=[zeros(p,1);1]; a=[1;zeros(p-1,1);-a]; y3=filter(b,a,x); y4=fft(y3); %wave of time domain infinite single echo filter figure subplot(121);plot(y3) title(无限个回声滤波的时域波形); subplot(122);plot(abs(y4)) title(滤波的波形);%wave of frequency domain through infinite echo filter sound(y3,fs,bits); %clear all [x,fs,bits]=wavread(1.wav); x=x(:,1);p=200; a=0.6; b=[a;zeros(p,1);1]; a=[1;zeros(p,1);-a]; y5=filter(b,a,x); y6=fft(y5); %time wave through full-pass filter figure subplot(121);plot(y5) title(全通结构的混响器的时域波形’)%time wave through full-pass filter subplot(122);plot(abs(y6)) title(波形)%frequence wave through full-pass filter sound(y5,fs,bits); %clear all [x,fs,bit