| 查看: 2687 | 回复: 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 |
回复此楼» 猜你喜欢
0703化学
已经有26人回复
-
274求调剂求调剂
已经有11人回复
-
327求调剂
已经有10人回复
-
调剂
已经有19人回复
-
275 求调剂
已经有7人回复
-
求材料调剂，一志愿郑州大学289分
已经有21人回复
-
电子信息346
已经有4人回复
-
求助071001调剂！！！
已经有7人回复
-
一志愿深大085601材料工程专业（专硕）300分可以调剂去哪
已经有12人回复
-
一志愿华东理工085601材料工程303分求调剂
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- 偏微分方程Matlab求解
已经有7人回复
- Matlab求解微分方程(组)及偏微分方程(组)
已经有224人回复
- Matlab数值求解二阶常微分方程
已经有9人回复
- Matlab数值求解非线性常微分方程
已经有5人回复
- matlab求二阶偏导数
已经有3人回复
- 求助，matlab求解二元二阶的常微分方程组
已经有9人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- matlab解偏微分方程求助
已经有12人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 二阶非线性微分方程求解
已经有14人回复
- 求解两个简单非线性常微分方程~~
已经有4人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点
已经有5人回复
- matlab二阶微分方程求解程序
已经有5人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
4楼2014-08-30 15:50:18
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
5楼2014-08-30 16:08:02
Coco_hdu
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.3
- 帖子: 14
- 在线: 4.7小时
- 虫号: 2620473
- 注册: 2013-08-28
- 性别: GG
- 专业: 制造系统与自动化
从哪里看，求指导 6楼2014-09-03 16:19:39
Coco_hdu
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.3
- 帖子: 14
- 在线: 4.7小时
- 虫号: 2620473
- 注册: 2013-08-28
- 性别: GG
- 专业: 制造系统与自动化
7楼2014-09-03 16:20:59
Coco_hdu
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.3
- 帖子: 14
- 在线: 4.7小时
- 虫号: 2620473
- 注册: 2013-08-28
- 性别: GG
- 专业: 制造系统与自动化
8楼2014-09-03 16:22:54
Coco_hdu
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.3
- 帖子: 14
- 在线: 4.7小时
- 虫号: 2620473
- 注册: 2013-08-28
- 性别: GG
- 专业: 制造系统与自动化
9楼2014-09-03 16:25:45
Coco_hdu
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.3
- 帖子: 14
- 在线: 4.7小时
- 虫号: 2620473
- 注册: 2013-08-28
- 性别: GG
- 专业: 制造系统与自动化
10楼2014-09-03 16:26:46
