Plotting a polynomial problem

Question

Josie 2014-7-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/144245-plotting-a-polynomial-problem

回答： Star Strider 2014-7-31

Trying to plot this polynomial for different values of A but no luck so far.

M=linspace(-10,10,100);
g=10;
for A=0.1:0.1:0.9
    D=-(A^2-A+1/3)*(((A-1/2)*M*g-0.25*A^2)/(2*M*((A-1/2)^2+1)+0.5*(A^2-A+1/3)))-0.5*A^2
    plot(M,D)
    hold all
end
legend('A=0.1','A=0.2','A=0.3','A=0.4','A=0.5','A=0.6','A=0.7','A=0.8','A=0.9')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Image Analyst 2014-7-31

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/144245-plotting-a-polynomial-problem#answer_147023

在 MATLAB Online 中打开

You forgot to use .* instead of *.

Corrected code:

workspace;
clc;
M=linspace(-10,10,100);
g=10;
lineColors = lines(9);
counter = 1;
for A=0.1:0.1:0.9
  D = -(A^2-A+1/3) .* (((A-1/2)*M*g-0.25.*A^2) ./ (2*M*((A-1/2)^2+1)+0.5*(A^2-A+1/3)))-0.5*A^2
  plot(M, D, 'LineWidth', 2, 'Color', lineColors(counter, :));
  hold on
  counter = counter + 1;
end
legend('A=0.1','A=0.2','A=0.3','A=0.4','A=0.5','A=0.6','A=0.7','A=0.8','A=0.9')
grid on;
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Star Strider 2014-7-31

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/144245-plotting-a-polynomial-problem#answer_147024

在 MATLAB Online 中打开

I don’t understand what ‘no luck so far’ means (I can only assume that a plot of several parallel lines of dots wasn’t your desired result), but after ‘vectorising’ your code (the line that assigns ‘D’) it produces a much different plot than the one your code originally produced:

M=linspace(-10,10,100);
g=10;
for A=0.1:0.1:0.9
    D=-(A^2-A+1/3).*(((A-1/2).*M.*g-0.25.*A.^2)./(2*M.*((A-1/2).^2+1)+0.5.*(A.^2-A+1/3)))-0.5*A.^2;
    plot(M,D)
    hold all
end
legend('A=0.1','A=0.2','A=0.3','A=0.4','A=0.5','A=0.6','A=0.7','A=0.8','A=0.9')

‘Vectorising’ produces element-by-element operations, rather than MATLAB’s default matrix operations. It requires using the ‘dot operator’ in front of multiplication, division, and exponentiation operators: (.*), (./), (.^).

(The only exception is the transpose operator ('), where (.') produces a simple transpose rather than the default complex conjugate transpose. I mention that for completeness. It doesn’t relate to your code.)