xiazaide.m

x=2*sin(4*pi*t)+5*cos(8*pi*t) 对其进行快速傅立叶变换 加入高斯噪声 再进行频谱分析

N=length(x); %x为电压信号数据
y=fft(x)/sqrt(N);
f=(1:N/2)*fs/N; %fs为采样频率
PSD=abs(y).^2;
plot(f, 2*PSD);