odeFunction matlab,MATLAB 如何在function裡使用ode45輸出值

最新推荐文章于 2023-02-18 20:18:21 发布
最新推荐文章于 2023-02-18 20:18:21 发布
阅读量312 收藏
点赞数
文章标签： odeFunction matlab

clear;

M=0.024;

a=1/1.5;

b=1;

alpha=0;

beta=0;

maxiter=100;

for i=1:3

C=1;

if i==1

C=C-10^-3;

elseif i==2

C=C+10^-3;

elseif i==3

C=C

end

aa=0

epsilon=1e-3;

x0=0;

h=(b-a)/maxiter

t0=a

s(1)=(beta-alpha)/(b-a);

t=zeros(maxiter+1,1);

t=a:h:b;

x1(1)=0

x2(1)=s(1)

z1(1)=0

z2(1)=1

for ii=1:maxiter

H=0.25;

m=1.2;

Y=x1(ii);

Y1=x2(ii);

tt=t(ii);

sigmay=((Y/tt)^2-Y*Y1/tt+(Y1)^2)^0.5;

if sigmay>1

aa=(sigmay.^m-1)/(H*sigmay)

elseif sigmay<=1

aa=0

end

f=@(t,x1,x2)x2  %f代表y

g=@(t,x1,x2)1/(t*(1+aa))*C+1/(t.^2)*x1-1/t.*x2  %g代表y'

F=@(t,z1,z2)z2;%F代表z

G=@(t,z1,z2)1/(t.^2)*z1-1/t*z2;%G代表z'

f0=x1(ii);

g0=x2(ii);

F0=z1(ii);

G0=z2(ii);

f1=f(t0,f0,g0)

g1=g(t0,f0,g0)

F1=F(t0,F0,G0)

G1=G(t0,F0,G0)

f2=f(t0+h/2,f0+h/2*f1,g0+h/2*g1)

g2=g(t0+h/2,f0+h/2*f1,g0+h/2*g1)

F2=F(t0+h/2,F0+h/2*F1,G0+h/2*G1)

G2=G(t0+h/2,F0+h/2*F1,G0+h/2*G1)

f3=f(t0+h/2,f0+h/2*f2,g0+h/2*g2)

g3=g(t0+h/2,f0+h/2*f2,g0+h/2*g2)

F3=F(t0+h/2,F0+h/2*F2,G0+h/2*G2)

G3=G(t0+h/2,F0+h/2*F2,G0+h/2*G2)

f4=f(t0+h,f0+h*f3,g0+h*g3)

g4=g(t0+h,f0+h*f3,g0+h*g3)

F4=F(t0+h,F0+h*F3,G0+h*G3)

G4=G(t0+h,F0+h*F3,G0+h*G3)

x1(ii+1)=x1(ii)+h/6*(f1+2*f2+2*f3+f4)

x2(ii+1)=x2(ii)+h/6*(g1+2*g2+2*g3+g4)    %龙格库塔法迭代

z1(ii+1)=z1(ii)+h/6*(F1+2*F2+2*F3+F4)

z2(ii+1)=z2(ii)+h/6*(G1+2*G2+2*G3+G4)

end

% iterations:

%

flag=0;

n=length(x1(1,:));

diff=abs(x1(1,n)-beta);

if diff

flag=1;

end

k=1;

while flag==0,

s(k+1)=s(k)+(beta-x1(1,n))/z1(1,n);

clear t;

clear x1;

clear x2;

x1(1)=0

x2(1)=s(k+1)

z1(1)=0

z2(1)=1

t=zeros(maxiter+1,1);

t=a:h:b;

for ii=1:maxiter

H=0.25;

m=1.2;

Y=x1(ii);

Y1=x2(ii);

tt=t(ii);

sigmay=((Y/tt)^2-Y*Y1/tt+(Y1)^2)^0.5;

if sigmay>1

aa=(sigmay.^m-1)/(H*sigmay)

elseif sigmay<=1

aa=0

end

f=@(t,x1,x2)x2  %f代表y

g=@(t,x1,x2)1/(t*(1+aa))*C+1/(t.^2)*x1-1/t.*x2  %g代表y'

F=@(t,z1,z2)z2;

G=@(t,z1,z2)1/(t.^2)*z1-1/t*z2;

f0=x1(ii);

g0=x2(ii);

F0=z1(ii);

G0=z2(ii);

f1=f(t0,f0,g0)

g1=g(t0,f0,g0)

F1=F(t0,F0,G0)

G1=G(t0,F0,G0)

f2=f(t0+h/2,f0+h/2*f1,g0+h/2*g1)

g2=g(t0+h/2,f0+h/2*f1,g0+h/2*g1)

F2=F(t0+h/2,F0+h/2*F1,G0+h/2*G1)

G2=G(t0+h/2,F0+h/2*F1,G0+h/2*G1)

f3=f(t0+h/2,f0+h/2*f2,g0+h/2*g2)

g3=g(t0+h/2,f0+h/2*f2,g0+h/2*g2)

F3=F(t0+h/2,F0+h/2*F2,G0+h/2*G2)

G3=G(t0+h/2,F0+h/2*F2,G0+h/2*G2)

f4=f(t0+h,f0+h*f3,g0+h*g3)

g4=g(t0+h,f0+h*f3,g0+h*g3)

F4=F(t0+h,F0+h*F3,G0+h*G3)

G4=G(t0+h,F0+h*F3,G0+h*G3)

x1(ii+1)=x1(ii)+h/6*(f1+2*f2+2*f3+f4)

x2(ii+1)=x2(ii)+h/6*(g1+2*g2+2*g3+g4)    %龙格库塔法迭代

z1(ii+1)=z1(ii)+h/6*(F1+2*F2+2*F3+F4)

z2(ii+1)=z2(ii)+h/6*(G1+2*G2+2*G3+G4)

end

n=length(x1(1,:));

diff=abs(x1(1,n)-beta);

if diff

flag=1;

end

k=k+1;

end

n=length(x1(1,:));

for iii =1:n

cc(1,iii)=x2(1,iii).*t(1,iii);

end

ff(i)=trapz(t,cc)+M;

end

while abs(ff(3))>0.001

CC=C-2*10^-3*ff(3)/(ff(2)-ff(1));

for i=1:3

if i==1

C=CC-10^-3;

elseif i==2

C=CC+10^-3;

elseif i==3

C=CC;

end

for ii=1:maxiter

H=0.25;

m=1.2;

Y=x1(ii);

Y1=x2(ii);

tt=t(ii);

sigmay=((Y/tt)^2-Y*Y1/tt+(Y1)^2)^0.5;

if sigmay>1

aa=(sigmay.^m-1)/(H*sigmay)

elseif sigmay<=1

aa=0

end

f=@(t,x1,x2)x2  %f代表y

g=@(t,x1,x2)1/(t*(1+aa))*C+1/(t.^2)*x1-1/t.*x2  %g代表y'

F=@(t,z1,z2)z2;

G=@(t,z1,z2)1/(t.^2)*z1-1/t*z2;

f0=x1(ii);

g0=x2(ii);

F0=z1(ii);

G0=z2(ii);

f1=f(t0,f0,g0)

g1=g(t0,f0,g0)

F1=F(t0,F0,G0)

G1=G(t0,F0,G0)

f2=f(t0+h/2,f0+h/2*f1,g0+h/2*g1)

g2=g(t0+h/2,f0+h/2*f1,g0+h/2*g1)

F2=F(t0+h/2,F0+h/2*F1,G0+h/2*G1)

G2=G(t0+h/2,F0+h/2*F1,G0+h/2*G1)

f3=f(t0+h/2,f0+h/2*f2,g0+h/2*g2)

g3=g(t0+h/2,f0+h/2*f2,g0+h/2*g2)

F3=F(t0+h/2,F0+h/2*F2,G0+h/2*G2)

G3=G(t0+h/2,F0+h/2*F2,G0+h/2*G2)

f4=f(t0+h,f0+h*f3,g0+h*g3)

g4=g(t0+h,f0+h*f3,g0+h*g3)

F4=F(t0+h,F0+h*F3,G0+h*G3)

G4=G(t0+h,F0+h*F3,G0+h*G3)

x1(ii+1)=x1(ii)+h/6*(f1+2*f2+2*f3+f4)

x2(ii+1)=x2(ii)+h/6*(g1+2*g2+2*g3+g4)    %龙格库塔法迭代

z1(ii+1)=z1(ii)+h/6*(F1+2*F2+2*F3+F4)

z2(ii+1)=z2(ii)+h/6*(G1+2*G2+2*G3+G4)

end

% iterations:

%

flag=0;

n=length(x1(1,:));

diff=abs(x1(1,n)-beta);

if diff

flag=1;

end

k=1;

while flag==0,

s(k+1)=s(k)+(beta-x1(1,n))/z1(1,n);

clear t;

clear x1;

clear x2;

x1(1)=0

x2(1)=s(k+1)

z1(1)=0

z2(1)=1

t=zeros(maxiter+1,1);

t=a:h:b;

for ii=1:maxiter

H=0.25;

m=1.2;

Y=x1(ii);

Y1=x2(ii);

tt=t(ii);

sigmay=((Y/tt)^2-Y*Y1/tt+(Y1)^2)^0.5;

if sigmay>1

aa=(sigmay.^m-1)/(H*sigmay)

elseif sigmay<=1

aa=0

end

f=@(t,x1,x2)x2  %f代表y

g=@(t,x1,x2)1/(t*(1+aa))*C+1/(t.^2)*x1-1/t.*x2  %g代表y'

F=@(t,z1,z2)z2;

G=@(t,z1,z2)1/(t.^2)*z1-1/t*z2;

f0=x1(ii);

g0=x2(ii);

F0=z1(ii);

G0=z2(ii);

f1=f(t0,f0,g0)

g1=g(t0,f0,g0)

F1=F(t0,F0,G0)

G1=G(t0,F0,G0)

f2=f(t0+h/2,f0+h/2*f1,g0+h/2*g1)

g2=g(t0+h/2,f0+h/2*f1,g0+h/2*g1)

F2=F(t0+h/2,F0+h/2*F1,G0+h/2*G1)

G2=G(t0+h/2,F0+h/2*F1,G0+h/2*G1)

f3=f(t0+h/2,f0+h/2*f2,g0+h/2*g2)

g3=g(t0+h/2,f0+h/2*f2,g0+h/2*g2)

F3=F(t0+h/2,F0+h/2*F2,G0+h/2*G2)

G3=G(t0+h/2,F0+h/2*F2,G0+h/2*G2)

f4=f(t0+h,f0+h*f3,g0+h*g3)

g4=g(t0+h,f0+h*f3,g0+h*g3)

F4=F(t0+h,F0+h*F3,G0+h*G3)

G4=G(t0+h,F0+h*F3,G0+h*G3)

x1(ii+1)=x1(ii)+h/6*(f1+2*f2+2*f3+f4)

x2(ii+1)=x2(ii)+h/6*(g1+2*g2+2*g3+g4)    %龙格库塔法迭代

z1(ii+1)=z1(ii)+h/6*(F1+2*F2+2*F3+F4)

z2(ii+1)=z2(ii)+h/6*(G1+2*G2+2*G3+G4)

end

n=length(x1(1,:));

diff=abs(x1(1,n)-beta);

if diff

flag=1;

end

k=k+1;

end

n=length(x1(1,:));

for iiii =1:n

cc(1,iiii)=x2(1,iiii).*t(1,iiii);

end

ff(i)=trapz(t,cc)+M;

end

end

for iiiii=1:n

T(1,iiiii)=x1(1,iiiii)/t(1,iiiii);

end

plot(t/a,T,t/a,x2)

孟海丽
关注
matlab ode45使用,ODE45函数的使用——翻译
weixin_36078604的博客
04-01 5680
《ODE45函数的使用——翻译》由会员分享，可在线阅读，更多相关《ODE45函数的使用——翻译(7页珍藏版)》请在人人文库网上搜索。1、在Matlab中使用ode45简介Matlab中常微分方程常用的函数是ODE45，这个函数能够利用-龙哥库塔法-有效求解带时间变量步长的计算。Ode45用于求解如下的一般问题：(1)其中，时间t是独立变量，x为时间相关矢量，是时间t和x的函数。当(1)右边的是固定...
ODEfunction.m
06-21
分子动力学 拟合
matlab中ode45怎么用
最新发布
weixin_42610010的博客
02-18 716
ode45是MATLAB中用于解决常微分方程(ODE)的函数，它是一种根据Runge-Kutta方法对ODE进行数值求解的函数。它可以自动根据输入的微分方程，选择最佳的步长来解决问题。使用ode45函数，需要编写一个问题函数，它接受t和y作为输入，并计算dy/dt，然后将该函数作为参数传递给ode45函数。 ...
用MATLAB函数绘制系统的,用matlab绘制odes系统定义的函数
weixin_33514300的博客
03-17 222
我试图用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_...
matlab中odefun,matlab 中ode45的源代码如何看懂
weixin_39672296的博客
03-23 310
该楼层疑似违规已被系统折叠隐藏此楼查看此楼求常微分方程的数值解 ode45方法的源代码怎么看懂呢？四百多行 如何理解这些代码的核心思想 以方便未来自己使用呢？求大神指点迷津 感激不尽function varargout = ode45(ode,tspan,y0,options,varargin)%ODE45 Solve non-stiff differential equations, ...
1.15 学习matlab里面自己编写函数以及使用ode函数
weixin_54910455的博客
01-15 3745
格式 子函数（.m）文件 function [outarg1,outare2,...]=fname(inarg1,inarg2,...) %function [输出量1，输出量2,...]=函数名(输入量1，输入量2,...) [coding] end 调用格式 [a,b,...]=fname(c,d,...); %c,d,...是在主程序里面已经确定好的（也可以递归）输入的量 %然后把叫fname的函数看成一个黑箱子，向其中输入值后它会输出值，把输出值赋给a,b,... 一个简单的
matlab.rar_ODE45_matlab ode45_matlab-ode45
09-24
在MATLAB中，`ode45`是求解常微分方程（Ordinary Differential Equations，ODE）的标准函数，适用于非 stiff 方程组。它采用四阶Runge-Kutta方法，这是一种数值积分方法，能够高效且准确地逼近微分方程的解。本文件...
matlab_滑模控制简单ode45例程及simulink仿真框图，控制对象为直流变换器
05-30
滑模控制简单ode45例程及simulink仿真框图，控制对象为直流变换器 Sliding mode control is a simple routine and Simulink simulation block diagram. The control object is DC-DC converter
ode45的笔记_ODE45_matlabode45_
10-03
若需要在特定时间点输出结果，可以使用` ode45(@dydt,tspan,y0,[],[],tout)`，其中`tout`是一个包含期望输出时间点的向量。 6. **事件检测：** ode45支持事件检测功能，可以在满足特定条件时停止求解或记录额外信息...
matlab ode45使用方法
10-10
在MATLAB中，`ode45`是求解非线性微分方程初值问题（ODEs）的标准函数，适用于常微分方程（ Ordinary Differential Equations）。它基于四阶Runge-Kutta方法，是一种数值积分算法，能有效地处理各种类型的微分方程。...
odefunc_ode45齿轮_ODE45_齿轮odefunc_齿轮_齿轮ode45_源码
10-04
行星齿轮动力学分析微分方程的matlab中ode45求解
odeFunction matlab,Imitate ode45 function from MATLAB in Python
weixin_35866747的博客
03-18 236
The interface of integrate.ode is not as intuitive as of a simpler method odeint which, however, does not support choosing an ODE integrator. The main difference is that ode does not run a loop for yo...
【MATLAB】关于ode45的一部分用法
热门推荐
Loggsy的博客
06-24 10万+
1 简介       ode45，常微分方程的数值求解。MATLAB提供了求常微分方程数值解的函数。当难以求得微分方程的解析解时，可以求其数值解（解析解就是给出解的具体函数形式，从解的表达式中就可以算出任何对应值；数值解就是用数值方法求出近似解，给出一系列对应的自变量和解）。         Matlab中求微分方程数值解的函数有七个：ode45，ode23，ode113，ode15s，ode...
matlab ode45
weixin_43238322的博客
09-18 9342
matlab ode45 求解微分方程 举例老师的例子： 写一个function： ode的公式 odefunction就是刚刚写的那个function tspan时输入变量的区间或者点 y0就是选择刚开始时候的y0的值 比如说下面这个图 [0,3]就是x在零至三的区间 -1就是起始点y=0 plot一下 ...
matlab保存ode45计算中间值并绘图
anbuqi的博客
03-24 3374
matlab保存ode45计算中间值并绘图
matlab 定义一个有自变量的方程_Matlab求微分方程组(非刚性)
weixin_39810901的博客
11-28 1738
Matlab中ode系列函数的用法 Matlab的ode系列函数（统称odesolver）专门用于求解常微分方程（组）的初值问题。ode是常微分方程（ordinary differential equation）的英文缩写。 Matlab中采用ode系列函数求解常微分方程（组）的初值问题一般需要两个文件：主程序文件（如：main.m）和定义函数odefun的文件（即odefun.m...
matlab中ode45用法,ode45(ode45用法举例)
weixin_42567279的博客
03-19 1万+
ode45是用4阶方法提供候选解，5阶方法控制误差，是一种自适应步长的方法。而我们平时用的4阶和5阶龙格库塔法的公式中步长是给定的。具体算法和原理你可以看.ode45的初始条件是否必须是在x=0处没有必要 只要是选取的初值对应的t等于tspan(1) 我们看下面一个例子，对于dx/dt=2t. %注意下面的t0必须等于2，因为我们的选择的初值条件是t=2,ode45默认将t0作为初值.请问这个要怎...
多变量微分方程组带事件控制的ODE45函数写法
Kasen‘s experience on reserch & study
12-27 8726
1.来源 在CUHK学习advanced robotics时做hopping robot仿真时多次使用自己的判断条件终止ODE45方程组，遂整理在此 2.ODE45函数调用 [t,x,te,xe，ie] = ode45(@(t,x) myStanceODEFct(t,x,modelPara,Tau(t,x)), t_span, x0, option); [t,y,te,ye,ie] = ode...
ode中函数参数传递
youhuakongzhi的博客
06-17 4936
转自：http://ma-yue.net/blog/archives/89在用odesolver(ode45, ode15s, …)来解微分方程的时候，最基本的用法是：[t, y] = odesolver(odefun, tspan, y0);这里的odefun是待求的微分方程。那么odefun中一般会含有多个系统参数，通常要通过改变参数来观察系统动态的变化。那么如何在调用odesolver的时候...
MATLAB中ODE45函数的使用方法和详解
MATLAB中函数ODE45的使用方法和详解 MATLAB中函数ODE45是用于解决常微分方程的初值问题的求解器之一。与其他求解器（如ode23、ode113、ode15s、ode23s、ode23t、ode23tb）相比，ODE45是MATLAB官方推荐的优先选择。 ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值