24小时热门版块排行榜    

查看: 1732  |  回复: 6

amethyst356

木虫 (小有名气)

[求助] 请教一个Mathematica积分的问题(重金求解)

问题大致是这样的:现有100个txt文件,文件当中是二维的数组;这100个文件的数据是从1秒到100秒所测得。现在所需要做的是需要对1秒时候所测得数据进行积分;然后把1秒作为x1,所求的积分为y1;然后对对2秒时候所测得数据进行积分;然后把2秒作为x2,所求的积分为y2;依次类推,得到x3,y3.......x100,y100,然后画出x-y曲线。

谢谢大家!
回复此楼

» 猜你喜欢

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

弱水三千,我一瓢一瓢饮。。。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zl030091126

金虫 (正式写手)

【答案】应助回帖

★ ★
amethyst356(金币+10): 2011-11-09 23:14:30
余泽成(金币+2): 谢谢参与应助!欢迎常来程序语言版! 2011-11-10 21:41:52
楼主没有说清楚一个问题,对1秒积分的时候,是从0秒开始到1秒结束;但是对2秒积分时候,是从0秒还是从1秒开始再到2秒结束?积分的区间没有弄清楚。我的理解是楼主想从1到2秒积分,再3到4秒,依次类推。我下面的程序楼主可以做个参考(随机生成的y)。
x=1:1:100;%1到100秒
y=rand(100,1);%1到100对应的y值
m(:,1)=[0;y(1)];%t=0时,y我设的值是0,不确定楼主是多少
for i=2:1:100;
    m(:,i)=[y(i-1);y(i)];%构造了个新的m矩阵
    Y(i)=trapz(m(:,i))%用梯形求积分公式积分
end
plot(x,Y)%作图
A--B,虚线变成实线则是你的轨迹
2楼2011-11-09 04:13:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amethyst356

木虫 (小有名气)

引用回帖:
2楼: Originally posted by zl030091126 at 2011-11-08 13:13:24:
楼主没有说清楚一个问题,对1秒积分的时候,是从0秒开始到1秒结束;但是对2秒积分时候,是从0秒还是从1秒开始再到2秒结束?积分的区间没有弄清楚。我的理解是楼主想从1到2秒积分,再3到4秒,依次类推。我下面的程 ...

谢谢你的回复。

我把我的问题说具体点:
1秒的时候,测得1组数据,即一个txt文件,文件里面是一个二维数组,1列是波长,1列是波的强度,现在需要积分出整个波长区域的总强度,所以积分变量不是时间,时间只是作为一个计数点。

更具体来说,1秒时,有五个波长,对应5个波强,假若数据如下:({300,1},{400,2},{500,3},{600,4},{700,5})这样的一个txt文件,第一列是波长,第二列对应的是波强,现在积出1秒时总的波强,然后把1秒当作x,积出的总的波强为y。

有100个这样的txt文件,因为是从1秒测到100秒测了100次得到的。

希望用mathematica求解。

谢谢。
弱水三千,我一瓢一瓢饮。。。
3楼2011-11-09 10:04:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)

【答案】应助回帖

★ ★
余泽成(金币+2): 谢谢参与应助! 2011-11-10 21:42:07
我想,你的提问里积分不是主要问题,简单的数值积分即可,甚至直接将第二列数据求和都可以正比于积分值了。
你应该是想用mathematica操作这100个文件,不用手动去一个一个文件打开并计算,最后能给出图就ok。
是这样吧?
如果这样,那就简单了。如果文件名是有规律的,比如1.txt,2.txt,...那么做一个循环,循环变量 i 转换成文件名,从而提取相应文件中的数据作积分,得yi,将{i,yi}赋值到表(Table)中,然后listplot这个表。
4楼2011-11-09 10:31:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amethyst356

木虫 (小有名气)

引用回帖:
4楼: Originally posted by mshwangg at 2011-11-08 19:31:56:
我想,你的提问里积分不是主要问题,简单的数值积分即可,甚至直接将第二列数据求和都可以正比于积分值了。
你应该是想用mathematica操作这100个文件,不用手动去一个一个文件打开并计算,最后能给出图就ok。
是 ...

谢谢你的回复,很受启发,我给你发短消息了。
弱水三千,我一瓢一瓢饮。。。
5楼2011-11-09 13:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)

【答案】应助回帖

★ ★ ★
xzhdty(金币+2): 欢迎常来程序语言 2011-11-09 17:06:50
amethyst356(金币+80): 2011-11-09 23:14:07
余泽成(金币+1): 专家考核存档! 2011-11-10 21:42:27
引用回帖:
5楼: Originally posted by amethyst356 at 2011-11-09 13:06:12:
谢谢你的回复,很受启发,我给你发短消息了。


写了个简单的小程序,主要想说文件的顺序操作,所以这里积分用求和代替了,循环直接用建表的那个,数据是随机生成的如下图。希望对你有所帮助
6楼2011-11-09 14:47:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

markh

银虫 (初入文坛)


余泽成(金币+1): 谢谢参与应助! 2011-11-13 21:02:14
引用回帖:
6楼: Originally posted by mshwangg at 2011-11-09 14:47:42:

写了个简单的小程序,主要想说文件的顺序操作,所以这里积分用求和代替了,循环直接用建表的那个,数据是随机生成的如下图。希望对你有所帮助
[eimg]50/b6/5767 ...

提取出各个文件的“date”之后, 积分部分可以如下
 
yeah, 终于升级成为淫虫。
7楼2011-11-11 04:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 amethyst356 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
信息提示
请填处理意见