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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 新手求助matlab读入文件格式设置及画三维图
91/1返回列表
查看: 1134  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

godchenkun

银虫 (正式写手)

[求助] 新手求助matlab读入文件格式设置及画三维图 已有1人参与

matlab菜鸟一枚，求教高手指导：
从如下两种格式的txt文件读入文件：
格式一：
s axis  -2.50000E+01 to -2.45000E+01
t axis  -2.50000E+01 to -2.45000E+01
            1.00537E-06   0.4723

s axis  -2.50000E+01 to -2.45000E+01
t axis  -2.45000E+01 to -2.40000E+01
            1.03832E-06   0.4831
...省略很多行

格式二：
direct ray-trace source image flux
s axis  -2.50000E+01 to -2.45000E+01
t axis  -2.50000E+01 to -2.45000E+01
            1.00537E-06   0.4723

direct ray-trace source image flux
s axis  -2.50000E+01 to -2.45000E+01
t axis  -2.45000E+01 to -2.40000E+01
            1.03832E-06   0.4831
...省略很多行

目标：s axis行的两个数值取平均作为x坐标，t axis行两数值取平均作为y坐标，t axis行下一行第一个数据作为z坐标，画三维曲面。
          我的想法是将x、y、z坐标存入三个数组，然后mesh或surface，但不知道怎么写读入文件的格式来达到上述目的。
         两个文件样本在附件，请高手教我，谢谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-04-20 10:37:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

godchenkun

银虫 (正式写手)
纠正一下，是将xyz存入三个矩阵，不是数组。
赞一下 回复此楼
2楼2015-04-20 10:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
godchenkun: 金币+5, ★★★★★最佳答案 2015-04-20 14:20:12
假定第一个文件名Data1.txt，第二个文件名Data2.txt，如下的代码分别读取这两个文件中的数字，存在MyData1和MyData2中，每一组6个数，接下来的你自己处理一下吧。

fid1 = fopen('Data1.txt');
MyData1 = textscan(fid1,'%*s %*s %f %*s %f \r\n %*s %*s %f %*s %f \r\n %f %f \r\n');
fclose(fid1);

fid2 = fopen('Data2.txt');
N = 1;
while ~feof(fid2)
    fgetl(fid2);
    MyData2{N} = textscan(fid1,'%*s %*s %f %*s %f \r\n %*s %*s %f %*s %f \r\n %f %f \r\n');
    N = N + 1;
end
fclose(fid2);
赞一下 回复此楼
知其然，知其所以然。
3楼2015-04-20 12:59:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

godchenkun

银虫 (正式写手)
引用回帖:
3楼: Originally posted by mygt_hit at 2015-04-20 12:59:25
假定第一个文件名Data1.txt，第二个文件名Data2.txt，如下的代码分别读取这两个文件中的数字，存在MyData1和MyData2中，每一组6个数，接下来的你自己处理一下吧。

fid1 = fopen('Data1.txt');
MyData1 = textsc ...

你好，我用你的方法试过了，确实能读取文件中的数据，但是我用下列语句提取列向量时遇到错误：

x1=MyData1(:,1);
x2=MyData1(:,2);
y1=MyData1(:,3);
y2=MyData1(:,4);
z=MyData1(:,5);
x=(x1+x2)/2;y=(y1+y2)/2;

到这里，提示加法不能用于cell的运算，也就是提取出来的列还是矩阵格式不是向量格式？不知道是不是这个意思，有什么办法可以解决呢。。。
赞一下 回复此楼
4楼2015-04-20 16:28:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
引用回帖:
4楼: Originally posted by godchenkun at 2015-04-20 16:28:11
你好，我用你的方法试过了，确实能读取文件中的数据，但是我用下列语句提取列向量时遇到错误：

x1=MyData1(:,1);
x2=MyData1(:,2);
y1=MyData1(:,3);
y2=MyData1(:,4);
z=MyData1(:,5);
x=(x1+x2)/2;y=(y1 ...

将MyData1后面的圆括号改为花括号就可以了。圆括号提出来的是cell格式，花括号是提取cell中的内容，具体到这里就是向量或矩阵，就可以直接运算了。
赞一下 回复此楼
知其然，知其所以然。
5楼2015-04-20 18:30:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
引用回帖:
4楼: Originally posted by godchenkun at 2015-04-20 16:28:11
你好，我用你的方法试过了，确实能读取文件中的数据，但是我用下列语句提取列向量时遇到错误：

x1=MyData1(:,1);
x2=MyData1(:,2);
y1=MyData1(:,3);
y2=MyData1(:,4);
z=MyData1(:,5);
x=(x1+x2)/2;y=(y1 ...

fid1 = fopen('Data1.txt');
MyData1 = textscan(fid1,'%*s %*s %f %*s %f \r\n %*s %*s %f %*s %f \r\n %f %f \r\n');
fclose(fid1);
x1 = MyData1{:,1};
x2 = MyData1{:,2};
y1 = MyData1{:,3};
y2 = MyData1{:,4};
z = MyData1{:,5};
x = (x1+x2)/2;
y = (y1+y2)/2;

但是这样得到的数据不能直接用mesh画曲面图，因为mesh需要网格化数据，配合meshgrid使用。所谓网格化数据就是得到x-y平面内某个区域内网格状点及对应的z值，假定该区域内点为m*n个，即m*n个数据对(x,y)，此时再加一个m*n矩阵z，对应每个点(x,y)的z取值，此时用mesh(x,y,z)才能画图。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

godchenkun
知其然，知其所以然。
6楼2015-04-20 18:44:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

godchenkun

银虫 (正式写手)
送红花一朵
引用回帖:
6楼: Originally posted by mygt_hit at 2015-04-20 18:44:21
fid1 = fopen('Data1.txt');
MyData1 = textscan(fid1,'%*s %*s %f %*s %f \r\n %*s %*s %f %*s %f \r\n %f %f \r\n');
fclose(fid1);
x1 = MyData1{:,1};
x2 = MyData1{:,2};
y1 = MyData1{:,3};
y2 = MyDa ...

太感谢了，花括号我摸索出来了，就是mesh搞不清楚啊，现在懂了，赞！
赞一下 回复此楼
7楼2015-04-21 08:05:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

godchenkun

银虫 (正式写手)
引用回帖:
6楼: Originally posted by mygt_hit at 2015-04-20 18:44:21
fid1 = fopen('Data1.txt');
MyData1 = textscan(fid1,'%*s %*s %f %*s %f \r\n %*s %*s %f %*s %f \r\n %f %f \r\n');
fclose(fid1);
x1 = MyData1{:,1};
x2 = MyData1{:,2};
y1 = MyData1{:,3};
y2 = MyDa ...

mesh也解决了，因为xy在网格上都是等差均匀分布的，所以我干脆自己生成xy坐标，分别为m个和n个，然后将向量z reshape成m*n矩阵，就能画了^_^
赞一下 回复此楼
8楼2015-04-21 08:14:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
引用回帖:
8楼: Originally posted by godchenkun at 2015-04-21 08:14:08
mesh也解决了，因为xy在网格上都是等差均匀分布的，所以我干脆自己生成xy坐标，分别为m个和n个，然后将向量z reshape成m*n矩阵，就能画了^_^...

恩，挺好。只要保证（x,y）->z的映射正确就可以。
赞一下 回复此楼
知其然，知其所以然。
9楼2015-04-21 08:35:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 godchenkun 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分人工智能专业请求调剂愿意跨考！ +3 lemonzzn 2026-03-17 3/150 2026-03-19 23:58 by 23Postgrad
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 12/600 2026-03-19 16:17 by 脱颖而出
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 一志愿北京化工大学0703化学318分，有科研经历，求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 317求调剂 +3 申子申申 2026-03-19 6/300 2026-03-19 14:16 by 申子申申
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
信息提示
关闭
请填处理意见
关闭