在水一方801

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線: 12.8小時(shí)
蟲號(hào): 2796876
注冊(cè): 2013-11-12
專業(yè): 光學(xué)

[求助] 求助：matlab解超越方程組得到數(shù)值解，要對(duì)此數(shù)值解的求定積分，得到一個(gè)確定的數(shù)值。已有1人參與

求助：matlab解3個(gè)超越方程組得到數(shù)值解，解完畫出圖后，要對(duì)其中一個(gè)函數(shù)圖形求定積分，此時(shí)該用什么函數(shù)；如果沒有合適的函數(shù)，利用定積分的定義寫成for循環(huán)，這種方法是否靠譜，請(qǐng)大神賜教！
例如寫成下面的循環(huán)求定積分。求y(1)所對(duì)應(yīng)的定積分，t為自變量，M為積分值。
for i=0:1e-14:7e-8;
t=i;
M = M+y(1)*1e-14;
end

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有225人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求解決matlab解非線性方程組出現(xiàn)的問題已經(jīng)有6人回復(fù)
求助matlab高手指點(diǎn)：大型非線性方程組（牛頓拉佛森法）已經(jīng)有7人回復(fù)
大神！這個(gè)求定積分不等式的證明怎么做呢？已經(jīng)有3人回復(fù)
如何用最小二乘法及MATLAB解這個(gè)方程組？已經(jīng)有6人回復(fù)
MATLAB解非線性方程組已經(jīng)有9人回復(fù)
用含有定積分的隱函數(shù)擬合數(shù)據(jù)，matlab或者origin 已經(jīng)有18人回復(fù)
求高手解方程組（應(yīng)該不難）已經(jīng)有22人回復(fù)
1stopt/MATLAB 解非線性方程組問題已經(jīng)有4人回復(fù)
求大神解決matlab算積分方程組！急已經(jīng)有6人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
matlab解方程組求助攻已經(jīng)有4人回復(fù)
matlab非線性方程組該如何編程。求助已經(jīng)有10人回復(fù)
用matlab解非線性方程組已經(jīng)有4人回復(fù)
求助用matlab或者（1st0pt)編程解方程組已經(jīng)有3人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
求助，用matlab解微分方程組，希望高手能給予指點(diǎn) 已經(jīng)有17人回復(fù)
這個(gè)方程組怎么解啊已經(jīng)有14人回復(fù)
Matlab數(shù)值積分結(jié)果是虛數(shù) 已經(jīng)有6人回復(fù)
有償求助：matlab怎么解聯(lián)立方程組？已經(jīng)有11人回復(fù)
matlab求解非線性方程組，錯(cuò)誤提示怎么解決已經(jīng)有5人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
matlab 積分計(jì)算算不出結(jié)果，求幫助�。�！已經(jīng)有10人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

1樓 2015-05-29 16:11:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

在水一方801

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線: 12.8小時(shí)
蟲號(hào): 2796876
注冊(cè): 2013-11-12
專業(yè): 光學(xué)

自己頂一下，我先問如何對(duì)數(shù)值解擬合的函數(shù)求積分

贊一下

回復(fù)此樓

2樓2015-06-08 10:36:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
在水一方801: 金幣+5, ★★★★★最佳答案, 非常感謝，先給分，后試驗(yàn)。 2015-06-09 13:48:19

引用回帖:

2樓: Originally posted by 在水一方801 at 2015-06-08 10:36:40
自己頂一下，我先問如何對(duì)數(shù)值解擬合的函數(shù)求積分

擬合出來的方程，int 函數(shù)求積分即可。
只有散點(diǎn)，trapz函數(shù)求積分即可。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

在水一方801

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

3樓2015-06-08 16:41:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

在水一方801

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線: 12.8小時(shí)
蟲號(hào): 2796876
注冊(cè): 2013-11-12
專業(yè): 光學(xué)

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2015-06-08 16:41:54
擬合出來的方程，int 函數(shù)求積分即可。
只有散點(diǎn)，trapz函數(shù)求積分即可。...

我還是沒有試出來，大神能幫忙寫寫程序嗎。我附上我的程序。
對(duì)程序中y(1)的圖中的曲線求積分，如何求，希望能寫出詳細(xì)的程序！
麻煩了！

贊一下

回復(fù)此樓

4樓2015-06-24 10:23:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

在水一方801

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線: 12.8小時(shí)
蟲號(hào): 2796876
注冊(cè): 2013-11-12
專業(yè): 光學(xué)

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2015-06-08 16:41:54
擬合出來的方程，int 函數(shù)求積分即可。
只有散點(diǎn)，trapz函數(shù)求積分即可。...

附上程序

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

在水一方801

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線: 12.8小時(shí)
蟲號(hào): 2796876
注冊(cè): 2013-11-12
專業(yè): 光學(xué)

送紅花一朵

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2015-06-08 16:41:54
擬合出來的方程，int 函數(shù)求積分即可。
只有散點(diǎn)，trapz函數(shù)求積分即可。...

麻煩您幫忙看一下程序，幫忙對(duì)y（1）進(jìn)行積分，麻煩了，我嘗試了很久沒能成功。

贊一下

回復(fù)此樓

6樓2015-06-25 11:27:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

引用回帖:

6樓: Originally posted by 在水一方801 at 2015-06-25 11:27:01
麻煩您幫忙看一下程序，幫忙對(duì)y（1）進(jìn)行積分，麻煩了，我嘗試了很久沒能成功。...

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對(duì)t的積分：
jifen1 =

4.2736e+003

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

7樓2015-06-25 13:42:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

在水一方801

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線: 12.8小時(shí)
蟲號(hào): 2796876
注冊(cè): 2013-11-12
專業(yè): 光學(xué)

引用回帖:

7樓: Originally posted by 月只藍(lán) 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，非常感謝！

贊一下

回復(fù)此樓

8樓2015-06-25 15:00:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	14/700	2026-03-28 15:48 by Micky11223
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5	niko- 2026-03-27	5/250	2026-03-28 14:25 by 唐沐兒
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 083000學(xué)碩274求調(diào)劑 +7	Li李魚 2026-03-26	7/350	2026-03-28 08:01 by baoball
[考研] 0856材料化工調(diào)劑總分330 +10	zhubinhao 2026-03-27	10/500	2026-03-28 03:34 by fmesaito
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 266求調(diào)劑 +11	陽(yáng)陽(yáng)哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 316求調(diào)劑 +5	江辭666 2026-03-26	5/250	2026-03-27 08:08 by hypershenger
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 085602化學(xué)工程求調(diào)劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 332求調(diào)劑 +6	032500 2026-03-25	6/300	2026-03-25 22:45 by 418490947
[考研] 296求調(diào)劑 +4	汪！？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪�。�！
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518