自己看吧,但愿你能理解,而不是COPY完事
-------------------------------------------------------
clear;
close all;
N=50;
t=1:N;
x=zeros(1,50);
x(1)=1; %单位冲击
figure;
subplot(2,1,1);
stem(x);
X=fft(x);
magX=abs(X);
subplot(2,1,2);
plot(t(1:N/2),magX(1:N/2));
N=1024; %信号的抽样点数
fs=100; %抽样频率,高于信号中最高频率的2倍
t=(0:N-1)/fs;
f=(0:N-1)*fs/N;
x0=sign(t); %阶跃
x1=sign(sign(8-t)+1); %单个方波,宽度为8个采样周期
x2=square(2*pi*5*t,50); %周期方波,5Hz,50%占空比
x3=exp(-t.^2);
y0=abs(fft(x0)); %FFT后求模
figure;
subplot(2,1,1)
plot(t,x0);
subplot(2,1,2)
plot(f,y0);
y1=abs(fft(x1)); %FFT后求模
figure;
subplot(2,1,1)
plot(t,x1);
subplot(2,1,2)
plot(f,y1);
y2=abs(fft(x2));
figure;
subplot(2,1,1)
plot(t,x2);
subplot(2,1,2)
plot(f,y2);
y3=abs(fft(x3));
figure;
subplot(2,1,1)
plot(t,x3);
subplot(2,1,2)
plot(f,y3);