fft_ideal_signal_sample_pulse.m

MATLAB数字信号处理的入门学习例程及实验结果

n=0:50; %定义序列的长度是50
A=444.128; %设置信号有关的参数
a=50*sqrt(2.0)*pi;
T=0.001; %采样率
w0=50*sqrt(2.0)*pi;
x=A*exp(-a*n*T).*sin(w0*n*T); %pi是MATLAB定义的π，信号乘可采用“.*”
close all %清除已经绘制的x(n)图形
subplot(3,1,1);stem(x); %绘制x(n)的图形
title('理想信号采样序列');

%下面用两种方法得出fft
% 1.使用数学计算得出的fft变换
k=0:50;
Wn=exp(-i*pi/25);     %旋转因子
X=x*(Wn.^(k'*n));       %x的fft变换－－X(N)
amgX=abs(X);            %求X的幅度
angleX=angle(X);        %求X的相位
subplot(3,1,2);stem(amgX);
title('方法1：x的fft变换的幅度谱');
subplot(3,1,3);stem(angleX);
title('方法1：x的fft变换的相位谱');

% % 2.使用fft函数得出的fft变换
% X=fft(x,50);
% amgX=abs(X);
% angleX=angle(X);
% subplot(3,1,2);stem(amgX);
% title('方法2：x的fft变换的幅度谱');
% subplot(3,1,3);stem(angleX);
% title('方法2：x的fft变换的相位谱');