matlab求解dy dx=-ay b

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

方程两边求关x的导数ddx(xy)＝(y+xdydx)；     ddxex+y＝ex+y(1+dydx)；所以有  （y+xdy

由微分方程dydx=2xy，得dyy＝2xdx（y≠0）两边积分得：ln|y|=x2+C1即y＝Cex2（C为任意常数）

可以的,假设行向量X=（x1,x2,.,xn）,Y=（y1,y2,.,yn）,那么x+ty=(x1+ty1,x2+ty2,...,xn+tyn);||x+ty||^2=(x+ty)*(x+ty)'=c

fun=(64*pi*(27889/10000-n^2)^(1/2))/31-atan((2479*(n^2-9/4)^(1/2))/(2000*(27889/10000-n^2)^(1/2)))-a

你确定你用symsasolve('arcsin(1/a)+1/a*sqrt(1-(1/a)^2)=pi/4')这段代码解出的值是负数?可我用你的代码求出的怎么不是负数呢,而且还是你说的2.5左右当然,

y='abs(exp(-x).*cos(10*x))-0.05';ts=fzero(y,2)结果ts=2.0040这里面只能用x,用t就会出错,至于原因我也不知道.

A=[-816;4-41;44-7];b=[5;1;2];x=A\b%直接利用matlab中函数即可.还可分析A是否可逆等.当然也可自己编写程序求解.

二次方程组用fsolve求解吧,用solve可能接不出来呢,因为solve是解线性方程组的.再问：fsolve也试过了，依然无法求解呀clearP0=[1e26;1e26;1e26;1e26;1e26

1、在文件编辑区建立待求方程组文件并保存：functiony=fun(x)y=[0.56-1.1018*x(1)*(exp(-0.1855014*x(2))-exp(-2.007944*x(2))),

exp前面那个*改为.*（点乘）保险的方法,从头到尾检查凡是乘和除都改成点乘和点除.混着用容易漏掉.

没错啊,解不出来很正常,solve是符号工具箱的,用来求解析解,很弱的你这种竖直问题去看fsolve函数吧

functionf=fun(x,y)f=[y(2);1/x-y(1)/x^2-y(2)/x];在commandwindowtype:y0=[1,2];[x,y]=ode45[@fun,[0,8],y0

相关系数就用命令corrcoefmin(min(corrcoef(x1,x2)))就是x1,x2之间的相关系数.比如t=(1:0.1:100)';w=2*pi;x1=sin(w*t)+randn(si

小问题1似乎是特征分解.[V,D]=eig(K);这样就可以得矩阵V和对角阵D,满足K*V=V*D再问：恩。。这样特征值对角阵的确可以求出来，变化向量P怎么求了呢再答：P不就是V么。。。。V是单位正交

1.首先,这里有8个方程,3个未知量,当然你可以用其中任意三个线性无关的方程求出a,b,c的一组值,但是这样我们就浪费了很多数据.2.如果同时求解8个方程肯定无解,但是我们想要找到一个向量Y使得AY与

x应当是一维行向量,x'共轭转置后,变成一维列向量,尤其是实数时.length(x)是求出x的元素个数ones(length(x),1)是构造一个矩阵,length(x)行1列.[x'ones(len

这是个超越方程,答案含有朗伯W函数clearall;clc;symsay=solve('log10(x)-426.5*x=a','x');y=exp(-1.*lambertw(-982.0525421

dydx要是等式才行吧.如果是的话,这句话就是求这个等式的根,用r表示x.

Z=[0.50.40.30.250.20.150.1];X=[630349213544288022481433964];Y=[980738574500365267150];M=[X;Y]';z=Z';