MATLAB，BVP4C问题非线性微分方程问题
functiondydx=uodeinit(x,y)globalrh;dydx=[y(2)-h*y(1)/r+y(1)^2/r];functionres=ubc(ya,y...
function dydx=uodeinit(x,y)
global r h;
dydx=[y(2)
-h*y(1)/r+y(1)^2/r];
function res = ubc(ya,yb)
res = [ya(2)
yb(2) ];
global r h;
r=0.1;h=0.5;
solinit = bvpinit(linspace(0,1,21),[1 0]); %bvp4c所需的猜测网格
usolutioninit = bvp4c(@uodeinit,@ubc,solinit);%求出u1
x=usolutioninit.x
t=usolutioninit.y;
length(x)
怎么算出来x的列向量数是length(x)=61，很神奇啊，我给的差分是分成21个点的.为什么会这样，是不是哪里有问题，而且，r=0.3,h=0.5算出来length(x)=41 展开
global r h;
dydx=[y(2)
-h*y(1)/r+y(1)^2/r];
function res = ubc(ya,yb)
res = [ya(2)
yb(2) ];
global r h;
r=0.1;h=0.5;
solinit = bvpinit(linspace(0,1,21),[1 0]); %bvp4c所需的猜测网格
usolutioninit = bvp4c(@uodeinit,@ubc,solinit);%求出u1
x=usolutioninit.x
t=usolutioninit.y;
length(x)
怎么算出来x的列向量数是length(x)=61，很神奇啊，我给的差分是分成21个点的.为什么会这样，是不是哪里有问题，而且，r=0.3,h=0.5算出来length(x)=41 展开
2个回答
展开全部
正常现象，微分方程的求解在算法上要讨论求解精度的问题，有可能和参数的数值有些关系。
比如说有时我们指定求解步长是1e-3，而实际上的求解过程可能不会按此步长计算，只是在更高精度求解结束后，以1e-3的精度给你返回求解结果。
比如说有时我们指定求解步长是1e-3，而实际上的求解过程可能不会按此步长计算，只是在更高精度求解结束后，以1e-3的精度给你返回求解结果。
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
瑞地测控
2024-08-12 广告
在苏州瑞地测控技术有限公司，我们深知频率同步与相位同步的重要性。频率同步确保两个或多个设备的时钟频率变化一致，但相位（即时间点）可保持相对固定差值。而相位同步，即时间同步，要求不仅频率一致，相位也必须完全一致，即时间差恒定为零。相位同步的精...
点击进入详情页
本回答由瑞地测控提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
