| 查看: 376 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】傅里叶变换
|
|||
|
我们用的是采集卡采集连续信号，采样率25MHZ，采集数据10000个，采集完后保存到文本文档，这10000个点只有100个点是信号，其余的都是幅值0左右的噪声，我想做这个信号的傅里叶变换，请教这个程序哪里有问题 数据已经读取到h1，上面的程序不再写了 for i=1:dtnum x2=h1(i,[1000:2000]); k=10000; x2k=fft(x2,20000); y=abs(x2k); f=(0:length(x2k)-1)'*2.5*10^8/length(x2k); subplot(2,1,2) plot(f,y); title('频谱') axis([0, 1*10^7, 0, 2*10^4]) xlabel('Hz') i=i+1; end |
回复此楼» 猜你喜欢
面上本子正文33页，违规吗？会被低分嘛？
已经有13人回复
-
评审有感
已经有23人回复
-
提交了我也来说说感想
已经有6人回复
-
Sci. Bull. 悲剧经验
已经有7人回复
-
今年审到国自然15份，谈谈感受
已经有20人回复
-
国自然评分
已经有4人回复
-
反应很差，大量原料没有反应
已经有6人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
qianbh
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 28061.8
- 红花: 1
- 帖子: 2525
- 在线: 238.5小时
- 虫号: 245509
- 注册: 2006-04-23
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★ ★
loulou85(金币+2,VIP+0):没用过vb呀 最好是matlab的 3-29 09:55
loulou85(金币+6,VIP+0):333333 4-7 20:36
loulou85(金币+2,VIP+0):没用过vb呀 最好是matlab的 3-29 09:55
loulou85(金币+6,VIP+0):333333 4-7 20:36
|
提供一个本人用过的VB程序 '☆供分析调用的程序——FOUR1的过程☆ Public Sub FOUR1(DATA() As Double, NN As Integer, ISIGN As Integer) N = 2 * NN j = 1 For I = 1 To N Step 2 If j > I Then TempR = DATA(j) TempI = DATA(j + 1) DATA(j) = DATA(I) DATA(j + 1) = DATA(I + 1) DATA(I) = TempR DATA(I + 1) = TempI End If M = N / 2 While M >= 2 And j > M j = j - M M = M / 2 Wend j = j + M Next I Mmax = 2 While N > Mmax Istep = 2 * Mmax Theta = 6.28318530717959 / (ISIGN * Mmax) Wpr = -2# * Sin(0.5 * Theta) ^ 2 Wpi = Sin(Theta) Wr = 1# Wi = 0# For M = 1 To Mmax Step 2 For I = M To N Step Istep j = I + Mmax TempR = CSng(Wr) * DATA(j) - CSng(Wi) * DATA(j + 1) TempI = CSng(Wr) * DATA(j + 1) + CSng(Wi) * DATA(j) DATA(j) = DATA(I) - TempR DATA(j + 1) = DATA(I + 1) - TempI DATA(I) = DATA(I) + TempR DATA(I + 1) = DATA(I + 1) + TempI Next I Wtemp = Wr Wr = Wr * Wpr - Wi * Wpi + Wr Wi = Wi * Wpr + Wtemp * Wpi + Wi Next M Mmax = Istep Wend End Sub |
2楼2009-03-29 08:59:20
3楼2009-03-29 10:04:05
4楼2009-04-03 15:59:53
5