此微分方程组可用ode45()函数求得其数值解 。
首先,建立自定义函数
function dy = rigid(t,y)
dy = zeros(3,1);
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
end
其二,用ode45()函数求出x、y值
[x,y] = ode45(@rigid,[0 12],[0 1 1]);
其三,用plot()函数绘出x—y,x—dy/dx,x—d2y/dx2曲线图
plot(x,y(:,1),'-',x,y(:,2),'-.',x,y(:,3),'.')
运行结果
文章插图
【写出求解下列微分方程模型的matlab程序。】
推荐阅读
- 清朝小说家吴敬梓为什么写出了儒林外史 儒林外史作者是谁
- 如何写出热门短视频标题?热门短视频标题模板
- 二阶矩阵求逆,用伴随矩阵方法。有点懵,求解答
- 广汽汇理汽车金融怎么样靠谱吗?急求解答!!
- 求解析 高二地理会考
- 天珠变中都哪几章有推倒??求大神求解。
- 中政行测怎么样啊 ,,求解。。。
- 沙姜与生姜有什么不同?求解
- 安装CAD2012出现C++ 2005错误 求解
- 句号打不出来了但是书名号却能打出来 这是怎么回事啊 求解啊同样都是一个按键 按键也没坏啊! 求大神帮忙