[转载]信号的采样与恢复 matlab模拟

2016-09-27 13:38阅读：

http://blog.sina.cn/dpool/blog/u/1238979861

原文作者：fpcsoc

以上为对直线的采样与恢复不同采样周期效果不同。
采样周期越小，信号恢复符合度越高。

对函数 sin(w0*m*T)+sin(w0*(m*T+300))的采样与恢复;
function [s]=Super_Command()
%%调用函数处理恢复采样 Super_Command
clf;
clear;
w0=0.03; %信号频率
T=5; %采样周期
M=500/T;
m=-M:M; %求和项数范围越多越精确 M=500;
% 待采样函数（可修改）采样点数与求和项数一致
Xat=sin(w0*(m*T+300)); % w0越小越好 w0=0.01;
%w0=2*pi/T0; T0=200*pi
subplot 211;
stem(m*T,Xat); %采样信号
title('采样信号信号频率w0=0.03 采样周期 T=5 ');
a=1;
for t=-M*T:.1:M*T %计算各恢复点值
s(a)=Super_Sign_Rebuild(w0,M,T,t);
a=a+1;
end;
subplot 212;
t=-M*T:.1:M*T ; %恢复信号作图
plot(m*T,Xat,'red',t,s,'k'); %采样信号与恢复信号对比
title('采样信号与恢复信号对比');

function [S]=Super_Sign_Rebuild(w0,N,T,t)
%%信号重建 xa(t)=sum( x(nT)sin[(pi/T)(t-nT)]/(pi/T)(t-nT) ）
%%采样周期T
n=-N:N; %求和项数越多越精确
Xat=sin(w0*(n*T+300)); %待恢复函数 w0越小越好
Samp= sin((pi/T)*(t-n*T))./((pi/T)*(t-n*T));
S=sum(Xat.*Samp);