文档介绍:该【matlab信号抽样与恢复 】是由【淘气小宇】上传分享,文档一共【3】页,该文档可以免费在线阅读,需要了解更多关于【matlab信号抽样与恢复 】的内容,可以使用淘豆网的站内搜索功能,选择自己适合的文档,以下文字是截取该文章内的部分文字,如需要获得完整电子版,请下载此文档到您的设备,方便您编辑和打印。实验一信号抽样与恢复
一、实验目的
学会用MATLAB实现连续信号的采样和重建
二、实验原理
若f(t)是带限信号,带宽为m,f(t)经采样后的频谱Fs( )就是将f(t)的频谱
F( )在频率轴上以采样频率s为间隔进行周期延拓。因此,当sm时,不会发生频率混
叠;而当s<m时将发生频率混叠。
经采样后获取信号fs(t)经理想低通h(t)则可获取重建信号f(t),即:
f(t)=fs(t)*h(t)
其中:fs(t)=f(t)(tnTs)=f(nTs)(tnTs)
因此:
f(t)=fs(t)*h(t)=f(nTs)(tnTs)*TscSa(ct)
c
f(nTs)Sa[c(t
nTs)]
=Ts
上式表示,连续信号可以张开成抽样函数的无量级数。
利用MATLAB中的sinc(t)
sin(t)来表示Sa(t),有Sa(t)
sinc(t),因此可以得
t
到在MATLAB
中信号由f(nTs)重建f(t)的表达式以下:
f(t)=Tscf(nTs)sinc[c(tnTs)]
我们采用信号f(t)=Sa(t)作为被采样信号,当采样频率s=2m时,称为临界采样。我
们取理想低通的截止频率c=m。下面程序实现对信号f(t)=Sa(t)的采样及由该采样信号恢
复重建Sa(t):
三、上机实验内容
;
(t)=*(1+cost)*(u(t+pi)-u(t-pi))
,由于不是严格的频带有限信号,
但其频谱大部分集中在
[0,2]之间,带宽
wm
可依照必然的精度要求做一些近似。试依照以下两种情况用
MATLAB
实
现由
f(t)的抽样信号
fs(t)重建
f(t)
并求两者误差,解析两种情况下的结果。
wm=2,wc=,Ts=1;
wm=2,wc=2,Ts=
低通采样定理以程序实现
,详尽参数参照框图内参数。
例1-1Sa(t)的临界采样及信号重构;
wm=1;%信号带宽
wc=wm;%滤波器截止频率
Ts=pi/wm;%采样间隔
ws=2*pi/Ts;%采样角频率
n=-100:100;%时域采样电数
nTs=n*Ts%时域采样点
f=sinc(nTs/pi);
Dt=;t=-15:Dt:15;
fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));
t1=-15::15;
%信号重构
f1=sinc(t1/pi);
subplot(211);
stem(t1,f1);
xlabel('kTs');
ylabel('f(kTs)');
title('sa(t)=sinc(t/pi)的临界采样信号');
subplot(212);
plot(t,fa)
xlabel('t');
ylabel('fa(t)');
title('由sa(t)=sinc(t/pi)的临界采样信号重构sa(t)');
grid;
例1-2Sa(t)的过采样及信号重构和绝对误差解析
程序和例4-1近似,将采样间隔改成Ts=*pi/wm,滤波器截止频率该成wc=*wm,
增加一个误差函数
wm=1;
wc=*wm;
Ts=*pi/wm;
ws=2*pi/Ts;
n=-100:100;
nTs=n*Ts
f=sinc(nTs/pi);
Dt=;t=-15:Dt:15;
fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));
error=abs(fa-sinc(t/pi));%重构信号与原信号误差
t1=-15::15;
f1=sinc(t1/pi);
subplot(311);
stem(t1,f1);
xlabel('kTs');
ylabel('f(kTs)');
title('sa(t)=sinc(t/pi)的采样信号');
subplot(312);
plot(t,fa)
xlabel('t');
ylabel('fa(t)');
title('由
sa(t)=sinc(t/pi)
的过采样信号重构
sa(t)');
grid;
subplot(313);
plot(t,error);
xlabel('t');
ylabel('error(t)');
title('过采样信号与原信号的误差error(t)');
例1-3Sa(t)的欠采样及信号重构和绝对误差解析
程序和例
4-2
近似,将采样间隔改成
Ts=*pi/wm,
滤波器截止频率该成
wc=wm=1