Hello, I'm trying to code this expression, this is my attempt but I'm not getting the same answer/value, not sure what I'm doing wrong. Any help/guidance will be appreciated.
% Given values
%AInfinity = 1;
%AInlet = 6.5;
MInlet = 0.5;
MInfinity = 2;
gamma = 1.4;
a = (1 + (gamma-1/2)* MInfinity^2)
b = (1 + (gamma-1/2)* MInlet^2)
b1 = (MInlet/MInfinity)
r = b1*(a/b)^(gamma + 1/(2*gamma -1))

 채택된 답변

VBBV
VBBV 2024년 3월 2일
편집: VBBV 2024년 3월 2일

0 개 추천

% Given values
%AInfinity = 1;
%AInlet = 6.5;
MInlet = 0.5;
MInfinity = 2;
gamma = 1.4;
a = (1 + ((gamma-1)/2)* MInfinity^2) % for numerator
a = 1.8000
b = (1 + ((gamma-1)/2)* MInlet^2) % for numerator
b = 1.0500
b1 = (MInlet/MInfinity)
b1 = 0.2500
r = b1*(a/b)^((gamma + 1)/(2*(gamma -1))) % for num and den
r = 1.2595

댓글 수: 3

VBBV
VBBV 2024년 3월 2일
use the parenthesis for num & den in the expression
Thank you, I just would like to know if these are well coded? Thanks
p02 = pa * (1 + eta_d * (gamma - 1) / 2 * Ma^2) ^ (gamma / (gamma - 1))
p2 = p02 * (1 + (gamma - 1) / 2 * M2^2) ^ (gamma / (gamma - 1))
Torsten
Torsten 2024년 3월 17일
I just would like to know if these are well coded?
Yes, the two expressions are correctly coded.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Gamma Functions에 대해 자세히 알아보기

질문:

2024년 3월 2일

댓글:

2024년 3월 17일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by