zhuqi.m

来自「fft做谱分析各种连续信号 用MATLAB实现」· M 代码 · 共 10 行

% 连续周期时间信号的频谱 x（t）＝cos(5t)   
clc;
clear;
N=input('N=');
T=0.1;n=1:N;%原始数据
D=2*pi/(N*T);             %频率分辨率 
xa=sin(5*n*T);            %采样 
Xa=T*fftshift(fft (xa));Xa(1)%求h(n)的FFT，移到对称位置
k=floor(-(N-1)/2);        %奈奎斯特频率下标向量
plot(k*D,abs(Xa))