matlab-傅立叶

fs=100;
x=-2:1/fs:2;
y=sin(3*pi*x);
z=rectpuls(x);
figure;plot(x,y,x,z,':r');

my=abs(fft(y));
mz=abs(fft(z));
my=my/max(my); %归一化
mz=mz/max(mz); %归一化

f=(0:1/length(x):1)*fs;
figure;plot(f(1:fs/2),my(1:fs/2),f(1:fs/2),mz(1:fs/2),':r');
%为什么这么做？
%还记得Nyquist采样定理吗？哎～～你对了，就是这个原因。
%频谱图同理。

my=20*log10(my+eps);mz=20*log10(mz+eps);
figure;plot(f(1:fs/2),my(1:fs/2),f(1:fs/2),mz(1:fs/2),':r');

 

 

A=1HI=pi/3; Freq=1.5; %信号的参数
DT=0.1; %信号采样间隔
N=1024; %信号采样点数
TV=[0:N-1]*DT; %采样时刻数组;

SignalVec=A*cos(Freq*2*pi*TV+PHI); %生成信号数据

FV=fft(SignalVec); %调用FFT计算傅立叶谱;

subplot(3,1,1);%转到最上面一个图
plot([0:N-1]/(N*DT), abs(FV)); %幅值谱

subplot(3,1,2);%在中间那张图上画
plot([0:N-1]/(N*DT), real(FV)); %实部

subplot(3,1,3);%在下面的那张图上画
plot([0:N-1]/(N*DT), imag(FV)); %实部