| 查看: 663 | 回复: 4 | |||
[交流]
微分方程组就解，求大神帮忙 已有2人参与
|
|
function dy=rigid(t,y) dy=zeros(4,1); ey=11.65e-25;Ny=2.01e+27;k1=5e-22; k2=5e-22;ry=1000;c=3e+8;n1=1.5;ee=8e-25; re=125;C=6e-25;n2=1.75;l1=0.0031;l2=0.0019;L=0.005;egsa=3.5e-25;eesa=1e-26;R=0.85; Ns0=1.29e24;l=(n1-1).*l1+(n2-1).*l2+L; tr=2*l/(c*(log(1/R)+0.04)); ta=3.4e-7; Pin=200;wp=0.1e-3;h=6.62e-34; vp=c/(1535e-9); Wp=2.*Pin.*exp(-280.*0.0031)./(pi.*h.*vp.*(wp).^2); dy(1)=ey.*Ny.*Wp-k1.*y(1).*(9.88e25-y(2))-k2.*y(1).*y(2)-y(1).*ry; dy(2)=k1.*y(1).*(9.88e25-y(2))-(2*y(2)-9.88e25).*(c./n1).*ee.*y(3)-y(2).*re-C.*(y(2)).^2; dy(3)=y(3)./tr.* (2.*ee.*(2.*y(2)-9.88e+25).*l1-2.*egsa.*y(4).*l2-2.*eesa.*(Ns0-y(4)).*l2-log(1./R)-0.04); dy(4)=-egsa.*c./n2.*y(3).*y(4)+(Ns0-y(4))./ta; end [T,Y] = ode45(@rigid,[0 5e-08],[0 0 0 1.29e24] ,5e-09) 结果中y（3）y(4).值一直不变，和理论不一致。图片是结果。上边是源程序。求大神看看有什么问题，谢谢了~  发自小木虫Android客户端 |
回复此楼» 猜你喜欢
354求调剂
已经有9人回复
-
化学工程321分求调剂
已经有22人回复
-
考研调剂
已经有3人回复
-
326求调剂
已经有8人回复
-
333求调剂
已经有7人回复
-
一志愿东华大学控制学硕320求调剂
已经有3人回复
-
考研化学学硕调剂，一志愿985
已经有7人回复
-
0703化学调剂
已经有15人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有14人回复
-
326求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程的数值解求法
已经有3人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2016-09-20 09:40:05
3楼2016-09-20 14:11:45
4楼2016-09-20 20:22:09
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
5楼2016-09-21 15:37:40
