用matlab求dy/dx=y+x/y-x的通解

 我来答
Miafifi
2015-06-20 · TA获得超过451个赞
知道小有建树答主
回答量:277
采纳率:93%
帮助的人:171万
展开全部

Hi your question is quite difficult to solve out. 


Since the common method combined with [dsolve()] + [diff] doesn't work at all.

e.g if you type the following command to generate the general solution of this ODE, the IDE should warn you that "Explicit solution could not be found."

syms y(x)
Sol = dsolve( diff(y,x) == y+x/y-x );

Therefore, my current solution is just alternatively get the graphic answer of your equation. But I will dig further in the field and give a better solution if I would find in future.

Generally, I manually set a x domain at  [-8 8], and inital y value at 2, and then I use ode45 command to help me generate the equation graph. Hope you could like it.

Anyhow, here is my code.

clc; clear all; clf

%  dy          x
% ---- == y + --- - x
%  dx          y
ode1 = @(x,y)(y+x/y-x);

% assume -8 < x < 8 and initial y0 = 2
[x,y] = ode45(ode1, [-8:.01:8], 2.0);

plot(x, y, 'r', 'linewidth', 4)
xlabel('x'), ylabel('y'), grid on
title('solution to ODE dy/dx = y+x/y-x', 'fontsize', 14)
ylim([min(y)-1E7 max(y)+1E7])

line([0 0], ylim, 'color', 'k', 'linewidth', 2);
line(xlim, [0 0], 'color', 'k', 'linewidth', 2);

The graph should look like

推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式