matlab 矩阵微分方程
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/12 23:58:53
x=dsolve('Dx=r*(1-x/xm)*x','x(0)=x0','t')x=xm/(1+exp(-r*t)*(xm-x0)/x0)
结果:代码:clearallclcf=@(x,y)([y(2); 0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]
dsolve函数默认自变量是t,试试soln_1=dsolve('Dy=3*x^2','y(2)=0.5',’x‘)如果还是不行就不知道了可以了我运行过了>>soln_1=dsolve('Dy=3*x
(1)B矩阵需要预定义(2)你所谓的8*8是这个意思吧32*32\x0932*32\x0932*32\x0932*32\x0932*32\x0932*32\x0932*32\x0932*3232*32
一维导热问题,应该采用有限差分法求解.具体方法随便找个数值计算方法书都有
设缉私艇航速为v1,走私船航速为v2,航行时间为t,以缉私艇出发点为原点,正东方向为x轴,正北方向为y轴,缉私艇经t时间行驶到(x,y)点.经过的航迹为y=f(x);y'=(v2*t-y)/(c-x)
matlab中导数是diff.help一下,就什么都有了
朋友,要根据初值积分对dM/dt积分求得M,才能求解M(t)=0或t(M)=0.solver(积分函数,积分时间,初值,设置)也是这样的数学方法.[时间,解]=solver(积分函数,积分时间,初值)
请查看dsolve命令,输入helpdsolve回车即可!
这次看看,程序通了.functionhhh[t,x]=ode45(@xprim2,[0,20],[30;20]);plot(t,x);xlabel('timet0=0,tt=20');ylabel('
恩你没有给初值不要解答呀记住加分哦%bydynamic%2009.2.16%%我这里假设初值全部为1%下面是解析解root=dsolve('D2x=2*x*y','D2y=y^2-x^2-20^2',
前两天刚回答了一个类似的问题:再问:非常感谢你,那个问题也是我发的,但我处理二阶就不会了,失败了。听同学说没有解析解,想问下怎么写二阶形式的求解,只要ODE45的方法就行了,别的不麻烦您了再答:>>o
这只是求偏导数而已,你别搞笑哈.给你代码symsxyzz=(x^2+y^2)*exp((x^2+y^2)/x/y);z_x=diff(z,x)z_xx=diff(z,x,2)z_xy=diff(dif
最常用的就是广义特征向量基础矩阵解方法.你要一个思路,我给一个2维情况的例子,其中特解x(t0)=x0的理解和如何使用都有,你看看是否够用.. Matlab下二维的例子:再问:嗯,这个不错,
这样可以吗?t=1:6m=[556569727885]qi=[151418192120]ppi=[212434]v=[202225273133]plot(t,[m;v],'o-')vsp=spline
建立函数文件-(文字部分自己替换为表达式)functiondifxy=fun(t,y)difxy=[x对t的导数表达式;y对t的导数表达式]调用ode45指令--绘图区间可换[t,y]=ode45('
这个方程是解不出解析解的. 我用ODE45给您算了一下. 先
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functiondy=zhi
clearallclcf=@(t,x)([x(2);-x(2)+100*x(1)+1+200*cos(2.5*t)]);[t,X]=ode45(f,[01],[142.510604]);plot(X(
functionhahaha[t,r]=ode45(@myfun,[01],0)functiondrdt=myfun(r,t)drdt=-i*(50.24-235.5*t)*r+3.14*(1-r^2