matlab 解四元四次方程

用solve命令解

function[x,y,z]=sp(A,B,C)symsxyz[x,y,z]=solve(y*y-z*z-x*x,y+z-A,x*x-B*x-C);

MATLAB固然强大,但不是万能的但要想充分发挥其威力,需要使用者具备数学知识,不能硬算,否则难以得到相要的结果,或者即使算出个数值结果来,也很丑陋这里应当先做三角函数的万能替换设t=tgβsin&#

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

中间还有几个乘号也加点试试呗x=-35:35在函数中用x（1,：）表示是可以的ps.你的程序可读性差了点吧……再问：程序可读性差什么意思啊？再答：就是别人看不懂……还有，提醒你一下，MATLAB的si

初值条件看不太懂

x2的表达式修改如下：X2=solve(['(x-1)*(x-2)*(x-3)*(x-4)*(x-5)*(x-6)*(x-7)*(x-8)*(x-9)*(x-10)*(x-11)*(x-12)*(x-

y=@(x)-((16-x).^2*(576-(24-x).^2)./(24-x).^2);[x,fval]=fminunc(y,8)Optimizationterminated:relativein

%bydynamic%2009.2.%直接将数据导入到Matlab中,再使用循环语句求解ab=load('data');%每行一组a、bforii=1:size(ab,1)root(ii)=Newto

solve('5*(-0.0000003)*X^4+4*(0.0000308)*X^3+3*(-0.00105)*X^2+2*0.01209*X-0.000021=0','0

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

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

a=1;b=2;symsxy;[xy]=solve('x+y=b','x*y=a');X=eval(x),Y=eval(y)按上面的语句试试

fun=inline('-504.4-0.131.*u.^2+14.297.*(-19.313+295.27.*(0.0421.*u)-165.44.*(0.0421.*u).^2+40.874.*(

设该四次方程为a0x^4+a1x^3+a2x^2+a3x+a4=0输入roots([a0a1a2a3a4])即可再问：怎么算出解析式

eqn=@(x)4*x^4-4*x^2;fzero(eqn,0)fzero(eqn,-1.1)fzero(eqn,-1)fzero(eqn,-0.9)fzero(eqn,0.0)fzero(eqn,1

不带参数：solve('x^4-x^3+x^2=0')单引号内式子可以任意改变,但形式要与例子一致.带参数：symsabcx;solve('a*x^4-b*x^3+c*x^2=0',x),要解变量a就

设该四次方程为a0x^4+a1x^3+a2x^2+a3x+a4=0输入roots([a0a1a2a3a4])即可

symsaxbcy=a*x^4+b*x+c;Y=solve(y,'x')结果如下：ans=1/12*6^(1/2)*(12^(1/3)*(((9*b^2+(-768*c^3*a+81*b^4)^(1/

用Mathematica这样求解Solve[d^4-1.2410^-4d-1.29610^-5==0,d]答案是：{{d->-0.0508017},{d->-0.00860384-0.0606388I