matlab如何求解4次方程

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

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

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

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])

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

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

%%定义xx为最终解x=solve('x*(10+30*sqrt(1+(20*x)^2))-y');i=1;foryy=0:0.01:1xx(:,i)=subs(x,'y',yy);i=i+1;end

这类方程可以用fsove（）函数来求解.myfun=@(x)2*sin(0.8*x)*sinh(0.8*x)-88.85*x^3*(sin(0.8*x)*cosh(0.8*x)-cos(0.8*x)*

clcsymsabcthetay=b*b*c*c-a*a*(cos(theta)^2*c^2-sin(theta)^2*b^2);theta=solve(y,theta)大概就是这个命令

symsxkeq=x^3+x^2+0.25*x-(k+1);x=solve(eq,x)%你说的那种形式不可能

举个例子,比如被积函数是x^2,上限已知为5,定积分的值是10,求下限a.f=@(x)x.^2;%被积函数b=5;%上限val=10;%定积分值g=@(a)quadl(f,a,b)-val;a=fze

在较低版本的matlab中,solve只能用来求解解析解,当没有解析解时返回一个空矩阵较高版本中好想也能部分地数值求解了也就是说当solve的对象内含有字符的话默认进行符号求解,若没有除变量外的字符则

这个是没有实数解的……我用Mathematics帮你算了下复数解

为什么要解方程呢.直接写出解析解当函数就可以了啊.

subs(SA)可以把工作区的数值带入表达式求函数值,solve(SA)是用于求解SA=0时的自变量的值

%%下面方法比较笨,但是可行.%%先看一下为什么你解不出来%%把你的k*1.57e-6看成一个整体x%%变成求解tan(x)./x-1/0.2159==0x=-2:0.0001:2;y=tan(x).

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

1、y=solve('9.8*t^2+t-2')2、y=solve('98*t^2+10*t-20=0')3、a=[9.81-2];roots(a)自己试下第一种和第二种的区别!

方程在0附近高度跌宕,没法用solve方程求解.你可以画图看一下就知道了,解都集中在x=-1.6到0.2之间,其他区域无解.图就不上了,你自己画吧.