| 查看: 1483 | 回复: 7 | ||||
[求助]
求助：matlab解超越方程组得到数值解，要对此数值解的求定积分，得到一个确定的数值。 已有1人参与
|
|
求助：matlab解3个超越方程组得到数值解，解完画出图后，要对其中一个函数图形求定积分，此时该用什么函数；如果没有合适的函数，利用定积分的定义写成for循环，这种方法是否靠谱，请大神赐教！ 例如写成下面的循环求定积分。 求y(1)所对应的定积分，t为自变量，M为积分值。 for i=0:1e-14:7e-8; t=i; M = M+y(1)*1e-14; end |
回复此楼» 猜你喜欢
哈尔滨理工大学物理系招收物理学考研调剂
已经有6人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
物理学I论文润色/翻译怎么收费?
已经有202人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有30人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有65人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解决matlab解非线性方程组出现的问题
已经有6人回复
- 求助matlab高手指点：大型非线性方程组（牛顿拉佛森法）
已经有7人回复
- 大神！这个求定积分不等式的证明怎么做呢？
已经有3人回复
- 如何用最小二乘法及MATLAB解这个方程组？
已经有6人回复
- MATLAB解非线性方程组
已经有9人回复
- 用含有定积分的隐函数拟合数据，matlab或者origin
已经有18人回复
- 求高手解方程组（应该不难）
已经有22人回复
- 1stopt/MATLAB 解非线性方程组问题
已经有4人回复
- 求大神解决matlab算积分方程组！急
已经有6人回复
- 利用matlab求一阶变系数微分方程组的数值解。
已经有7人回复
- matlab解方程组求助攻
已经有4人回复
- matlab非线性方程组该如何编程。求助
已经有10人回复
- 用matlab解非线性方程组
已经有4人回复
- 求助用matlab或者（1st0pt)编程解方程组
已经有3人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 这个方程组怎么解啊
已经有14人回复
- Matlab数值积分结果是虚数
已经有6人回复
- 有偿求助：matlab怎么解联立方程组？
已经有11人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
2楼2015-06-08 10:36:40
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
在水一方801
3楼2015-06-08 16:41:54
4楼2015-06-24 10:23:42
5楼2015-06-24 10:24:52
6楼2015-06-25 11:27:01
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
|
clc clear close all T0 = 0.938; R = 0.9; Rp = 4.7463e27; y0 = [1.0e7;5.3676e24;1.4e21]; tspan=[0 6e-7]; tic [t,y] = ode23('rate_eq',tspan,y0,[],Rp,T0,R); toc figure subplot(3,1,1); plot(t,y(:,1)); xlabel('t(s)'); ylabel('y1'); y1=y(:,1); jifen1=trapz(t,y1) subplot(3,1,2); plot(t,y(:,2)); xlabel('t(s)'); ylabel('y2'); subplot(3,1,3); plot(t,y(:,3)); xlabel('t(s)'); ylabel('y3'); 加了两行： y1=y(:,1); jifen1=trapz(t,y1) 即得y1对t的积分： jifen1 = 4.2736e+003 |
7楼2015-06-25 13:42:07
8楼2015-06-25 15:00:31
