如何解决耦合的二阶2依赖和1独立变量方程?

1视图(30天)

显示旧的评论

Atharva骡 2021年9月6日

0
链接

这个问题直接联系

https://it.mathworks.com/matlabcentral/answers/1447579-how-to-solve-coupled-2nd-order-2-dependent-and-1-independent-variable-equation

评论道: Atharva骡2021年9月15日

eqn.png

我新MATLAB,想解下列方程(形象)。

这里c1和ck是变量的函数z . v1, D1, vk、Dk, k, Kd和c_cmc都是常数。

边界条件- > c1 (0) = 0, c1 (0.1) = 0.025, ck (0) = 0.125, ck (0.1) = 0

我怎么解这个方程?请帮助。

甚至与其他值的常量,我怎么解决这种类型的方程?

v1 = 2.03 * 10 ^ -14;

vk = 1.3 * 10 ^ -12;

D1 = 9.8 * 10 ^ -10;

Dk = 1.2 * 10 ^ -10;

k = 3;

Kd = 10 ^ 5;

c_cmc = 8.3334;

我= imread (“eqn.png”);

imshow(我)

2的评论
显示1年长的评论藏1年长的评论

Atharva骡 2021年9月15日

                                v1 = 1
                               
                                vk = 5
                               
                                D1 = 3
                               
                                Dk = 1
                               
                                k = 2
                               
                                Kd = 10
                               
                                c_cmc = 8
                               
                                信谊c1 (z) ck (z)
                               
                                eq1 = diff (c1、2) = = (1 / D1) * (Kd * c_cmc ^ (1 k) * k * c1 ^ k - Kd * k * ck - v1 * diff (c1))
                               
                                eqk = diff (ck, 2) = = (1 / Dk) * (Kd * ck - Kd * c_cmc ^ (1 k) * c1 ^ k - vk * diff (ck))
                               
                                eqn = [eq1 eqk]
                               
                                [V S] = odeToVectorField (eqn)
                               
                                M = matlabFunction (V,“var”,{“t”,“Y”})
                               
                                间隔= 0.1 [0]
                               
                                yInit = (2 0 1 1)
                               
                                ySol =数值(M,间隔,yInit)
                               
                                ySol.y
                               
                                zValues = linspace (0、0.1、10);
                               
                                yValues =德瓦尔(ySol zValues 1);
                               
                                yValues =德瓦尔(ySol zValues 2);
                               
                                yValues =德瓦尔(ySol zValues 3);
                               
                                yValues =德瓦尔(ySol zValues 4);
                               
                                情节(zValues yValues)

谁能请让我知道什么yValues =德瓦尔(1、2、3和4)表示吗?1个人ck值或diff (ck)价值观?也请让我知道yInit =[2 0 1 1]表示初始值的变量(c1 diff (c1)或差异(c1、2) ?

@Bill格林

登录置评。

在回答这个问题。

答案(0)

在回答这个问题。

类别

数学和优化偏微分方程的工具箱几何和网格

找到更多的在几何和网格在帮助中心和文件交换

下载188bet金宝搏

释放

R2020b

社区寻宝

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

开始狩猎!

翻译的