Matlab求解方程组概述.docx

Matlab解方程组解法1符号解法，即使用solve函数解法步骤：第一步:定义变量syms?x?y?z?...;?第二步:求解[x,y,z,...]=solve(eqn1,eqn2,...,eqnN,var1,var2,...varN);?第三步:求出n位有效数字的数值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);看一个简单的例子[x ,y]=solve(y=2*x+3,y=3*x-7);写成[x y]=solve(y-2*x-3,y-3*x+7)也可以网上可能有些代码在之前加了一句 syms x y,表示定义符号变量x，y 不加也是可以的。但是代码如果写成这样[x y]=solve(y-2*x-3,y-3*x+7)（无单引号括起来，则必须先说明x，y是符号变量）syms xy[x y]=solve(y-2*x-3,y-3*x+7)解二（多）元二（高）次方程组[x,y]=solve(x^2+3*y+1=0,y^2+4*x+1=0);?x=vpa(x,4)%vpa(x,n)用来控制运算精度，这里表示保留4位有效数字y=vpa(y,4)x = -0.2824 -2.988 1.635 - 3.029*i 1.635 + 3.029*iy = -0.3599 -3.309 1.835 + 3.303*i 1.835 - 3.303*i解法2采用左除运算解方程组x=A\b原方程可以变为2*x-y=-33*x-y=7代码：A=[2,-1;3,-1]; 系数矩阵B=[-3;7];X=A\B 解法3采用求逆运算解方程组x=inv(A)*b?从这段话，可以看出虽然可以使用法3来求解。但是从速度和精度来说，法2更好例1：syms A B P[A,B,P]= solve(12.56*(A-311)=-20.9*(B-311),B/311=P^0.71,933*P=A*20+B)输出参数顺序solve函数求解出的参数按字母表顺序排序syms x y z t[x,y,z,t]=solve(2*x+3*y-z+t-2,5*x+y+z-t-13,x-y+2*z+2*t-3,3*x+2*y+2*z+9*t+3)正确写法应该是： solve函数会求解出xyzt的值，并把他们按字母表顺序排序，所以第一个是t 后三个一次为xyzsyms x y z t[t,x,y,z]=solve(2*x+3*y-z+t-2,5*x+y+z-t-13,x-y+2*z+2*t-3,3*x+2*y+2*z+9*t+3)已知参数求方程组例如求解该方程组[x,y]=solve(x^2+y^2=1,x+y=1)x =1 y=0或x=0 y=1修改该方程此时方程加入了已知参数 [x,y]=solve(a*x^2+b*y^2=c,d*x+e*y=f)此时求解出来的变量用已知参数表示。如果已知参数值a=1;b=1;c=1;d=1;e=1;f=1;然后用subs函数来替换已知参数中的值，即计算变量的值。完整代码：[x,y]=solve(a*x^2+b*y^2=c,d*x+e*y=f);a=1;b=1;c=1;d=1;e=1;f=1;subs([x y])(a-x)^2+(b-y)^2=e^2(C-x)^2+(D-y)^2=v^2已知a,b,c,d,e,v 值求解 X,Y 请问用 matlab 如何写，就是求2个园的交点问题。仿真程序为：global a b c d e v;a=1;b=0;c=-1;d=0;e=1.5;v=1.5;%参数的值%可以不用展开[x,y]=solve(x^2+y^2-2*a*x-2*b*y=e^2-a^2-b^2,x^2+y^2-2*c*x-2*d*y=v^2-c^2-d^2)运行结果subs(x)subs(y)带入参数数值后的结果