[交流] Matlab: 实现任意个方程之间求解。

10楼: Originally posted by csgt0 at 2012-07-10 13:37:49
就是解方程组啊，肯定可以，但是有的不一定能求出解析解。

10楼: Originally posted by csgt0 at 2012-07-10 13:37:49
就是解方程组啊，肯定可以，但是有的不一定能求出解析解。

11楼: Originally posted by daicong at 2012-07-11 11:28:42
当然是解方程，问题是如何通过编程自动解决。
我只要输入方程即可。
然后给出一系列的可行解。...

19楼: Originally posted by hhlxg at 2012-07-11 00:08:05
求解有限个的方程做过，像你描述的这个有难度，关键是程序能自己识别的有限。

18楼: Originally posted by baobiao007 at 2012-07-10 23:25:03
这属于某种反演问题了，属于专门的研究课题了。不是三下五除二能搞定的。

20楼: Originally posted by daicong at 2012-07-11 22:44:03
目前我只做线性方程，比如：
A*x=b,方程组用矩阵表示。A,b为已知量，x是需要求的。
例如(举一个例子)：
5*x1+3*x2-4*x3=5;
6*x1+3*x2-2*x3=3;
13*x1+2*x2-9*x3=13;
76*x1+43*x2-24*x3=32;
24*x1+43*x2-23*x ...

24楼: Originally posted by hhlxg at 2012-07-11 10:43:03
如果是你描述的是线性方程组的话
用矩阵求解就可以了
A*x=b
x='b
'为矩阵A的逆矩阵
x=
>> A=

A =

     5     3    -4
     6     6    -2
    13     2    -9

>> b=

b =

      ...

26楼: Originally posted by wanghao9999 at 2012-07-12 03:51:58
应该能解出，但是有的不一定能求出解析解，能否具体描述下问题

28楼: Originally posted by userhung at 2012-07-12 09:53:04
解方程组哦~~~~~~~~~~~~~~~~~~~~`

27楼: Originally posted by daicong at 2012-07-13 00:03:19
目前我只做线性方程，比如：
A*x=b,方程组用矩阵表示。A,b为已知量，x是需要求的。
例如(举一个例子)：
5*x1+3*x2-4*x3=5;
6*x1+3*x2-2*x3=3;
13*x1+2*x2-9*x3=13;
76*x1+43*x2-24*x3=32;
24*x1+43*x2-23*x ...

27楼: Originally posted by daicong at 2012-07-13 00:03:19
目前我只做线性方程，比如：
A*x=b,方程组用矩阵表示。A,b为已知量，x是需要求的。
例如(举一个例子)：
5*x1+3*x2-4*x3=5;
6*x1+3*x2-2*x3=3;
13*x1+2*x2-9*x3=13;
76*x1+43*x2-24*x3=32;
24*x1+43*x2-23*x ...

31楼: Originally posted by wangww2011 at 2012-07-13 03:23:44
你这个问题不需要什么算法，其实就是组合呗，matlab里有个函数combnk
http://www.mathworks.co.uk/help/toolbox/stats/combnk.html
然后好像你的问题就能解决了...

37楼: Originally posted by csgt0 at 2012-07-15 17:52:20
大哥啊，最简单的方法你可以给每个方程一个序号，然后从这些序号中随机选就可以了。更方便点的，你构成一个m*n+1的矩阵，从其中的m行中任选k行进行行列式计算就可以了。...

38楼: Originally posted by daicong at 2012-07-17 01:28:37
这个函数确实对我很有帮助，不过即使这样对整个程序来说还需要很复杂的过程。...

A=[5 3 -4;
    6 3 -2;
    13 2 -9;
    76 43 -24;
    24 43 -23];

b=[5;3;13;32;43];


indices=combnk(1:5,3);

for i=1:length(indices)
    index=indices(i,:);
    B=A(index,:)
    if abs(det(B))>1e-12
        res=B\b(index)
    else
        det(B)
    end
   
end

40楼: Originally posted by wangww2011 at 2012-07-16 11:17:42
忍不住给你写了

A=;

b=;


indices=combnk(1:5,3);

for i=1:length(indices)
    index=indices(i,;
    B=A(index,
    if abs(det(B))>1e-12
        res=B\b(index)
    else
        ...

41楼: Originally posted by daicong at 2012-07-17 22:24:26
谢谢你写的代码。但我运行之后发现结果与我想象的不一样：按道理每三个方程联合求解可以得到一个解，那么这五个方程就应该有十个，当然也有可能有的三个方程没有解；但运行之后的结果只有一个解。请问这个是否有问 ...

42楼: Originally posted by wangww2011 at 2012-07-17 03:22:13
额 不会吧 明明有10个解的啊 至少在我电脑上就是10个解啊 matlab不可能给用户提供不一样的用户体验啊...

42楼: Originally posted by wangww2011 at 2012-07-17 03:22:13
额 不会吧 明明有10个解的啊 至少在我电脑上就是10个解啊 matlab不可能给用户提供不一样的用户体验啊...

44楼: Originally posted by daicong at 2012-07-18 02:18:13
对了，我还有个很长的代码。是计算顶点的，和这个有异曲同工之处，但那个更加复杂。我给你个链接，看你能否帮忙看看。谢谢。
http://muchong.com/bbs/viewthread.php?tid=4737207&pid=9&page=1#pid9...