m06_lowsample.m

matlab的经典讲义

t0=10;          %定义时间长度
ts=0.005;        %采样周期
fs=1/ts;
df=0.5;
t=[-t0/2:ts:t0/2];%采样时间序列
x=sin(200*t)./(200*t);%计算除零时产生的误差点
w=t0/(2*ts)+1;    %修正零点误差
x(w)=1;
x=x.*x;x=50.*x;
[M,mn,df]=fft_seq(x,ts,df);%傅立叶变换
M=M/fs;
M1=[M,M,M,M,M,M,M];  %延拓到多个周期
f=[-4*df*length(mn):df:3*df*length(mn)-df]-fs/2;
subplot(2,1,1);
stem(t,x);title('采样信号');
axis([-0.15,0.15,-1,50]);
subplot(2,1,2);
plot(f,abs(fftshift(M1)));title('采样信号频谱')