谁会解这个矩阵方程组
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/05 22:01:04
谁会解这个矩阵方程组
谁会解这个矩阵方程组
谁会解这个矩阵方程组
(Ps Pf PS PF)= 0.653222 0.00434825 0.340162 0.00226754
采用MATLAB解,很简单,把第二个方程合并和第一个矩阵合并,直调用线性代数方程求解函数就可以了
A =
-0.0023 0.0313 0.0040 0
0.0002 -0.0333 0 0.0040
0.0021 0 -0.0044 0.0600
0 0.0021 0.0004 -0.0640
1.0000 1.0000 1.0000 1.0000
C=[0 0 0 0 1]
C=conj(C')
C =
0
0
0
0
1
>> P=linsolve(A,C)
P =
0.6532
0.0043
0.3402
0.0023
>> vpa(P,6)
ans =
.653222
.434825e-2
.340162
.226754e-2
很明显嘛,都等于0.25
观察一下很容易出来的
问题很简单,计算量大些而已,用计算机会比较方便些
通过MATLAB随机产生p1,p2,p3,p4来求解符合条的P向量,附上程序
clear all;
clc;
A=[-0.002291 0.000208 0.002083 0;
0.031250 -0.033333 0 0.002083;
0.004000 0 -0.004400 0.000400;
0 0.004000 0.060000 -0.06...
全部展开
通过MATLAB随机产生p1,p2,p3,p4来求解符合条的P向量,附上程序
clear all;
clc;
A=[-0.002291 0.000208 0.002083 0;
0.031250 -0.033333 0 0.002083;
0.004000 0 -0.004400 0.000400;
0 0.004000 0.060000 -0.064000];
%随机产生一组和为1的P
n=4;s=1;
while s>=1
u=rand(n-1,1)/(n-1);
s=sum(u);
end
u(n)=1-sum(u);
P=[u(1) u(2) u(3) u(4)];
C=P*A;
K=1000;
if C>0
K=K-1;
elseif C<0
K=K-1;
else
P
end
结果不唯一,任举一例
[0.092832739622349 0.182293839734995 0.319168945144766 0.405704475497890]
收起