matlab如何使用龙格库塔,怎么用龙格库塔法

初值给一下。

在Matlab下输入:edit，然后将下面两行百分号之间的内容，复制进去，保存

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function dxdt=ode_Miss_ghost(t,x)

%分别用x(1),x(2),x(3),x(4)代替N1,P1,N2,P2

N1=x(1);

P1=x(2);

N2=x(3);

P2=x(4);

K=2;

tau_c=3e-9;

tan_p=6e-12;

beta =5e-5;

delta=0。

692;

eta =0。0001;

fm =8e6;

Ith =26e-3;

Ib =1。5*Ith;

Im =0。

3*Ith;

I1=Ib Im*sin(2*pi*fm*t) K*P2;

I2=Ib Im*sin(2*pi*fm*t) K*P1;

dxdt=[

(I1/Ith-N1-(N1-delta)/(1-delta)*P1)/tau_e;

((N1-delta)/(1-delta)*(1-eta*P1)*P1-P1 beta*N1)/tau_p;

(I2/Ith-N2-(N2-delta)/(1-delta)*P2)/tau_e;

((N2-delta)/(1-delta)*(1-eta*P2)*P2-P2 beta*N2)/tau_p;

];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

在Matlab下面输入:

t_start=0;

t_end=2e-9;

y0=[1e-3;1e-4;0;0]; %初值

[x,y]=ode15s('ode_Miss_ghost',[0,t_end],y0);

plot(x,y);

legend('N1','P1','N2','P2');

xlabel('x');

如果对您有帮助，请记得采纳为满意答案，谢谢!祝您生活愉快!

vaela。

全部

matlab如何使用龙格库塔,怎么用龙格库塔法相关推荐

matlab：二阶龙格库塔求解欧拉方程
%书籍:常用数值算法及其matlab实现 %第10章常微分方程初值问题的数值解法 %二阶龙格库塔方法 function S = heunsec(fun, x0, xn, y0, h) %fun:微分 ...
龙格库塔法matlab求解微分方程组,微分方程组的龙格库塔公式求解matlab版.pdf
微分方程组的龙格库塔公式求解matlab版微分方程组的龙格-库塔公式求解matlab版南京大学王寻 1. 一阶常微分方程组考虑方程组     y'f x,y,z , y x y ...
matlab 数值计算课二阶微分方程-龙格库塔方法 ODE45
详见mathworks 龙格库塔方法写成矩阵(状态方程)的形式更简洁一点(其实这两种方法结果是一样的,如果C是[1,0,0]的话,就很明显了) 例如:求系统在0-5s内的单位阶跃相应,已知传递函数: ...
比例导引 matlab,单步龙格库塔比例导引弹道计算matlab源程序
单步龙格库塔比例导引弹道计算matlab源程序本程序特别适合于弹道计算等方面使用.比ODE45函数的速度加快了很多,且程序充分展示了Matlab向量运算的强大功能,以及编程的简单快捷. 但愿本程序, ...
二阶龙格库塔公式推导_连续系统数值仿真方法——龙格库塔法
在决策理论与方法最后一课中,老师讲了连续系统建模与仿真,事实上在本科就已经接触过许多,在这里我把数值仿真方法中的龙格库塔法给详细介绍一下,这是个比较有趣,同时也比较实用的方法.在求解常微分方程初值问题 ...
2021-01-07 matlab数值分析常微分方程初边值问题数值解标准龙格库塔四阶四段公式欧拉法
matlab数值分析常微分方程初边值问题数值解标准龙格库塔四阶四段公式欧拉法 1.标准龙格库塔四阶四段公式 function y=rk4(fun,a,b,y0,n) h=(b-a)/n; y(1 ...
matlab:使用4阶龙格库塔方法求解常微分方程组
%书籍:常用数值算法及其matlab实现 %第10章常微分方程初值问题的数值解法,例10.14使用 %四阶龙格库塔方法 function [t,z] = rk4symeq(fun, t0, tf, ...
matlab 龙格-库塔法求解常微分方程
最近学习分室模型,里面碰到了用matlab 龙格-库塔法求解常微分方程研究了一阵子终于明白到底怎么实现了: 1. matlab 新建.m文件,编写龙格-库塔法求解函数 function [x,y] ...
四阶龙格库塔方程（Rungekutta）解二阶常微分方程组并计算船舶在迎浪下的纵摇埀荡耦合运动方程-附Matlab代码
今年年初的时候给师姐做了DDPG算法的船舶减横摇控制算法,师姐还有想法要让我把纵摇-埀荡两个自由度的减摇也做出来,这个任务归我了.实际上不管是多少个自由度的减摇,其实都需要解运动方程,当初做单自由度横 ...

matlab如何使用龙格库塔,怎么用龙格库塔法

matlab如何使用龙格库塔,怎么用龙格库塔法相关推荐

最新文章

热门文章