初值为1,X'=-0.75x(t-1)时滞区间为-1到0的一阶时滞微分方程的MATLAB三维图代码
时间: 2024-03-25 10:35:25 浏览: 61
假设时间步长为0.01,时滞区间为-1到0,则MATLAB代码如下:
```
t = -1:0.01:10; % 定义时间范围
x = zeros(size(t)); % 初始化x值
x(t>=0) = 1; % 给时刻t>=0的x赋值为1
for i = 2:length(t)
x(i) = x(i-1) - 0.75*x(max(i-101,1));
end
plot3(t,x,zeros(size(t))); % 画出三维图
xlabel('t');
ylabel('x');
zlabel('x''');
```
其中,`max(i-101,1)`表示i-101和1中的较大值,这是因为时滞区间为-1到0,所以最早需要取到t=0时刻的x值,而时滞长度为100(即100/0.01=10000个时间步长),因此需要用到i-101。运行后可以得到一个以时间t为横轴,x为纵轴,x'为纵轴的三维图,表示该时滞微分方程的解。
相关问题
初值为y(n)=sin(10n)其中n为-1到0,X'=-X+tanh(10X(t-1))的一阶时滞微分方程的MATLAB的二维图
以下是MATLAB代码和图形:
```matlab
% 定义初值和时间步长
n = -1:0;
y = sin(10*n);
h = 0.01;
% 定义时滞微分方程
f = @(t, x) -x + tanh(10*x(t-1));
% 使用ode45求解微分方程
[t, x] = ode45(f, n, y);
% 绘制图形
plot(n, y, 'o-', t, x, '-')
xlabel('n')
ylabel('y(n) and x(n)')
legend('y(n)', 'x(n)')
```
图形如下:
![二维图](https://img-blog.csdnimg.cn/20210925132523627.png)
阅读全文
相关推荐
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![pdf](https://img-home.csdnimg.cn/images/20241231044930.png)
![doc](https://img-home.csdnimg.cn/images/20241231044833.png)
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![doc](https://img-home.csdnimg.cn/images/20241231044833.png)
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![doc](https://img-home.csdnimg.cn/images/20241231044833.png)
![.zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![pdf](https://img-home.csdnimg.cn/images/20241231044930.png)
![zip](https://img-home.csdnimg.cn/images/20241231045053.png)
![doc](https://img-home.csdnimg.cn/images/20241231044833.png)