matlab解常微分方程出现[ empty sym ]

你这个方程应该是在初条件v(0)=0时解的咯,没有初条件matlab解不出来.语句如下（我现在没条件在软件中试,如果不行你告诉我再改~）G=527.436*0.4536*9.8;F=470.327*0

%欧拉法解一阶常微分方程% y'=xy^(1/3)f = inline('x*y^(1/3)','x','y');

一维导热问题,应该采用有限差分法求解.具体方法随便找个数值计算方法书都有

设缉私艇航速为v1,走私船航速为v2,航行时间为t,以缉私艇出发点为原点,正东方向为x轴,正北方向为y轴,缉私艇经t时间行驶到（x,y）点.经过的航迹为y=f（x）；y'=(v2*t-y）/(c-x)

请查看dsolve命令,输入helpdsolve回车即可!

>> y=dsolve('Dy=1-y^2','y(0)=0') y = tanh(t)>> 

你这个属于边值问题.得用打靶算法.如果想简单,你令一个变量为x=y-2001,就变成正常的微分方程.

这个你没给初值条件啊我假设是x0=[0,0];x(1)是A,x(2）是R程序如下首先定义这个方程functiondx=funct(t,x)dx=zeros(2,1);dx(1)=0.00875+7.5

这次看看,程序通了.functionhhh[t,x]=ode45(@xprim2,[0,20],[30;20]);plot(t,x);xlabel('timet0=0,tt=20');ylabel('

第一题y=dsolve('Dy=y-2*y^3*sin(t)','y(0)=1')[T,Y]=ode45(@(t,y)y-2*y^3*sin(t),[0,3],1);plot(T,Y)第二题（1）[T

恩你没有给初值不要解答呀记住加分哦%bydynamic%2009.2.16%%我这里假设初值全部为1%下面是解析解root=dsolve('D2x=2*x*y','D2y=y^2-x^2-20^2',

首先建立vdp1.m文件functiondx=vdp1(t,x)dx=zeros(2,1);dx(1)=x(2);dx(2)=(10000*x+(96.321*exp(-t/378.47108)+28

x(:,1)意思是访问x矩阵的任意行的第一列举个例子>>x=randn(3,3)x=0.7254-0.20501.4090-0.0631-0.12411.41720.71471.48970.6715>

这只是求偏导数而已,你别搞笑哈.给你代码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

x = 0:0.2:6;y = 0:0.1:2;[x,y] = meshgrid(x,y);fxy = y.*(1-y)

欧拉方法的matlab先定义函数eulerfunction[x,y]=euler(fun,x0,xfinal,y0,n);ifnargin

这样可以吗?t=1:6m=[556569727885]qi=[151418192120]ppi=[212434]v=[202225273133]plot(t,[m;v],'o-')vsp=spline

    这个方程是解不出解析解的.    我用ODE45给您算了一下.    先

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

functionhahaha[t,r]=ode45(@myfun,[01],0)functiondrdt=myfun(r,t)drdt=-i*(50.24-235.5*t)*r+3.14*(1-r^2