Error :Nonscalar arrays of function handles are not allowed; use cell arrays instead.

2 次查看（过去 30 天）
Shimon Katzman
Shimon Katzman 2019-11-30
回答： Star Strider 2019-11-30
Hi everyone,
trying to plot the graph like in the picture but got Error using this code:
Es=200000; %Mpa
Esh=8500; %Mpa
fy=500; %Mpa
fsu=750; %Mpa
epssh=0.009;
epssu=0.075;
P=Esh*((epssu-epssh)/(fsu-fy));
epsy=fy/Es;
epscmv = linspace(0.1, 10, 500)*1E-3;
for i=1:numel(epscmv);
epscm = epscmv(i);
sigmaSteel(i)=@(epscm) Es*epscm .* (epscm<=epsy) + fy .* (epscm>epsy & epscm<=epssh) + fsu+(fy-fsu)*abs((epssu-epscm)./(epssu-epssh))^(1/P) .* (epscm>epssh & epscm<=epssu) + 0 .* (epscm>epsu);
end
plot(epscmv, sigmaSteel)
grid on
Thank you very much
  3 个评论
显示 1更早的评论
Shimon Katzman
Shimon Katzman 2019-11-30
Dear Sthephen.
1.Im new to Matlab.
2.Later I will need to use this function, this is just the beggining of the code that i am writing.
Shimon Katzman
Shimon Katzman 2019-11-30
Es=200000; %Mpa
Esh=8500; %Mpa
fy=500; %Mpa
fsu=750; %Mpa
epssh=0.009;
epssu=0.075;
P=Esh*((epssu-epssh)/(fsu-fy));
epsy=fy/Es;
epscmv = linspace(0, 10, 5000)*1E-3;
for i=1:numel(epscmv);
epscm = epscmv(i);
sigmaSteel(i)= Es*epscm .* (epscm<=epsy) +fy .* (epscm>epsy & epscm<=epssh) + fsu+(fy-fsu)*abs((epssu-epscm)./(epssu-epssh))^(1/P) .* (epscm>epssh & epscm<=epssu) + 0 .* (epscm>epssu);
end
plot(epscmv, sigmaSteel)
grid on

请先登录，再进行评论。

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

采纳的回答

Star Strider
Star Strider 2019-11-30
The loop is not necessary.
This required a few corrections (a typographical error and to vectorise the exponentiation), and now appears to work:
Es=200000; %Mpa
Esh=8500; %Mpa
fy=500; %Mpa
fsu=750; %Mpa
epssh=0.009;
epssu=0.075;
P=Esh*((epssu-epssh)/(fsu-fy));
epsy=fy/Es;
epscmv = linspace(0.1, 10, 500)*1E-3;
sigmaSteel=@(epscm) Es*epscm .* (epscm<=epsy) + fy .* (epscm>epsy & epscm<=epssh) + fsu+(fy-fsu)*abs((epssu-epscm)./(epssu-epssh)).^(1/P) .* (epscm>epssh & epscm<=epssu) + 0 .* (epscm>epssu);
figure
plot(epscmv, sigmaSteel(epscmv))
grid on
It does not exactly reproduce the plot in the 1.jpg attachment, although it appears to be reasonably close.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Graphics Object Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by