应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab txt文件数据提取转换
111/2返回列表12下一页
查看: 2079  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
一勤天下无难事，处世无难忍最高。
1楼 2012-11-02 14:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

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

jiangdonal

铁虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢建议 2012-11-04 19:23:14
mei3014: 金币+1, ★有帮助, 谢谢回帖 2012-11-07 13:07:53
用fgetl这种比较底层一点的函数，一行一行读取文件，搜索到你需要的就可以了
赞一下 回复此楼
6楼2012-11-03 11:25:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑
引用回帖:
6楼: Originally posted by jiangdonal at 2012-11-03 11:25:35
用fgetl这种比较底层一点的函数，一行一行读取文件，搜索到你需要的就可以了

一行行读？50多M，太费时间！
赞一下 回复此楼
一勤天下无难事，处世无难忍最高。
7楼2012-11-03 12:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★
xiegangmai: 金币+1, 谢谢参与 2012-11-07 20:30:35
你还是用excel弄一下更快
赞一下 回复此楼
showmethemoney
8楼2012-11-06 10:43:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shcnxjy

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
mei3014: 金币+5, ★有帮助, 谢谢回帖，不过你这textread用的不对，在我这读不出来！ 2012-11-07 13:04:45
xiegangmai: 金币+2, 谢谢参与 2012-11-07 20:30:46
textread函数的功能是十分强大的，如果没记错的话（具体用法请参考matlab帮助文件）
[data1,data2,data3,data4]=textread('X1_X2_45cm.txt','%n%n%n%n',N1,'headerlines',N2);
就可以跳过开头的N2行，而N1表示重复该读取动作N1次，
这样差不多就可以解决你的文件格式问题了！
赞一下 回复此楼
9楼2012-11-06 16:36:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

satosi

银虫 (初入文坛)

【答案】应助回帖

★ ★ ★
mei3014: 金币+2, ★有帮助, 谢谢回帖，最好能有具体程序，用法还是比较复杂的！ 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函数的功能是十分强大的，如果没记错的话（具体用法请参考matlab帮助文件）
=textread('X1_X2_45cm.txt','%n%n%n%n',N1,'headerlines',N2);
就可以跳过开头的N2行，而N1表示重复该读取动作N1次，
这样差 ...

textread可以试试。还有textscan，这个函数比textread要有一定的优势，尤其是在处理大文件时。楼主不妨都试试。具体用法可以去参考help。
赞一下 回复此楼
10楼2012-11-07 10:48:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mei3014 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭