MATLAB的答案

BVP4C无法解决方程——雅可比矩阵奇异的搭配

33视图(30天)
路加福音马兰
路加福音马兰 2021年4月25日
评论道: MOSLI卡里姆2022年10月8日
你好,我收到这个错误与BVP4c试图解决一个二阶的颂歌。
我解决的内部温度分布使用ode核燃料棒:
边界条件:
在R = 0.2米,外径吗
如果你有什么事,请让我知道。
问和k是常数
r是一个半径变量
%这是我设置的功能。
函数就要= rod_BV (r, y)%的细节需要解决的颂歌
就要= 0 (2,1);
它们(1)= y (2)
它们(2)= (q / k) - ((1 / r) * y (1))%这个方程是无效的在r = 0
结束
函数res = rod_BC (ya,)%细节边界条件
res = 0 (2, 1);
res(1) =丫(2);
res (2) = yb (1) - T_s;
结束
% BVP4C的这是我实现
solinit = bvpinit (linspace (0, R, 4), [0]);%我试过设置这些0值最小正浮点数
选择= bvpset (“RelTol”10的军医,“AbsTol”10 e);
索尔= bvp4c (@rod_BV @rod_BC、solinit选项)

在回答这个问题。

答案(2)

Divija Aleti
Divija Aleti 2021年6月24日
你好卢克,
单一的雅可比矩阵表明,最初的猜测引起分歧的解决方案。BVP4C函数发现解决方案通过求解非线性代数方程组。非线性动力学是一样有效的初始猜他们首先,所以开始改变你猜会有帮助。同时,BVP4C负责提供猜测为解决方案从一个迭代到另一个地方。
详细的解释和可能的解决方案对于这个错误给出了下面的链接:金宝搏官方网站
MATLAB的答案
同时,基于方程的图片你解决,我相信 rod_BV 函数应该写成:
函数就要= rod_BV (r, y)%的细节需要解决的颂歌
就要= 0 (2,1);
它们(1)= y (2)
它们(2)= (q / r) - ((1 / r) * y (2))%这个方程是无效的在r = 0
结束
这是因为 T = y (1) , T ' = y (2) 这意味着, T”= - (q / r) - (1 / r) * y (2) 。
希望这可以帮助!
问候,
Divija
MOSLI卡里姆
MOSLI卡里姆 2022年7月31日
编辑:Torsten 2022年7月31日
跑:
你好,这是我对你提出的解决方案
%这是我设置的功能。
R = 0.2;
q = 1;
k = 1;
每股收益= 1.0000 e-06;
r = linspace (0, r, 4);
solinit = bvpinit (linspace (0, R, 4), [0]);%我试过设置这些0值最小正浮点数
选择= bvpset (“RelTol”10的军医,“AbsTol”10 e);
索尔= bvp4c (@ (r, y) rod_BV (r, y, q, k), @rod_BC, solinit,选项);
y =德瓦尔(溶胶,r);
情节(r, y (:))
函数就要= rod_BV (r, y, q, k)%的细节需要解决的颂歌
就要= 0 (2,1);
它们(1)= y (2);
它们(2)= (q / k) - ((1 / (r + eps)) * y (1));%这个方程是无效的在r = 0
结束
函数res = rod_BC (ya,)%细节边界条件
% res = 0 (2, 1);
% res(1) =丫(2);
% res (2) = yb (1) - 200;
res = [ya; (2) (1) -200];
结束
% BVP4C的这是我实现
% solinit = bvpinit (linspace (0, R, 4), [0]);%我试过设置这些0值最小正浮点数
%选项= bvpset (RelTol, 10的军医,AbsTol, 10 e);
%索尔= bvp4c (@rod_BV、@rod_BC solinit,选项)
2的评论
显示1年长的评论

登录置评。

类别

找到更多的在边值问题在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!