清空筛选条件
清空筛选条件

Differential equation error.

2 次查看（过去 30 天）
Reji G
Reji G 2022-10-18
评论： Walter Roberson 2022-10-19
The following equation need to be implemented and the graph should be obtained like this. I tried all the way, but couldn't.
clc;
close all;
clear all;
t=0:0.001:10;
I = @(t)100*sin(100*pi*t);
r0=0.0002;
Lp=1;
U0=(0.95+(0.005./I(t)));
Ust=(U0+(r0*abs(I(t))))*Lp;
G=(abs(I(t)))./Ust;
g0 = 0.1;
fun = @(t,g) (1/0.00001)*(G-g);
[T,g] = ode45(@(t,g) fun(t,g,I),0.01:0.001:5,g0);
plot(T,g);

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

回答（1 个）

Walter Roberson
Walter Roberson 2022-10-18
U0=@(t)(0.95+(0.005./I(t)));
Ust=@(t)(U0(t)+(r0*abs(I(t))))*Lp;
G=@(t)(abs(I(t)))./Ust(t);
g0 = 0.1;
fun = @(t,g) (1/0.00001)*(G(t)-g);
However using abs() violates the restriction that the function you use has continuous second derivatives.
Please see the ballode example.
  3 个评论
显示 1更早的评论
Reji G
Reji G 2022-10-19
I'm getting 5 errors. I don't know how to rectify it. But I need the graph as plotted above. Please edit the code to run and get the plot.
Walter Roberson
Walter Roberson 2022-10-19
The code needs bigger changes. The function violates the continuity requirements of all ode*() solvers

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 2-D and 3-D Plots 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by