-物电-刘石山-数字信号处理实验报告.docVIP

数字信号处理实验 姓 名： 刘石山 专 业: 物理电子学 学 号： 2011211011 任课教师： 赵建军 日 期： 2012年1月3日 目录 实验一：用matlab对音频文件进行处理······················3 实验二：画sin函数···········································5 实验三：画方波················································6 实验四 验证定理2.1··········································9 实验五.连续时间信号的傅里叶级数（CTFS）···············10 实验六：相位的重要性·······································12 实验七 用MATLAB画欧拉公式····························14 实验八对 shift（）的理解··································15 实验九：离散信号的傅里叶变换·····························16 实验十 如何选取N值·········································18 Program 2.2··················································19 Example 2.8 DTFS系数的幅度与相位·······················19 Example 2.9··················································21 Example 2.12 CTFS from DTFS····························23 Example2.13 CTFS周期窗函数的FFT··················24 实验一：用matlab对音频文件进行处理 声卡对声音的处理质量可以用三个基本参数来衡量即声道数S、采样频率Fs和采样位数bit。用matlab获取声音文件时可以用wavread函数。音频信号图如图（1.1）. 程序：(1) 获取声音信号 [s,Fs,bit]=wavread(liuvoice.wav);%读取声音信号 a=s(:,1); %作单声道处理 a1=a(100:100100); %截取100000个信号采样点 b=length(a1); t1=(0:b-1)/Fs; figure;plot(t1,a1); xlabel(Time(second)); ylabel(Amplitude); Title(Time signal) （1.1） （2）对上述声音信号进行傅里叶变换，并用matlab画出其频谱图，如图（1.2）. 程序： x=wavread(liushishan.wav); T=1/8000; N=length(x); x=fft(x); D=2*pi/(T*N); m=ceil(-N/2):ceil(N/2-1); plot(m*D/(2*pi),abs(fftshift(x))); xlabel(Frequency(Hz)); ylabel(Magnitude spectrum); title(Frequency content); （1.2） 下图是上图只取正频率的情况下，所呈现的图。如图（1.3）. 程序： x=wavread(liuvoice.wav); T=1/8000; x=fft(x); N=length(x); D=2*pi/(T*N); m=0:ceil(N/2); plot(m*D,abs(x(m+1))); xlabel(Frequency(rad/s)); ylabel(Magnitude spectrum); title(Frequency content in the positive frequency range)； （1.3） 实验二：画