24小时热门版块排行榜

返回列表

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
am99lx: 金币+80, ★★★很有帮助 2016-01-27 14:30:11

function dTdt
clear all;clc
format long

[t T]=ode15s(@odes,[0 1e-12],[300 300])
figure(1)
plot(t,T(:,1),'b',t,T(:,2),'r'),axis([0 max(t) min(min(T))-10 300+10]);

function f=odes(t,y)

Te=y(1);
T1=y(2);
Ce=88*Te;
C1=2.7e16;

J=44;
R=0.6;
b=12.5e-9;
tp=0.1e-12;
G=2.7e16;

Q=-sqrt(4*log(2)/pi)*(1-R)*J/b/tp*exp( -4*log(2)*((t-tp)/tp)^2 );

f1=-G*(Te-T1)+Q;
f2=G*(Te-T1);

f=[f1/Ce;f2/C1];

附图1.png

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

11楼2016-01-27 13:54:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 am99lx 的主题更新

返回列表