实验三 信号采样与重建(参考答案)

1、实验3:信号采样与重构，讲师指导2013.5.5，1，1。实验目的，1。进一步研究MATLAB函数及其表示。2.掌握并验证信号的SHANNON采样定理。3.从采样序列中重构并恢复原始信号。2，2，实验内容，1。对于连续时间信号y(t)sin(24t) sin(40t)，它具有12Hz和20Hz的两个等幅分量。奈奎斯特频率2fmax用MATLAB绘图法计算。% AnalogSignalDt=0.005t=-1: dt :1；%采样频率200赫兹=正弦(24 * t)正弦(40 * t)；子情节(3，1，1)；曲线图(t*4，y)；xlabel(tin1/4sec。)；伊拉贝尔(y(t)；标题(模

2、拟信号)；子情节(3，1，2)；图(t*12，sin(24 * pi * t)；xlabel(tin1/12sec。)；伊拉贝尔(sin(24 * pi * t)；% f=子图(3，1，3)；图(t*20，sin(40 * pi * t)；XL Abel(tin1/20秒。)；伊拉贝尔(sin(40 * pi * t)；%f=20Hz，复合时看不到频率，分为两个分量，一般分为许多正弦和余弦，最高频率为两倍。3、另一种解决方案：频谱法(适用于许多正弦分量)t=033300.013331% t 0.01秒，采样频率100赫兹=正弦(24 * t)正弦(40 * t)；Y=fft(y)。f=(0:

3、length(y)-1)/长度(y)* 100；各点对应的频率图(1)；曲线图(f，绝对值(Y)；网格；%振幅图(2)；绘图(f，展开(角度(Y)* 180/pi)；网格；相位曲线不是折叠的，而是完整绘制的。两个频率分量，采样频率200，4、2:让连续信号x(t)=exp(1000|t|) a，找到连续傅里叶变换X(jw)(具有非常密集的点)，所以：满足采样定理，所以x(t)可以被xan替换，范围从0.005到0.005 0相位，并且在DC的最大相对幅度是0.002，5、2:当连续信号X(t)=exp(1000 | t |)a时，计算连续总共201点Wmax=2 * pi * 2000只要频率

4、为2000赫兹，M=500m=0:1:MW=m*Wmax/M。总共获得了501个钨(1)钨(501)的光谱点。xa=x * exp(j * t * W)* 0.00005；计算点Xa(1)Xa(501)。Xa=光谱501的实数(Xa)；Xa=fliplr(Xa)，Xa(2:501)；%左折叠，但Xa(1)不移动。频率点也保持折叠。图(1)；图(t*1000，x，r)；网格；% ms图(2)；图(W/(2*pi*1000)，Xa*1000，r)；振幅放大为1000，6、频谱和波形，振幅放大为1000倍后，7、2:当连续信号x(t)=exp(1000 | t |)b、Fs=5000和1000Hz时

5、，首先采样，并计算离散DFT频谱。t=-0.00533600.00005:0.005；%原始信号条件不变x1=exp (1000 *绝对值(t)；Ts=0.0002。n=25:1:25。%Fs=5000，x=exp(1000 * ABS(n * Ts)；M=500。m=0:1:Mw=pi * M/M；频率点Xa=x * exp(j * n * W)；计算离散频谱Xa=实数(Xa)；Xa=fliplr(Xa)，Xa(2: m1)；W=fliplr(W)，W(2: m1)；图(1)；图(t*1000，x1，r)；坚持；阀杆(n*Ts*1000，x)；拒绝；图(2)；绘图(W/pi，Xa)；8、光谱

6、。序列，9，C，降低Fs=1000Hz，首先采样，然后计算离散DFT频谱。t=-0.00533600.00005:0.005；%连续条件不变x1=exp (1000 *绝对值(t)；Ts=0.001。n=5:1:5。%Fs=1000，TS1mSx=exp(1000 * ABS(n * Ts)；n * Ts=0.0050.005M=500m=0:1:Mw=pi * M/M；频率点Xa=x * exp(j * n * W)；计算离散频谱Xa=实数(Xa)；Xa=fliplr(Xa)，Xa(2: m1)；W=fliplr(W)，W(2: m1)；图(1)；图(t*1000，x1，r)；坚持；阀杆(n

7、*Ts*1000，x)；拒绝；图(2)；绘图(W/pi，Xa)；10、降低采样速率和频谱，并注意两侧的衰减不接近0。时域采样序列。11、D信号的重建和恢复。Ts=0.0002n=25:1:25。Fs=5000。t=-0.005:0.00005:0.005；nTsn * Tsx=exp(1000 * ABS(NTs)；序列xaspline(nTs，x，t)；三次样条函数的重构。error 1 max(ABS(xaexp(1000 * ABS(t)；图(1)；图(t*1000，xa，r)；text(2，0.5，strcat(error=，num 2 str(error 1)；坚持；已绘制的曲线将被

8、保留。阀杆(n*Ts*1000，x)；拒绝；12、最大误差，以及13、D信号的重建，是不断减少和逐点重建的，不能恢复。Ts=0.001n=5:1:5。Fs=1000。t=-0.005:0.00005:0.005；nTsn * Tsx=exp(1000 * ABS(NTs)；序列xaspline(nTs，x，t)；三次样条函数法重构。error 2 max(ABS(xaexp(1000 * ABS(t)；图(1)；图(t*1000，xa，r)；text(2，0.5，strcat(error=，num 2 str(error 2)；坚持；阀杆(n*Ts*1000，x)；拒绝；14，最大误差，大失真，15，3，思考问题：连续时间信号的傅里叶变换matlab解，这里使用的近似公