24小时热门版块排行榜

返回列表

liuyalei820

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 2159.6
散金: 598
红花: 1
帖子: 877
在线: 236.2小时
虫号: 996558
注册: 2010-04-14
性别: GG
专业: 导航、制导与传感技术

[交流] 【求助】急。。。，FFt变换已有2人参与

对x中数据的第一、二、三列数据有FFT，采样频率Fs=20Hz，发现MATLAB编程后所求频率不对，希望高手帮忙看看，程序如下：
clear
fs=20;
x=[-1.84955E-03,-.2829812,-.1720082
-.2237956,-.36991,.0554865
-.1979019,-.1405658,-4.253966E-02
-.1202208,-.238592,-.147964
-.2718839,-7.398201E-03,-8.322976E-02
-.1035748,-.2367424,-.2237956
-.1442649,-.2866803,-1.664595E-02
-.1794064,-.2163974,-.2404415
-.1942028,-.3273704,2.774325E-02
-.1202208,-.2607866,-.1202208
-.2237956,-.4142992,2.404415E-02
-3.144235E-02,-.147964,-3.884055E-02
-.1868046,-.3551136,-.036991
-.2367424,-.2829812,-.1923532
-9.802616E-02,-.1794064,-.2829812
-.2330433,-1.84955E-03,.0221946
-.1091235,-5.54865E-03,-.1424154
-.147964,-.3070253,-8.138021E-02
-.1535127,-1.664595E-02,-.2570875
-.147964,-.4697857,-.1091235
-.1794064,-7.953066E-02,-.1350172
-8.322976E-02,-.3273704,-.2367424
-.16461,-.1350172,-.0110973
-.1183712,-.2145478,-.295928
-.2367424,-7.768111E-02,-.0628847
-.1350172,-.3717596,-.2237956
-.16461,-.4753344,-.1091235
-.1183712,-4.623875E-02,-.2607866
-.1997514,-.2533884,-.1794064
-6.103516E-02,-3.144235E-02,-9.062796E-02
-.1035748,-.3273704,5.54865E-03
-.1146721,-.2811316,-.0036991
-8.692886E-02,-.386556,-.238592
-.1498136,.1035748,-.1923532
-.1183712,-.1442649,-.3828569
-1.84955E-03,-.1794064,-.1202208
-.2367424,-.1627604,-.3051758
-.1035748,-.4383434,-5.918561E-02
-.201601,-.3273704,-.1498136
-9.247751E-03,-.2533884,-.275583
-.2089992,-.1794064,-.1812559
-.2700343,-.2829812,-.1720082
-.1979019,-.2071496,.0258937
-.2237956,-.3181226,-.1461145
-.1405658,-.2071496,-2.034505E-02
-7.583155E-02,-.2903794,-.1553622
-.1775568,-.1720082,-.2034505
-.3199722,4.993786E-02,-.1923532
-.1035748,-.4401929,-5.178741E-02
-.221946,-.2071496,-3.144235E-02
-6.103516E-02,-.3606623,-.1294685
-.1331676,-.3421668,-.1350172
-.127619,-9.432706E-02,-4.623875E-02
-3.144235E-02,-.3847064,-.1942028
-.184955,0,-.295928
-.238592,-.4457416,-.1627604
-.1997514,-8.322976E-02,-.2681848
-9.062796E-02,-.5863074,-.2126983
-.2977776,-.1257694,-.2607866
-.2274947,-.3569632,-.2311938
-.1202208,3.144235E-02,-4.993786E-02
-.1868046,-.3939542,-.1035748
-.1202208,-.3921046,-.1165217
-.1202208,-7.583155E-02,-.110973
-.1997514,-.2903794,-7.583155E-02
-.2774325,-.1183712,-.1960523
-.2200965,-.221946,-.2811316
-8.507931E-02,-.147964,-6.473426E-02
-.1350172,-.3810073,-.2644857
-.2459902,-.4309452,-7.398201E-02
-9.617661E-02,-.1794064,-.16461
-.1202208,-5.178741E-02,-8.138021E-02
-.0554865,.1627604,-.1054244
-.1794064,-.6917318,-1.84955E-03
-.2515388,-.3255208,-6.843336E-02
-.1572118,-.4161488,-.1424154
-.1775568,-1.294685E-02,-.1794064
-.1997514,-.2330433,-.1738577
-.1942028,-.2663352,-1.84955E-03
-1.664595E-02,-8.322976E-02,-.2108487
-.2293442,-.4827326,5.178741E-02
-4.253966E-02,.1035748,-.2515388
-.1701586,-.1794064,-.1905037
-.2515388,-.1294685,-.1350172
-.1202208,-.4364938,-.1664595
-.127619,-.5178741,-4.623875E-02
-.2404415,-.4531398,-.3310695
-.1627604,-8.507931E-02,-.1072739
-.2108487,-.238592,-.2829812
-8.138021E-02,-.386556,-1.84955E-03
-1.294685E-02,-.1498136,-9.247751E-02
-.2681848,-.2977776,-3.514145E-02
-.1313181,-.0147964,-.1794064
-.2607866,-.2422911,-.2848307
-.2163974,-5.178741E-02,-7.213246E-02
-.2478397,-1.664595E-02,-.238592
-.2256451,.0110973,-.0443892
-.2459902,-.1405658,-.2459902
-.238592,-.1905037,-1.84955E-03
-.1683091,-.3310695,-.1535127
-9.617661E-02,-.201601,-.2163974
-6.843336E-02,-.4161488,-8.322976E-02
-.2367424,-.2811316,-.184955
-.1757073,-.5326704,-.0258937
-.2367424,-.2422911,-6.103516E-02
-.127619,-.5197236,-9.802616E-02
-5.363696E-02,-6.103516E-02,-8.322976E-02
-.2274947,-.3514145,3.884055E-02
-.1072739,-.0221946,-9.062796E-02
-.1331676,-.1775568,-.1923532
];
N=length(x);
n=0:N-1;t=n/fs;
y=fft(x(:,1),N);
mag=abs(y);
f=n*fs/N;
subplot(211),plot(t,x(:,1));xlabel('t/s');ylabel('x')
subplot(212),plot(f(1:N/2),mag(1:N/2)*2/N);xlabel('w/Hz');ylabel('x')

回复此楼

» 猜你喜欢

Cas 72-43-5需要30g，定制合成，能接单的留言已经有7人回复
磺酰氟产物，毕不了业了！已经有5人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有5人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2010-06-26 15:47:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

robert2020

银虫 (著名写手)

仿真EPI: 2
应助: 0 (幼儿园)
贵宾: 1.307
金币: 292.3
散金: 9143
红花: 19
沙发: 6
帖子: 1950
在线: 620.1小时
虫号: 861676
注册: 2009-10-02

★ ★
nono2009(金币+1):辛苦了。 2010-06-27 07:15:59
liuyalei820(金币+1): 2010-06-27 14:01:40
zzuwangshilei(金币+1):积极参与 2010-06-27 14:13:33

波形很杂，不知楼主采样的时候是否正确？
程序如下：

引用回帖:

m=2^12;
y=fft(x(:,1),m);
mag=abs(y);
f=(0:m-1)*fs/m;
plot(f(1:m/2),mag(1:m/2)*2/m);

引用回帖:

Originally posted by liuyalei820 at 2010-06-26 15:47:59:
对x中数据的第一、二、三列数据有FFT，采样频率Fs=20Hz，发现MATLAB编程后所求频率不对，希望高手帮忙看看，程序如下：
clear
fs=20;
x=[-1.84955E-03,-.2829812,-.1720082
-.2237956,-.36991,.0554865
-.19 ...