matlab数值法求解多元二次方程组

还是回答你的那个例子,你需要用subs函数将符号运算转化成常数运算,subs(z,{x,y},{2,4})关于subs的具体信息可以参考matlab的帮助文件SUBS  Symbo

clc;clear;y=sym('a+b*x+c*x^2+d*x^3+e*x^4+f*x^5')%是不是c*x^2啊y1=diff(y)y2=diff(y,2)y3=diff(y,3)eq1=subs

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

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

solve('a*x^2+b*x+c=0')ans=-1/2*(b-(b^2-4*a*c)^(1/2))/a-1/2*(b+(b^2-4*a*c)^(1/2))/a所以你如果带入直接的数字的话,出来的

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

……一定要用matlab吗,只会mathematica的路过.而且啊,你一面说要求方程的数值解,一面又只给一个全是字母系数的方程,你是要怎样……再问：呃，abcdef都是常量，可以随便假定的，就赋1、

在Matlab下输入：editzhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

对于线性方程组,可以建立矩阵来求解.把ABCDEF可以分别看成x（1:6）；具体的指令如下：cleara=[000110;101010;010sqrt(3)0-sqrt(3);00sqrt(3)1-s

解方程用solve()函数,然后注意自己定义的变量即可,下面我们举一个例子来讲x+2*y=5,3*x-y=7这是一个二元一次线性方程组所以在matlab中按照下面输入：>>[xy]=solve('x+

s=solve('c1=0','c2=0','c3=0','a1','a2','a3');改为s=solve(c1,c2,c3,'a1','a2','a3')

clear;clc[T,Dxyz]=ode45('sdfd',[012],[0;1;1]);plot(T,Dxyz)--------------------------functiondxyz=sdf

你去看下fsolve的帮助吧简单来说,你先把方程化成F(x)=0的形式,编个function,保存在同一文件夹下,然后[x,fval]=fsolve(F,x0).x是解,x0是迭代初值,fval是最后

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

试了一下,发现在一定的边界条件下,可以如：①设定每个变量的初值边界条件M文件functionres=ivbc(ya,yb)q=5;res=[ya(1)-qya(2)-2ya(3)-3ya(4)];②设

eq1eq2eq3未知?再问：都等于0，且x小于等于28654，y小于等于18449，z小于等于308，x,y,z都是大于等于0的。再答：前几天忙，现在才帮你做。sorry了。m文件如下：f=@(x)

因为你x=0时2/x是无穷大呀,然后y'又是0,然后(2/x)y'就是nan了,所以后面算的全是nan了.

有结果啊：>>Tab=(sqrt(3)-1)/340;>>Tac=(sqrt(3)-3/2)/340;>>v=340;>>SA=1/2*(-4+4*v^2*Tab^2+4*v^2*Tac^2)*(-4

fsolve(@myfun,x0,options)

这题就是为了说明定理条件"偏导连续"不能换成"偏导存在"啊就是给了一个具体的函数z,然后说明对z来说那个公式不成立.因为你没有具体说你是哪里不明白,所以我只能如上大概地说一下了.再问：我的理偏导数存在