matlab矩阵求求解微分方程组

x=dsolve('Dx=r*(1-x/xm)*x','x(0)=x0','t')x=xm/(1+exp(-r*t)*(xm-x0)/x0)

functionmyCalcclear;clc;closeall;opt=odeset('RelTol',1E-5,'AbsTol',[1E-61E-6]);[TY]=ode45(@getdY,[01

结果：代码：clearallclcf=@(x,y)([y(2);   0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]

请参考ode45函数的用法将所求值设为终值即可

首先,得有微分方程,然后在Simulink中用几个模块就搭建出仿真模型了,比较简单.你试试吧再问：我试过了，不过没成功，能截个图吗？

V=3*ones(1000,1);V1=ones(999,1);A=diag(V,0)+diag(V1,1)+diag(V1,-1);B=[1:1000]';X=A\B;得出X就是想求得值.

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

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.

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

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

估计没人会

前两天刚回答了一个类似的问题：再问：非常感谢你，那个问题也是我发的，但我处理二阶就不会了，失败了。听同学说没有解析解，想问下怎么写二阶形式的求解，只要ODE45的方法就行了，别的不麻烦您了再答：>>o

先建立一个M文件：functiony=t(x)y=[1xx.^2x.^3];然后在MATLAB内执行：>>Y=[t(1);t(2);t(3)]Y=1111124813927这样子便可以得到你想要的矩阵

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

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

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

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

好吧,再修改成|K-w^2*M|=0实际上就是个三次多项式的求解程序：symsvM=[2100;1410;0141;0012]K=[1-100;-12-10;0-12-1;00-11]f=det(K-

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)