马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear;clc; Fs=1000;%采样频率 N=1000;%采样点数 Fs/N=t(采样总时间) %频率分辨率=Fs/N需满足关注频率为频率分辨率的整数倍 n=0:1/Fs:(Fs/N-1/Fs);%不是应该是n=0:1/Fs:(N-1/Fs)吗 xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+5*cos(2*pi*150*n)+rand(size(n));%数据 ff=n(1:N/2)*Fs;%不是应该是ff=n*Fs/N吗 subplot(421),plot(n,xn),title('时域图'); xlabel('时间'); ylabel('幅值'); %%fft y=fft(xn,N); y1=abs(y)*2/N; plot_y=y1(1:(N/2)); %ploty2=plot_y/sqrt(2); subplot(422),plot(ff,plot_y),title('频谱图(FFT)'); %hold on;stem(ff,plot_y,'r'); xlabel('频率/Hz'); ylabel('幅值');
|