T0=[0.002 0.0015 0.001 0.005]; %四种采样间隔
LT=0.05; %信号记录长度（秒）
for i=1:4
T=T0(i);
N=LT/T+1; %信号记录点数
n=0:N-1;
F=1/LT; %频率分辨率
f=(floor(-(N-1)/2):floor((N-1)/2))*F;
x=cos(1400*pi*n*T); %对模拟信号进行抽样
X=T*fftshift(fft(x)); %求傅立叶变换
subplot(2,2,i),plot(f,abs(X));xlabel('模拟角频率（Hz）');
axis([min(f) max(f) 0 inf]); %坐标限制
str=['T=' num2str(T) '; fs=' num2str(1/T)];title(str); %标题显示采样间隔和采
样频率
end
%
T0=[0.09 0.08 0.05 0.03]; %四种采样间隔
LT=10; %信号记录长度（秒）
for i=1:4
T=T0(i);
N=LT/T+1; %信号记录点数
n=0:N-1;
F=1/LT; %频率分辨率
f=(floor(-(N-1)/2):floor((N-1)/2))*F;
x=exp(-14*n*T); %对模拟信号进行抽样
X=T*fftshift(fft(x)); %求傅立叶变换
subplot(2,2,i),plot(f,abs(X));xlabel('模拟角频率（Hz）');
axis([min(f) max(f) 0 inf]); %坐标限制
str=['T=' num2str(T) '; fs=' num2str(1/T)];title(str); %标题显示采样间隔和采
样频率
end
%
T0=[0.08 0.05 0.0015 0.001]; %四种采样间隔
