Matlab求解常微分方程组

求解这个常微分方程组。

$\left\{ \begin{array}{l} \frac{d}{{dt}}(P{R^5}) = a{P^2}{R^5}\\ \\ \frac{{{d^2}R}}{{d{t^2}}} = {R^2P} \end{array} \right.$

初始条件为

$R(0) = \frac{1}{\varepsilon }$ $\frac { d R } { d t } ( 0 ) = - 1$ $\frac { d ^ { 2 } R } { d t ^ { 2 } } ( 0 ) = \varepsilon ^ { 3 }$ $P ( 0 ) = \varepsilon ^ { 5 } \$

其中ε取0.01，a是有上限的参数，求解方程的目的其实是找出a的临界值。

syms y(t)
for i = [0:0.5:1.5,1.763]   %1.763是临界值eqn = diff(diff(y,2)*y^3) == i*diff(y,2)^2*y;%化三阶方程为一阶V = odeToVectorField(eqn);M = matlabFunction(V,'vars',{'t','Y'});e = 0.01;% 坐标范围和初始条件interval = [0 2/e];yInit = [1/e -1 e^3];ySol = ode45(M,interval,yInit);tValues =  linspace(0,2/e,1000);yValues = deval(ySol,tValues,1);pValues = deval(ySol,tValues,3)./yValues.^2;vValues = yValues.^3;plot(tValues,yValues)hold on
endylim([-10 100])
xlabel('$\hat{t}$','interpreter','latex');
ylabel('$\hat{R}$','interpreter','latex');
text(60,80,'\epsilon=0.01','FontName','Times New Roman','FontSize',15)
lgd = legend(num2str(0),num2str(0.5),num2str(1),num2str(1.5),num2str(1.763));
title(lgd,'\xi_n_o_\alpha')
set(gca,'LineWidth',2)
set(findall(gcf,'type','line'),'linewidth',2)
set(gca,'FontName','Times New Roman','FontSize',15)

Matlab求解常微分方程组相关推荐

matlab求解一阶微分方程组
在matlab里运用ode计算一阶微分方程组可以根据需求选择合适的ode 举个例程序为:function dy = odefun(t,y) dy = zeros(3,1); % 一个列向量 dy( ...
matlab:使用4阶龙格库塔方法求解常微分方程组
%书籍:常用数值算法及其matlab实现 %第10章常微分方程初值问题的数值解法,例10.14使用 %四阶龙格库塔方法 function [t,z] = rk4symeq(fun, t0, tf, ...
matlab求解常微分方程,matlab 求解常微分方程式
MATLAB解常微分方程式的语法是dsolve('equation','condition'),其中equation代表常微分方程式即y'=g(x,y),且须以Dy代表一阶微分项y' D2y代表二阶微 ...
matlab解二阶微分方程组ode,MATLAB解含参数方程、矩阵方程、二阶微分方程组
1.如下公式: 其中x.y.m.n为参数,a.b为未知数,利用MATLAB求解方程: syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y= ...
matlab常系数线性矩阵微分方程组,基于Matlab常系数线性微分方程组的求解
·基础数学· 基于 Matlab 常系数线性微分方程组的求解* 严水仙 (赣南师范大学数学与计算机科学学院,江西赣州 341000) 摘要: 在常微分方程课程教学中,常系数线性微分方程组可以通过 ...
Matlab：Matlab编程语言应用之数学计算(求极限/渐近线求导数常微分方程求解求微分方程组的解求临界阻尼系数的解)的简介、案例实现之详细攻略
Matlab:Matlab编程语言应用之数学计算(求极限/渐近线&求导数&常微分方程求解&求微分方程组的解&求临界阻尼系数的解)的简介.案例实现之详细攻略目录三.极 ...
8 精通matlab入门求微分方程组的通解特解数值解看这篇文章就够了
%% 学习目标:使用MATLAB求解常微分方程通解特解数值解 %% 求微分方程的通解 clear all; f1=dsolve('Dy-y=sin(x)') %默认情况下以t为自变量求解求出的 ...
一阶微分方程的物理意义_MIT—微分方程笔记24 一阶常微分方程组
18.03 微分方程 Differential Equations 第四单元一阶常微分方程组 Unit 4 First-Order Systems 第24讲一阶常微分方程组第25讲常系数齐次线 ...
高等工程数学--求解非齐次线性微分方程组
今天终于弄懂了怎么求解非齐次线性微分方程组,大家都要学会哦,非常重要的,不懂的留言哦! 嘻嘻第一张图片的标题写错了,实在抱歉,是非齐次些许抽象,来看例题吧刚开始卡在了积分那里,后来灵光一现,会了 ...

Matlab求解常微分方程组

Matlab求解常微分方程组相关推荐

最新文章

热门文章