ode45中odeset函数

这个方程要解出来还要初始条件吧不妨设(x(0),y(0))=(A,B);还有用ode45只能得到数值解,不能给出解析式表达那么先建立一个函数文件test_fun.mfunctiondx=tsst_fu

我去,你这个……肯定得先把方程化成一阶线型方程组哇,二阶导数怎么算!百度文库找个实例看看

初值条件看不太懂

函数FreeOcillation(t,x,dummy,zeta)写错了,应该把那个dummy删掉.该参数在函数中并未使用,而且调用的时候也没有提供相应的数据.命令窗口的空矩阵代表ode求解器的选项.各

把方程和初值条件都列出来这一大堆东西几个人看得懂?再问：公式（1）（2）就是方程啊，我写的那两个就是对应的方程。倒数第二张图片有实验数据。还需要什么参数吗？非常抱歉，我学的专业不太扎实，问题可能有点弱

示例如下：odefun=@(t,y)(y+3*t)/t^2;%定义函数tspan=[02.3];!y0=-2;%初值[t,y]=ode45(odefun,tspan,y0);其中求解区间为0~2.3,

在Matlab下输入：edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functiondx=myf

functiontest()[t,y]=ode45(@func1,[0,1],[0;0;1;2;2;2]);figure(1);clf;plot(t,y);legend('x','y','z','dx

％先建立一个M文件,然后都复制进去,运行试一试functionMyexamp1cleartspan=[0:0.01:20];y0=[0;0.25];[tt,yy]=ode45(@vandepol,ts

建立如下m文件：function xdot=fun(t,x)xdot=zeros(2,1); xdot(1)=-2*x(1)+x(2); xdot(2)=-3*x(1)+

zs是多少?再问：T=2;%调用函数的时候改成100*T；fs=1000;noise1=randn(1,fs*T);N=6;%滤波器的阶数wn=[0.2,0.4];%截止频率[b,a]=butter(

步长可以设定的,比如设为5t=0:5:100,就可以了再问：ode本身是变步长的？这个要怎么设定啊再答：ode45('函数名',[0:5:100],0)

不能这么干的呢ode45()都是要有函数调用的再说了eq是什么啊,这点eq应该是函数名应该把函数单独保存一个m文件,另外的求解再单独保存一个m文件,运行求解的m文件就能调用相应的函数的了哦

dy(1)是对y(1)求导数的意思后面一样了

1、你贴出来的报错信息和代码对不上号：前面显示错误的那行代码和你贴出来的完全不一样；而后面的错误（Inputargument'u1'isundefined）也不可能是目前的代码所导致的

我来解答你的关键问题：将之前主程序计算出来的系数用globalcof定义,再在ode45的子函数里同样定义globalcof,这样就能把值计算出来并且传递进去了．

没有必要只要是选取的初值对应的t等于tspan(1)我们看下面一个例子,对于dx/dt=2t来说吧我们知道精确解是x=t^2下面我们给出数值解法,编写代码%%选取状态变量%x1=x%%则各个状态变量的

1、后两个方程都是dH?2、最后一个dH的表达式有误（括号不匹配）.3、st的初值?

因为你x=0时2/x是无穷大呀,然后y'又是0,然后(2/x)y'就是nan了,所以后面算的全是nan了.

参考代码： %参数定义m1=1;m2=2;m3=3;k1=4;k2=5;k3=6;k4=7;P0=8;w=9;%微分方程A=[k1+k2-k10;-k2k2+k3-k3;0-k3k3+k4]