0

我需要帮助绘制一个微分方程......它总是很时髦，而且图表不是它应该看起来的样子。

function [dydt] = diff(y,t)

dydt = (-3*y)+(t*(exp(-3*t)));

end

tI = 0;
yI = -0.1;
tEnd = 5;
dt = 0.5;

t = tI:dt:tEnd;
y = zeros(size(t));
y(1) = yI;

for k = 2:numel(y)
    yPrime = diff(t(k-1),y(k-1));
    y(k) = y(k-1) + dt*yPrime;
end

plot(t,y)
grid on
title('Engr')
xlabel('Time')
ylabel('y(t)')
legend(['dt = ' num2str(dt)])

那是我的代码，但图表与它应该的样子完全不同。我是否缺少for声明索引之类的东西？

编辑

我收到一个错误：

Error using diff
Difference order N must be a positive integer scalar.

Error in diff3 (line 12)
    yPrime = diff(t(k-1),y(k-1));
4

1 回答 1

1

在修复了 Danil Asotsky 和 ​​horchler 在评论中指出的错误之后：

  1. 避免与内置函数 'diff' 的名称冲突
  2. 将参数的顺序更改为t,y.
  3. 将时间步长减小dt到 0.1
  4. 将 ODE 右侧转换为匿名函数

（并删除函数定义中不必要的括号），您的代码可能如下所示：

F = @(t,y) -3*y+t*exp(-3*t);

tI = 0;
yI = -0.1;
tEnd = 5;
dt = 0.1;

t = tI:dt:tEnd;
y = zeros(size(t));
y(1) = yI;

for k = 2:numel(y)
    yPrime = F(t(k-1),y(k-1));
    y(k) = y(k-1) + dt*yPrime;
end

plot(t,y)
grid on
title('Engr')
xlabel('Time')
ylabel('y(t)')
legend(['dt = ' num2str(dt)])

它按预期执行：

输出

于 2015-01-05T21:13:19.937 回答