[经验]

使用Matlab对噪声信号进行FFT变换

李郝荫

2019-3-22 09:45:34

5131 FFT matlab

Fs = 1000; % 采样频率

复制代码

先用上述代码生成一个1Hz，1.5s的采样信号。
然后我们加入一个正常的幅值0.7频率120的正弦波信号

S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);

复制代码

为该信号加入零均值白噪声，方差为4

X = S + 2*randn(size(t));

复制代码

画出时域上的噪声信号

plot(1000*t(1:50),X(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('t (milliseconds)')
ylabel('X(t)')

复制代码

计算该信号的FFT

Y = fft(X);

复制代码

计算双边谱P2，然后计算单边谱P1和均值信号长度L

P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);

复制代码

定义时域f，并画出P1

f = Fs*(0:(L/2))/L;
plot(f,P1)
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')

复制代码

对原信号进行傅里叶变换

Y = fft(S);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
plot(f,P1)
title('Single-Sided Amplitude Spectrum of S(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')

复制代码

只看该作者 淘帖0 举报

2 条评论

只看该作者

李郝荫 2019-3-22 15:50:55

沙发

谢谢分享

赞回复

林昕宇 2020-9-20 19:04:32

板凳

非常不错的资料感谢楼主分享

赞回复

声明：本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人，不代表电子发烧友网立场。文章及其配图仅供工程师学习之用，如有内容图片侵权或者其他问题，请联系本站作侵删。侵权投诉

发经验

精选推荐

怎么修改matlab与comsol联合的代码啊

1415 浏览 1 评论
200圆！求助大佬给一份VSG并网和离网模式的simulink仿真

1713 浏览 0 评论
MATLAB(3)--矩阵的引用(sub2ind、ind2sub、reshape函数使用)

2524 浏览 0 评论
求问MATLAB仿真KK接收器出现的问题？

2863 浏览 0 评论
浅析MATLAB和Simulink嵌入式视觉应用

4089 浏览 1 评论

快速回复 返回顶部 返回列表

关注微信公众号

电子发烧友网

电子发烧友论坛

社区合作: 刘勇; 联系电话：15994832713; 邮箱地址：liuyong@huaqiu.com

社区管理: elecfans短短; 微信：elecfans_666; 邮箱：users@huaqiu.com

返回 Matlab论坛

搜索历史

使用Matlab对噪声信号进行FFT变换

相关推荐

2 条评论

评论

精选推荐

站长推荐 /7