问题详见知乎:https://www.zhihu.com/question/68000713

我的问题:有5个方程,6个变量,其实我是想求出来de1=(系数)*dS1的形式,系数有Cij组成,Cij为常数。怎样求?具体方程下面代码可以看出(sys1里时方程,var1里时变量,b是一个大于等于0小于等于1常数)

我是用mathematica求的,参考http://reference.wolfram.com/language/tutorial/SolvingEquations.html

Eliminate[{de1-C11*dS1-C12*dS2-C13*dS3==0,de2-C21*dS1-C22*dS2-C23*dS3==0,de3-C31*dS1-C32*dS2-C33*dS3==0,de1+de2+de3==0,(dS1-dS3)*b-(dS2-dS3)==0},{de2,de3,dS2,dS3}]

Eliminate[{de1 - C11*dS1 - C12*dS2 - C13*dS3 == 0,
de2 - C21*dS1 - C22*dS2 - C23*dS3 == 0,
de3 - C31*dS1 - C32*dS2 - C33*dS3 == 0, (de1 + de2 + de3)/de1 ==
0.5, (dS1 - dS3)*b - (dS2 - dS3) == 0}, {de2, de3, dS2, dS3}] Solve[C11 (2.` C22 - 2.` b C22 + 2.` C23 + 2.` C32 - 2.` b C32 +
2.` C33) dS1 ==
C12 de1 - 1.` b C12 de1 + C13 de1 + 2.` C22 de1 - 2.` b C22 de1 +
2.` C23 de1 + 2.` C32 de1 - 2.` b C32 de1 + 2.` C33 de1 +
2.` C12 C21 dS1 - 2.` b C12 C21 dS1 + 2.` C13 C21 dS1 +
2.` b C13 C22 dS1 - 2.` b C12 C23 dS1 + 2.` C12 C31 dS1 -
2.` b C12 C31 dS1 + 2.` C13 C31 dS1 + 2.` b C13 C32 dS1 -
2.` b C12 C33 dS1, {de1, dS1}]

{{dS1 -> 0. - (0.5 (-1. C12 + 1. b C12 - 1. C13 - 2. C22 + 2. b C22 - 2. C23 - 2. C32 + 2. b C32 - 2. C33) de1)/(-1. C12 C21 + 1. b C12 C21 - 1. C13 C21 + 1. C11 C22 - 1. b C11 C22 - 1. b C13 C22 + 1. C11 C23 + 1. b C12 C23 - 1. C12 C31 + 1. b C12 C31 - 1. C13 C31 + 1. C11 C32 - 1. b C11 C32 - 1. b C13 C32 + 1. C11 C33 + 1. b C12 C33)}}

用mathematica求六元一次方程组且方程个数比变量个数少一个的更多相关文章

  1. 【Java例题】4.3 3. 使用Gauss消元法求解n元一次方程组的根,

    3. 使用Gauss消元法求解n元一次方程组的根,举例,三元一次方程组:0.729x1+0.81x2+0.9x3=0.6867x1+x2+x3=0.83381.331x1+1.21x2+1.1x3=1 ...

  2. math & 三元一次方程组的解法

    math & 三元一次方程组的解法 class 6 math 例题 问题: 1. 已经做好的与没有做好的比例是 5 比 7; 2 再做好51,完成总数的 70%; 3. 问,一共要做多少朵花? ...

  3. 求一个n元一次方程的解,Gauss消元

    求一个n元一次方程的解,Gauss消元 const Matrix=require('./Matrix.js') /*Gauss 消元 传入一个矩阵,传出结果 */ function Gauss(mat ...

  4. Java基础_循环嵌套_打印乘法口诀、菱形,各种图形,计算二元一次和三元一次方程组_7

    循环嵌套 打印乘法口诀 for(int j=1;j<=9;j++){ for(int i=1;i<=j;i++){ System.out.print(i+"*"+j+& ...

  5. Java开发中经典的小实例-( 鸡蛋0.1元一个,鸭蛋3元一个,鹅蛋6元一个。求一百元买一百个蛋。)

    public class Test24 {    public static void main(String[] args) {        // 鸡蛋0.1元一个,鸭蛋3元一个,鹅蛋6元一个.求 ...

  6. 算法导论(第三版) Exercises4.2(求最大和子数组的算法优化过程)

    4.1-1 如所有元素都为负,则返回所有元素中最大的负数. 4.1-2(暴力法求最大和子数组) struct subarray { int start, end, sum; }; void brute ...

  7. XJTUOJ wmq的队伍(树状数组求 K 元逆序对)

    题目链接:http://oj.xjtuacm.com/problem/14/[分析]二元的逆序对应该都会求,可以用树状数组.这个题要求K元,我们可以看成二元的.我们先从后往前求二元逆序对数, 然后对于 ...

  8. 求数组的子数组之和的最大值II

    这次在求数组的子数组之和的最大值的条件下又增加了新的约束:  1.要求数组从文件读取.      2.如果输入的数组很大,  并且有很多大的数字,  就会产生比较大的结果 (考虑一下数的溢出), 请保 ...

  9. 【c语言】实现一个函数,求字符串的长度,不同意创建第三方变量

    // 实现一个函数,求字符串的长度.不同意创建第三方变量. #include <stdio.h> #include <assert.h> int my_strlen_no(ch ...

随机推荐

  1. mybatis插入数据并获取主键值

    有时候我们的主键是自增的,但是我们想要在插入一条数据以后获取这条数据的主键值,而我们知道,mybatis执行完插入操作以后返回的是生效的记录数.那如何才能获取这个主键值呢. 1.在配置文件mapper ...

  2. XML文件的小结

    1.Layer-list <?xml version="1.0" encoding="utf-8"?><layer-list xmlns:an ...

  3. 举个通俗易懂的例子告诉你IAAS,SAAS,PAAS的区别

    作者:何足道链接:https://www.zhihu.com/question/21641778/answer/62523535来源:知乎著作权归作者所有.商业转载请联系作者获得授权,非商业转载请注明 ...

  4. python之virtualenv

    一 virtualenv简介 virtualenv------用来建立一个虚拟的python环境,一个专属于项目的python环境.用virtualenv 来保持一个干净的环境非常有用 在开发Pyth ...

  5. oracle 11g用exp无法导出空表的解决方案

    racle 11g中有个新特性,当表无数据时,不分配segment,以节省空间,当我们用exp导出空表时,无法导出. 解决方法是两个方面, 一是处理现有的空表,让其能导出: 二是设置参数,让后续的新的 ...

  6. 组委会正在为美团点评CodeM大赛的决赛设计新赛制

    比赛有 n 个人参加(其中 n 为2的幂),每个参赛者根据资格赛和预赛.复赛的成绩,会有不同的积分.比赛采取锦标赛赛制,分轮次进行,设某一轮有 m 个人参加,那么参赛者会被分为 m/2 组,每组恰好 ...

  7. 解决JFinal多文件上传时只获取到第一个文件名

    我的思路: 用户生成时随即生成一串随机字符作为该用户的文件上传目录,并保存该字符串到用户的某一字段.需要显示上传的附件时,遍历这个文件夹.上传时可把文件名设置为上传时间. 1.生成上传路径 可写在注册 ...

  8. 为Linux虚拟机设置网络

    安装虚拟机的时候为了使用方便我们除了需要设置静态ip为了能够让虚拟机也能够上网我们需要设置虚拟机网络 当然也可以使用虚拟机和主机共享上网,这个比较简单,这里就不说了,现在我们来通过桥接的方式为虚拟机设 ...

  9. 自己动手开启QUIC(转载)

    源:https://www.bennythink.com/quic.html#title-0 今天在推上偶然发现 Google 在自己的服务器上启用了 QUIC,QUIC 这东西嘛(发音同 quick ...

  10. DEDE 首页调用指定栏目链接的代码

    {dede:type typeid='6'} <a href='[field:typelink /]' target="_blank" >更多</a> {/ ...