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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：matlab解超越方程组得到数值解，要对此数值解的求定积分，得到一个确定的数值。
81/1返回列表
查看: 1495  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

在水一方801

新虫 (初入文坛)

[求助] 求助：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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-05-29 16:11:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

在水一方801

新虫 (初入文坛)
自己顶一下，我先问如何对数值解拟合的函数求积分
赞一下 回复此楼
2楼2015-06-08 10:36:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
在水一方801: 金币+5, ★★★★★最佳答案, 非常感谢，先给分，后试验。 2015-06-09 13:48:19
引用回帖:
2楼: Originally posted by 在水一方801 at 2015-06-08 10:36:40
自己顶一下，我先问如何对数值解拟合的函数求积分

拟合出来的方程，int 函数求积分即可。
只有散点，trapz函数求积分即可。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

在水一方801
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2015-06-08 16:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

在水一方801

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by 月只蓝 at 2015-06-08 16:41:54
拟合出来的方程，int 函数求积分即可。
只有散点，trapz函数求积分即可。...

我还是没有试出来，大神能帮忙写写程序吗。我附上我的程序。
对程序中y(1)的图中的曲线求积分，如何求，希望能写出详细的程序！
麻烦了！
赞一下 回复此楼
4楼2015-06-24 10:23:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

在水一方801

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by 月只蓝 at 2015-06-08 16:41:54
拟合出来的方程，int 函数求积分即可。
只有散点，trapz函数求积分即可。...

附上程序
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : Qswitch.m
    • 2015-06-24 10:24:30, 385 bytes
  • 附件 2 : rate_eq.m
    • 2015-06-24 10:24:30, 582 bytes
5楼2015-06-24 10:24:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

在水一方801

新虫 (初入文坛)
送红花一朵
引用回帖:
3楼: Originally posted by 月只蓝 at 2015-06-08 16:41:54
拟合出来的方程，int 函数求积分即可。
只有散点，trapz函数求积分即可。...

麻烦您帮忙看一下程序，帮忙对y（1）进行积分，麻烦了，我尝试了很久没能成功。
赞一下 回复此楼
6楼2015-06-25 11:27:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
6楼: Originally posted by 在水一方801 at 2015-06-25 11:27:01
麻烦您帮忙看一下程序，帮忙对y（1）进行积分，麻烦了，我尝试了很久没能成功。...

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
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
7楼2015-06-25 13:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

在水一方801

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by 月只蓝 at 2015-06-25 13:42:07
clc
clear
close all

T0 = 0.938;
R = 0.9;
Rp = 4.7463e27;

y0 = ;
tspan=;
tic
= ode23('rate_eq',tspan,y0,[],Rp,T0,R);
toc

figure
subplot(3,1,1);
plot(t,y(:,1));
xlabel('t(s)');
y ...

原来要这样，get it，非常感谢！
赞一下 回复此楼
8楼2015-06-25 15:00:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 在水一方801 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中科技大学，080502，354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 354求调剂 +5 Tyoumou 2026-03-18 8/400 2026-03-21 00:35 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by ！本暗一次！
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林，给自己招师弟中～ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
信息提示
关闭
请填处理意见
关闭