| 查看: 2447 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助，我现在要求离散点的定积分，请问可以通过1 stopt或者其他软件实现吗？ 已有1人参与
|
|||
数据已经附上，第一列就是r，第二列是σz，求这个定积分的结果。求方法，因为有很多组数据。谢谢大神~~~ LZ$$5[0J[Y0EW}PIQ[00(GC.jpg |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有71人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有52人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
lijf517
木虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 3391.9
- 红花: 14
- 帖子: 293
- 在线: 390.9小时
- 虫号: 2249203
- 注册: 2013-01-18
- 性别: GG
- 专业: 力学
|
嗯，matlab的语法以及调用的函数名称与我上面贴的python代码略有差异。为尊重他人版权，我现在主要在linux下做模拟，使用的都是开源免费的软件，所以不再用matlab了。你其实也可以用matlab中的函数把文件中的数据读进数组，不需要一个个输进代码，容易出错。datFmt = '%f %f'; sizeDat = [2 Inf]; fileID = fopen('data.txt', 'r'); data = fscanf(fileID, dataFmt, sizeDat); data = data'; 这样就把文件data.txt中的数据读进数组data了 发自小木虫Android客户端 |
10楼2016-06-06 15:32:10
lijf517
木虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 3391.9
- 红花: 14
- 帖子: 293
- 在线: 390.9小时
- 虫号: 2249203
- 注册: 2013-01-18
- 性别: GG
- 专业: 力学
|
自己写程序吧。最简单的，可以用梯形积分法，就是计算一系列梯形面积加起来。搞计算模拟的，没这点最基本的数值分析技能，在科研路上会受到很多牵制。随便找本数值分析的书都有这方面内容 发自小木虫Android客户端 |
3楼2016-06-03 14:49:16
我不是学这个的，只是需要有些数据，需要量化表达出来，文献里使用积分表达的。
4楼2016-06-03 15:40:47
lijf517
木虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 3391.9
- 红花: 14
- 帖子: 293
- 在线: 390.9小时
- 虫号: 2249203
- 注册: 2013-01-18
- 性别: GG
- 专业: 力学
【答案】应助回帖
感谢参与，应助指数 +1
|
我看你还在拟合J-C本构的参数，说明我们专业一致啊。这点数值模拟的基本功还是得学学吧。下面是我用python写的代码，python是免费软件，matlab也可以做类似工作。 from scipy import integrate import numpy as np data = np.loadtxt("data.txt"  rr = data[:, 0] sigZ = data[:, 1] fr = rr*np.abs(sigZ) intgl = 2.0*np.pi*integrate.trapz(fr, rr) print "The integral is ", intgl 运行结果显示，最终积分值为129534.146566 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
拐点5楼2016-06-03 16:41:36
投票: