1、1因为指数的系数较大,所以所取的时间间隔应当足够小t=-0.01:0.00001:0.01;然后表示出信号x(t) x=exp(-1000*abs(t);调用plot函数画出模拟的信号图形。2x(t)的频谱首先要考虑角频率的取值范围,选取w= -104:104;然后根据公式写出频谱的表达式X=2000/(10002+W.2);然后调用plot函数画出w和X的图形。3. 以抽样频率3000Hz对x(t)抽样,所以抽样的时间间隔为T=1/3000;离散化以后的时域为t=-0.01:1/3000:然后表示出离散序列xn=exp(-1000*abs(t);调用stem(t,xn);画出离散序列图形。4
2、.先调用子函数dtft function DTFT(xn,N,M);然后根据上一个程序表示出xn。再求出离散后的角频率w:给定长度为N的序列x(n)的离散时间傅里叶变换(DTFT)的公式为 。可以看出x(n)的DTFT仍然是一个连续函数,所以需要将数字角频域w离散化,设一个频率周期内离散点有M个,则第k个点所代表的数字角频率这样x(n)的DTFT变成:所以表示出来Xw=xn*exp(-j*(t*w)。最后调用plot函数画出幅频响应和相频响应。5. 将(1)中求出的原信号x(t)的频谱的估计与2中求出的x(t)的真实频谱在相同的频率点处的值做均方误差。计算时,均方误差是指两者对应点之差的平方和
3、的平均值的平方根。原信号x(t)的频谱的估计为,因为关于对称,所以应当取w的有效部分信号的重建(理想内插公式):其中x(mT)是某一个抽样值,t是任意的一个时刻。它表示重建信号的任意时刻的函数值是所有抽样值与该时刻插值函数的乘积的总和。求出原信号的估计y(t)后,就可以求出均方误差,从而可以从连续时间域上考察信号的恢复程度。6将抽样频率3000HZ换成800HZ重做。.三、程序脚本,并注释:1. t=-0.01: %取t的取值范围x=exp(-1000*abs(t); %表示出信号xplot(t,abs(x); %画出信号波形xlabel(t);ylabel(xtitle(原信号xgrid;2
4、. w=-104: %规定模拟角频率w的取值范围X=2000./(10002+w.2); %信号的频谱 figure;plot(w,X);频谱3T=1/3000; %抽样时间间隔t=-0.01: %规定时间t取值范围xn=exp(-1000*abs(t); %离散后的序列xnstem(t,xn);xn离散信号xn 4. function DTFT(xn,N,M);M=500;k=-M:M;w=2*pi*1000/M*k;Xw=xn*exp(-j*(t*w);subplot(1,2,1);plot(w,abs(Xw); grid;wrad/s幅度abs幅频特性subplot(1,2,2);plo
5、t(w,phase(Xw);相位phase相频特性5. fs=3000; %抽样频率T=1/fs;T: %规定t的取值范围 %离散系列xn 采样信号y=-0.01:0.001:n=-0.01/T:0.01/T;M=ones(length(n),1)*y-n*T*ones(1,length(y);fs=xn*sinc(fs*M);%把不是0的地方用sinc函数来计算,隐含着为零的那一个是为1的.plot(y,fs,r N=length(t); %t的长度i=1; %设置一个累加变量for tn=ty(i)=x_n*(sin(pi/T*(tn-n*T)+eps)./(pi/T*(tn-n*T)+e
6、ps); %内插重建t的数值 i=i+1;endfigure,plot(t,y); %画出时域内插重建函数,内插恢复时间t内插函数值 y3000Hz时域内插重建函数 %原信号diff=sqrt(sum(x-y).2)/N); %均方误差6,将频率改为800HZ重做T=1/800;f=-0.01:1/800:xn=exp(-1000*abs(f);stem(xn);f HZ离散信号xn7. function DTFT(xn,N,M);w rad/s8. fs=800;重建信号内插函数值y800Hz时域内插函数值 四仿真结果、图形:1. 2. 3.抽样频率为3000HZ时 4.抽样频率为3000H
7、Z时5抽样频率3000hz 6,.抽样频率为800HZ时. 7.抽样频率为800HZ时8. 抽样频率800HZ五、结果分析和结论(1)3000Hz抽样时模拟频域上频谱的估计与真实频谱的均方误差比800Hz抽样时小,3000Hz连续时域上信号的估计与真实信号的均方误差也较小。(2)从实验的结果可以看出3000Hz抽样时抽样信号更能够反映原来真实的信号。800Hz由于抽样频率较小,所得结果有些失真。(3)从抽样后频谱图可以看出3000Hz时比800Hz时的频谱更接近真实频谱;从内插重建信号图形中也可以看出3000Hz时复原的信号也更接近真实信号。六、遇到的问题、解决方法及收获(1)时间t的取值区间必须足够小,并且时间间隔足够小才能得出信号波形,当间隔较大时看不到波形。(2)模拟信号的角频率的取值范围很难选取,因为信号指数系数较大,角频率选择范围不合适时,不能得出函数图形。(3)DTFT子函数的调用,之前已经练过很多次了,但是还不能熟练掌握,仍须多次调试。(4)对最后的两个问题内插函数以及函数恢复问题,难度很大,程序有参考同学的部分,仍有部分疑点大家还未解决,希望老师可以再详细讲讲。(5)本次实验的难度较大,花了较长的时间,但是从实验中验证了一些理论,加深了对理论的理解。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1