数字信号处理实验6.doc

实用标准文案 精彩文档 数字信号处理实验六 张达成 201000121241 通信二班 给定信号 。 信号图像编程： n=0:100; for n1=0:13 x(n1+1)=n1+1; end for n2=14:26 x(n2+1)=27-n2; end for n3=27:100 x(n3+1)=0; end stem(n,x); xlabel(n); ylabel(x(n)); title(给定信号序列); 编程图像： 图0 给定信号图像 1．利用DTFT计算信号的频谱，一个周期内角频率离散为M=1000点，画出频谱图，标明坐标轴。 n=0:26; for n1=0:13 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=1000; %设定一个周期内1000个点 k=0:M-1; w=2*pi*k/M; X=dtft(xn,M); plot(w,abs(X)); xlabel(w/rad); ylabel(|X(ejw)|); title(信号频谱图 M=1000); 编程图像： 图1 给定信号的DTFT 2．分别对信号的频谱在区间上等间隔抽样16点和32点，得到和。离散傅里叶反变换后得到时域信号和。 编程思路：只需要将代码中M的值分别更改为16和32即可完成对于DTFT的16点和32点抽样。 编程代码： n=0:100; for n1=0:13 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=16; %M值即为抽样点数 k=0:M-1; w=2*pi*k/M; X=dtft(xn,M); stem(w,abs(X)); xlabel(k); ylabel(X16(k)); title(Xejw16点抽样); 32点抽样代码不再重复赘述 编程图像： 图2-1 16点DTFT抽样图像 图2-2 32点DTFT抽样图像 3．画出信号和的图形，计算与和的均方误差。从时域角度上进行对比和分析，验证频域抽样定理。 编程思路：抽样后套用IDFT函数对抽样进行还原，还原后的序列分别是从0到15以及0到31 编程代码： for n1=0:13 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=16; %频率抽样16个点，故还原为时域信号也应包含16个点 X=dtft(xn,M); xn16=idft(X,M); nn=0:15; stem(nn,xn16); xlabel(n); ylabel(x16(n)); title(16点频率抽样信号还原); for n1=0:13 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=32; X=dtft(xn,M); xn32=idft(X,M); nn=0:31; %频率抽样32个点，故还原为时域信号也应包含32个点 stem(nn,xn32); xlabel(n); ylabel(x32(n)); title(32点频率抽样信号还原); 编程图像: 图3-1 16点DTFT抽样还原时域信号 图3-1 32点DTFT抽样还原时域信号 计算与和的均方误差： 由于16点抽样还原信号产生混叠，故和原信号计算方差时应根据周期性变化补齐剩下的点；而计算32点抽样还原信号时，可以忽略27至31这5个零点。 16点计算方差编程代码： for n1=0:13 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=16; X=dtft(xn,M); xn16=idft(X,M); for c=17:27 xn16(c)=12; %根据周期性将剩余的点补齐，值均为12 end sum=0; for n=1:27 sum=sum+(abs(xn(n))-abs(xn16(n)))^2; end d=sum/27 matlab计算得到16点抽样还原信号与原信号的方差为37.4815。 32点计算方差编程代码： for n1=0:13 xn(n1+1)=n1+1; end for n2=14:26 xn(n2+1)=27-n2; end M=32; X=dtft(xn,M); xn32=idft(X,M); sum=0; for n=1:27 sum=sum+(abs(xn(n))-abs(xn32(n)))^2; %计算方差时舍去27至31这5个零点 end d=sum/27 matlab计算得到3