24小时热门版块排行榜

返回列表

木西666

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 11.5
帖子: 3
在线: 5小时
虫号: 5121000
注册: 2016-10-14
专业: 机械工程

[求助] 非线性微分方程组的边值问题已有1人参与

我现在是一个六阶的非线性微分方程，但是只有三个初始值，其他都是边界条件。我现在的难处就是想把边界条件转换成初始条件。
clear
clc
close all
global T;
for T = 1:20000
solinit = bvpinit(linspace(0.00000001,20,200),[10^(-5) 0 0 0.2 0 0.9]);% 解的初始估计
bcfun = @(ya,yb)[ya(1)-10^(-5); yb(1)-2.25;ya(2);yb(2)-10;ya(3);yb(3)-18];
sol = bvp4c(@ringbal1,bcfun,solinit);
s=linspace(0.00000001,20,2000);
y=deval(sol,s);
end
X=10.*y(:,1).*cos(y(:,2));
Z=-10.*y(:,3);
plot3(T,X,Z);
程序报错为：
??? Error using ==> bvp4c at 252
Unable to solve the collocation equations -- a singular Jacobian encountered

Error in ==> dabafa_ at 8
sol = bvp4c(@ringbal1,bcfun,solinit);
求大神指点，急急急

回复此楼