马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我碰到了两个问题,已经多天查找,没有结果。请大家帮我解决这两个难题。
第一个问题是这样:
根据图一的方程(其中T0 = 0.1, epsilon = 0.1 a = 0.1 0 <= b <= 4.0; 0 <= c <= 100),求出图三中的相图,它的坐标为b, c, 对应的稳定点区域和非稳定的相图区域边界。I:单稳极限环,II单稳,III:双稳
我的想法是先求jacobian,其行列式为0, 然后求含有b, c 参数的方程 (左边微分用0代替),最后将三个方程联合,得到bc的一个函数,然后ezplot就可以得到最后的图。
但是我不会用solve解三个方程,四个参数的非线性方程,方程中不包含x(1), x(2). 下面是我的代码,但是在solve的地方我不会了。- %Condition:
- %(1) F1 = F2 = 0 (2) Jacobian = 0
- %
- a = 0.1; epslion = 0.1; T0 = 5.0;
- %use the symbolic calcution
- syms x1 x2 y b c
- F = [((epslion ^ 2 + x1 ^ 2 ) / ( 1 + x1 ^ 2) ) / (1 + x2) - a * x1;
- b / T0 - x2 / ( T0 * (1 + c * x1 ^ 2) )];
- V = [x1, x2];
- %calculate the jacobian matrix
- J = jacobian(F, V);
- J_det = det(J);
- %function fc = g(x)
- y(1) = ((epslion ^ 2 + x1 ^ 2 ) / ( 1 + x1 ^ 2) ) / (1 + x2) - a * x1;
- y(2) = b / T0 - x2 / ( T0 * (1 + c * x1 ^ 2) );
- y = [J_det y(1) y(2)];
- result = solve(y, 'b', 'c');
复制代码 第二个问题,是图2的偏微分。其中D= 0, 右边只有前两项。matlab中没有这种类型的方程,它只能解elliptic, parobalic 之类的,然后PDE tool双只能解两阶,我这个是一阶三参数方程(t, x(1), x(2))这三个变量。
它的边值条件是反射边值条件。J(x, t) = FP − D· 偏P/偏x = 0 (t在任一时候都是)。这个好像第三类边界条件。
|