MATLAB的答案

帮我把这个方程

3视图(30天)
夏玛穆罕默德
夏玛穆罕默德 2023年4月4日
评论道: 夏玛穆罕默德2023年4月4日
我转换方程,得到这个数字。
我忘了怎样但我失去了我写代码的文件上面的图。
下面我想rewrtie,但它不工作。指数函数为零或无穷大。
编辑!
这是R_e原始方程。
我检查我的报告和他们的参数是相同的。
也许这里有一个问题吗?
我做了你问到最后的分母,仍然相同的结果。
clc
Contact_angle = 2.4783666516568;
体积= 5;
ST_L = 72;
密度= 0.973;
重力= 9.807;
Shape_factor = 37.1;
T_zero = 0;
Dynamic_viscosity = 8.9的军医;
R_e =(4 *体积/(π* Contact_angle)) ^ 1/3;
T = 0:0.1:6;
R = R_e * (1 - exp (- (2 * ST_L / R_e ^ 12 +重力密度* / (9 * R_e ^ 10)) * (24 * Shape_factor *体积^ 4 * (T + T_zero) /(π^ 2 * Dynamic_viscosity))))。^ (1/6);
情节(T, R);
谢谢大家的帮助
8的评论
显示 7年长的评论
夏玛穆罕默德
夏玛穆罕默德 2023年4月4日
是的对不起它不是原始方程。检查更新的问题

登录置评。

在回答这个问题。

答案(3)

山姆翟
山姆翟 2023年4月4日
编辑:山姆翟 2023年4月4日
跑:
嗨 @shammas穆罕默德
编辑: 有两个校正方程。谢谢 @Walter罗伯森 指出。
%的参数
Contact_angle = 2.4783666516568;
体积= 5;
ST_L = 72;
密度= 0.973;
重力= 9.807;
Shape_factor = 37.1;
T_zero = 0;
Dynamic_viscosity = 8.9的军医;
R_e =(4 *体积/(π* Contact_angle)) ^ 1/3;
%条款建议分手
T = linspace (0 6 e-12 601);
f5 = 24 * Shape_factor *体积^ 4 * (T + T_zero) /((π^ 2)* Dynamic_viscosity);% < - - - - - -修正
f4 =(密度*重力)/ (9 * (R_e ^ 10));
f3 = (2 * ST_L) / (R_e ^ 12);
f2 = f3 + f4;
f1 = - f2 * f5;
R = R_e * (1 - exp (f1)) ^ (1/6);% < - - - - - -修正
情节(T, R)、网格在,包含(“T”)
1评论
夏玛穆罕默德
夏玛穆罕默德 2023年4月4日
谢谢你的更正。它给了我一个好的开始的初始速率半径会呈几何级数增长。然而仍有问题,函数并不增加。(指数的部分是0)。
我不明白这是为什么发生。我编辑的问题,如果你需要更多的信息。

登录置评。

沃尔特·罗伯森
沃尔特·罗伯森 2023年4月4日
x ^ 1/3意味着(x ^ 1)除以3 x / 3。如果你想提高1/3的力量需要x ^ (1/3)
当心^ 1/6同样的问题。
1评论
Dyuman Joshi
Dyuman Joshi 2023年4月4日
沃尔特,我在代码中实现这个在我的评论,它仍然没有给OP实现任何相似的结果。

登录置评。

克里斯·拉皮埃尔
克里斯·拉皮埃尔 2023年4月4日
编辑:克里斯·拉皮埃尔 2023年4月4日
跑:
写,你需要包括期末分母在括号,或者除以 Dynamic_Viscocity 。
我们需要你给其余的参数是正确的。一个潜在的错误调查——你的角度应该在弧度而不是度吗?
Contact_angle = 142;
体积= 5;
ST_L = 72;
密度= 0.973;
重力= 9.807;
Shape_factor = 37.1;
T_zero = 0;
Dynamic_viscosity = 8.9的军医;
R_e =(体积400 * /(π* Contact_angle)) ^ 1/3;
T = 0:300;
R = R_e * (1 - exp (- (2 * ST_L / R_e ^ 12 +(密度*重力)/ (9 * R_e ^ 10)) * 24 * Shape_factor *体积^ 4 * (T + T_zero) /(π^ 2 * Dynamic_viscosity)))。^ (1/6);
情节(T, R);
6个评论
显示 5年长的评论
夏玛穆罕默德
夏玛穆罕默德 2023年4月4日
我明白了谢谢

登录置评。

类别

找到更多的在表面和网格图在帮助中心和文件交换

标签

s manbetx 845


释放

R2022b

社区寻宝

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

开始狩猎!