Matlab 快速傅里叶变换

来源:百度知道 编辑:UC知道 时间:2024/05/24 11:29:39
如何使用matlab中的fft函数对x(t)=sin(wt+pi/6)进行快速傅里叶变换,采样频率是原频率的8倍,截断长度N=16

t=0:15
w=10;
x=sin(w*t+pi/6);
X=fft(x,16)
Pyy = X.* conj(X)/16;
plot(Pyy)

直接调用MATLAB中的fft函数
假设你的基频为50Hz,自己可以设定;
clc;
clear all;
clf;
w=50;
fs=8*w;
t=0:1/fs:15/fs;
x=sin(w*t+pi/6);
y=fft(x,16);
n=0:length(y)-1;
f=fs*n/length(y);
plot(f,abs(y));
就可以看到频谱图了。
初次回答问题,请多关照

i=0:15;
thita=2*pi/8;
x=sin(thita+pi/6);
y=fft(x,16)
py=abs(y);
bar(py);