我怎样才能使我歌唱图有意义吗?

1视图(30天)

显示旧的评论

伊恩•萨缪尔森 2021年10月18日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1566551-how-can-i-make-my-ode-graph-make-sense

评论道: 明星黾 2021年10月18日

答:接受明星黾

我不知道为什么我的图是这样的!我试着画出微分和解决歌唱但我的图不合理。

我想做一个图的颂歌球滴油,所有值都是随机的,情节的

清晰;clc;

clearvars;

信谊t v (t) aa g p p1

aa = 5;

u = 3;

g = 10

p1 = 9;

p = 3;

v0 = 0;

气孔导度= (0)= = v0

eqn = diff (v, t) = = (1-p1 / p) * g - (9 * u / aa ^ 2 * 2 * p) * v

vSol (t) =简化(dsolve (eqn电导率))

vSol fplot (t (t), 10000年[0]);抓住;网格;

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

登录置评。

在回答这个问题。

接受的答案

明星黾 2021年10月18日

1
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1566551-how-can-i-make-my-ode-graph-make-sense answer_811161

的 fplot 电话不正确,功能下降得如此之快,这是不可能的时间间隔从看得清楚一些 0 来 1 e + 5 。

正确的和它的工作原理!

信谊t v (t) aa g p p1

aa = 5;

u = 3;

g = 10

p1 = 9;

p = 3;

v0 = 0;

气孔导度= (0)= = v0

气孔导度=

eqn = diff (v, t) = = (1-p1 / p) * g - (9 * u / aa ^ 2 * 2 * p) * v

eqn (t) =

vSol (t) =简化(dsolve (eqn电导率))

vSol (t) =

图

fplot (vSol (t)[0 1依照]);持有在;网格在;

。

8的评论
显示7年长的评论隐藏7年长的评论

伊恩•萨缪尔森 2021年10月18日

编辑:伊恩•萨缪尔森 2021年10月18日

@Star黾抱歉打扰你,你帮助我很多,但是现在我想阴谋vetorial营的流箭袋,但我不知道我做错了什么

信谊t v (t) aa g p p1

%α= r

aa = 5;

% u = mi

u = 3;

%的重力

g = 10

%密度'

p1 = 9;

%密度

p = 3;

v0 = 0;

气孔导度= (0)= = v0

eqn = diff (v, t) = = (1-p1 / p) * g - (9 * u / aa ^ 2 * 2 * p) * v

vSol (t) =简化(dsolve (eqn电导率))

fplot (vSol (t)[0 1依照]);抓住;网格;

ang =每股(eqn (V, T));

xx = linspace (5、5、30);yy = xx;

[V, T] = meshgrid (xx和yy);

U = cos (ang);

V =罪(ang);

箭袋(X, Y, U, V);抓住;网格;

%颤抖(t, v (t) [0 10 100]);

我试图让箭头和箭袋,但是我不知道这是一个失败的原因,我在xx和yy单独的最主要的事情,我不知道问题是一切都对t。

明星黾 2021年10月18日

没烦恼!

的箭袋功能想要的数值,而不是象征性的价值观,和那些很容易提供。

我不确定什么是期望的结果,不过我创建了一些结果,可能接近于它。如果他们不是,请详细描述期望的结果是什么,我将尽我所能帮助提供。

信谊t v (t) aa g p p1

%α= r

aa = 5;

% u = mi

u = 3;

%的重力

g = 10

%密度'

p1 = 9;

%密度

p = 3;

v0 = 0;

气孔导度= (0)= = v0

气孔导度=

eqn = diff (v, t) = = (1-p1 / p) * g - (9 * u / aa ^ 2 * 2 * p) * v

eqn (t) =

vSol (t) =简化(dsolve (eqn电导率))

vSol (t) =

fplot (vSol (t)[0 1依照]);持有在;网格在;

甘氨胆酸Ax =;

%得到(Ax)

T = Ax.Children.XData;%从情节获得数据

V = Ax.Children.YData;%从情节获得数据

ang =量化(V, T);%使用四象限反正切

% xx = linspace (5、5、30);yy = xx;

% (V, T) = meshgrid (xx和yy);

% Uq = cos (ang);

% Vq =罪(ang);

Uq =梯度(T);%的数值微分

Vq =梯度(V);%的数值微分

箭袋(T V Uq Vq);持有在;网格在;

%颤抖(t, v (t) [0 10 100]);

Experimment得到期望的结果。

。

伊恩•萨缪尔森 2021年10月18日

男人我apreciate你的帮助,但是看看我有什么,我让我的颤抖,meshgrid告诉我梯度,现在我得到了我想要的,我试着调整规模。

信谊t v (t) aa g p p1

%α= r

aa = 3;

% u = mi

u = 1;

%的重力

g = 10

%流体密度

p1 = 2;

%密度

p = 1.5;

v0 = 0;

气孔导度= (0)= = v0

eqn = diff (v, t) = = (1-p1 / p) * g - (9 * u / aa ^ 2 * 2 * p) * v

vSol (t) =简化(dsolve (eqn电导率))

fplot (vSol (t)[0 1依照]);抓住;网格;

xx = linspace(0,1的军医,10);

yy = vSol (xx);

(X, Y) = meshgrid (xx和yy);

% ang =每股(eqn (X, Y));

% U = cos (ang);

% V =罪(ang);

箭袋(X, Y);抓住;网格;

%颤抖(t, v (t) [0 10 100]);

明星黾 2021年10月18日

的持有条件可以连接在或从取决于其他数据绘制同一轴,或不是。

登录置评。

类别

MATLAB 数学数值积分和微分方程

找到更多的在数值积分和微分方程在帮助中心和文件交换

社区寻宝

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

开始狩猎!