数字信号处理实验六报告.docVIP

实验六 频域抽样定理和音频信号的处理实验报告 （一）频域抽样定理 给定信号 1．利用DTFT计算信号的频谱，一个周期内角频率离散为M=1024点，画出频谱图，标明坐标轴。 n=0:100; %设定n及其取值范围 for n1=0:13 %对于n处于不同的取值范围将n代入不同的表达式 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end for n3=27:100 xn(n3+1)=0; end M=1024; %设定抽样离散点的个数 k=0:M-1; %设定k的取值范围 w=2*pi*k/M; %定义数字角频率 [X,w] = dtft2( xn,n, M ) %调用dtft2子程序求频谱 plot(w,abs(X)); %画出幅度值的连续图像 xlabel(w/rad);ylabel(|X(exp(jw))|);title( M=1024时的信号频谱图像); %标明图像的横纵坐标和图像标题 function [X,w] = dtft2(xn, n, M ) %定义x(n)的DTFT函数 w=0:2*pi/M:2*pi-2*pi/M; %将数字角频率w离散化 L=length(n); %设定L为序列n的长度 for (k=1:M) %外层循环，w循环M次 sum=0; %每确定一个w值，将sum赋初值为零 for (m=1:L) %内层循环,对n求和，循环次数为n的长度 sum=sum+xn(m)*exp(-j*w(k)*n(m)); %求和 X(k)=sum; %把每一次各x(n)的和的总值赋给X，然后开始对下一个w的求和过程 end%内层循环结束 end%外层循环结束 M=1024时的信号频谱图像如图1-1所示： 图1-1 M=1024时的信号频谱图像 2．分别对信号的频谱在区间上等间隔抽样16点和32点，得到和。离散傅里叶反变换后得到时域信号和。 M=16时： n=0:100;%设定n及其取值范围 for n1=0:13 %对于n处于不同的取值范围将n代入不同的表达式 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end for n3=27:100 xn(n3+1)=0; end M=16;%设定抽样点的个数 k=0:M-1; %设定k的取值范围 w=2*pi*k/M; %定义数字角频率 [X,w] = dtft2( xn,n, M ) %调用dtft2子程序求频谱 stem(w,abs(X),.); %画出幅度值的抽样离散图像 xlabel(k);ylabel(X16(k));title(X（exp(jw))的16点抽样);%标明横纵坐标和图像标题 M=32时只需将程序中的M值改变即可，子函数程序同实验一。 16点抽样的图形如图1-2所示： 图1-2 的16点抽样点抽样的图形如图1-2所示： 图1-2-2 的32点抽样 3．画出信号和的图形，计算与和的均方误差。从时域角度上进行对比和分析，验证频域抽样定理。 M=16时：function xn=idft(Xk,N) N=length(Xk); for n=0:N-1 sum=0; for k=0:N-1 sum=sum+(1/N)*(Xk(k+1)*exp(j*2*pi*n*k/N)); end xn(n+1)=sum; end n=0:26; %设定n及其取值范围 for n1=0:13 %对于n处于不同的取值范围将n代入不同的表达式 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=16; %设定抽样点的个数 k=0:M-1; %设定k的取值范围 w=2*pi*k/M; %定义数字角频率 [X,w] = dtft2( xn,n, M ) %调用dtft2子程序求频谱 xn=idft(X,M) %调用idft求时域信号xn m=0:15; %限定m的取值范围 stem(m,xn);%画出xn16的离散图像 xlabel(m);ylabel(xn16);title( M=16时的时域信号图像); %标明图像的横纵坐标和图像标题 图1-3-1 M=16时的时域信号图像 M=32时： n=0:26; %设定n及其取值范围 for n1=0:13 %对于n处于不同的取值范围