matlab求解jacobian矩阵

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

可以的,假设行向量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

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

使用数值解法吧x=1%给x赋值fun=@(w)exp(i*x*cos(w))quadl(fun,0,2*pi)

你的b不该是3*1的向量吗?返回的结果c是3*1的向量,这样a(3*3)*c(3*1)-b(3,1)才有意义啊.

你要放在M文件中,在命令窗口不能定义函数

functionhhX0=[17,195,1500,4,62,200];TSPAN=[04];[t,X]=ode45(@myfun,TSPAN,X0)plot(t,X)legend('x','y','

S和e分别是什么?再问：哦，S也已知，是693.5，e上面给过了，是10再答：可以用fsolve求解。参考代码（我的理解，21140.25应该是指σ²）：a=50;b=1;c=29;e=10

func1=@(x)[log(x)-cos(x)]root=fzero(func1,[pi/42*pi])

那是因为你没有指定输出宗量的顺序!你的u,v,w,p并没有指定表示那个变量的值.再问：能不能说下怎么指定？谢谢再答：[uvwp]=solve(a,b,c,d,'x','y','z','l').前面和后

请问：①你是得到了数值解了吗?②laplace方程只有在特殊边界条件下存在真解,若要求的真解,请具体给出计算域的边界条件.③使用surf函数并不难,需要得到二维数组Z（若laplace方程是二维的,而

symsabcdx%定义符号变量abcdxs=a*(b*log(x)/x)^0.5-c;solve(s,'x')