24小时热门版块排行榜    

查看: 2146  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

mei3014

金虫 (正式写手)

木有头屑

[求助] matlab txt文件数据提取转换

假设我有一个1.txt的数据文本,里面有这么一段:
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
我需要的数据是
4001         4.0000000000e+003        -2.25194e-001        -5.07143e-001
换言之就是4.0000000000e+003这个频率下的-2.25194e-001-5.07143e-001i 这个声压数据
然后重复刚才的过程,一个txt文件里需要提取60个4.0000000000e+003这个频率下对应的声压数据。再把这60个复数数据导入另一个txt文本,每个复数数据之间留个空格,那么就大功告成了。
请问各位大神,这matlab程序该怎么编写才能一步完成以上操作,谁帮下忙,非常感谢,如果做不到一步完成,提供些思路也可以。可以追加分。谢谢!
回复此楼

» 收录本帖的淘帖专辑推荐

matlab

» 猜你喜欢

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

一勤天下无难事,处世无难忍最高。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

xiegangmai: 帖中代码部分可使用“插入代码”功能 2012-11-02 19:56:34
code:
[a(:,1),a(:,2),a(:,3),a(:,4)]=textread('ddt.txt');  %读入数据,假设你的数据文件名为ddt.txt
ind=find(a(:,1)==4000);
b=a(ind,; % 此时已把ddt.txt中所有含4000的行都提取出来单独组成了矩阵b
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
4楼2012-11-02 17:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

mei3014

金虫 (正式写手)

木有头屑

X1——X2 50cm.txt(56.41MB)
http://kuai.xunlei.com/d/GOYWKZJKEDUD?p=130497
这是原始数据文件,内容较多
一勤天下无难事,处世无难忍最高。
2楼2012-11-02 15:35:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢应助 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');  %读入数据,假设你的数据文件名为ddt.txt
ind=find(a(:,1)==4000);
b=a(ind,; % 此时已把ddt.txt中所有含4000的行都提取出来单独组成了矩阵b

写到文件里,应该很容易了吧,自己弄吧
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
3楼2012-11-02 17:38:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑

引用回帖:
3楼: Originally posted by baobiao007 at 2012-11-02 17:38:09
=textread('ddt.txt');  %读入数据,假设你的数据文件名为ddt.txt
ind=find(a(:,1)==4000);
b=a(ind,; % 此时已把ddt.txt中所有含4000的行都提取出来单独组成了矩阵b

写到文件里,应该很容易了吧,自己弄吧...

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
程序比较复杂,在一楼所示程序内容这段上面还有一段是不规则的数据,如上所示,而且txt文本中内容很多,txt文本大小就有56M多了,可想而知。
另外我用你的程序运行语句:
[a(:,1),a(:,2),a(:,3),a(:,4)]=textread('X1_X2_45cm.txt');  %读入数据,假设你的数据文件名为ddt.txt。
matlab运行不下去了,提示错误为:
??? 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{:});
所以你编的程序不能应用于我的含有复杂内容的txt文件数据提取,只适合含有规则数据内容的数据提取,不过还是谢谢你的回答!
一勤天下无难事,处世无难忍最高。
5楼2012-11-03 09:32:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 3/150 2026-03-18 00:13 by PY-Hzb
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[考研] 材料与化工304求B区调剂 +7 邱gl 2026-03-11 8/400 2026-03-17 09:36 by 努力学习赚彩礼
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
请填处理意见