Graphing a variable function with boundaries problem

조회 수: 1 (최근 30일)
Mohsen Nashel
Mohsen Nashel 2019년 9월 24일
댓글: Mohsen Nashel 2019년 9월 24일
i am trying to plot the function that is shown in the picture. Matlab keep giving me the error "Error: File: Q4.m Line: 11 Column: 64
Invalid expression. When calling a function or indexing a variable, use parentheses. Otherwise, check for
mismatched delimiters."
Refers to P_0_1 in the loop Screen Shot 2019-09-24 at 12.52.30 PM.png
%Part a)
R=287; % Specific gas constant in J/kg.k
M=[2:0.1:6];
gamma=1.4;
P_0_1=101325; % Pressure at sea level in Pa unit
for n=1:length(M);
P_0(n,1)=((1+((gamma-1))/2).*M(n).^2).^(gamma/(gamma-1)).*P_0_1);
end
P_0_ratio=P_0./P_0_1;
figure (1)
plot (M, P_0_ratio)
xlabel('Moch Number')
ylabel('Total stagnation pressure ratio P_0/P_0_1')
title('Total stagnation pressure ratio P_0/P_0_1 vs. Moch Number')
%Part b)
for n=1:length(M);
delta_s_r(n,1)=-log(P_0_ratio(n));
end
figure(2)
plot (M,delta_s_r)
xlabel('Moch Number')
ylabel('Dimensionless entropy change \delta s')
title('Dimensionless entropy change \delta s vs. Moch Number')

채택된 답변

Ruger28
Ruger28 2019년 9월 24일
편집: Ruger28 2019년 9월 24일
Removing the extra " ) " at the end of that line, yields plots and no errors.
%Part a)
R=287; % Specific gas constant in J/kg.k
M=[2:0.1:6];
gamma=1.4;
P_0_1=101325; % Pressure at sea level in Pa unit
for n=1:length(M);
P_0(n,1)=((1+((gamma-1))/2).*M(n).^2).^(gamma/(gamma-1)).*P_0_1;
end
P_0_ratio=P_0./P_0_1;
figure (1)
plot (M, P_0_ratio)
xlabel('Moch Number')
ylabel('Total stagnation pressure ratio P_0/P_0_1')
title('Total stagnation pressure ratio P_0/P_0_1 vs. Moch Number')
%Part b)
for n=1:length(M);
delta_s_r(n,1)=-log(P_0_ratio(n));
end
figure(2)
plot (M,delta_s_r)
xlabel('Moch Number')
ylabel('Dimensionless entropy change \delta s')
title('Dimensionless entropy change \delta s vs. Moch Number')
  댓글 수: 1
Mohsen Nashel
Mohsen Nashel 2019년 9월 24일
Thank you! Rushing to finish homework can't figure out anything! You are a scholar!!

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

추가 답변 (0개)

Community Treasure Hunt

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

Start Hunting!

Translated by