| 查看: 1757 | 回复: 6 | ||
amethyst356木虫 (小有名气)
|
[求助]
请教一个Mathematica积分的问题（重金求解）
|
|
问题大致是这样的：现有100个txt文件，文件当中是二维的数组；这100个文件的数据是从1秒到100秒所测得。现在所需要做的是需要对1秒时候所测得数据进行积分；然后把1秒作为x1，所求的积分为y1；然后对对2秒时候所测得数据进行积分；然后把2秒作为x2，所求的积分为y2；依次类推，得到x3,y3.......x100,y100，然后画出x-y曲线。 谢谢大家！ |
回复此楼» 猜你喜欢
本9 一志愿西工大085601 324求调剂
已经有4人回复
-
一志愿东北大学控制工程085406数二英二385，求调剂
已经有6人回复
-
272分材料子求调剂
已经有22人回复
-
0856专硕求调剂 希望是a区院校
已经有23人回复
-
346，工科0854求调剂，专硕
已经有5人回复
-
337求调剂
已经有3人回复
-
266求调剂
已经有29人回复
-
309求调剂
已经有13人回复
-
一志愿0703化学招61最终排名62化学求调剂
已经有28人回复
-
286求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个求弧长的问题！
已经有4人回复
- 一个很不错的mathematica教程
已经有507人回复
- Mathematica关于绘制隐函数方程组的问题，跪求，快疯了
已经有11人回复
- Mathematica问题求助
已经有7人回复
- 请教一个图论的题
已经有5人回复
- 请教接一个接收函的问题
已经有12人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- 请教一个积分用留数定理计算的问题
已经有15人回复
- 请教瞬时相位的求解？
已经有7人回复
- FFT算法的一个问题请教
已经有18人回复
- mathematica请教计算结果问题
已经有4人回复
- 请教一个while loop程序的问题
已经有5人回复
- mathematica里的运算
已经有4人回复
- 请教一个矩阵求逆问题
已经有4人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】请教一个关于Mathematica中计算小量的问题
已经有4人回复
- 【求助】想请教一个如何求解旋转矩阵的问题
已经有6人回复
- 【求助】【求助】mathematica中的一种数值积分问题
已经有5人回复
- 【求助】请教一个关于解方程的问题
已经有8人回复
- 【求助】请教COmsol求解问题
已经有10人回复
- 【求助】一个经典的数模问题求解！！！！急
已经有7人回复
- 【求助】紧急求助有关MATHEMATICA的问题
已经有9人回复
zl030091126
金虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 303.8
- 散金: 605
- 红花: 2
- 帖子: 524
- 在线: 141.5小时
- 虫号: 1444718
- 注册: 2011-10-16
- 性别: GG
- 专业: 化学反应工程
【答案】应助回帖
★ ★
amethyst356(金币+10): 2011-11-09 23:14:30
余泽成(金币+2): 谢谢参与应助！欢迎常来程序语言版！ 2011-11-10 21:41:52
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)%作图 |
2楼2011-11-09 04:13:24
amethyst356
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1986
- 散金: 3555
- 红花: 3
- 帖子: 275
- 在线: 830小时
- 虫号: 474004
- 注册: 2007-12-05
- 专业: 基础物理学
|
谢谢你的回复。 我把我的问题说具体点： 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
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
4楼2011-11-09 10:31:56
amethyst356
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1986
- 散金: 3555
- 红花: 3
- 帖子: 275
- 在线: 830小时
- 虫号: 474004
- 注册: 2007-12-05
- 专业: 基础物理学
5楼2011-11-09 13:06:12
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
6楼2011-11-09 14:47:42
markh
银虫 (初入文坛)
- 应助: 7 (幼儿园)
- 金币: 765.1
- 帖子: 35
- 在线: 138.6小时
- 虫号: 366637
- 注册: 2007-05-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
![](\"https://muchongimg.xmcimg.com/data/edu/cc/7b/366637_1320956237_346.png#opennewwindow\")
7楼2011-11-11 04:19:29
5