| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1140 人关注过本帖
标题:MATLAB解微分方程
只看楼主 加入收藏
E个自己
Rank: 1
等 级:新手上路
威 望:1
帖 子:3
专家分:0
注 册:2009-10-13
结帖率:0
收藏
已结贴  问题点数:20 回复次数:1 
MATLAB解微分方程
Uc=0.9292;
 Ac=3.8513e-004;
 Tc=0.8237;
 Cc=0.9843;
 y0=0.9486;
 f=inline('Tc*(-y+Cc*((Uc-Ac/(y.^2)).^0.5)','Tc','Cc','Uc','Ac','x','y');%微分方程的右边项
 dx=0.004; %x方向步长
xleft=0.029000;   %区域的左边界
xright=1.020991;  %区域的右边界
xx=xleft:dx:xright; %一系列离散的点
n=length(xx); %点的个数
%%(2)龙格--库塔法
RK=y0;
for i=2:n
 k1=f(Tc,Cc,Uc,Ac,xx(i-1),RK(i-1));
 k2=f(Tc,Cc,Uc,Ac,xx(i-1)+dx/2,RK(i-1)+k1*dx/2);
 k3=f(Tc,Cc,Uc,Ac,xx(i-1)+dx/2,RK(i-1)+k2*dx/2);
 k4=f(Tc,Cc,Uc,Ac,xx(i-1)+dx,RK(i-1)+k3*dx);
 RK(i)=RK(i-1)+dx*(k1+2*k2+2*k3+k4)/6;
end
运行后老是提示错误为
??? Error using ==> inlineeval
Error in inline expression ==> Tc*(-y+Cc*((Uc-Ac/(y.^2)).^0.5)
??? Error: Incomplete or misformed expression or statement.

Error in ==> inline.subsref at 25
    INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);

Error in ==> gyddj2 at 69
 k1=f(Tc,Cc,Uc,Ac,xx(i-1),RK(i-1));
         请高手指点,谢谢!
搜索更多相关主题的帖子: 解微分 方程 MATLAB 
2009-10-14 14:35
leeqee
Rank: 2
等 级:论坛游民
帖 子:1
专家分:20
注 册:2009-10-11
收藏
得分:20 
您看的是哪一本书的亚 我还是个菜鸟
2009-10-21 22:29
快速回复:MATLAB解微分方程
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.026677 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved