求用matlab 解方程组的方法.
\u6c42\u52a9\uff01\uff01\uff01 \u7528MATLAB\u6c42\u4e0b\u9762\u7684\u7684\u65b9\u7a0b\u7ec4\u3002~~\u6025~~A = [7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13];
b = [4,7,-1,0];
inv(A)*b'
ans =
0.4979
0.1445
0.0629
-0.0813
clear
clc
syms iq ir p1 q1 u1 r x
f1=((p1-(p1^2+q1^2)*r/u1^2)/sqrt((u1+(p1*r+q1*x)/u1)^2+(p1*x-q1*r)^2/u1^2))-iq;
f2=((q1-(p1^2+q1^2)*x/u1^2)/sqrt((u1+(p1*r+q1*x)/u1)^2+(p1*x-q1*r)^2/u1^2))-ir;
[p1,q1]= solve(f1,f2,p1,q1)
\u4f60\u8981\u7684\u89e3\u6ca1\u6709\u610f\u4e49\u3002
\u4e24\u4e2a\u89e3\u7684\u89e3\u6790\u8868\u8fbe\u5f0f\u5728word\u91cc\u5360\u4e86a4\u7eb8100\u591a\u9875\u3002
\u5982\u679c\u4ee4iq=ir=0,\u89e3\u5c31\u7b80\u5355\u7684\u591a\u3002
B=[-3;-6;0;2]
X=A\B
%或
X=inv(A)*B
%验算
b=A*X
结果:
A =
2 -2 0 1
2 3 1 -3
3 4 -1 2
1 3 1 0
B =
-3
-6
0
2
X =
-1.8052
0.7273
1.6234
2.0649
X =
-1.8052
0.7273
1.6234
2.0649
b =
-3.0000
-6.0000
0.0000
2.0000
用矩阵写,就是
[2 -2 0 1
2 3 1 -3
3 4 -1 2
1 3 1 0]×[x1 x2 x3 x4]'=[-3 -6 0 2]'
记为A×[x1 x2 x3 x4]'=C
则
[x1 x2 x3 x4]'=inv(A)*C
ans =
-1.8052
0.7273
1.6234
2.0649
扩展阅读:matlab jacobi解方程组 ... matlab解复杂方程fsolve ... matlab解x y方程组 ... matlab solve解方程组 ... matlab求解多元方程组 ... 怎么用matlab解多元方程 ... matlab求方程组的唯一解 ... matlab 解微分方程 ... matlab求方程的全部根 ...