信号与系统课程设计2.docxVIP

《信号与系统》课程设计报告姓名学号年级专业及班级指导教师设计起止日期成绩指导教师签名1.课程设计目的和要求：32.课程设计原理：33.课程设计所用设备：44.课程设计内容和步骤：45.设计过程、实现过程、实现结果：56.设计结果分析；117.设计总结：12参考文献121.课程设计目的和要求：1)熟悉抽样定理、信号的抽样过程；2)通过实验观察欠采样时信号频谱的混叠现象；3)掌握抽样前后信号的频谱的变化，加深对抽样定理的理解；4)掌握抽样频率的确定方法。2.课程设计原理：离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。抽样信号fs(t)可以看成连续信号f (t)和一组开关函数s(t)的乘积。s(t)是一组周期性窄脉冲，Ts称为抽样周期，其倒数fs = 1/Ts 称抽样频率。对抽样信号进行傅里叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的原信号频率。平移的频率等于抽样频率fs 及其谐波频率2 f s 、3 f s ······。当抽样信号是周期性窄脉冲时，平移后的频率幅度按sinxx规律衰减。抽样信号的频谱是原信号频谱周期的延拓，它占有的频带要比原信号频谱宽得多。但原信号得以恢复的条件是fs = 2B，其中fs为抽样频率，B为原信号占有的频带宽度。而Fmin = 2B为最低抽样频率又称“奈奎斯特抽样率”。当fs 2B时，抽样信号的频谱会发生混迭，从发生混迭后的频谱中我们无法用低通滤波器获得原信号频谱的全部内容。在实际使用中，仅包含有限频率的信号是极少的。因此即使fs = 2B，恢复后的信号失真还是难免的。3.课程设计所用设备：MATLAB软件。4.课程设计内容和步骤：1.许多实际工程信号不满足带限条件，因此要进行抗混叠滤波。比较抗混叠滤波前后歌曲信号频谱并试听。2．混叠误差与截断误差比较。3.不同抽样频率的语音信号效果比较。3.根据时域抽样定理，对连续时间信号进行抽样时，只需抽样率 fs = 2fm。在工程应用中，抽样率常设为 fs =(3~5) fm，为什么？5.设计过程、实现过程、实现结果：1.[x,fs,bits]=wavread(chuanqi.wav);sound(y,fs,bits); %sound(x,fs,bits); %对原歌曲的听取N=5;wc=0.05;[b,a]=butter(N,wc,low);X=fft(x);X=fftshift(X);k=fs/length(x);subplot(2,1,1);plot([-fs/2:k: fs/2-k],abs(X));title(滤波前频谱); %滤波前的频谱y=filter(b,a,x);sound(y,fs,bits);%sound(y,fs,bits); %对进行抗混叠滤波后的听取Y=fft(y);Y=fftshift(Y);k=fs/length(y);subplot(2,1,2);plot([-fs/2:k: fs/2-k],abs(Y));title(滤波后频谱); %滤波后的频谱我选的是音质很好的音频文件，在抗混叠滤波处理之后音质有所下降，因为低通滤波器滤去了一部分高频信号，造成了信号的缺失；2．混叠误差：用5512HZ频率抽样[y,fs,nbits]=wavread(chuanqi.wav);y=y(1:1000000);y1=y(1:8:end);yy1=fft(y);yy1=fftshift(yy1);k=fs/length(y);plot([-fs/2:k: fs/2-k],abs(yy1));title(抽样前信号频谱);figure;yy2=fft(y1);yy2=fftshift(yy2);k=5512/length(y1);plot([-5512/2:k: 5512/2-k],abs(yy2));title(抽样前信号频谱);wavwrite(y1,5512,chuanqi_5512.wav);由图可以看出，以频率抽样后信号发生了混叠误差，信号有失真，因此音频听起来非常多的杂音。截断误差先进行抗混叠滤波然再进行抽样。[x,fs,bits]=wavread(chuanqi.wav);%sound(x,fs,bits);N=5;wc=0.05;[b,a]=butter(N,wc,low);y=filter(b,a,x);y=y(1:1000000);y1=y(1:8:end);yy1=fft(y);yy1=fftshift(yy1);k=fs/length(y);plot([-fs/2:k: fs/2-k],abs(yy1));title(滤波后抽样前频谱);figure;yy2=fft(y1);yy2=fftshift(yy2);k=5512/length(y1);plot([-5512