Matlab求解方程組總結(jié)

Matlab解方程組解法1符號解法，即使用solve函數(shù)解法步驟：第一步:定義變量syms

x

y

z

...;

第二步:求解[x,y,z,...]=solve('eqn1','eqn2',...,'eqnN','var1','var2',...'varN');

第三步:求出n位有效數(shù)字的數(shù)值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);看一個簡單的例子2x+3=y[x,y]=solve('y=2*x+3','y=3*x-7');寫成[xy]=solve('y-2*x-3','y-3*x+7')也可以網(wǎng)上可能有些代碼在之前加了一句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=0[x,y]=solve('x^2+3*y+1=0','y^2+4*x+1=0');

x=vpa(x,4)%vpa(x,n)用來控制運(yùn)算精度，這里表示保留4位有效數(shù)字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采用左除運(yùn)算解方程組x=A\b2x+3=y原方程可以變?yōu)?*x-y=-33*x-y=7代碼：A=[2,-1;3,-1];系數(shù)矩陣B=[-3;7];X=A\B解法3采用求逆運(yùn)算解方程組x=inv(A)*b

從這段話，可以看出雖然可以使用法3來求解。但是從速度和精度來說，法2更好例1：12.56×symsABP[A,B,P]=solve('12.56*(A-311)=-20.9*(B-311)','B/311=P^0.71','933*P=A*20+B')輸出參數(shù)順序solve函數(shù)求解出的參數(shù)按字母表順序排序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)正確寫法應(yīng)該是：solve函數(shù)會求解出xyzt的值，并把他們按字母表順序排序，所以第一個是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)已知參數(shù)求方程組例如求解該方程組x[x,y]=solve('x^2+y^2=1','x+y=1')x=1y=0或x=0y=1修改該方程ax此時方程加入了已知參數(shù)[x,y]=solve('a*x^2+b*y^2=c','d*x+e*y=f')此時求解出來的變量用已知參數(shù)表示。如果已知參數(shù)值a=1;b=1;c=1;d=1;e=1;f=1;然后用subs函數(shù)來替換已知參數(shù)中的值，即計(jì)算變量的值。完整代碼：[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個園的交點(diǎn)問題。仿真程序?yàn)椋篻lobalabcdev;a=1;b=0;c=-1;d=0;e=1.5;v=1.5;%參數(shù)的值%可以不用展開[x,y]=solve('x^2+y^2-2*