版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

mei3014

金蟲 (正式寫手)

木有頭屑

應(yīng)助: 14 (小學(xué)生)
金幣: 2957.5
散金: 274
帖子: 541
在線: 335小時
蟲號: 1333109
注冊: 2011-06-28
性別: GG
專業(yè): 信號理論與信號處理

[求助] matlab txt文件數(shù)據(jù)提取轉(zhuǎn)換

假設(shè)我有一個1.txt的數(shù)據(jù)文本，里面有這么一段：
3996    3.9950000000e+003 -1.46929e-005    1.39538e-005
3997    3.9960000000e+003 -1.65129e-005    6.65343e-005
3998    3.9970000000e+003    4.05215e-005 -2.06176e-004
3999    3.9980000000e+003 -2.65027e-005    1.81863e-004
4000    3.9990000000e+003    1.12635e-001    2.53541e-001
4001    4.0000000000e+003 -2.25194e-001 -5.07143e-001
4002    4.0010000000e+003    1.12585e-001    2.53582e-001
4003    4.0020000000e+003 -1.29527e-004 -7.78130e-005
4004    4.0030000000e+003    1.62932e-004    7.54291e-005
我需要的數(shù)據(jù)是
4001    4.0000000000e+003 -2.25194e-001 -5.07143e-001
換言之就是4.0000000000e+003這個頻率下的-2.25194e-001-5.07143e-001i 這個聲壓數(shù)據(jù)
然后重復(fù)剛才的過程，一個txt文件里需要提取60個4.0000000000e+003這個頻率下對應(yīng)的聲壓數(shù)據(jù)。再把這60個復(fù)數(shù)數(shù)據(jù)導(dǎo)入另一個txt文本，每個復(fù)數(shù)數(shù)據(jù)之間留個空格，那么就大功告成了。
請問各位大神，這matlab程序該怎么編寫才能一步完成以上操作，誰幫下忙，非常感謝，如果做不到一步完成，提供些思路也可以。可以追加分。謝謝！

回復(fù)此樓

一勤天下無難事，處世無難忍最高。

1樓 2012-11-02 14:17:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mei3014

金蟲 (正式寫手)

木有頭屑

應(yīng)助: 14 (小學(xué)生)
金幣: 2957.5
散金: 274
帖子: 541
在線: 335小時
蟲號: 1333109
注冊: 2011-06-28
性別: GG
專業(yè): 信號理論與信號處理

X1——X2 50cm.txt(56.41MB)
http://kuai.xunlei.com/d/GOYWKZJKEDUD?p=130497
這是原始數(shù)據(jù)文件，內(nèi)容較多

贊一下

回復(fù)此樓

一勤天下無難事，處世無難忍最高。

2樓2012-11-02 15:35:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時
蟲號: 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝應(yīng)助 2012-11-02 19:56:25
mei3014: 金幣+5, ★有幫助, 謝謝回帖！ 2012-11-03 09:33:51

[a(:,1),a(:,2),a(:,3),a(:,4)]=textread('ddt.txt'); %讀入數(shù)據(jù)，假設(shè)你的數(shù)據(jù)文件名為ddt.txt
ind=find(a(:,1)==4000);
b=a(ind,

; % 此時已把ddt.txt中所有含4000的行都提取出來單獨組成了矩陣b

寫到文件里，應(yīng)該很容易了吧，自己弄吧

贊一下

回復(fù)此樓

我同意叔本華的觀點，人們投身藝術(shù)和科學(xué)領(lǐng)域的強(qiáng)烈愿望之一就是逃離痛苦、殘酷和枯燥無味的現(xiàn)實生活，逃離自己飄忽不定的七情六欲的桎梏。--愛因斯坦

3樓2012-11-02 17:38:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時
蟲號: 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

【答案】應(yīng)助回帖

xiegangmai: 帖中代碼部分可使用“插入代碼”功能 2012-11-02 19:56:34

code:
[a(:,1),a(:,2),a(:,3),a(:,4)]=textread('ddt.txt'); %讀入數(shù)據(jù)，假設(shè)你的數(shù)據(jù)文件名為ddt.txt
ind=find(a(:,1)==4000);
b=a(ind,; % 此時已把ddt.txt中所有含4000的行都提取出來單獨組成了矩陣b

贊一下

回復(fù)此樓

4樓2012-11-02 17:39:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mei3014

金蟲 (正式寫手)

木有頭屑

應(yīng)助: 14 (小學(xué)生)
金幣: 2957.5
散金: 274
帖子: 541
在線: 335小時
蟲號: 1333109
注冊: 2011-06-28
性別: GG
專業(yè): 信號理論與信號處理

引用回帖:

3樓: Originally posted by baobiao007 at 2012-11-02 17:38:09
=textread('ddt.txt'); %讀入數(shù)據(jù)，假設(shè)你的數(shù)據(jù)文件名為ddt.txt
ind=find(a(:,1)==4000);
b=a(ind,

; % 此時已把ddt.txt中所有含4000的行都提取出來單獨組成了矩陣b

寫到文件里，應(yīng)該很容易了吧，自己弄吧...

Header Size:    79
Pulse Version:    80
Running Pulse Version: PULSE LabShop v. 15.1.0.15 - 2010-11-08
Decimal Symbol: .
Date Format: 'yyyy-MM-dd'
Time Format: 'HH:mm:ss.mmm'
Data Type: Complex
Slice: False
Z-Axis type: Linear
Z-Axis size: 1
Z-Axis unit:
Z-Axis first value:    0.0000000000e+000
Z-Axis delta:    1.0000000000e+000
X-Axis type: Linear
X-Axis size: 6401
X-Axis unit: Hz
X-Axis first value:    0.0000000000e+000
X-Axis delta:    1.0000000000e+000

AmplitudeUnit: Pa
Analyzer: FFT_Spectrum_Averaging
AnalyzerName: FFT Analyzer


Bipolar: True
CenterFrequency:    3.2000000000e+003




dBReference:    2.0000000000e-005



Domain: 1
Function: Fourier Spectrum

InputRange:    8.5304914740e+002


MultiBufferName: Input
NBW:    1.5000000000e+000

OverlapFailed: False

Power: False
PowerUnit: Pa^2
Ratio: False
RecordLength:    1.0000000000e+000

Ref_Signal: Signal 1





Signal: Signal 1

SignalUnit: Pa



Title: Fourier Spectrum(Signal 1) - Input
Title1: Working : Input : Input : FFT Analyzer


Weight: 0
ZoomOn: False

TemplateName: Working
MeasurementName: Input
MeasurementSlice: False



EnvelopeOn: False



Z-index:      1
Date:      '2012-11-01'
Time:      '14:25:27.624'
Relative time:         0.00000e+000
Z-axis:         0.0000000000e+000
1    0.0000000000e+000 -2.69714e-005    0.00000e+000
程序比較復(fù)雜，在一樓所示程序內(nèi)容這段上面還有一段是不規(guī)則的數(shù)據(jù)，如上所示，而且txt文本中內(nèi)容很多，txt文本大小就有56M多了，可想而知。
另外我用你的程序運(yùn)行語句：
[a(:,1),a(:,2),a(:,3),a(:,4)]=textread('X1_X2_45cm.txt');  %讀入數(shù)據(jù)，假設(shè)你的數(shù)據(jù)文件名為ddt.txt。
matlab運(yùn)行不下去了，提示錯誤為：
??? Error using ==> dataread
Trouble reading number from file (row 1, field 1) ==> Header Size:    79         \n

Error in ==> textread at 176
[varargout{1:nlhs}]=dataread('file',varargin{:});
所以你編的程序不能應(yīng)用于我的含有復(fù)雜內(nèi)容的txt文件數(shù)據(jù)提取，只適合含有規(guī)則數(shù)據(jù)內(nèi)容的數(shù)據(jù)提取，不過還是謝謝你的回答！

贊一下

回復(fù)此樓

一勤天下無難事，處世無難忍最高。

5樓2012-11-03 09:32:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiangdonal

鐵蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 13.6
帖子: 15
在線: 12.3小時
蟲號: 1556791
注冊: 2011-12-30
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝建議 2012-11-04 19:23:14
mei3014: 金幣+1, ★有幫助, 謝謝回帖 2012-11-07 13:07:53

用fgetl這種比較底層一點的函數(shù)，一行一行讀取文件，搜索到你需要的就可以了

贊一下

回復(fù)此樓

6樓2012-11-03 11:25:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mei3014

金蟲 (正式寫手)

木有頭屑

應(yīng)助: 14 (小學(xué)生)
金幣: 2957.5
散金: 274
帖子: 541
在線: 335小時
蟲號: 1333109
注冊: 2011-06-28
性別: GG
專業(yè): 信號理論與信號處理

引用回帖:

6樓: Originally posted by jiangdonal at 2012-11-03 11:25:35
用fgetl這種比較底層一點的函數(shù)，一行一行讀取文件，搜索到你需要的就可以了

一行行讀？50多M，太費(fèi)時間！

贊一下

回復(fù)此樓

一勤天下無難事，處世無難忍最高。

7樓2012-11-03 12:23:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

★
xiegangmai: 金幣+1, 謝謝參與 2012-11-07 20:30:35

你還是用excel弄一下更快

贊一下

回復(fù)此樓

showmethemoney

8樓2012-11-06 10:43:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shcnxjy

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 561.1
帖子: 288
在線: 13.5小時
蟲號: 1783333
注冊: 2012-04-27
專業(yè): 信號理論與信號處理

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
mei3014: 金幣+5, ★有幫助, 謝謝回帖，不過你這textread用的不對，在我這讀不出來！ 2012-11-07 13:04:45
xiegangmai: 金幣+2, 謝謝參與 2012-11-07 20:30:46

textread函數(shù)的功能是十分強(qiáng)大的，如果沒記錯的話（具體用法請參考matlab幫助文件）
[data1,data2,data3,data4]=textread('X1_X2_45cm.txt','%n%n%n%n',N1,'headerlines',N2);
就可以跳過開頭的N2行，而N1表示重復(fù)該讀取動作N1次，
這樣差不多就可以解決你的文件格式問題了！

贊一下

回復(fù)此樓

9樓2012-11-06 16:36:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

satosi

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 443.6
紅花: 1
帖子: 38
在線: 17.8小時
蟲號: 1567064
注冊: 2012-01-06
專業(yè): 通信理論與系統(tǒng)

【答案】應(yīng)助回帖

★ ★ ★
mei3014: 金幣+2, ★有幫助, 謝謝回帖，最好能有具體程序，用法還是比較復(fù)雜的！ 2012-11-07 13:06:01
xiegangmai: 金幣+1, 謝謝參與 2012-11-07 20:30:53

引用回帖:

9樓: Originally posted by shcnxjy at 2012-11-06 16:36:40
textread函數(shù)的功能是十分強(qiáng)大的，如果沒記錯的話（具體用法請參考matlab幫助文件）
=textread('X1_X2_45cm.txt','%n%n%n%n',N1,'headerlines',N2);
就可以跳過開頭的N2行，而N1表示重復(fù)該讀取動作N1次，
這樣差 ...

textread可以試試。還有textscan，這個函數(shù)比textread要有一定的優(yōu)勢，尤其是在處理大文件時。樓主不妨都試試。具體用法可以去參考help。

贊一下

回復(fù)此樓

10樓2012-11-07 10:48:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 mei3014 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 284求調(diào)劑 +10	Zhao anqi 2026-03-22	10/500	2026-03-24 00:08 by Equinoxhua
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 材料調(diào)劑 +3	匹克i 2026-03-23	3/150	2026-03-23 23:18 by peike
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +4	小羊睡著了? 2026-03-23	4/200	2026-03-23 21:42 by hyzs6688
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5	葵梓衛(wèi)隊 2026-03-18	7/350	2026-03-23 16:26 by lingjue
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-17	10/500	2026-03-23 13:05 by 醉在風(fēng)里
[考研] 276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷 +8	材料學(xué)257求調(diào)劑 2026-03-23	9/450	2026-03-23 13:01 by ztnimte
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 307求調(diào)劑 +11	冷笙123 2026-03-17	11/550	2026-03-22 20:16 by edmund7
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 305分求調(diào)劑（食品工程） +4	Sxy112 2026-03-21	6/300	2026-03-22 15:26 by 無懈可擊111
[考研] 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑 +3	Grand777 2026-03-21	3/150	2026-03-21 19:23 by 簡之-
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky

24小時熱門版塊排行榜

[求助] matlab txt文件數(shù)據(jù)提取轉(zhuǎn)換

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖