用mathematica求六元一次方程组且方程个数比变量个数少一个
问题详见知乎: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求六元一次方程组且方程个数比变量个数少一个的更多相关文章
- 【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 ...
- math & 三元一次方程组的解法
math & 三元一次方程组的解法 class 6 math 例题 问题: 1. 已经做好的与没有做好的比例是 5 比 7; 2 再做好51,完成总数的 70%; 3. 问,一共要做多少朵花? ...
- 求一个n元一次方程的解,Gauss消元
求一个n元一次方程的解,Gauss消元 const Matrix=require('./Matrix.js') /*Gauss 消元 传入一个矩阵,传出结果 */ function Gauss(mat ...
- Java基础_循环嵌套_打印乘法口诀、菱形,各种图形,计算二元一次和三元一次方程组_7
循环嵌套 打印乘法口诀 for(int j=1;j<=9;j++){ for(int i=1;i<=j;i++){ System.out.print(i+"*"+j+& ...
- Java开发中经典的小实例-( 鸡蛋0.1元一个,鸭蛋3元一个,鹅蛋6元一个。求一百元买一百个蛋。)
public class Test24 { public static void main(String[] args) { // 鸡蛋0.1元一个,鸭蛋3元一个,鹅蛋6元一个.求 ...
- 算法导论(第三版) Exercises4.2(求最大和子数组的算法优化过程)
4.1-1 如所有元素都为负,则返回所有元素中最大的负数. 4.1-2(暴力法求最大和子数组) struct subarray { int start, end, sum; }; void brute ...
- XJTUOJ wmq的队伍(树状数组求 K 元逆序对)
题目链接:http://oj.xjtuacm.com/problem/14/[分析]二元的逆序对应该都会求,可以用树状数组.这个题要求K元,我们可以看成二元的.我们先从后往前求二元逆序对数, 然后对于 ...
- 求数组的子数组之和的最大值II
这次在求数组的子数组之和的最大值的条件下又增加了新的约束: 1.要求数组从文件读取. 2.如果输入的数组很大, 并且有很多大的数字, 就会产生比较大的结果 (考虑一下数的溢出), 请保 ...
- 【c语言】实现一个函数,求字符串的长度,不同意创建第三方变量
// 实现一个函数,求字符串的长度.不同意创建第三方变量. #include <stdio.h> #include <assert.h> int my_strlen_no(ch ...
随机推荐
- Linux 线程】线程同步《四》
1.信号量 (1)概念 信号量和互斥锁(mutex)的区别:互斥锁只允许一个线程进入临界区,而信号量允许多个线程同时进入临界区. 不多做解释,要使用信号量同步,需要包含头文件semaphore.h. ...
- tableView与导航栏的偏移问题
//方法1.关闭自动布局 self.automaticallyAdjustsScrollViewInsets = NO; //方法2.改变相对于scrollview的位置 //contentInset ...
- 26-算法训练 Torry的困惑(基本型) 素数打表
算法训练 Torry的困惑(基本型) 时间限制:1.0s 内存限制:512.0MB 问题描述 Torry从小喜爱数学.一天,老师告诉他,像2.3.5.7……这样的数叫做质数.To ...
- css常见问题一
[1]禁止换行.class {word-break:keep-all;white-space:nowrap;}[2]强制换行.class{word-break:break-all;}普通容器中(Div ...
- C#实现发送给QQ邮件
最近在做一个通过点击忘记密码往用户邮箱中发邮件(邮件内容是一个超链接)点击进行修改的功能,发送原理,我们只是把邮件发送给smtp服务器,然后再由smtp服务器发送到邮箱,发送之前要校验一下. 1.微软 ...
- 一行代码轻松搞定各种IE兼容问题,IE6,IE7,IE8,IE9,IE10
在网站开发中不免因为各种兼容问题苦恼,针对兼容问题,其实IE给出了解决方案Google也给出了解决方案百度也应用了这种方案去解决IE的兼容问题? 百度源代码如下: <!Doctype html& ...
- XStream将XML转javaben,出现多余的tag,导致出错
今天在测试银联无卡快捷支付的案例时,多了一个多tag兼容性测试,它是指银联的XML报文中会出现多余的tag,如果我们用XStream解析的时候,没有Javabean的字段可以对应上,就会报错!提示: ...
- nested exception is java.lang.NoClassDefFoundError: org/springframework/dao/support/PersistenceExceptionTranslator
该问题是少了一个spring-tx-的jar包,把该包加入到buildpath中就行了. 参考链接:http://blog.csdn.net/Rongbo_J/article/details/4666 ...
- java中 this 关键字的三种用法
Java中this的三种用法 调用属性 (1)this可以调用本类中的任何成员变量 调用方法(可省略) (2)this调用本类中的成员方法(在main方法里面没有办法通过this调用) 调用构造方法 ...
- PC 上的 LVM 灾难修复
LVM 介绍 LVM 简介 LVM 是逻辑盘卷管理(Logical Volume Manager)的简称,最早是 IBM 为 AIX 研发的存储管理机制.LVM 通过在硬盘和分区之间建立一个逻辑层,可 ...