0

我目前正在尝试计算机翼上的速度分布以解决工作中的问题，并且一直坚持如何正确地将其输入 Matlab。我知道 x' 是一个虚拟弦坐标，我应该在积分中使用它来将它与 x 区分开来，但我很困惑如何将它传递给 Matlab。如何区分 x' 和 x，然后将所有这些值相加以获得 Matlab 中弦上的完整速度分布？

我的同事建议使用陷阱，但我不确定如何使用 x'。吨

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 插入正在积分的方程中

需要积分的方程

dT/dx 插入正在积分的方程中

4

0 回答 0