ms8r22.m

《数字信号处理实验(MATLAB版)》这

dt=0.1;
t=-4:dt:4;
N=length(t);            %求时间轴上采样点数
f=sinc(t);              %建立原连续信号
Tm=1;fm=1/Tm;           %由原连续信号周期,计算最高频率
wm=2*pi*fm;
k=1:N;
w1=k*wm/N;              %在频率轴上生成N个采样频率点
F1=f*exp(-j*t'*w1)*dt;  %傅里叶变换
subplot(4,1,1),plot(w1/(2*pi),abs(F1),'k');grid
axis([0 max(3*fm) 1.1*min(F1) 1.1*max(F1)]);
%生成fs>2fm.fs=2fm.fs<2fm三种抽样信号的振幅频谱
for i=1:3;
  if i<=2 c=0,else c=0.2,end
  fs=(4-i+c)*fm;Ts=1/fs;  %确定采样频率和周期      
  n=-4:Ts:4;      
  f=sinc(n);            %生成抽样信号
  N=length(n);          %求时间轴上采样点数
  wm=2*pi*fs;
  k=1:N;
  w=k*wm/N;
  F=f*exp(-j*n'*w)*Ts;
  subplot(4,1,i+1),plot(w/(2*pi),abs(F),'k');grid
  axis([0 max(3*fm) 1.1*min(F1) 1.1*max(F1)]);
end