24小时热门版块排行榜    

查看: 1669  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 中青基了要发朋友圈吗? +6 349506619 2026-08-28 6/300 2026-08-31 08:43 by kingkocxr
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 08:28 by xiangy672
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见