查看: 393|回复: 2|关注: 0
打印 上一主题 下一主题

快乐扑克遗漏数据: [未答复] 建立行列式等于0的矩阵方程组求解

[复制链接]

新手

5 麦片

财富积分


快乐扑克3如何跟对子 www.vqxik.tw 050


1

主题

3

帖子

0

最佳答案
跳转到指定楼层
1#
发表于 2019-10-9 12:16:27 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
先感谢各位大佬帮忙解答。

我想用matlab写一个方程组,目的是另矩阵行列式等于零,求解矩阵中未知数


我想求解的方程为  det(A+M-lambda(1,i)*E(9) )=0
也就是另A+M的和,减去lambda中第i个元素为系数乘以单位矩阵,的行列式等于零
以此可以列出9个方程组,再求解方程
我在建立方程组的一步就出现的报错,另外想请教求解此方程我应该用fslove函数吗

(1)其中A矩阵为9*9已知矩阵
     A =

    0.6658    1.1747   -1.9918   -0.1047   -0.0040    0.0506    2.0060   -1.1561    0.3213
   -0.5632    4.4128   -3.7787   -0.2729   -0.0358    0.1784    3.8157   -3.4099    0.5803
   -0.7129    4.6949   -3.6406   -0.3889   -0.0621    0.2774    4.6721   -4.6744    0.7428
   -0.2815    2.0190   -3.1366    0.7842   -0.0116    0.1174    3.1608   -1.9977    0.2894
    0.2151    0.1911   -1.9218   -0.0277    1.0029   -0.0116    1.9243   -0.1674   -0.2328
   -0.1915    2.0166   -3.2150   -0.1810   -0.0074    1.0873    3.2400   -1.9984    0.1927
   -0.7178    4.7048   -4.6400   -0.3902   -0.0621    0.2783    5.6717   -4.6841    0.7479
   -0.5689    3.4643   -3.7867   -0.2747   -0.0364    0.1801    3.8241   -2.4618    0.5863
   -0.3004    1.1419   -2.0133   -0.1037   -0.0033    0.0490    2.0271   -1.1227    1.2873

(2) M矩阵为对角元素为m11~m99 共九个未知数的对角矩阵

(3) lambda为1*9,共9个元素的矩阵

     lambda =

    0.9600    0.9600    0.9680    0.9720    0.9810    0.9880    0.9900    0.9980    0.9980



以下为我的程序

lambda=[0.96 0.96 0.968 0.972 0.981 0.988 0.99 0.998 0.998]%扰动后;

syms m11 m22 m33 m44 m55 m66 m77 m88 m99;
n=9;
for i=1:n
    for j=1:n
        if i==j
            cmd = sprintf('sym(''m%i%i'')',i,j);
            M(i,j) = eval(cmd);
        end
    end
end %建立对角元素为未知数的M矩阵

A=xlsread('A矩阵.xlsx')%插入A矩阵,A矩阵在附件中

SUM=A+M

for i=1:n
    det(SUM-lambda(1,i).*eye(n))=0
end %建立方程组



报错
错误使用 sym/subsindex (line 845)
Invalid indexing or function definition. Indexing must follow MATLAB indexing. Function arguments must be symbolic variables, and function
body must be sym expression.


出错 findM (line 20)
    det(SUM-lambda(1,i).*eye(n))=0



新手

5 麦片

财富积分


050


1

主题

3

帖子

0

最佳答案
2#
 楼主| 发表于 2019-10-9 13:07:18 | 只看该作者
方程截图

A矩阵.xlsx

9.31 KB, 下载次数: 0

新手

5 麦片

财富积分


050


1

主题

3

帖子

0

最佳答案
3#
 楼主| 发表于 2019-10-9 13:08:22 | 只看该作者

截图中是4*4矩阵,我的是9*9矩阵,一共九个方程
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 快乐扑克3如何跟对子 返回列表