数字信号处理实验作业(2)

subplot(4,1,1);plot(t,x);axis([min(t),max(t),1.1*min(x),1.1*max(x)]); title('ó?ê±óò?í?y???¨3é?ùD?o?'); for i=1:3;

fs=i*fm;Ts=1/fs; n=0:(3*T0)/Ts; t1=-7*T0:Ts:7*T0; x1=sinc(f0*n/fs);

T_N=ones(length(n),1)*t1-n'*Ts*ones(1,length(t1)); xa=x1*sinc(fs*pi*T_N); subplot(4,1,i+1);plot(t1,xa);

axis([min(t1),max(t1),1.1*min(xa),1.1*max(xa)]); end

用时域卷积重建抽样信号10.50-610.50-610.50-610.50-6-4-20246-4-20246-4-20246-4-20246 3、已知一个时间序列的频谱为：X(e)=jωn=-??x(n)e?-jωn=2+4e-jω+6e-j2ω+4e-j3ω+2e-j4ω

分别取频域抽样点数N为3、5和10，用IFFT计算并求出其时间序列x(n)，绘图显示个时间序列。由此讨论由频域抽样不失真地恢复原时域信号的条件。

Ts=1;N0=[3,5,10]; for r=1:3; N=N0(r); D=2*pi/(Ts*N);

kn=floor(-(N-1)/2:-1/2); kp=floor(0:(N-1)/2); w=[kp,kn]*D;

X=2+4*exp(-j*w)+6*exp(-j*2*w)+4*exp(-j*3*w)+2*exp(-j*4*w); n=0:N-1; x=ifft(X,N)

subplot(1,3,r);stem(n*Ts,abs(x)); box end

76665554433222341110012002400510 4、已知一个频率范围在[-6.28,6.28]rad/s的频谱，在模拟频率|Ωc|=3.14处幅度为1，其它范围幅度为0.计算其连续信号xa(t)，并绘图显示信号曲线。

wc=3.14;Tmax=0.1*pi/wc; Ts=input('(TsNmin)N='); D=2*pi/(Ts*N); M=floor(wc/D);

Xa=[zeros(1,M),1,zeros(1,N-2*M-1),1,zeros(1,M-1)]; n=-(N-1)/2:(N-1)/2;

xa=abs(fftshift(ifft(Xa/Ts))); plot(n*Ts,xa);

输入Ts=0.1 N=100

0.250.20.150.10.050-5-4-3-2-1012345 四、实验预习： 1、认真阅读实验原理部分，明确实验目的，读懂例题程序，了解实验方法。 2、根据实验任务预先编写实验程序。 3、预习思考题：

① 什么是内插公式？在MATLAB中内插公式可用什么函数来编写？ ② 从频域抽样序列不失真地恢复离散时域信号的条件是什么？ 五、实验报告：

1、列写调试通过的实验程序，打印实验程序产生的曲线图形。 2、给出预习思考题答案。 3、思考题：

① 试归纳用IFFT数值计算方法从频谱恢复离散时间序列的方法和步骤。 ② 从频谱恢复连续时间信号与恢复离散时间序列有何不同？