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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 请 月只蓝 进来看一下
71/1返回列表
查看: 1314  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wahrrysu

木虫 (正式写手)

[求助] 请 月只蓝 进来看一下

你好，上次请您帮忙计算的方程（http://muchong.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
非常感谢！
回复此楼

» 本帖@通知

@月只蓝

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-08-15 20:27:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
ben_ladeng: 专家考核 2013-08-16 18:55:20
现在的方程是：
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函数求解。
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2013-08-16 12:58:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wahrrysu

木虫 (正式写手)
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-08-16 12:58:42
现在的方程是：
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函数求解。

请问这个方程有其它方法可以解么？
赞一下 回复此楼
3楼2013-08-16 14:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
ben_ladeng: 专家考核 2013-08-16 18:55:26
引用回帖:
3楼: Originally posted by wahrrysu at 2013-08-16 14:04:47
请问这个方程有其它方法可以解么？...

比如直接差分法，用差商替代导数，把原方程在要求解的区间上离散成代数方程组，解非线性代数方程组。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2013-08-16 18:41:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wahrrysu

木虫 (正式写手)
引用回帖:
4楼: Originally posted by 月只蓝 at 2013-08-16 18:41:00
比如直接差分法，用差商替代导数，把原方程在要求解的区间上离散成代数方程组，解非线性代数方程组。...

感谢您耐心的回复，听上去很复杂啊，您有时间的话能帮我写一个matlab的代码么？像上次那样的就行，现在我把金币提高了，因为除了这个，我也没有可以报答您的了，非常感谢！！
赞一下 回复此楼
5楼2013-08-16 19:31:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

成cheng
6楼2013-08-16 20:01:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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的代码么？像上次那样的就行，现在我把金币提高了，因为除了这个，我也没有可以报答您的了，非常感谢！！...

不必客气。我之前很少涉及解隐式的常微分方程，其实MATLAB有求解隐式常微分方程的函数，程序如下，结果见附图：
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小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
7楼2013-08-16 21:35:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wahrrysu 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 23:45 by vZfe6xYu34yj
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:21 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:12 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:24 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
关闭
请填处理意见
关闭