How to plot a (which is changing from 0-1 in 0.01 increments) vs x(2) (using a for loop and fsolve to find the solution of a nonlinear equation containing x(s) sol based on a

Question

Victor Jimenez Carrillo 2021-9-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1454239-how-to-plot-a-which-is-changing-from-0-1-in-0-01-increments-vs-x-2-using-a-for-loop-and-fsolve-t

评论： Alan Stevens 2021-9-16

V=1000; Q=50; Ca0=1; k=1;

for a=0:1:0.1

f=@(x) [Q*Ca0-Q*x(1)-k*x(1)^2*(a*V); Q*x(1)-Q*x(2)-k*x(2)^2*(1-a)*V];

fsolve(f,[0.5,0.5])

end

plot(a,x(2))

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

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2021-9-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1454239-how-to-plot-a-which-is-changing-from-0-1-in-0-01-increments-vs-x-2-using-a-for-loop-and-fsolve-t#answer_788554

在 MATLAB Online 中打开

Your first equation is a simple quadratic in x(1); your second is a quadratic in x(2) that depends on x(1), so, assuming you are only interested in the positive roots, these can be solved as follows:

V=1000; Q=50; Ca0=1; k=1;

a = 0:0.01:1;

x1 = zeros(1,numel(a));

x2 = zeros(1,numel(a));

for i=1:numel(a)

% assuming you want positive values of x1 and x2

if a(i) == 0

A = k*V;

x1(i) = Ca0;

x2(i) = (-Q + sqrt(Q^2 + 4*A*Q*x1(i)))/(2*A);

elseif a(i) == 1

A = k*V;

x1(i) = (-Q + sqrt(Q^2 + 4*A*Q*Ca0))/(2*A);

x2(i) = x1(i);

else

A1 = k*a(i)*V;

A2 = k*(1-a(i))*V;

x1(i) = (-Q + sqrt(Q^2 + 4*A1*Q*Ca0))/(2*A1);

x2(i) = (-Q + sqrt(Q^2 + 4*A2*Q*x1(i)))/(2*A2);

end

subplot(2,1,1)

plot(a,x1),grid

xlabel('a'),ylabel('x1')

subplot(2,1,2)

plot(a,x2),grid

xlabel('a'),ylabel('x2')

% The two equations can be expressed as:

% k*a*V*x1^2 + Q*x1 - Q*Ca0 = 0

% k*(1-a)*V*x2^2 + Q*x2 - Q*x1 = 0

2 个评论
显示无隐藏无

Victor Jimenez Carrillo 2021-9-16

This code is amazing, but I was wondering if I could use Fsolve becaue after plotting the graphs for the 2 simple ^2 values. I'd like to change the values to ^3. I was wondering if I could still use this code. That was the whole reason I wanted to use Fsolve. I aprreciate your help sir.