| 查看: 664 | 回复: 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客户端 |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有23人回复
-
求调剂院校信息
已经有4人回复
-
085600材料与化工306
已经有4人回复
-
286求调剂
已经有10人回复
-
328求调剂，英语六级551，有科研经历
已经有9人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有4人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有8人回复
-
资源与环境 调剂申请(333分)
已经有5人回复
-
280求调剂
已经有12人回复
-
269专硕求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程的数值解求法
已经有3人回复
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
