Question

求大神帮忙解决一道关于matlab求解微分方程组的边值问题

clear;clc;
upt=pi;%*(2*rand(1,151)-1);
ita=0.5;%rand(1,151)';
tspan=linspace(0,479,3);
yinit=[1700 100 57 ;0 0.5 0.6014;pi/4 0.7 0;19.03 19.31 19.51;28 33 44.12;15000 4000 3000];
solinit=bvpinit(tspan,yinit);
sol=bvp4c(@rigid,@bcfunc,solinit);
结果无论怎么改总是出现这个，有大神知道是怎么回事吗？？？？？？

注：这是6个微分方程组成的方程组，@bcfunc为边值条件函数，yint为一组猜测值

• 你的回答被采纳后将获得：
• 系统奖励15（财富值+成长值）+难题奖励30（财富值+成长值）+提问者悬赏50（财富值+成长值）

Answer 1

yinit 应该是向量（代表对微分方程几个变量的估计值），你定义成矩阵是什么意思？

追问

意思是说只能对每个变量分别估计一次吗？

追答

对。bvpinit会把你的估计值按照指定的时间范围进行复制，然后再用bvp4c求解。

追问

还有啊，这个bcfunc边界条件应该怎么写啊？我写的是[ya(1) yb(1) ya(2) yb(2)]，这样不行啊怎么？

追答

6个微分方程，应该对应6个边界条件才对，否则方程不封闭。你现在的写法只有4个条件。