matlab解五元方程 matlab解三次方程

你说的那个可不是四元一次方程组哦,这个非线性方程组,一般只能用数值解法,solve函数就不行了.用fsolve取适当的初值是可以解出的

clear all; clc; syms x [S]=solve('121*x^5-1210*x^3+1705*x^2-440*x-109=0'); for i=1:length(S) disp(strcat(sprintf('第%d个解:',i),sprintf('%.4f',double(S(i))))) end; 看结果:第1个解:1.9844 第2个解:1.2787 第3个解:0.6301 第4个解:-0.1505 第5个解:-3.7426

首先要把这个二阶的方程组转换为四个一阶的方程,然后用Romberg方法求解,也就是ode45这个函数.

X=[2.5 25 62.5 6.25 625 21.252.5 20 50 6.25 400 16.752.5 50/3 125/3 6.25 2500/9 13.252.5 100/7 250/7 6.25 10000/49 10.752 25 50 4 625 20.752 20 40 4 400 162 50/3 .

假如你的方程从最高次项到低次项的系数依次为a b c d e f ,则构建一个向量假设为p,则程序如下:p=[a b c d e f]roots(p),执行这两句之后你想要的结果就会显示在命令框里.

^syms x f =(x+1)*(x+2)*(x-3)*(x-4)*(x-6)>> simple(f) simplify:(x+1)*(x+2)*(x-3)*(x-4)*(x-6) radsimp:(x+1)*(x+2)*(x-3)*(x-4)*(x-6) combine(trig):x^5-10*x^4+17*x^3+64*x^2-108*x-144 方程为(x+1)*(x+2)*(x-3)*(x-4)*(x-6)=0 x^5-10*x^4+17*x^3+64*x^2-108*x-144=0

首先要明确,simulink的作用为求解常微分方程(组)!且这是他唯一的作用!也就是说偏微分方程在simulink中是无法求解的,需要其他工具或软件作为接口.dx=diff(a,n,dim)计算数值微分用dsolve(''表达式1','表达是',.,'表达式n','初始值','独立变量');也可以用ode23或者ode45函数计算.要查函数调用格式可以在matlab workspace输入>>help ode23即可.

syms a b c d s=solve( (50+2*a)^2==c^2+d^2, a^2 + b^2 == c^2, b^2==(50+a)*a, a*d==b*c) s.a,s.b,s.c,s.d % 解很多.全0就是一个解!其他一些(部分) [ 0, 0, 0, 50] [ 0, 0, 0, -50] [ -75, 25*3^(1/2), -50*3^(1/2), 50] [ -75, -25*3^(1/2), 50*3^(1/2), 50] [ -75, 25*3^(1/2), 50*3^(1/2), -50] [ -75, -25*3^(1/2), -50*3^(1/2), -50]

在t=solve('cos1=cos2')中,'cos1=cos2'被认为是一个表达式,里面的变量由solve函数自己负责解析出来,和变量cos1,cos2的值没有关系,而且solve函数在内部决定哪一个.

先编一个m文件:function y = myfun_1003(x);temp = 1;for n = 1:40,temp = temp + x^n/prod(1:n);endy = x^40/prod(1:40)*temp - 0.01;然后主窗口运行:[x,feval,flag] = fsolve(@myfun_1003, 10, optimset('Display', 'off'))x就是解,feval是y在x情况下的取值(理想应为零),flag为返回值,表示有解与否.搞定~