| 查看: 1197 | 回复: 11 | ||
wodeai12345银虫 (小有名气)
|
[求助]
请大神帮我看看，这个微分方程组要怎么编程计算 已有1人参与
|
|
 B=2.84 ρl/ρg=890 Fr=5650000 初始条件为X=B 时 y=0, dy/dx=0 想求出X Y的曲线的坐标点（自己来画图） 自己编了一个 发现不对 求帮助 建立m文件jisuan1 function dy=eq2(x,y); dy=zeros(3,1); dy(1)=y(2); dy(2)=1/(5650000*2.84^2)*(1+y(2)^2)^(1/2)*x^2*dy(3); dy(3)=2.84/int(sqrt(1+y(2)^2),x,0,x)*[dy(3)+(1-dy(3))*890]^(1/2); 计算程序 x0=2.84,xf=100; [x,y]=ode45('jisuan1',[x0,xf],[0,0,1]); plot(x,y); |
回复此楼» 猜你喜欢
299求调剂
已经有6人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有5人回复
-
279求调剂
已经有3人回复
-
336求调剂
已经有3人回复
-
286求调剂
已经有3人回复
-
一志愿华中科技大学071000，求调剂
已经有3人回复
-
279分求调剂 一志愿211
已经有14人回复
-
332求调剂
已经有3人回复
-
0856材料专硕353求调剂
已经有3人回复
-
材料学学硕080502 337求调剂-一志愿华中科技大学
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 小弟是土木的，请大神们帮我看看这三张红外光谱图，万分感谢！
已经有19人回复
- 请大神帮我看看，这个正常吗？
已经有13人回复
- 请大神帮小弟翻译这个句子，我看了好长时间都不知道该咋翻译，谢谢
已经有7人回复
- 有两个小问题，请大神帮看看。
已经有3人回复
- 四区sci，一个月后修改意见，major revise，请大神帮我看看希望大不大啊
已经有32人回复
- 请大神帮我看看这个分子式中的斜线代表啥意思，有没有类似的文献也是这么写的
已经有6人回复
- 请大神帮我看看这个XPS的数据，新手求助
已经有0人回复
- 请大神帮帮忙，我实在是没有金币啊，可以给看看这个循环伏安曲线是怎么回事？
已经有1人回复
- 请大神帮我看看这个图片，能不能写一下它的立体构型
已经有0人回复
- 请大神帮我看一下这个是怎么翻译的，自己没看懂怎么做 ，谢谢
已经有2人回复
- 请大神帮我看看，环糊精能实现我的这个用途吗？
已经有0人回复
- 请大神们看下我这个反应，有可能得到不？
已经有8人回复
- 请大神门帮我看看这个错误，在线等待答案
已经有11人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
wodeai12345: 金币+88, ★有帮助 2016-07-14 09:22:42
感谢参与，应助指数 +1
wodeai12345: 金币+88, ★有帮助 2016-07-14 09:22:42
 本帖仅楼主可见 |
2楼2016-07-14 06:33:05
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
wodeai12345: 金币+212, ★★★★★最佳答案, 太感谢你了 我要给你送花。。。 2016-07-14 12:38:17
月只蓝: 金币+10, 感谢应助 2016-07-14 13:51:50
wodeai12345: 金币+212, ★★★★★最佳答案, 太感谢你了 我要给你送花。。。 2016-07-14 12:38:17
月只蓝: 金币+10, 感谢应助 2016-07-14 13:51:50
|
本帖仅楼主可见 |
6楼2016-07-14 12:28:04
wodeai12345
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 369.6
- 散金: 10
- 帖子: 82
- 在线: 225.9小时
- 虫号: 3310922
- 注册: 2014-07-07
- 专业: 工程热物理相关交叉领域
3楼2016-07-14 08:20:26
wodeai12345
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 369.6
- 散金: 10
- 帖子: 82
- 在线: 225.9小时
- 虫号: 3310922
- 注册: 2014-07-07
- 专业: 工程热物理相关交叉领域
4楼2016-07-14 08:22:55
wodeai12345
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 369.6
- 散金: 10
- 帖子: 82
- 在线: 225.9小时
- 虫号: 3310922
- 注册: 2014-07-07
- 专业: 工程热物理相关交叉领域
|
后面又编了一个 还是不对 简直要崩溃了 X=B时 C确实应该等于1 dy(2)=1/(Fr*B^2)*(1+dy(1)^2)^(1/2)*x^2*c; c=B/int((1+dy(1)^2)^(1/2),x,0,x)*[c+(1-c)*(pj/pg)]^(1/2); 把第一个式子转换成C= XXX的形式 代入第二个式子里 消去C dy(2)=1/(Fr*B^2)*(1+dy(1)^2)^(1/2)*x^2*B*[dy(2)*Fr*B^2/(1+dy(1)^2)^(1/2)/x^2+(1-dy(2)*Fr*B^2/(1+dy(1)^2)^(1/2)/x^2)*(pj/pg)]^(1/2)/int((1+dy(1)^2)^(1/2),x,0,x) 编程如下 function dy=rigid(x,y); dy=zeros(2,1); dy(1)=y(2); dy(2)=1/(5650000*2.84^2)*(1+dy(1)^2)^(1/2)*x^2*2.84*[dy(2)*5650000*2.84^2/(1+dy(1)^2)^(1/2)/x^2+(1-dy(2)*5650000*2.84^2/(1+dy(1)^2)^(1/2)/x^2)*890]^(1/2)/int((1+dy(1)^2)^(1/2),x,0,x) 主程序 [x,y]=ode45('rigid',[2.84 30],[0 0]); plot(x,y); 提示错误 ??? Undefined function or method 'int' for input arguments of type 'double'. Error in ==> rigid at 4 dy(2)=1/(5650000*2.84^2)*(1+dy(1)^2)^(1/2)*x^2*2.84*[dy(2)*5650000*2.84^2/(1+dy(1)^2)^(1/2)/x^2+(1-dy(2)*5650000*2.84^2/(1+dy(1)^2)^(1/2)/x^2)*890]^(1/2)/int((1+dy(1)^2)^(1/2),x,0,x) Error in ==> odearguments at 109 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in ==> jisuanzhuchengxu at 1 [x,y]=ode45('rigid',[2.84 30],[0 0]); |
5楼2016-07-14 08:35:27
wodeai12345
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 369.6
- 散金: 10
- 帖子: 82
- 在线: 225.9小时
- 虫号: 3310922
- 注册: 2014-07-07
- 专业: 工程热物理相关交叉领域
7楼2016-07-14 12:45:05
8楼2016-07-14 21:22:10
wodeai12345
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 369.6
- 散金: 10
- 帖子: 82
- 在线: 225.9小时
- 虫号: 3310922
- 注册: 2014-07-07
- 专业: 工程热物理相关交叉领域
9楼2016-07-14 22:52:43
wodeai12345
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 369.6
- 散金: 10
- 帖子: 82
- 在线: 225.9小时
- 虫号: 3310922
- 注册: 2014-07-07
- 专业: 工程热物理相关交叉领域
10楼2016-07-14 22:57:40
