How to iterate with fsolve to find Mach number
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Hailey
el 10 de Sept. de 2019
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](https://www.mathworks.com/matlabcentral/answers/uploaded_files/237712/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 comentarios
Torsten
el 10 de Sept. de 2019
They are stored in the fsolve return parameter (pt2pinf_test in your case).
Respuesta aceptada
Torsten
el 10 de Sept. de 2019
Editada: Torsten
el 10 de Sept. de 2019
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 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Gas Dynamics en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!