我目前正在尝试针对工作中的问题计算机翼上的速度分布,并坚持如何正确地将其输入到Matlab中。我知道x'是我应该在积分中使用的伪弦向坐标,以将其与x区别开来,但是对于如何将其中继到Matlab感到困惑。我如何区分x'和x',然后将所有这些值求和以获得Matlab中弦的完整速度分布?
我的同事建议使用陷阱,但是我不确定如何使用x'。 T
Vt = zeros(length(XC),1); % XC corresponds to the x value used to calculate dT/dxfor i = 1:length(XC)
Prime = XC(i); % x'
dT = dT(i); % dT/dx'
X = XC;
X(X==XPrime) = 0; % remove the x = x' singularity to allow proper integration
int_expression = (1/2*pi)*dT*(X-XPrime);
Vt(i) =-trapz(X,int_expression);
end
积分在链接中并且很简单,并且dT / dx也作为照片发布。 dT / dx插入正在积分的方程式