matlab微分方程的加速度图像
来源:学生作业帮助网 编辑:作业帮 时间:2024/10/05 15:24:59
请参考pdepe函数的用法
x=dsolve('Dx=r*(1-x/xm)*x','x(0)=x0','t')x=xm/(1+exp(-r*t)*(xm-x0)/x0)
首先,得有微分方程,然后在Simulink中用几个模块就搭建出仿真模型了,比较简单.你试试吧再问:我试过了,不过没成功,能截个图吗?
dsolve函数默认自变量是t,试试soln_1=dsolve('Dy=3*x^2','y(2)=0.5',’x‘)如果还是不行就不知道了可以了我运行过了>>soln_1=dsolve('Dy=3*x
一维导热问题,应该采用有限差分法求解.具体方法随便找个数值计算方法书都有
可以用simulink做啊,很简单的,只需要托模块就可以了
matlab中导数是diff.help一下,就什么都有了
这里面fo(t)是什么假设fo(t)=f[c1,c2]=dsolve('Dc1=-(k12+k13)*c1+V2/V1*k21*c+f/V1','Dc2=V1/V2*k12*c1-k12*c2','t
你可以用mathematic试试,不懂的可以通过帮助函数查找,但是势都有边界条件,比其编程加入,个人感觉还是手算方便的多.
看了看,运行了一下确实出问题 原因是在用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
请查看dsolve命令,输入helpdsolve回车即可!
先定义函数,再用ode15s或ode45解.建一.m文件,在里面输入:functiondy=xiao(x,y)dy=zeros(1,1);dy(1)=5.42*exp(-5.42)-2.5*y(1);
你最好先给出数学模型再问:我就是数学模型没做对,我的方程不是收敛的,是发散的。再答:发散可能有两种原因:(1)数学模型正确,但m1、m2、k1、k2、c以及小车初始位置的数值给的不对,引起发散,建议换
symsxa=5;b=a*x;w=dsolve('Dw=b','x')subs(w)%将a,b代入w的表达式.
这个你没给初值条件啊我假设是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('
恩你没有给初值不要解答呀记住加分哦%bydynamic%2009.2.16%%我这里假设初值全部为1%下面是解析解root=dsolve('D2x=2*x*y','D2y=y^2-x^2-20^2',