24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1531  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

在水一方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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 8 个回答

在水一方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朵)

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 280求调剂 +12 李rien 2026-04-04 12/600 2026-04-10 21:19 by zhouxiaoyu
[考研] 计算机类求调剂,22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 346,工科0854求调剂,专硕 +5 moser233 2026-04-10 6/300 2026-04-10 18:04 by 周邵阳
[考研] 314求调剂 +23 wakeluofu 2026-04-09 24/1200 2026-04-10 15:31 by MOF_Catal
[考研] 调剂 +19 不逢春 2026-04-05 20/1000 2026-04-10 10:15 by may_新宇
[考研] 277求调剂 +19 倪建设 2026-04-06 19/950 2026-04-10 09:24 by guosr9609
[考研] 生物与医药调剂 +5 十七sa 2026-04-05 5/250 2026-04-10 08:14 by kangsm
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 考研调剂 +13 冰冰,,, 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 材料专硕322 +14 哈哈哈吼吼吼哈 2026-04-05 14/700 2026-04-09 13:25 by 5268321
[考研] 086000调剂 +4 十七sa 2026-04-07 4/200 2026-04-09 09:27 by 壹往無前
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 材料科学与工程320求调剂,080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 材料专硕(0856) 339分求调剂 +16 哈哈哈鹅哈哈哈 2026-04-05 16/800 2026-04-08 16:02 by luoyongfeng
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 298求调剂 +4 残荷新柳 2026-04-07 4/200 2026-04-07 23:02 by lbsjt
[考研] 297分083200求助 +9 aekx 2026-04-05 9/450 2026-04-06 20:57 by flysky1234
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
信息提示
请填处理意见