频谱分析

我是刚子仔 · 发表于 2014-4-4 00:08:30

求大神指导：帮我看看我的这个频谱分析的程序，其中振幅谱对着，为啥相位谱不对？

我是刚子仔 · 发表于 2014-4-11 09:35:42

chuxf 发表于 2014-4-4 18:33
你的代码没有错误，这也不是一个 Fortran 问题，而是算法问题。

对于采样率不是信号频率整数倍，或者整个 ...

好的，我试试，谢谢！

我是刚子仔 · 发表于 2014-4-12 21:01:45

chuxf 发表于 2014-4-4 18:33
你的代码没有错误，这也不是一个 Fortran 问题，而是算法问题。

对于采样率不是信号频率整数倍，或者整个 ...

您好，我按您说的改了一下程序，结果发现相位谱还是不太对。还有，如何将振幅谱中的振幅大小和实际的振幅大小一样？比如我这个程序中的振幅是1，结果在振幅谱中它的大小却不是1.
还望大神帮我把这程序改改，跪谢啊！

我是刚子仔 · 发表于 2014-4-12 21:03:43

chuxf 发表于 2014-4-4 18:33
你的代码没有错误，这也不是一个 Fortran 问题，而是算法问题。

对于采样率不是信号频率整数倍，或者整个 ...

哦，这是我修改后的程序

我是刚子仔 · 发表于 2014-4-16 22:28:35

chuxf 发表于 2014-4-12 23:57
第一：建议你把圆周率定义为一个常数，要保持有一定的精度。（3.14 太初略了）
real , parameter :: PI = 3 ...

这就是修改程序后做的相位谱，还是差得好远频率 50Hz

[通用算法] 频谱分析