如何编写解决这类方程组?

1视图(30天)
NIRUPAM SAHOO
NIRUPAM SAHOO 2021年9月11日
评论道: NIRUPAM SAHOO2021年9月12日
1评论
NIRUPAM SAHOO
NIRUPAM SAHOO 2021年9月12日
请人解决这个问题。这是u和v r的函数。

登录置评。

答案(2)

Wan霁
Wan霁 2021年9月12日
嘿,朋友
只剩下扩大项目的两个方程,提取u和v”,那么一个数值解算器。

沃尔特·罗伯森
沃尔特·罗伯森 2021年9月12日
编辑:沃尔特·罗伯森 2021年9月12日
我没能找出被提高到10/9。我使用 乱涂乱画 代替。
信谊u (r)(右)
信谊N潦草地写真实的
假设(r,“真实”的)
du = diff (u);
dv = diff (v);
left1 = diff (r ^ (n - 1) * du ^ 3)
left1 (r) =
right1 = r ^ (n - 1) * sqrt (u) *倍根号(v) / (3 * r ^ (2/3) * sqrt(1 + 9 *花体^ (10/9)/ (10 * (3 * 2)^ (1/3))))
right1 (r) =
候补委员= diff (r ^ (n - 1) * dv ^ 3)
候补委员(r) =
right2 = r ^ (n - 1) *你* v / (3 * r ^ (2/3)
right2 (r) =
eqn1 = left1 = = right1
eqn1 (r) =
eqn2 =候补委员= = right2
eqn2 (r) =
ic = [u (0) = = 1, v (0) = = 1, du (0) = = 0, dv (0) = = 0)
ic =
索尔= dsolve ([eqn1 eqn2, ic])
警告:无法找到具有象征意义的解决方案。
索尔=(空符号)
字符串(eqn1)
ans =“3 * r ^ (N - 1) * diff (u (r), r) ^ 2 * diff (u (r), r, r) + r ^ (N - 2) * (N - 1) * diff (u (r), r) ^ 3 = = (r ^ (N - 1) * u (r) ^ (1/2) * v (r) ^ (1/2)) / (3 * r ^(2/3) *((9 *花体^ (10/9))/ (10 * (3 * N - 2) ^ (1/3)) + 1) ^ (1/2))”
字符串(eqn2)
ans =“3 * r ^ (N - 1) * diff (v (r), r) ^ 2 * diff (v (r), r, r) + r ^ (N - 2) * (N - 1) * diff (v (r), r) ^ 3 = = (r ^ (N - 1) * u (r) * v (r)) / (3 * r ^ (2/3))”
字符串(集成电路)
ans =1×4字符串数组
“u (0) = = 1”“v(0) = = 1”“潜艇(diff (u (r), r), r, 0) = = 0”“潜艇(diff (v (r), r), r, 0) = = 0 "
缺乏一个象征性的数字解决方案解决方案意味着你将不得不做,但你不能做一个数字解决方案到正无穷,你肯定不会得到一个公式。金宝搏官方网站
5个评论
NIRUPAM SAHOO
NIRUPAM SAHOO 2021年9月12日
信谊p (t) m (t)t Y
(diff命令= ((t ^ (100 - 1)) * (diff (p (t), t))) = = (t ^ (100 - 1)) * (t - 1) * exp (t) * p (t) * m (t);
diff ((t ^ (100 - 1)) * (diff (t m (t)))) = = (t ^ (100 - 1)) * (t - 1) * exp (t) * p (t) ^ (1/2) * m (t) ^ (1/2)
[谜底,潜艇]= odeToVectorField(命令);
DEFcn = matlabFunction(谜底,“var”,{t、Y});
tspan = [0100];
y0 = [0 0 0 0];
[t、Y] =数值(DEFcn tspan, y0);
情节(t, Y)
# #我写但图不显示。请帮我

登录置评。

标签

下载188bet金宝搏


释放

R2021a

社区寻宝

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

开始狩猎!