fsint.m

我做了很久的

Fs=1000
N=1000
n=0:N-2
t=n/Fs
x=sin(2*pi*50*t)
subplot(221)
plot(t,x)
axis([0 0.1 -1 1])
xlabel('t(s)')
ylabel('x')
title('sin(2*pi*50*t)的时域图')
Nfft=1024
window=hamming(512)
noverlap=256
dflag='none'
[Pxx,f]=psd(x,Nfft,Fs,window,noverlap,dflag)
subplot(222)
plot(f,Pxx)
xlabel('f(HZ)')
ylabel('信号的幅频特性图')
grid