continuoussignaldft.m

来自「连续信号离散傅立叶变换」· M 代码 · 共 23 行

% ContinuousSignalDFT


clear;
disp('用DFT对连续信号进行谱分析');
Repeat = 1;
while Repeat == 1
    disp('输入采样间隔Ts，采样点数N，信号x(t)');
    Ts = input('Ts=');
    fs = 1/Ts;
    N = input('N=');
    n = 0:N-1;
    t = n*Ts;
    x = input('t=n*Ts,x(t)=');
    Xk = dft(x,N);
    k = 0:N-1;
    
    figure;
    subplot(1,2,1),plot(t,x),title(sprintf('x(t),Ts=%f',Ts));
    subplot(1,2,2),stem(k,abs(Xk)),title(sprintf('X(k) N=%d',N));
    
    Repeat = input('Still to continue, Yes=1 No=0? Ans=');
end