资源描述：
《matlab求解方程组总结》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、Matlab解方程组解法1符号解法，即使用solve函数解法步骤：l第一步:定义变量syms x y z ...; l第二步:求解[x,y,z,...]=solve('eqn1','eqn2',...,'eqnN','var1','var2',...'varN'); l第三步:求出n位有效数字的数值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);看一个简单的例子2x+3=y3x-7=y[x,y]=solve('y=2*x+3','y=3*x-7');写成[xy]=solve('y-2*x-3',
2、'y-3*x+7')也可以网上可能有些代码在之前加了一句symsxy,表示定义符号变量x，y不加也是可以的。但是代码如果写成这样[xy]=solve(y-2*x-3,y-3*x+7)（无单引号括起来，则必须先说明x，y是符号变量）symsxy[xy]=solve(y-2*x-3,y-3*x+7)解二（多）元二（高）次方程组2+3y+1=0y2+4x+1=0[x,y]=solve('x^2+3*y+1=0','y^2+4*x+1=0'); x=vpa(x,4)%vpa(x,n)用来控制运算精度，这里表示保留4位有效数
3、字y=vpa(y,4)x=-0.2824-2.9881.635-3.029*i1.635+3.029*iy=-0.3599-3.3091.835+3.303*i1.835-3.303*i解法2采用左除运算解方程组x=Ab2x+3=y3x-7=y原方程可以变为2*x-y=-33*x-y=7代码：A=[2,-1;3,-1];系数矩阵B=[-3;7];X=AB解法3采用求逆运算解方程组x=inv(A)*b 从这段话，可以看出虽然可以使用法3来求解。但是从速度和精度来说，法2更好例1：12.56×A-311=-20.9
4、×B-311B311=P0.71933×P=A×20+BsymsABP[A,B,P]=solve('12.56*(A-311)=-20.9*(B-311)','B/311=P^0.71','933*P=A*20+B')输出参数顺序lsolve函数求解出的参数按字母表顺序排序symsxyzt[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的值，并把他们按字母表顺序排序，所
5、以第一个是t后三个一次为xyzsymsxyzt[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)已知参数求方程组例如求解该方程组x2+y2=1x+y=1[x,y]=solve('x^2+y^2=1','x+y=1')x=1y=0或x=0y=1修改该方程ax2+by2=cdx+fy=g此时方程加入了已知参数[x,y]=solve('a*x^2+b*y^2=c','d*x+e*y=f')此时求解出来的变量用已知参数表示。
6、如果已知参数值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([xy])(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个园的交点问题。仿真程序为：globalabcdev;a=1;b=0;c=-1;d=0;e=1
7、.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)带入参数数值后的结果
