用matlab解微分方程时出现dsolve未定义是怎么回事?
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/23 07:23:21
![用matlab解微分方程时出现dsolve未定义是怎么回事?](/uploads/image/f/6250307-59-7.jpg?t=%E7%94%A8matlab%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E6%97%B6%E5%87%BA%E7%8E%B0dsolve%E6%9C%AA%E5%AE%9A%E4%B9%89%E6%98%AF%E6%80%8E%E4%B9%88%E5%9B%9E%E4%BA%8B%3F)
结果:代码:clearallclcf=@(x,y)([y(2); 0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]
%给你参考一下我当年做的一个类似的程序%cleardt=0.004;t=0:dt:5;a=-pi:0.2:pi;l=100;r=8;d=10;%导弹尺度参数:长度、半径、火焰尺度L=l;D=l/3;c
一维导热问题,应该采用有限差分法求解.具体方法随便找个数值计算方法书都有
设缉私艇航速为v1,走私船航速为v2,航行时间为t,以缉私艇出发点为原点,正东方向为x轴,正北方向为y轴,缉私艇经t时间行驶到(x,y)点.经过的航迹为y=f(x);y'=(v2*t-y)/(c-x)
请查看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
dy(1)是对y(1)求导数的意思后面一样了
A=[-53;100-301];[VD]=eig(A);lambda=diag(D);C=V\[50;100];t=linspace(0,1,1000)';y1=C(1)*V(1,1)*exp
x(:,1)意思是访问x矩阵的任意行的第一列举个例子>>x=randn(3,3)x=0.7254-0.20501.4090-0.0631-0.12411.41720.71471.48970.6715>
1、把G=1/3*((5*Pp-2*P1)/(P1-2*Pp)-P1*c1^2/Pp*cp^2);改成G=1/3*((5*Pp-2*P1)/(P1+2*Pp)-P1*c1^2/(Pp*cp^2));表
这只是求偏导数而已,你别搞笑哈.给你代码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
参考代码:a=0.5; b=0.6;dIdt=@(t,I)a*(1-I)*I-b*I;ode45(dIdt,[0 10],0)存在的问题是,如果初值为0,可以计算出导数dI/dt=
这样可以吗?t=1:6m=[556569727885]qi=[151418192120]ppi=[212434]v=[202225273133]plot(t,[m;v],'o-')vsp=spline
这个方程是解不出解析解的. 我用ODE45给您算了一下. 先
先定义地字符,再求eg:symsxy;dsolve(Dy=x)simiplify(y)或直接dsolve('Dy=x');simiplify(y)
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functiondy=zhi
functionf=ffun(t,x)f=zeros(3,1);f(1)=-x(1)+0.5*(x(3)+x(2))*(x(3)+x(2)+sin(x(3)));f(2)=-0.5*x(1)*(x(3
functionhahaha[t,r]=ode45(@myfun,[01],0)functiondrdt=myfun(r,t)drdt=-i*(50.24-235.5*t)*r+3.14*(1-r^2