24小时热门版块排行榜    

查看: 378  |  回复: 2

liuyalei820

金虫 (正式写手)

[交流] 【求助】急。。。,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')
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

robert2020

银虫 (著名写手)

★ ★
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 ...

Truthfulness,Compassion,Tolerance,Beauty
2楼2010-06-26 23:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyalei820

金虫 (正式写手)


zzuwangshilei(金币+1):鼓励讨论 2010-06-27 14:13:43
引用回帖:
Originally posted by robert2020 at 2010-06-26 23:34:50:
波形很杂,不知楼主采样的时候是否正确?
程序如下:





你算的和我自己算的差不多,我想大概是采样有问题,也许采样频率有点低,导致采样数据不能正确反映系统的特性。。。。。
不过谢谢你。。。
3楼2010-06-27 14:01:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuyalei820 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +7 不知道叫什么! 2026-03-15 8/400 2026-03-15 20:11 by 棒棒球手
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 化学工程321分求调剂(南京工业,浙江工业) +3 大米饭! 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 求调剂,一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见