用MATLAB函数绘制系统的,用matlab绘制odes系统定义的函数

我试图用Matlab求解一阶ODE系统,

$$\left\{

\begin{array}{l}

x_{1}^{\prime }=-\frac{1}{t+1}x_{1}+x_{2} \\

x_{2}^{\prime }=-(1+e^{-2t})x_{1}-\frac{1}{t+1}x_{2}+\frac{e^{-3t}}{t+1}x_{3}

x_{3}^{\prime }=-\frac{1}{t+1}x_{3}+x_{4} \\

x_{4}^{\prime }=e^{-3t}\left( t+1\right) x_{1}-\left( 1+e^{-2t}\right) x_{3}-%

\frac{1}{t+1}x_{4}-\frac{1}{t+1}x_{3}^{2}%

\end{array}%

\right. $$

我已经定义了函数:

function dzdt=odefun(t,z)

dzdt=zeros(4,1);

dzdt(1)=-(1/(t+1))*z(1)+z(2);

dzdt(2)=-(1+exp(-2*t))*z(1)-(1/(t+1))*z(2)+(exp(-3*t))/(t+1)*z(3);

dzdt(3)=z(4)-(1/(t+1))*z(3);

dzdt(4)=(exp(-3*t))*(t+1)*z(1)-(1+exp(-2*t))*z(3)-(1/(t+1))*z(4)-(1/(t+1))*z(3)^2;

end

时间间隔是

[0,100]

初始条件是

z0 = [0.01 0.01 0.01 0.01]

与

ode45

解算器,我使用了以下命令:

>> tspan = [0 100];

>> z0 = [0.01 0.01 0.01 0.01];

>> [t,z] = ode45(@(t,z) odefun(t,z), tspan, z0);

>> plot(t,z(:,1),'r')

我很容易得到

z(1)=x_1

但我想画出函数

f(t)=(t+1)*x_1(t)

t\in [0,100]

在哪里

x_1=z(1)

是第一个未知的系统。我怎么能这样做?

用MATLAB函数绘制系统的,用matlab绘制odes系统定义的函数相关推荐

python中用函数编写程序_python编写程序,在程序中定义一个函数,计算1+1/2+1/3+1/4+……+1/n...
python语言编写的程序需要运行库.Python(英语发音:/?pa?θ?n/),是一种面向对象.解释型计算机程序设计语言,由Guido van Rossum于1989年发明,第一个公开发行版发行于 ...
matlab中m文件是什么,MATLAB中M文件的使用
运行方法: 使用matlab的时候,可以在"Command Window"内直接书写matlab代码,也可以将代码保存到M文件中,然后运行该文件.使用matlab主界面菜单&quo ...
在类内定义成员函数、在类外定义成员函数、计算长方体的体积【C++面向对象编程类的使用经典案例】
文章目录一.在类内定义成员函数二.在类外定义成员函数(使用符号::) 三.计算3个长方体的体积(class成员函数) 一.在类内定义成员函数 #include <iostream> u ...
python中函数的定义_Python函数是什么_如何定义和调用函数?
函数是什么? 本节Python培训教程是要大家认识一个"老朋友"--函数. 与数学中的函数不同,在Python中,函数不是看上去冰冷无聊的规则和公式,而是实打实的.有自己作用的代码 ...
mysql调用java函数_mysql概要（十五）存储过程(2)(java调用函数和存储过程）
使用存储过程的目的: 即将数据的函数处理放在数据库,这样主要是为了,多语言跨平台的使用(java,.net,php),不用重复的函数代码,直接调用存储过程. (对于使用存储过程还是代码来处理数据的效率 ...
C++如何定义一个函数指针
首先给出一个函数 int function(int a) {return a; } 一.定义一个函数指针 int (*fun)(int a); // 赋值 fun = function; // 使用 ...
matlab绘制系统函数频率特性,实验七用MATLAB进行系统频率特性曲线绘制
实验七用MATLAB进行系统频率特性曲线绘制自动控制原理课程验证性实验报告实验名称用MATLAB进行系统频率特性曲线绘制实验时间 2013年 05月30日学生姓名实验地点 070312 ...
matlab状态空间法算反馈阵,matlab中已知系统的状态方程怎样绘制系统阶跃响应曲线...
matlab中,怎样设置坐标系统的原点在左上角? axisijx=-pi:0.1:pi;y=sin(x);plot(x,y,'g:'); 谁能帮我做Matlab的作业:分析系统能否通过状态反馈实现镇定 ...
matlab 绘制符号函数,DAY8 MATLAB学习笔记—simulink入门、MATLAB符号函数的图形绘制...
如何打开simulink: 启动simulink: 先打开MATLAB软件界面第一步打开simulink 第二步在command windows输入 simulink然后enter,等待有很多模块 ...

用MATLAB函数绘制系统的,用matlab绘制odes系统定义的函数

用MATLAB函数绘制系统的,用matlab绘制odes系统定义的函数相关推荐

最新文章

热门文章