| 查看: 574 | 回复: 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客户端 |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
微分方程的数值解求法
已经有3人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.5小时
- 虫号: 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
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
5楼2016-09-21 15:37:40