审校颂歌的解决方案

4视图(30天)
塞吉奥Manzetti
塞吉奥Manzetti 2018年1月2日
评论道: Sana Syed2018年2月28日
嗨,MATLAB在线解决了这个:
如果真正的
%的代码
结束
信谊h x Y (x) g B E
eqn = (h ^ 2) * diff (Y, x, 2) +(2 *我* h * g) * diff (Y、x) + (g ^ 2 - E) * Y = = 0;
DY = diff (Y)
cond1 = Y (0) = = 1;
cond2 = DY (0) = = 0;
Y (x) = dsolve (eqn,气孔导度)
然而MATLAB计算机,2017,不解决它,我得到:
diff (Y (x), x)
没有足够的输入参数。
气孔导度误差(第24行)如果issparse (A)
错误PDE_sol(第6行)Y (x) = dsolve (eqn,气孔导度)
哪一个我应该“相信”吗?

接受的答案

捕鸟者
捕鸟者 2018年1月2日
更换线
Y (x) = dsolve (eqn,气孔导度)
Y (x) = dsolve (eqn [cond1 cond2])
9日评论
Sana Syed
Sana Syed 2018年2月28日
如果我有两个耦合的四阶非线性常微分方程和8 BCS它应该在8条件对吧?它不接受任何非线性的边界条件。
信谊uy (x) txd (x)是乌斯(x)用户体验(x);%梁维泰= 0.001,Tz = 0.01 L = 0.1;%常量基于excel的尺寸计算c21 = 42.5015, = 8.42 * 10 ^ 4, c3 = 374.81, = 100, d = 0.03;%装载条件归一化负载fxl = 0, fyl = 0,万立升= 0,手段= 0,mzl = 0, mxdl = 0.3125,万立升= 0;%纯扭
%微分方程输入ode1 = diff (uy x 4) -fxl * diff (uy x, 2) = =万立升* diff (txd x, 2) 2 *手段* diff (txd x) + (1 - x) * * txd手段;%为纯扭转,它减少了D4y / Dx4 = 0;ode2 = (c21 + c22 * fxl) * diff (txd x 2) c3 * diff (txd x 4)——(myl-fzl * (1 - x)) * (mzl + fyl * (1 - x) -fxl * (uy (L) uy)——(myl-fzl * (1 - x)) * txd) = = 0;% ode2 = (c21 + c22 * fxl) * diff (txd x 2) c3 * diff (txd x 4) = = 0;ode3 = diff (x是乌斯,2)= = txd * diff (uy x 2)——(myl-fzl * (1 - x)) /;% ode3 = diff (x是乌斯,2)= = txd * diff (uy x 2) - (1 /) * (myl-fzl * (1 - x));ode4 = d * (diff (ux, x) + 0.5 * (diff (uy x)) ^ 2 + 0.5 * (diff(乌斯x)) ^ 2) + (((+ 1) * (diff (txd x)) ^ 2) / 2) = = fxl;%的微分方程输入
%定义导数值的变量Duy = diff (uy x);D2uy = diff (uy x, 2);D3uy = diff (uy x 3);Dtxd = diff (txd x);D2txd = diff (txd x, 2);D3txd = diff (txd x 3);Duz = diff(乌斯x);D2uz = diff (x是乌斯,2);D3uz = diff(是乌斯x 3);%的导数值的变量
%的几何边界条件cond1 =用户体验(0)= = 0;cond2 = uy (0) = = 0;cond3 =是乌斯(0)= = 0;cond4 = txd (0) = = 0;cond5 = Duz (0) = = 0;cond6 = Duy (0) = = 0;cond7 = Dtxd (0) = = 0;cond8 = Dtxd (L) = = 0;%的几何边界条件
%加载边界条件
%这结束条件给错误:cond9 = -D3uy (L) - * (2 * txd (L) * Dtxd (L) * D2uy (L) + txd (L) ^ 2 * D3uy (L) -txd (L) * D3uz (L) -Dtxd (L) * D2txd (L)) + fxl * Duy (L) = = fyl;
%这个条件也给错误cond10 = D2uy (L) +一个* (txd (L) * D2uy (L) -D2uz (L)) * txd (L) = = mzl
cond11 = c3 * D3txd (L) = = mxdl;
%常微分方程= [ode1 ode2 ode3 ode4]
常微分方程= [ode1 ode2]
气孔导度= [cond2 cond4 cond6 cond7 cond8]
% (uySol (x) txSol (x) uzSol (x) uxSol (x) = dsolve(诗赋,气孔导度)
[uySol (x) txSol (x) = dsolve(诗赋,气孔导度)

登录置评。

更多的答案(0)

社区寻宝

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

开始狩猎!