loops

7 次查看（过去 30 天）
marvin corado
marvin corado 2011-3-11
ok i made my self a simpler code. what i am trying to accomplish is for MatLab to give me the M at which A and z are equal. Here is my code.
A=1.687;
gamma=1.4;
tol=1e4;
for M=1:1:5
a=(1/M^2);
b=(2/(gamma+1));
c=(1+(((gamma-1)/2)*M^2));
d=(gamma+1)/(gamma-1);
z(M)=sqrt(a*((b*c)^d));
end
z'
if abs(A-z)<tol
M=z(M)
end
i cannot seem to get it to work. i got z' and i know M=2 is pretty close to my A value. FOr M=2 my code returns 1.6875. the if statement prints M=25.00, but i want it to print 2 because thats the answer. Thank you

请先登录，再回答此问题。

回答（5 个）

Moes
Moes 2011-3-11
Try this code instead (no loops):
A=1.687;
gamma=1.4;
tol=1e4;
M = 1:5;
a=(1./M.^2);
b=(2/(gamma+1));
c=(1+(((gamma-1)/2).*M.^2));
d=(gamma+1)/(gamma-1);
z=sqrt(a.*((b.*c).^d));
diffrnc = abs(A*ones(size(z)) - z)
diffrnc(find(diffrnc < tol))
  1 个评论
Moes
Moes 2011-3-11
If you just want the index and not the actua;l number replace "diffrnc(find(diffrnc < tol))" with "find(diffrnc < tol)".

请先登录，再进行评论。

Walter Roberson
Walter Roberson 2011-3-11
z is an array, so A-z is an array and so abs(A-z) is an array. That array is compared to tol, which gives a logical array in return. When you use a logical array in an "if" statement, the "if" is only true if all of the logical values are true.
Replace the "if" with something like,
Mpos = find(abs(A-z)<tol);
M = z(Mpos);
That is a pretty large tol you have. Do you perhaps mean 1e-4 ?
marvin corado
marvin corado 2011-3-11
well i initially put 1e-4 but it did not give me any answer. so i decreased it and then increased it until i got something. I was just trying to get a result, but i will try the 1e-4 because that made sense to me. thankyou
marvin corado
marvin corado 2011-3-11
i tried and here is the result. Really what i am after is the M for which z is equal to A. That is really what i am after. i modified it as you suggested.
A=1.687;
gamma=1.4;
tol=1e-4;
for M=1:1:5
a=(1/M^2);
b=(2/(gamma+1));
c=(1+(((gamma-1)/2)*M^2));
d=(gamma+1)/(gamma-1);
z(M)=sqrt(a*((b*c)^d));
end
Mpos=find(abs(A-z)<tol);
M=z(Mpos)
here is the result M =
Empty matrix: 1-by-0
thank you again for helping
Matt Fig
Matt Fig 2011-3-11
Is this what you are looking for?
A=1.687;
gamma=1.4;
tol=1e-4;
IDX = 0:.0001:2;
z = zeros(1,length(IDX)); % Pre-allocate!
cnt = 0;
for M=IDX
a = (1/M^2);
b = (2/(gamma+1));
c = (1+(((gamma-1)/2)*M^2));
d = (gamma+1)/(gamma-1);
cnt = cnt + 1;
z(cnt) = sqrt(a*((b*c)^d));
if abs(A-z(cnt))<tol
M % Turns out there is more than one!
z(cnt)
% break
end
end
.
.
.
Here is another method of finding these two locations:
g = gamma; % For convenience...
f = @(M) sqrt((1./M.^2).*(((2./(g+1)).*(1+(((g-1)./2).*M.^2))).^((g+1)/(g-1)))) - A;
rt1 = fzero(f,2)
rt2 = fzero(f,.5)

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by