ms15_2.m

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

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