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 展开
展开全部
正常现象，微分方程的求解在算法上要讨论求解精度的问题，有可能和参数的数值有些关系。
比如说有时我们指定求解步长是1e-3，而实际上的求解过程可能不会按此步长计算，只是在更高精度求解结束后，以1e-3的精度给你返回求解结果。
比如说有时我们指定求解步长是1e-3，而实际上的求解过程可能不会按此步长计算，只是在更高精度求解结束后，以1e-3的精度给你返回求解结果。
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
Sievers分析仪
2025-02-09 广告
是的。传统上，对于符合要求的内毒素检测，最终用户必须从标准内毒素库存瓶中构建至少一式两份三点标准曲线；必须有重复的阴性控制；每个样品和PPC必须一式两份。有了Sievers Eclipse内毒素检测仪，这些步骤可以通过使用预嵌入的内毒素标准...
点击进入详情页
本回答由Sievers分析仪提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
