Plotting for various values of the same parameter on the same plot
显示 更早的评论
Asking here on how to plot "Y1(2,:)" for multiple "S" values on the same figure?. Consider values such as "S=0.1, 0.5, 1.1, 1.5" etc.
dYdX = @(X,Y) [Y(2); Y(3); S*Y(2).^2+Y(1).*Y(3)];
res = @(ya,yb) [ya(1)-1; ya(2); yb(2)];
SolYinit = bvpinit(linspace(0,5,50),[0; 0.2; 0.4]);
Fsol = bvp4c(dYdX, res, SolYinit);
X1 = Fsol.x;
Y1 = Fsol.y;
figure (1)
plot(X1, Y1(2,:), 'LineWidth', 2);
hold on
采纳的回答
Set ‘S’ as a vector, add ‘S’ as an additional parameter to the ‘dYdX’ function, then iterate through the values of ‘S’.
Try this —
Sv = [0.1, 0.5, 1.1, 1.5];
dYdX = @(X,Y,S) [Y(2); Y(3); S*Y(2).^2+Y(1).*Y(3)];
res = @(ya,yb) [ya(1)-1; ya(2); yb(2)];
SolYinit = bvpinit(linspace(0,5,50),[0; 0.2; 0.4]);
for k = 1:numel(Sv)
S = Sv(k);
Fsol = bvp4c(@(X,Y)dYdX(X,Y,S), res, SolYinit);
X1{k} = Fsol.x;
Y1{k} = Fsol.y;
end
figure (1)
Nrsp = numel(Sv);
for k = 1:Nrsp
subplot(Nrsp,1,k)
plot(X1{k}, Y1{k}(2,:), 'LineWidth', 2)
grid
title(sprintf('S = %.1f',Sv(k)))
end
I opted for the subplots because the scales vary so significantly on the solutions that the details in the solution for ‘S(1)’ gets lost if they are all plotted on the same axes.
.
4 个评论
Ke Le
2021-8-29
Star Strider
2021-8-29
The only reason I can think of that could have thrown that error is that you preallocated ‘X1’ and ‘Y1’ as matrices (rather than as cell arrays) prior to the loop, or have them defined as matrices (rather than as cell arrays) somewhere else earlier.
One way to avoid the error is to call them something else, such as ‘Xc’ and ‘Yc’ in the loop (and in the following plot loop), again providing that those variables have also not been previously defined as something else earlier in the code.
There could be other reasons, however that is the only situation that I have ever encountered a similar error in my experience. In any event, giving them unique names will likely solve it.
.
Ke Le
2021-8-30
Star Strider
2021-8-30
My pleasure!
If my Answer helped you solve your problem, please Accept it!
.
更多回答（1 个）
% Some data
x = 1:10;
% Create the figure window, use the hold command
% to prevent overwriting of axes
figure
hold on
% Plot line with different slope parameters
for s = 1:3
plot(x,s*x)
end
类别
在 帮助中心 和 File Exchange 中查找有关 Function Creation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
