用ode求解二阶非线性微分方程

输入参数数目不足求大神帮助

clear all;
close all;
clc;
f0=2*10^4;
R0=6*10^(-6);
P0=1.013*10^5;
time=1/f0;
Pa=1.2*10^5;
Ki=1.67;
rho=1000;
sigma=0.0725;
mu=0.001;
Pi=3.1415926;
c=1481;
a=R0/10;
w=2*Pi*f0;
P1=1*10^5;
P2=1*10^5;
P3=1*10^5;
P4=1*10^5;
step=time/30000;
ttime=1*time;
ts=0:step:ttime;
y0=[R0,0];
[t,y]=ode45(fun5,ts,y0);
plot(t/time,y(:,1)/R0)
function f=fun5(t,y)
f0=2*10^4;
R0=6*10^(-6);
P0=1.013*10^5;
time=1/f0;
Pa=1.2*10^5;
Ki=1.67;
rho=1000;
sigma=0.0725;
mu=0.001;
Pi=3.1415926;
c=1481;
a=R0/10;
w=2*Pi*f0;
f=[y(2);((P0+2*sigma/R0)*(R0/y(1))^(3*Ki)-2*sigma/y(1)-4*mu*y(2)/y(1)-(P0-Pa*sin(w*t)))/rho/y(1)-(3/2*(y(2)^2))/y(1)-3*(P0+2*sigma/R0)*Ki*R0^(3*Ki)*y(2)*y(1)^2/(rho*c*(y(1)^3-a^3)^(Ki+1))+(Pa*w*cos(w*t))/(rho*c)];
end

用ode求解二阶非线性微分方程相关推荐

matlab二阶非线性微分方程求解,二阶非线性常微分方程的打靶法.ppt
二阶非线性常微分方程的打靶法二阶非线性常微分方程的打靶法计算思路主要分为以下五步: 给定容许误差ε,迭代初始值γ1,对k=1,2,...做: (1)用四阶Runge-Kutta 方法求解初值问, ...
matlab求解f非线性微分方程数值解,非线性﹑微分方程数值求解.PPT
非线性﹑微分方程数值求解数学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 计算方法(B) 主 ...
Petsc求解一维非线性微分方程，SNES对象-2
微分方程引入 − u ′ ′ + ρ ( u ) = f , x ∈ ( 0 , 1 ) -u^{''}+\rho\sqrt(u)=f,x \in(0,1) −
Matlab通过ode求解微分方程
MATLAB有很多用于求解微分方程的内置功能.MATLAB包含了用于求解常微分方程(ODE)的函数高阶微分方程必须重新表述为一个一阶系统微分方程. 并不是所有的微分方程都可以用同样的方法求解,所以M ...
关于求解微分方程——初学Matlab里的 ODE求解器
学习背景最近想挖掘一下自己项目的理论深度,于是找到了老师.在老师的建议下,我们开始了漫长的研读老师的论文的旅程(论文名:Optimal Design of Adaptive Robust Contr ...
二阶边值问题的数值解matlab,二阶线性微分方程边值问题的MATLAB求解
2012 年 3 月第 26 卷第 1 期阴山学刊 YINSHAN ACADEMIC JOURNAL Mar． 2012 Vo1． 26 No． 1 二阶线性微分方程边值问题的 MATLAB ...
matlab解方程x 2-x-2=0,matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0...
问题描述: matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0 1个回答分类: 数学 2014-11-28 问题解答: 我来补答 ...
用MATLAB求解非线性微分方程
总结一下MATLAB中求解微分方程的思路和步骤.固然,网上很多关于此类的技术型文章,但往往一看下来发现,文章中的友情链接比文章字数还多,要了解这一篇文章,你要先了解那个:要了解那个,你又要了解那个那个 ...
二阶常系数线性微分方程
二阶常系数线性微分方程一般形式 y'' +p y' + qy = f(x) ① (下面用到r ...

用ode求解二阶非线性微分方程

用ode求解二阶非线性微分方程相关推荐

最新文章

热门文章