24小时热门版块排行榜

>论坛更新日志 (3228)
>虫友互识 (827)
>文献求助 (222)
>休闲灌水 (157)
>导师招生 (151)
>考博 (130)
>教师之家 (86)
>硕博家园 (80)
>博后之家 (77)
>找工作 (77)
>论文投稿 (74)
>论文道贺祈福 (67)
>基金申请 (65)
>考研 (57)
>公派出国 (51)
>招聘信息布告栏 (33)

返回列表

木西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);
求大神指点，急急急

回复此楼