shuzixinhao课程设计实验报告.docVIP

数字信号处理课程设计 专业：测控技术与仪器 班级：09级二班 姓名：李兵阵 学号：20092240213 指导老师：岳学东 1.【源程序如下】： [a,fs,bit]=wavread(I do片段); y1=a( : ,1); f1=fft(y1); n=length(f1); w=2/n*[0:n-1]; figure subplot(2,1,1);plot(a); grid on; title(音频); xlabel(时间); ylabel(幅度); subplot(2,1,2);plot(w,abs(f1)); grid on; title(频谱); xlabel(频率/pi); ylabel(幅度); 【程序运行结果】： 2【源程序如下】： [y,fs,bit]=wavread(I do片段) n1=length(y); D=2;%设置抽样间隔 y2=y(1:D:n1);%减抽样 n2=length(y2);%减抽样后信号长度 t2=(0:n2-1)/fs;%设置横坐标 figure subplot(2,1,1); plot(t2,y2); %绘制减抽样信号波形图 ylabel(幅度); xlabel(时间（s）); title(2:1减抽样信号波形); Y2=fft(y2); %对y2进行n2点fft谱分析 w2=2/n2*[0:n2-1]; subplot(2,1,2); plot(w2,abs(Y2));%绘制减抽样信号频谱图 title(2:1减抽样信号频谱); xlabel(数字角频率w); ylabel(幅度); grid on; sound(y2,fs/D); [y,fs,bit]=wavread(I do片段) n1=length(y); D=4;%设置抽样间隔 y2=y(1:D:n1);%减抽样 n2=length(y2);%减抽样后信号长度 t2=(0:n2-1)/fs;%设置横坐标 figure subplot(2,1,1); plot(t2,y2); %绘制减抽样信号波形图 ylabel(幅度); xlabel(时间（s）); title(4:1减抽样信号波形); Y2=fft(y2); %对y2进行n2点fft谱分析 w2=2/n2*[0:n2-1]; subplot(2,1,2); plot(w2,abs(Y2));%绘制减抽样信号频谱图 title(4:1减抽样信号频谱); xlabel(数字角频率w); ylabel(幅度); grid on; sound(y2,fs/D); 【运行结果】： 分析：采用较小的抽样间隔进行抽样时，频谱并未发生混叠，而采用较大的抽样频间隔进行抽样时，频谱发生严重混叠。这是因为，　奈奎斯特抽样定理：要从抽样信号中无失真地恢复原信号，抽样频率应大于2倍信号最高频率。 　　抽样频率小于2倍频谱最高频率时，信号的频谱有混叠。 　　抽样频率大于2倍频谱最高频率时，信号的频谱无混叠 3【源程序如下】： close all;clear all;clc [a,fs,bit]=wavread(I do片段); y1=a( : ,1); f1=fft(y1); n=length(f1); w=2/n*[0:n-1]; figure subplot(2,1,1);plot(a); grid on; title(音频); xlabel(时间); ylabel(幅度); subplot(2,1,2);plot(w,abs(f1)); grid on; title(频谱); xlabel(频率/pi); ylabel(幅度); t=0:n-1; y2=cos(pi*1/2*t); y3=y1.*y2; figure subplot(2,1,1);plot(y3); grid on; title(AM调制音频信号 低); xlabel(时间); ylabel(幅度); f3=fft(y3); n2=length(f3); w2=2/n2*[0:n2-1]; subplot(2,1,2);plot(w2,abs(f3)); grid on; title(AM调制频谱 低); xlabel(频率/pi); ylabel(幅度); y4=cos(pi*4/5*t); y5=y1.*y4; figure subplot(2,1,1);plot(y5); grid on; title(AM调制音频信号 高); xlabel(时间); ylabel(幅度); f5=fft(y5); n3=length(f5); w3=2/n2*[0:n