error occurs in the present code
显示 更早的评论
function main
clc
clear all
x=3;
options=optimset('Display','iter');
x1=fsolve(@solver,x,options);
function F=solver(x)
options=odeset('RelTol',le-8,'AbsTol',[le-8, le-8, le-8]);
[t,u]=ode45(@equation,[0,20],[4 -1 x],options);
s=length(t);
F=u(s,2);
figure(1)
plot(t,u(:,2))
hold on
end
end
function dy=equation(t,y)
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=y(3);
dy(3)=y(2)^2-y(1)*y(3);
end
After running the above code following error occurs:
Not enough input arguments.
Error in fsolve (line 230)
fuser = feval(funfcn{3},x,varargin{:});
Error in (line 10)
x1=fsolve(@solver,x,options);
Caused by:
Failure in initial objective function evaluation. FSOLVE cannot continue.
采纳的回答
Star Strider
2019-2-23
You have a typographical error in your odeset call. You typed ‘l’ (lower-case ‘L’) instead of the number 1.
Use this instead:
options=odeset('RelTol',1e-8,'AbsTol',1e-8);
and your code runs without error.
When I ran it,, the result was:
x1 =
3.732
7 个评论
MINATI
2019-2-23
Star Strider
2019-2-23
I am not certain what you want. If you want to plot a single curve (with the optimised value for the third initial condition), just plot the last one.
Try this:
x=3;
options=optimset('Display','iter');
x1=fsolve(@solver,x,options)
function F=solver(x)
options=odeset('RelTol',1e-8,'AbsTol',1e-8);
[t,u]=ode45(@equation,[0,20],[4 -1 x],options);
s=length(t);
F=u(s,2);
% figure(1)
% plot(t,u(:,2))
% hold on
% end
% end
function dy=equation(t,y)
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=y(3);
dy(3)=y(2)^2-y(1)*y(3);
end
end
function dy=equation(t,y)
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=y(3);
dy(3)=y(2)^2-y(1)*y(3);
end
[t,u]=ode45(@equation,[0,20],[4 -1 x1],options);
figure(1)
plot(t,u(:,2))
Esperiment to get the result you want.
MINATI
2019-2-23
Star Strider
2019-2-23
As always, my pleasure.
MINATI
2019-2-24
Star Strider
2019-2-24
It is one of the initial conditions for your ‘equation’ function differential equation system.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
