作业帮求大神帮忙解决一道关于matlab求解微分方程组的边值问题
来源：学生作业帮 编辑：拍题作业网作业帮 分类：数学作业 时间：2024/04/28 04:34:26
求大神帮忙解决一道关于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为一组猜测值
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为一组猜测值
yinit 应该是向量（代表对微分方程几个变量的估计值）,
再问： 意思是说只能对每个变量分别估计一次吗？
再答： 对。bvpinit会把你的估计值按照指定的时间范围进行复制，然后再用bvp4c求解。
再问： 还有啊，这个bcfunc边界条件应该怎么写啊？我写的是[ya(1) yb(1) ya(2) yb(2)]，这样不行啊怎么？
再答： 6个微分方程，应该对应6个边界条件才对，否则方程不封闭。你现在的写法只有4个条件。
再问： 意思是说只能对每个变量分别估计一次吗？
再答： 对。bvpinit会把你的估计值按照指定的时间范围进行复制，然后再用bvp4c求解。
再问： 还有啊，这个bcfunc边界条件应该怎么写啊？我写的是[ya(1) yb(1) ya(2) yb(2)]，这样不行啊怎么？
再答： 6个微分方程，应该对应6个边界条件才对，否则方程不封闭。你现在的写法只有4个条件。