ode45解微分方程是怎么降阶处理
来源:学生作业帮助网 编辑:作业帮 时间:2024/05/16 01:52:35
function z=krsx(t,y) &nb
我去,你这个……肯定得先把方程化成一阶线型方程组哇,二阶导数怎么算!百度文库找个实例看看
初值条件看不太懂
函数FreeOcillation(t,x,dummy,zeta)写错了,应该把那个dummy删掉.该参数在函数中并未使用,而且调用的时候也没有提供相应的数据.命令窗口的空矩阵代表ode求解器的选项.各
初值不够呀,应该有四个,题目能不能贴出来?
把方程和初值条件都列出来这一大堆东西几个人看得懂?再问:公式(1)(2)就是方程啊,我写的那两个就是对应的方程。倒数第二张图片有实验数据。还需要什么参数吗?非常抱歉,我学的专业不太扎实,问题可能有点弱
请查看dsolve命令,输入helpdsolve回车即可!
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
你把lgkt中的x改成x(1)试试看,我试过应该没问题的function xx=lgkt(t,x)xx=3*x(1)+cos(t)+sin(2*t);[t,x]=ode45(@lgkt,[0
新的matlab版本好像不鼓励采用global了.你的全局变量有点多了,哈哈.简单例子:m=2;[t,y]=ode45(@(t,x)f1(t,x,m),[0,10],[2])functiondy=f1
代入验证,可见不是y=x^2才是,你输错了!此时y'=2xxy'=2x^22y=2x^2再问:y''=1+y^2y=xe^2谢谢啊再答:应该是y=xe^x吧y'=(x+1)e^xy"=(x+2)e^x
function[x,y]=runge_kutta1(ufunc,y0,h,a,b)%参数表顺序依次是微分方程组的函数名称,初始值向量,步长,时间起点,时间终点,n=floor((b-a)/h);%求
dy(1)是对y(1)求导数的意思后面一样了
1、你贴出来的报错信息和代码对不上号:前面显示错误的那行代码和你贴出来的完全不一样;而后面的错误(Inputargument'u1'isundefined)也不可能是目前的代码所导致的
因为你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]
特征方程为r^2+ar=0其解为:r=0,r=-a(1)若a=0,则对应方程为y''=b积分两次得:y=C1+C2·x+b/2·x^2(2)若a≠0,则对应齐次线性方程的通解为Y=C1
functiondy=myfunc(x,y)a=b=c=dy=zeros(2,1)dy(1)=y(2)dy(2)=(b*dy(2)^2-c*sin(y(1)))/a这里的dy(1)是一阶导数,所以你的