单面的振幅傅里叶谱

                             D = dlmread (“Muhsin A00.txt”,' \ t '1 0);
                            
                             t = D (: 1);%时间(年代)
                            
                             = D (:, 2);%加速度(g)
                            
                             L =长度(t);
                            
                             Ts = t - t (1) (2);%采样间隔(秒)
                            
                             Fs = 1 / Ts;%采样频率(赫兹)
                            
                             Fn = f / 2;%奈奎斯特频率(赫兹)
                            
                             自由贸易协定= fft (a) / L;%傅里叶变换(比例)
                            
                             阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量
                            
                             4 = 1:长度(艘);%指数向量
                            
                             图(1)
                            
                             情节(阵线,abs (FTa (Iv)) * 2)%的振幅图
                            
                             包含(的频率(赫兹))
                            
                             ylabel (“振幅(g)”)
                            
                             网格

编辑 - - - - - - 改变了 “t” 计算容纳 ‘t (1) = 0.01 和 t(结束)= 0的。

19日的评论
显示18个年长的评论隐藏18年长的评论

明星黾 2017年10月12日

我的荣幸。

所不同的是,你的数据是在因为某些原因。我用它们作为你张贴他们。

试试这个(限制数据的非零值 “t” ):

                                  D = dlmread (“Muhsin A00.txt”,' \ t '1 0);
                                 
                                  t = D (: 1);%时间(年代)
                                 
                                  = D (:, 2);%加速度(g)
                                 
                                  t_end =找到(t = = 0, 1,“第一”);去年零' t ' %
                                 
                                  idx = 1: t_end-1;%指数为零' t '
                                 
                                  L =长度(t (idx));
                                 
                                  Ts = t - t (1) (2);%采样间隔(秒)
                                 
                                  Fs = 1 / Ts;%采样频率(赫兹)
                                 
                                  Fn = f / 2;%奈奎斯特频率(赫兹)
                                 
                                  自由贸易协定= fft ((idx)) / L;%傅里叶变换(比例)
                                 
                                  阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量
                                 
                                  4 = 1:长度(艘);%指数向量
                                 
                                  图(1)
                                 
                                  情节(阵线,abs (FTa (Iv)) * 2)%的振幅图
                                 
                                  包含(的频率(赫兹))
                                 
                                  ylabel (“振幅(g)”)
                                 
                                  网格

我不确定如果你想保留零或消除它们。我在这段代码中消除。注意,它也会影响频率向量。如果你想在你的信号提高频率分辨率,使用第二个 “NFFT” 参数 fft 用MATLAB函数:

                                  D = dlmread (“Muhsin A00.txt”,' \ t '1 0);
                                 
                                  t = D (: 1);%时间(年代)
                                 
                                  = D (:, 2);%加速度(g)
                                 
                                  t_end =找到(t = = 0, 1,“第一”);去年零' t ' %
                                 
                                  idx = 1: t_end-1;%指数为零' t '
                                 
                                  L =长度(t (idx));
                                 
                                  Ts = t - t (1) (2);%采样间隔(秒)
                                 
                                  Fs = 1 / Ts;%采样频率(赫兹)
                                 
                                  Fn = f / 2;%奈奎斯特频率(赫兹)
                                 
                                  NFFT = 2 ^ nextpow2 (t_end-1);
                                 
                                  自由贸易协定= fft ((idx) NFFT) / L;%傅里叶变换(比例)
                                 
                                  阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量
                                 
                                  4 = 1:长度(艘);%指数向量
                                 
                                  图(1)
                                 
                                  情节(阵线,abs (FTa (Iv)) * 2)%的振幅图
                                 
                                  包含(的频率(赫兹))
                                 
                                  ylabel (“振幅(g)”)
                                 
                                  网格

这是正确的计算方法和策划你的信号的傅里叶变换。奈奎斯特频率是 50 赫兹,最大可分解的频率的采样频率 One hundred. 赫兹。

Muhsin 2017年10月15日

编辑:Muhsin 2017年10月15日

请检查为什么下面的脚本不工作吗?它实际上是在行503 - 1903之间。我做错了什么。谢谢你！

D_main = dlmread (A00。txt ', ' \ t ', 1 0);

D = D_main (503:1903:);

t = D (: 1);%时间(年代)

= D (:, 2);%加速度(g)

t_end =找到(t = = 0, 1, '第一次');去年零' t ' %

idx = 1: t_end-1;%指数为零' t '

L =长度(t (idx));Ts = t - t (503) (504);%采样间隔(秒)

Fs = 1 / Ts;%采样频率(赫兹)

Fn = f / 2;%奈奎斯特频率(赫兹)

自由贸易协定= fft ((idx)) / L;%傅里叶变换(比例)

阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量

4 = 1:长度(艘);%指数向量图(1)

情节(阵线,abs (FTa (Iv)) * 2) %片面的振幅阴谋

包含(频率(赫兹))

ylabel(振幅(g))

网格

明星黾 2017年10月18日

的 smoothdata 函数可能会工作。(它了 R2017a ,因为我不知道你有什么版本我没有表明它。)

你会使用它的绝对值单面傅里叶转换数据。

登录置评。

答案(1)

图像分析 2017年10月12日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/360889-single-sided-amplitude-fourier-spectrum answer_285488

这不正是pwelch()呢?(信号处理工具箱pwelch所需)。

1评论
显示没有隐藏没有

图像分析 2017年10月12日

您可能还喜欢光谱图()函数。

或信号分析仪的应用程序的应用程序选项卡上带的工具。

登录置评。

在回答这个问题。

类别

信号处理信号处理工具箱光谱分析光谱测量

找到更多的在光谱测量在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

单面的振幅傅里叶谱

0评论
显示1年长的评论藏1年长的评论

接受的答案

19日的评论
显示18个年长的评论隐藏18年长的评论

答案(1)

1评论
显示没有隐藏没有

另请参阅

类别

标签

社区寻宝

单面的振幅傅里叶谱

0评论 显示1年长的评论藏1年长的评论

接受的答案

19日的评论 显示18个年长的评论隐藏18年长的评论

答案(1)

1评论 显示没有隐藏没有

另请参阅

类别

标签

社区寻宝

0评论
显示1年长的评论藏1年长的评论

19日的评论
显示18个年长的评论隐藏18年长的评论

1评论
显示没有隐藏没有