| 查看: 2701 | 回复: 8 | |||
[求助]
【求组】用matlab的fft求空间频谱
|
|
各位大神，小弟前些日子拿到一个片子，上面有微结构，想得到微结构的参数，然后就用白光干涉显微镜测了一下，得到了0.5mm*0.5mm视场下2500*2500的一个数组，数组里每个数据代表一个点的高度，也就是说把片子表面的形貌测出来了，但是由于镀的反射膜（白光干涉需要反射率较高）可能并不是非常平整，引入了误差，再加上片子保护的不算很好，可能也误差比较大，因此得到的数据并不能很清晰的看出来周期结构。因此想通过傅里叶分析来找一下其中的结构。取数组里的一行或者一列（记为列向量X）用matlab里的fft变换以后，得到一个相同大小的复数数组（记为列向量fre），取模以后可以得到一个实数数组（记为freamp）。 小弟的问题是： 1，画出freamp的图形以后的确可以看到一些连在一起的峰，那么这些峰对应的频率是多少？比如最高峰对应的坐标是50，那么对应的结构的周期（以毫米为单位）应该是多少？ 2，fft以后，得到的freamp，两端是低频，中间是高频是么？ 在网上看了不少东西，将的基本都是分析时间信号的，这种分析空间周期的没怎么看到，望各位大神答疑解惑，感激不尽！！！ |
回复此楼» 猜你喜欢
科研孤儿太难了
已经有9人回复
-
时间戳又变了
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教matlab傅里叶变换求解
已经有5人回复
- 求有关图像处理的一些建设性建议，多谢多谢多谢了！
已经有4人回复
- 用origin里的fft做了几张频谱图，请大神指教这些图都代表什么意思啊？
已经有10人回复
- comsol中的频谱图如何读取？？？
已经有6人回复
- 求助大虫，有做航迹关联matlab仿真的吗？
已经有5人回复
- 求助matlab/gui中画四个子图，出来的位置不对，求大虫解决
已经有8人回复
- 时间序列 振幅和周期提取
已经有5人回复
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有30人回复
- Matlab在数字信号处理中的应用
已经有137人回复
- FFT算法的一个问题请教
已经有18人回复
- 利用Matlab模拟空间滤波实验
已经有43人回复
- 求助随机信号处理方面的一个问题，20金币。
已经有10人回复
- matlab 傅里叶变换 频谱分析
已经有12人回复
- 关于MATLAB实现连续傅里叶变换绘图的问题
已经有3人回复
2楼2013-06-01 20:22:24
忆如沉
银虫 (小有名气)
- 应助: 20 (小学生)
- 金币: 460.5
- 帖子: 138
- 在线: 83.7小时
- 虫号: 2196842
- 注册: 2012-12-20
- 性别: GG
- 专业: 信号理论与信号处理
3楼2013-06-02 09:53:49
4楼2013-06-02 10:52:58
5楼2013-06-02 12:49:03
|
嗯，是的，我上面提到0.5mm得到了2500个数据，也就是说两个像素的实际距离是0.2um，也就是Fs为5000.我仔细看了fft的sample，里面的这句话，f = Fs/2*linspace(0,1,NFFT/2+1);也就是绘频率图时的横坐标，是归一化后乘以采样频率Fs的1/2，也就是说通过fft变化出来的频率谱的最高频率是采样频率的1/2，对么？（是不是依据奈奎斯特采样定律得到的？因为奈奎斯特定律说如果采样频率高于信号最高频率的2倍就能完整的还原出信号，所以反过来，对一个信号进行采样，能采样到的最高频率也就是采样频率的1/2？小弟对信号与系统方面的东西忘光了，见谅），感谢大神解释！ |
6楼2013-06-02 14:12:36
7楼2013-06-02 14:54:52
gzl9901
铁杆木虫 (文坛精英)
very good
- 应助: 3 (幼儿园)
- 金币: 6279.5
- 红花: 28
- 沙发: 2
- 帖子: 11263
- 在线: 1769.2小时
- 虫号: 1294464
- 注册: 2011-05-13
- 专业: 计算机软件
努力搞科研，祝早日出成果！ 8楼2013-06-02 15:36:48
9楼2013-06-02 21:58:28
40