matlab求解微分方程du dt=cos(sint),u(0)=1

请参考pdepe函数的用法

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]

朋友,要根据初值积分对dM/dt积分求得M,才能求解M(t)=0或t(M)=0.solver（积分函数,积分时间,初值,设置）也是这样的数学方法.[时间,解]=solver（积分函数,积分时间,初值）

用dsolve求不出解析解,用数值解法吧,最好能提供初始条件,这里假设t=0时,x0=1,y0=2,z0=3.　　functionhh　　[t,x]=ode23t(@fun,0:1,[123])　　p

你可以用mathematic试试,不懂的可以通过帮助函数查找,但是势都有边界条件,比其编程加入,个人感觉还是手算方便的多.

symstv=dsolve('Dv=(190.708-90.64*v^2)/47.27','v(0)=0','t');t=0:0.00001:0.002;v=eval(v);plot(t,v)使用这样

看了看,运行了一下确实出问题 原因是在用ode数值求解时,x并不是1:0.01:3均匀分散的 解决方法：1.在画解析解和欧拉解时横轴用x的转置；在画数值解时横轴用x,已运行成功2.

symsxyy=sqrt(2*pi/x)*sin(x);y1=diff(y,'x');y2=diff(y,'x',2);h=x^2*y2+x*y1+(x^2-1/2)*y;h=simplify(h)结

我运行的>>symsaknNzz1>>z=dsolve('Dx=a*x*(N-x)','t')z1=dsolve('Dx=a*x*(n-x)','t')结果：z=N0N/(exp(-N*(C3+a*t

在matlab中默认的自变量是t,所以你的结果中在自变量是t,x被当做一个参数了如果你的自变量是x,你在写方程的时候把x写出t就行了!例如y=dsolve('Dy=t')y=t^2/2+C18

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

Matlab偏微分方程工具箱应用简介1．概述本文只给出该工具箱的函数列表,读者应先具备偏微分方程的基本知识,然后根据本文列出的函数查阅Matlab的帮助,便可掌握该工具箱的使用.2．偏微分方程算法函数

由于方程比较复杂,解析解不能用初等函数表示只是要获得图像的话,用数值计算的办法可能更方便fun=@(t,y)1.44*(10^9)*(1-y).*exp(-109170./(8.314*t));[T,

是可以求解的,参看pdepe函数

i是虚数单位?那个1/3γ（u*v）中*是什么符号

解一阶偏微分方程,这是一道难题哦.你连初始条件都没有啊,没法做再问：x=0u=5000K,x=Lu=2000K,L=10m,t=0u=5000K.

最常用的就是广义特征向量基础矩阵解方法.你要一个思路,我给一个2维情况的例子,其中特解x(t0)=x0的理解和如何使用都有,你看看是否够用.. Matlab下二维的例子：再问：嗯，这个不错，

改动：functionS=history(t)S=ones(3,1);改为S=ones(4,1);OVER!PS.多翻翻help!

y=dsolve('2000*Dy-(0.08-y*0.08)','y(0)=0','t')y=1-exp(-1/25000*t)即：C(t)=1-exp(-1/25000*t)