|
|
【答案】应助回帖
aimforao(金币+100): 40 2011-10-21 08:56:00
方程可以转化成一下的线性形式：
![]()
计算步骤：
ie=[0.365 0.479 0.592 0.427 0.575 0.711 0.595 0.701 0.891]';
x1=[-7.022 -7.022 -7.022 -7.096 -7.096 -7.096 -6.648 -6.648 -6.648]';
x2=[-2.391 -2.391 -2.391 -2.524 -2.524 -2.524 -2.392 -2.392 -2.392]';
x3=[0.7963 0.7963 0.7963 1.9082 1.9082 1.9082 1.9277 1.9277 1.9277]';
x4=[106.444 106.444 106.444 98.851 98.851 98.851 135.257 135.257 135.257]';
x5=[10e-4 10e-3 10e-2 10e-4 10e-3 10e-2 10e-4 10e-3 10e-2]';
AA1 = x1.*x5;
AA2 = x2.*x5;
AA3 = x3.*x5;
AA4 = x4.*x5;
AA5 = x5;
AA = [AA1,AA2,AA3,AA4,AA5];
BB = 1./(1./ie-1);
x = AA\BB
计算结果：
Warning: Rank deficient, rank = 3, tol = 3.9865e-014.
> In ddddd at 15
x =
20.1351
0
19.7714
1.3244
0
>> [S,V,D]=svd(AA)
S =
-0.0053 -0.0032 0.0078 -0.9992 -0.0358 -0.0030 0.0033 0.0078 0.0071
-0.0534 -0.0319 0.0777 -0.0265 0.7190 0.0307 0.0563 0.3172 -0.6064
-0.5338 -0.3186 0.7769 0.0126 -0.0715 -0.0030 -0.0057 -0.0318 0.0606
-0.0050 -0.0062 -0.0060 0.0048 -0.1155 -0.9639 -0.0131 -0.0759 -0.2269
-0.0496 -0.0623 -0.0596 -0.0035 -0.2080 0.2283 -0.1128 -0.7088 -0.6161
-0.4960 -0.6232 -0.5963 0.0003 0.0220 -0.0132 0.0114 0.0716 0.0639
-0.0068 0.0071 -0.0018 0.0061 -0.1055 0.0196 0.9900 -0.0633 -0.0656
-0.0678 0.0707 -0.0176 0.0240 -0.6358 0.1306 -0.0609 0.6143 -0.4327
-0.6776 0.7072 -0.1755 -0.0025 0.0646 -0.0133 -0.0038 -0.0608 0.0439
V =
19.9909 0 0 0 0
0 0.2101 0 0 0
0 0 0.0871 0 0
0 0 0 0.0000 0
0 0 0 0 0.0000
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
D =
0.0595 0.9413 -0.0662 0.2992 -0.1282
0.0210 0.3092 0.0781 -0.6744 0.6656
-0.0135 -0.0383 -0.9947 -0.0690 0.0651
-0.9979 0.0641 0.0111 -0.0012 -0.0009
-0.0086 -0.1128 0.0059 0.6715 0.7323
>> format long
>> AA
AA =
-0.007022000000000 -0.002391000000000 0.000796300000000 0.106444000000000 0.001000000000000
-0.070220000000000 -0.023910000000000 0.007963000000000 1.064440000000000 0.010000000000000
-0.702200000000000 -0.239100000000000 0.079630000000000 10.644400000000001 0.100000000000000
-0.007096000000000 -0.002524000000000 0.001908200000000 0.098851000000000 0.001000000000000
-0.070960000000000 -0.025240000000000 0.019082000000000 0.988510000000000 0.010000000000000
-0.709600000000000 -0.252400000000000 0.190820000000000 9.885100000000001 0.100000000000000
-0.006648000000000 -0.002392000000000 0.001927700000000 0.135257000000000 0.001000000000000
-0.066480000000000 -0.023920000000000 0.019277000000000 1.352570000000000 0.010000000000000
-0.664800000000000 -0.239200000000000 0.192770000000000 13.525700000000001 0.100000000000000
可见系数矩阵的秩为3. |
|
回复此楼