24小时热门版块排行榜    

查看: 1728  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 8/400 2026-02-28 23:36 by xyx2012xyx
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 290求调剂 +5 材料专硕调剂; 2026-02-28 6/300 2026-02-28 21:40 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 材料类求调剂 +6 wana_kiko 2026-02-28 6/300 2026-02-28 21:20 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂;
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见