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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 请教一个Mathematica积分的问题（重金求解）
71/1返回列表
查看: 2073  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

amethyst356

木虫 (小有名气)

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

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

谢谢大家！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
弱水三千，我一瓢一瓢饮。。。
1楼 2011-11-09 01:53:06
已阅   回复此楼   关注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)%作图
赞一下(1人) 回复此楼
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这个表。
赞一下(1人) 回复此楼
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:
谢谢你的回复，很受启发，我给你发短消息了。


写了个简单的小程序，主要想说文件的顺序操作，所以这里积分用求和代替了，循环直接用建表的那个，数据是随机生成的如下图。希望对你有所帮助
赞一下(2人) 回复此楼
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”之后， 积分部分可以如下
 
赞一下(1人) 回复此楼
yeah， 终于升级成为淫虫。
7楼2011-11-11 04:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 amethyst356 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:33 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-30 01:58 by ZPa0EcMwuECS
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:03 by ZPa0EcMwuECS
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-29 23:14 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
关闭
请填处理意见
关闭