信号与系统实验五六报告.docVIP

实验五 采样 §5.1由欠采样引起的混叠 基本题 考虑正弦信号 若用频率采样，那么离散时间信号就等于 假定采样频率固定在。 1．假设并定义T=1/8192。创建向量n=[0:8192]，使得t=n*T包含了区间内8192个时间样本。创建向量x，它包含在t的时间样本上的样本。 实验思路：定义向量n，取8192个点，对t=n*T得到n个样本值 实验过程：clc syms t; w0=2000*pi; T=1/8192; x1=sin(w0*t); n=[0:8192]; xn=sin(w0*n*T) 2．用stem对展示前50个样本，用plot对采样时间展示的前50个样本。 实验思路：定义向量n，得到xn的离散表达式，用stem画图 实验过程： clc syms t; w0=2000*pi; T=1/8192; x1=sin(w0*t); n=[0:50] xn=sin(w0*n*T) subplot(2,1,2);stem(n,xn); subplot(2,1,1);plot(xn); 实验结果： 实验分析：结果如图所示 为了计算带限重建信号的连续时间傅立叶变换的样本，今用下列函数 function [X,w]=ctffts(x,T) N=length(x); X=fftshift(fft(x,N))*(2*pi/N); w=linspace(-1,1-1/N,N)/(2*T); 这个函数用fft计算重建信号的傅立叶变换。文件ctffts.m应该装在相应的MATLABPATH中。 3．用[X,w]=ctffts(x,T)计算重建信号的连续时间傅立叶变换。画出X对w的幅值图。X在合理的频率值上是非零吗？假定当X接近于零时，相位等于零，X的相位正确吗？ 实验思路：编辑ctffts函数，再用函数计算重建信号的连续时间傅立叶变换再画出幅度和相位图 实验过程：clc syms t; w0=2000*pi; T=1/8192; x1=sin(w0*t); n=[0:8192]; xn=sin(w0*n*T); [X,w]=ctffts(xn,T); xr=abs(X); plot(xr);title(幅度谱);figure xm=angle(X); plot(xm);title(相位谱) 实验结果： 实验分析：X在合理的频率值上是非零的，假定当X接近于零时，相位等于零，X的相位不正确，在那一点是突变点。 中等题 4．对正弦频率和重作1～3。X的幅值对于所预计的频率还是非零吗？X的相位正确吗？ 实验思路：按照1-3的思路来做 实验过程： clc syms t; w0=3000*pi; T=1/8192; x1=sin(w0*t); n=[0:50] xn=sin(w0*n*T) subplot(2,1,2);stem(n,xn); subplot(2,1,1);plot(xn); clc syms t; w0=3000*pi; T=1/8192; x1=sin(w0*t); n=[0:8192]; xn=sin(w0*n*T); [X,w]=ctffts(xn,T); plot(xr);title(幅度谱);figure xm=angle(X); plot(xm);title(相位谱); 实验结果： 实验分析：对和，可知X的幅值对于所预计的频率是非零的，X的相位不正确，。5．用sound(x,1/T)将4中创建的每个采样信号放出来。你听到的音调随频率的提高而提高吗？注意，和plot一样，sound函数也有内插的作用。 6．现在对正弦频率，，，和重作1和3。也用sound将每个采样信号放出来。你所听到的音调高度随每次频率的增加而提高吗？如果不是，你能解释这个现象吗？ 答：通过对多个频率信号释放，感觉音调高度并不是随频率的增加而增加，但趋势却是又是在增高的 §5.2 由样本重建信号 在下面的练习中，既用带限内插，又用线性内插从采样时刻得到的样本来重建下列信号： 基本题 1．用解析法证明，在采样时刻和都等于样本值。这样的内插器称为严格内插器，因为它们保留了原始信号在采样时刻的真正值。这个带限内插和线性内插滤波器时因果的吗？ 不是因果的，因为t0是，、不等于零。 2．和带限吗？若是，带宽是多少？ 由图可知，x1是带限的，大小为8*pi/5*2, x2t不是带限的。 3．创建一个向量ts，它包含在内的采样间隔。将对应于ts的和的样本存入向量xs1和xs2，用stem画出xs1和xs2对ts的图。 实验思路：创建向量ts，对T=0.5，t=nT,对信号进行