0

该表显示了cD作为雷诺数函数的球体阻力系数Re。在. cD_ Re = 5, 50, 500, 5000提示:使用对数-对数标度。

Re 0.2 2 20 200 2000 20000
cD 103 13.9 2.72 0.800 0.401 0.433

我不明白如何“使用对数刻度”来解决这个问题。我尝试了以下代码,但我不知道它是否正确。我如何在这里使用对数刻度?

您可以在下面的链接中找到问题的原件

https://d2vlcm61l7u1fs.cloudfront.net/media%2Faef%2Faef9a122-1a83-423d-9258-66fac0cac8e4%2FphpCF4Vwd.png

r = [0.2, 2, 20 ,200, 2000 ,20000];
c = [103, 13.9, 2.72 ,0.800 ,0.401 ,0.433];
rI = [5,50,500,5000];
cI = spline (r,c, rI);
4

1 回答 1

1

对于这个问题,您需要做的是计算X 和 Y 值的对数,然后才执行插值。如果你看一个球体的阻力系数

镉与稀土

您会看到此图表(和其他类似图表)是使用对数刻度绘制的。请注意,除了在 3E5 附近发生流动分离的过渡点之外,红色曲线非常平滑且表现良好。

要解决您的任务,您需要“使用红色曲线”执行插值,即在对数域中。这样做的原因是因为 X 和 Y 值的跨度非常大,并且原始域上的多项式或样条不能正确捕获行为。实际上——

r = [0.2, 2, 20 ,200, 2000 ,20000];
c = [103, 13.9, 2.72 ,0.800 ,0.401 ,0.433];
rI = [5,50,500,5000];
cI = exp( spline( log(r), log(c), log(rI)) ); % interpolation is performed for log(y) vs log(x)
%{
cI =
    6.9390    1.5843    0.5636    0.3717
%}

可以使用图表手动验证这些结果的正确性。

于 2019-05-06T09:41:26.770 回答