24小时热门版块排行榜    

查看: 1686  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 4/200 2026-02-07 22:24 by ioveiuleh5
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 4/200 2026-02-07 22:11 by ioveiuleh5
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 22:06 by ioveiuleh5
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 22:04 by ioveiuleh5
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 21:46 by ioveiuleh5
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:44 by ioveiuleh5
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:31 by ioveiuleh5
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:26 by ioveiuleh5
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:11 by ioveiuleh5
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 3/150 2026-02-07 18:31 by 6ntj4f1rhx
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[硕博家园] 博士延得我,科研能力直往上蹿 +7 偏振片 2026-02-02 7/350 2026-02-04 17:36 by 陈氏帝国
[基金申请] 面上基金申报没有其他的参与者成吗 +5 默默的小虫子 2026-01-31 5/250 2026-02-04 11:26 by 学员fOzRO9
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见