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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 请教一个Mathematica积分的问题（重金求解）
71/1返回列表
查看: 1612  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭