В настоящее время я пытаюсь рассчитать распределение скорости по профилю для проблемы на работе и застрял на том, как правильно ввести это в Matlab. Я знаю, что x - фиктивная хордовая координата, которую я должен использовать в интеграле, чтобы отличить guish от x, но я не совсем понимаю, как передать это Matlab. Как отличить guish 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, который включается в интегрируемое уравнение
