24小时热门版块排行榜    

查看: 2143  |  回复: 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

木虫 (职业作家)

中国特色

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖
查看全部 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

木虫 (职业作家)

中国特色

【答案】应助回帖

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268求调剂 +7 好运连绵不绝 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 289求调剂 +6 步川酷紫123 2026-03-11 6/300 2026-03-17 10:23 by Sammy2
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见