untitled2.m

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

n=0:50;
A=441.128;
a=50*sqrt(2.0)*pi;
T=0.001;    %采样率
w0=50*sqrt(2.0)*pi;
x=A*exp(-a*n*T)
close all
subplot(4,1,2);
stem(x);
title('理想采样信号序列');
t=0:0.001:0.05;
y=A*exp(-a*t)
subplot(4,1,1);
plot(t,y);
title('连续信号')
N=256;
f=-0.5:1/N:0.5-1/N;
X=fft(x,N);
subplot(4,1,3);
plot(f,fftshift(abs(X)));
title('fft频谱');
y=ifft(X);
subplot(4,1,4);
plot(real(x(1:N/8)));
title('ifft信号');