问题详见知乎：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. Django xadmin后台添加富文本编辑器UEditor的用法

    效果图: 步骤: 1.利用命令:pip install DjangoUeditor,安装DjangoUeditor,但由于DjangoUeditor没有python3版本的,从的Github上把修改好 ...

  2. php libev扩展使用

    在WorkerMan源码分析 - 实现最简单的原型文章中提到了libevent网络库,前者和libev都是事件驱动编程库高性能.简单说libev对libevent做了改进和精简.libevent使用全 ...

  3. python学习-(__new__方法和单例模式)

    class Dog(object): __instance = None __init_flag = False def __new__(cls, name): if cls.__instance = ...

  4. C#中DllImport用法

    http://blog.csdn.net/u011981242/article/details/52622923 http://www.jb51.net/article/46384.htm 读取身份证 ...

  5. c# 运行大运算程序主窗体卡掉的解决

    写了一个运算过滤大文本的程序, 其中方法里边使用了多线程,并行线程等方法.  但主窗体控件直接使用此方法时,页面卡顿.所以主线程被堵塞. 代码如下, splitfile 这个方法运行时页面卡顿,阻塞了 ...

  6. essential c++ 第一章 array及vector相关使用

    将对象初始化: 1.用等号(=)赋值运算符初始化,针对对象是内置类型或者对象可以单一值初始化 2.构造函数初始化,针对对象需要多个初始值的情况 单括号括住的字符表示字符常量(‘ ’): 第一个反斜线表 ...

  7. PAT 1087 有多少不同的值（20）（STL—set）

    1087 有多少不同的值(20 分) 当自然数 n 依次取 1.2.3.--.N 时,算式 ⌊n/2⌋+⌊n/3⌋+⌊n/5⌋ 有多少个不同的值?(注:⌊x⌋ 为取整函数,表示不超过 x 的最大自然数 ...

  8. OSGI引入Spring DM实现对服务对象的管理

    一.异同 熟悉Spring的应该也都了解它的IOC的功能,那么对于在OSGI开发环境下,在使用IOC功能时有什么不同呢?最重要的一点就是Spring上下文对象,每个Spring-Powered Bun ...

  9. R语言读取Hive数据表

    R通过RJDBC包连接Hive 目前Hive集群是可以通过跳板机来访问 HiveServer, 将Hive 中的批量数据读入R环境,并进行后续的模型和算法运算. 1. 登录跳板机后需要首先在Linux ...

  10. linux学习第四天 (Linux就该这么学)2018年11月16日

    今天主要讲了 管道符,重写向与环境变量 输入输出重写向 标准输出重写向 (标准,覆盖,错误) > 将标准输出重写向到一个文件中 >> 追加到文件 2>错误输出重定向 2> ...