Unable to meet integration tolerances without reducing the step size problem.

5 次查看（过去 30 天）
Md. Golam Zakaria
评论： Torsten 2023-2-7
Hello Everyone.
I am facing 'Warning: Failure at t=0.000000e+00. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (7.905050e-323) at time t, error for the following code. How can I fix it.
clc
clear all
A =1.4036e-18;
B = 3.4302e-18;
e=1.60217663e-19;
T_e=3;
f = @(t,x) [(sqrt(A/x(1)))*(-x(4)-(0.424)+(1.6977)*exp(-(e/T_e)*(x(1))));
(sqrt(A/x(1)))*(x(4)-(0.663)+(2.6540)*exp(-(e/T_e)*(x(2))));
-(2.994*10^5)*x(4);
((x(1))-(x(2))+x(3)+60*(cos(t)))-2.55*10^7*x(4)];
[t,xa] = ode15s(f,[0, 80],[0,0,-62,0]);
figure(1)
plot(t,xa(:,1),'Color','blue'),grid on;
title('xa1')
xlabel('t'), ylabel('Qsp (Powered Electrode Sheath Voltage)')

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

回答（1 个）

Torsten
Torsten 2023-2-7
You define an initial condition for x(1) as 0, but you divide by x(1) in equations 1 and 2.
This won't work.
  3 个评论
显示 1更早的评论
Md. Golam Zakaria
Md. Golam Zakaria 2023-2-7
@Steven Lord, If I chose nonzero initial conditions will it solve my problem. I didnt do it because It takes too much time to solve
Torsten
Torsten 2023-2-7
If I chose nonzero initial conditions will it solve my problem.
Nobody knows. Just try it.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by