wahrrysu

木蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 4221.3
散金: 499
紅花: 2
帖子: 346
在線: 206.8小時
蟲號: 889613
注冊: 2009-10-31
專業(yè): 金屬結(jié)構(gòu)材料

[求助] 請月只藍進來看一下

你好，上次請您幫忙計算的方程（http://www.gaoyang168.com/bbs/viewthread.php?tid=6227965），由于偏差較大，想在原來方程的基礎上增加一個變量w,其中w^2=lamda*delta*R*T/(2*Nv*Eb^2*DI)

原方程改為 v=M0*delta*Vm*(1/Rr-3fp/rp)(1+w^2*v^2)/(1+lamda*M0*Cs*w^2*v^2)
v=dr/dt
非常感謝！

回復此樓

» 本帖@通知

@月只藍

» 猜你喜歡

哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有7人回復
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有109人回復
求調(diào)劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調(diào)劑信息已經(jīng)有31人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有68人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經(jīng)有1人回復
調(diào)劑已經(jīng)有0人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

本人想考14年中財大的會計學研究生，急需知道的學長學姐進來看一下已經(jīng)有6人回復
目標函數(shù)出現(xiàn)錯誤已經(jīng)有3人回復
CASTEP SCF tolerance 和 convergence tolerance 有何區(qū)別已經(jīng)有4人回復
會origin的進來看一下啊已經(jīng)有18人回復
香港科技大學的待遇，進來看一下已經(jīng)有8人回復
關于二面角掃描的結(jié)果，大神進來幫忙看一下已經(jīng)有9人回復
MATLAB求解代數(shù)方程組問題，請幫忙已經(jīng)有9人回復
做水熱的進來看一下了大家目前都在做啥材料已經(jīng)有6人回復
用過氯丁烷和氯代正己烷的童鞋請進來看一下已經(jīng)有3人回復
用過乙腈作溶劑的童鞋請進來看一下撒~ 已經(jīng)有23人回復
急�。�！有組織過國際會議經(jīng)驗的蟲蟲進來看一下~~~~~~ 已經(jīng)有17人回復
求助循環(huán)伏安掃描，請各位蟲蟲進來看一下已經(jīng)有9人回復
學中藥的哥哥姐姐們進來看一下已經(jīng)有12人回復
【討論】無鉛壓電陶瓷進來看一下... 已經(jīng)有9人回復
【求助】投過《生態(tài)毒理學報》的蟲子進來看一下已經(jīng)有11人回復
【經(jīng)驗分享】山東大學的同學們進來看一下已經(jīng)有35人回復
【求助】涂料結(jié)合力問題，拜托進來看一下已經(jīng)有12人回復

1樓 2013-08-15 20:27:50

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

月只藍

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

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

ben_ladeng: 專家考核 2013-08-16 18:55:20

現(xiàn)在的方程是：
dr/dt=M0*delta*Vm*(1/Rr-3fp/rp)(1+w^2*v^2)/(1+lamda*M0*Cs*w^2*（dr/dt）^2)，可以看到，待求解的dr/dt的等號右端包含了一個dr/dt的平方，這樣dr/dt是隱式的，無法用ode函數(shù)求解。

贊一下(1人)

回復此樓

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

2樓2013-08-16 12:58:42

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

wahrrysu

木蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 4221.3
散金: 499
紅花: 2
帖子: 346
在線: 206.8小時
蟲號: 889613
注冊: 2009-10-31
專業(yè): 金屬結(jié)構(gòu)材料

引用回帖:

2樓: Originally posted by 月只藍 at 2013-08-16 12:58:42
現(xiàn)在的方程是：
dr/dt=M0*delta*Vm*(1/Rr-3fp/rp)(1+w^2*v^2)/(1+lamda*M0*Cs*w^2*（dr/dt）^2)，可以看到，待求解的dr/dt的等號右端包含了一個dr/dt的平方，這樣dr/dt是隱式的，無法用ode函數(shù)求解。

請問這個方程有其它方法可以解么？

贊一下

回復此樓

3樓2013-08-16 14:04:47

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

月只藍

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

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

ben_ladeng: 專家考核 2013-08-16 18:55:26

引用回帖:

3樓: Originally posted by wahrrysu at 2013-08-16 14:04:47
請問這個方程有其它方法可以解么？...

比如直接差分法，用差商替代導數(shù)，把原方程在要求解的區(qū)間上離散成代數(shù)方程組，解非線性代數(shù)方程組。

贊一下

回復此樓

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

4樓2013-08-16 18:41:00

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

wahrrysu

木蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 4221.3
散金: 499
紅花: 2
帖子: 346
在線: 206.8小時
蟲號: 889613
注冊: 2009-10-31
專業(yè): 金屬結(jié)構(gòu)材料

引用回帖:

4樓: Originally posted by 月只藍 at 2013-08-16 18:41:00
比如直接差分法，用差商替代導數(shù)，把原方程在要求解的區(qū)間上離散成代數(shù)方程組，解非線性代數(shù)方程組。...

感謝您耐心的回復，聽上去很復雜啊，您有時間的話能幫我寫一個matlab的代碼么？像上次那樣的就行，現(xiàn)在我把金幣提高了，因為除了這個，我也沒有可以報答您的了，非常感謝��！

贊一下

回復此樓

5樓2013-08-16 19:31:31

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

成cheng

6樓2013-08-16 20:01:02

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

月只藍

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

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應助指數(shù) +1
csgt0: 金幣+2, 謝謝 2013-08-17 09:42:30
ben_ladeng: 專家考核 2013-08-17 15:10:15
wahrrysu: 金幣+150, ★★★★★最佳答案 2013-08-18 19:02:21

引用回帖:

5樓: Originally posted by wahrrysu at 2013-08-16 19:31:31
感謝您耐心的回復，聽上去很復雜啊，您有時間的話能幫我寫一個matlab的代碼么？像上次那樣的就行，現(xiàn)在我把金幣提高了，因為除了這個，我也沒有可以報答您的了，非常感謝！！...

不必客氣。我之前很少涉及解隱式的常微分方程，其實MATLAB有求解隱式常微分方程的函數(shù)，程序如下，結(jié)果見附圖：

CODE:

function solve_eq

clear all;clc

format long

global T

T=1200;

Rr0=10e-9;

dRr0=0;

tspan=500:3500;

[t Rr]=ode15i(@dRrdt,tspan,Rr0,dRr0);

figure,plot(t,Rr,'b'),xlabel('t'),ylabel('Rr'),axis([450 3650 min(Rr)-(max(Rr)-min(Rr))/10 max(Rr)+(max(Rr)-min(Rr))/10]);

[t Rr];

function f=dRrdt(t,Rr,dRr)

global T

R=8.314;

phi=0.7;

DI=2.14*exp(-158840/R/T);

DNb=0.83*exp(-266500/R/T);

sigma=0.8;

delta=1e-9;

Vm=6.77e-6;

b=2.5e-10;

Eb=19.3e3;

Nv=1e15;

rp=1e-8;

M0=phi*DI*Vm*delta/(b^2*R*T);

lamda=2*delta*Nv*R^2*T^2*(sinh(Eb/R/T)-Eb/R/T)/(Eb*DI);

x0=[0.01 0.3];

x=fsolve(@nonfun,x0,optimset('Display','Off'));

Cs=x(1);

fp=(0.021-Cs)*0.0114;

%-----------------------------------------%

w2=lamda*delta*R*T/(2*Nv*Eb^2*DI);

f=dRr-M0*delta*Vm*(1/Rr-3*fp/rp)*(1+w2*dRr^2)/(1+lamda*M0*Cs*w2*dRr^2);

function ff=nonfun(x)

global T

Cs=x(1);

C=x(2);

ff(1)=Cs*C-10^(5.17-9810/T+(540.6/T-0.1651)*0.33);

ff(2)=(0.021-Cs)/(0.33-C)-92.9064/12.011;

附圖1.jpg

贊一下

回復此樓

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

7樓2013-08-16 21:35:29

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

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

[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 一志愿中國石油大學（華東）本科齊魯工業(yè)大學 +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 278求調(diào)劑 +6	煙火先于春 2026-03-17	6/300	2026-03-21 01:57 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 材料學求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 081700化工學碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 復試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 0703化學調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 334求調(diào)劑 +3	志存高遠意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 考研化學學碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客