T
Author
Em/Mình có viết một chương trình để giải 1 hệ phương trình vi phân. Nhưng Matlab báo lỗi sau đây:
Warning: Explicit solution could not be found.
> In solve at 160
Mình cũng mới học Matlab ko lâu, mong các anh chị và các bạn chia sẻ kinh nghiệm. Chương trình của mình như sau:
Cửa sổ Editor 1:
clear all; close all; clc;
init=bvpinit(linspace(0,4,100),[0 0 0 0 0 0 0 0 0])
sol=bvp4c(@rhs_bvp,@bc_bvp,init);
x=linspace(0,4,1000);
BS=deval(sol,x);
plot(x,BS);
Cửa sổ Editor 2:
function rhs=rhs_bvp(x,y)
rhs=[y(2)
y(3)
-3*y(1)*y(3)+2*y(2)*y(2)-y(5)+109.963*y(8)
y(6)
-3*5.076*y(1)*y(6)-0.0000007738*y(9)*y(6)-0.000002255*y(6)*y(6)
y(9)
-3*5.076*42770*y(1)*y(9)+2.914*3*5.076*y(1)*y(6)+2.914*0.0000007738*y(9)*y(6)+2.914*0.000002255*y(6)*y(6)];
Cửa số Editor 3:
function bc=bc_bvp(yl,yr)
bc=[yl(1)
yl(2)
yl(5)-1
yl(8)-1
yr(2)
yr(5)
yr(8)];
Warning: Explicit solution could not be found.
> In solve at 160
Mình cũng mới học Matlab ko lâu, mong các anh chị và các bạn chia sẻ kinh nghiệm. Chương trình của mình như sau:
Cửa sổ Editor 1:
clear all; close all; clc;
init=bvpinit(linspace(0,4,100),[0 0 0 0 0 0 0 0 0])
sol=bvp4c(@rhs_bvp,@bc_bvp,init);
x=linspace(0,4,1000);
BS=deval(sol,x);
plot(x,BS);
Cửa sổ Editor 2:
function rhs=rhs_bvp(x,y)
rhs=[y(2)
y(3)
-3*y(1)*y(3)+2*y(2)*y(2)-y(5)+109.963*y(8)
y(6)
-3*5.076*y(1)*y(6)-0.0000007738*y(9)*y(6)-0.000002255*y(6)*y(6)
y(9)
-3*5.076*42770*y(1)*y(9)+2.914*3*5.076*y(1)*y(6)+2.914*0.0000007738*y(9)*y(6)+2.914*0.000002255*y(6)*y(6)];
Cửa số Editor 3:
function bc=bc_bvp(yl,yr)
bc=[yl(1)
yl(2)
yl(5)-1
yl(8)-1
yr(2)
yr(5)
yr(8)];