| 查看: 2252 | 回复: 8 | ||
[求助]
【求组】用matlab的fft求空间频谱
|
|
各位大神，小弟前些日子拿到一个片子，上面有微结构，想得到微结构的参数，然后就用白光干涉显微镜测了一下，得到了0.5mm*0.5mm视场下2500*2500的一个数组，数组里每个数据代表一个点的高度，也就是说把片子表面的形貌测出来了，但是由于镀的反射膜（白光干涉需要反射率较高）可能并不是非常平整，引入了误差，再加上片子保护的不算很好，可能也误差比较大，因此得到的数据并不能很清晰的看出来周期结构。因此想通过傅里叶分析来找一下其中的结构。取数组里的一行或者一列（记为列向量X）用matlab里的fft变换以后，得到一个相同大小的复数数组（记为列向量fre），取模以后可以得到一个实数数组（记为freamp）。 小弟的问题是： 1，画出freamp的图形以后的确可以看到一些连在一起的峰，那么这些峰对应的频率是多少？比如最高峰对应的坐标是50，那么对应的结构的周期（以毫米为单位）应该是多少？ 2，fft以后，得到的freamp，两端是低频，中间是高频是么？ 在网上看了不少东西，将的基本都是分析时间信号的，这种分析空间周期的没怎么看到，望各位大神答疑解惑，感激不尽！！！ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教matlab傅里叶变换求解
已经有5人回复
- 求有关图像处理的一些建设性建议，多谢多谢多谢了！
已经有4人回复
- 用origin里的fft做了几张频谱图，请大神指教这些图都代表什么意思啊？
已经有10人回复
- comsol中的频谱图如何读取？？？
已经有6人回复
- 求助大虫，有做航迹关联matlab仿真的吗？
已经有5人回复
- 求助matlab/gui中画四个子图，出来的位置不对，求大虫解决
已经有8人回复
- 时间序列 振幅和周期提取
已经有5人回复
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有30人回复
- Matlab在数字信号处理中的应用
已经有136人回复
- 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
20