How to iterate with fsolve to find Mach number

10 次查看(过去 30 天)
I'm trying to interate static pressure ratio to find Mach number using fsolve, but I can't get the syntax correct. The Mathworks explaination didn't really help me. Could anyone out there help?
I have a set range for the pressure ratio, but the Mach number is unknown. Here is the equation:
Capture.PNG
This is what I have so far for the code:
pt2pinf = [1.9, 2.7, 3.5, 4.3, 5.1, 5.9, 6.7, 7.5, 27.5, 47.5, 67.5 ,87.5 ,107.5,127.5, 147.5, 200, 252.5, 305, 357.5, 410, 462.5 ,515 ,567.5, 620, 672.5, 725, 777.5, 830 ,882.5, 935, 987.5, 1000];
Me0=ones(1,32);
g=1.2;
pt2pinf_test = fsolve(@(Me)(((1+((g-1)/2)*Me.^2)^(g/(g-1)))/((((2*g)/(g+1))*Me.^2)-((g-1)/(g+1)))),Me0);
  8 个评论
Torsten
Torsten 2019-9-10
They are stored in the fsolve return parameter (pt2pinf_test in your case).
Hailey
Hailey 2019-9-10
编辑:Hailey 2019-9-10
The values are too small to be a Mach number after a shock wave.
Not sure what I did wrong.

请先登录,再进行评论。

采纳的回答

Torsten
Torsten 2019-9-10
编辑:Torsten 2019-9-10
function main
pt2pinf = [1.9, 2.7, 3.5, 4.3, 5.1, 5.9, 6.7, 7.5, 27.5, 47.5, 67.5 ,87.5 ,107.5,127.5, 147.5, 200, 252.5, 305, 357.5, 410, 462.5 ,515 ,567.5, 620, 672.5, 725, 777.5, 830 ,882.5, 935, 987.5, 1000];
Me0 = 5.0;
g = 1.2;
fun = @(Me,pquot)((((1+((g-1)/2)*Me.^2).^(g/(g-1)))./((((2*g)/(g+1))*Me.^2)-((g-1)/(g+1))))-pquot).^2
for i = 1:numel(pt2pinf)
pquot = pt2pinf(i);
machnumber(i) = fminsearch(@(Me)fun(Me,pquot),Me0,optimset('TolFun',1e-10));
end
machnumber
fun(machnumber,pt2pinf)
end

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Gas Dynamics 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by