| 查看: 2511 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
Coco_hdu新虫 (初入文坛)
|
[交流]
用Matlab求一个二阶微分方程，出错了，求指导（求助）已有4人参与
|
||
|
function dy=fun00050050(t,y)%方程定义 dy = zeros(3,1); %创建列向量 V=0.05e-6; %注射剂量 L=28e-3; %腔室长度 %Ap=V/L; %活塞面积 %D=sqrt(Ap/pi); %腔体内径 Ao=4.91e-8; %喷孔面积 Mp=0.12; %活塞质量 Sigama=1.31e-6; %运动粘度 %Re=y(2)*D/Sigama; %雷诺数 %landa=64/Re; %沿程摩擦系数 E=2.18e9; %集中模量 ro=1000; %注射液的密度 b=0.001; %密封圈接触宽度 Sita=0.2; %铝和密封圈之间的摩擦系数 Pc=3.5e6; %无压情况下安瓿内液体压强 I=6; %电流大小 dy(1)=((E+y(1))*y(2)-E*Ao/(V/L)*(2*y(1)/(64*Sigama*L/(4*V/L/pi)/y(2)+1)/ro)^(0.5))/(L-y(3)); %y(1)=p; dy(2)=(4*1000*y(3)-(0.01*I+0.06)*1000000*y(3)*y(3)+10*I)/Mp+(-2)*(V/L)*y(1)/Mp+(-1)*(Pc+y(1))*Sita*b*pi*sqrt(Ap/pi)/Mp; %y(2)=v dy(3)=y(2); %y(3)=xp 然后执行[t,y] = ode23s('fun00050050',[0 0.0137],[0 0 0]); 出来了： 错误使用 feval 未定义与 'double' 类型的输入参数相对应的函数 'fun00050050'。 出错 odearguments (line 87) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. 出错 ode23s (line 119) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 求指导哪里错了，我该怎么改。  1111111111.PNG |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有6人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
心脉受损
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
Coco_hdu
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.3
- 帖子: 14
- 在线: 4.7小时
- 虫号: 2620473
- 注册: 2013-08-28
- 性别: GG
- 专业: 制造系统与自动化
12楼2014-09-03 18:36:22
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 数学EPI: 2
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
2楼2014-08-25 18:09:06
小吴smallwu
木虫 (小有名气)
- 应助: 18 (小学生)
- 金币: 1627.6
- 散金: 200
- 红花: 2
- 帖子: 189
- 在线: 114.5小时
- 虫号: 3378771
- 注册: 2014-08-22
- 性别: GG
- 专业: 传热传质学
3楼2014-08-29 10:42:02
4楼2014-08-30 15:50:18