早教吧作业答案频道 -->数学-->
使用mathematica求解多元不等式整数解出错,tt={10.11,14.31,17.48,25.99};i=0;Do[m1=(tt[[i]]-d[[1]])*(tt[[i]]-d[[1]])/c[[1]]/c[[1]];Reduce[{m1-1≤39.1*k+16*o+14*n+10.8*b≤m1+1,k≥0,o≥0,n≥0,b≥0},{k,o,n
题目详情
使用mathematica求解多元不等式整数解出错,
tt = {10.11,14.31,17.48,25.99};
i = 0;
Do[m1 = (tt[[i]] - d[[1]])*(
tt[[i]] - d[[1]])/c[[1]]/c[[1]];
Reduce[{m1 - 1 ≤ 39.1*k + 16*o + 14*
n + 10.8*b ≤ m1 + 1,k ≥ 0,o ≥ 0,n ≥ 0,b ≥ 0},{k,
o,n,b},Integers],{i,1,4,1}]
这段代码是想分别用tt的五个元素,求解满足不等式的k,n,o,b的整数解,可是一执行语句就报错……为什么?
d[[1]]和c[[1]]是两个参数,在这里无关紧要,换成两个任意常数就行了,我这里是-2.846和3.294。tt这个列表里有四个元素,我的目的是在分别以这四个元素为参量的情况下求解我所给出的不等式的整数解,但是会报错。大家运行一下这段程序就明白了。提示信息的大意是mathematca无法精确求解。
话说这个其实昨天后来我也解决了。我在m1的外面加了个取整函数,又把整个不等式都乘了个10后就不会报出这个错了。(由此看来mathematica在整数规划方面并不是很强啊……)但现在产生了一个新的问题。
修改后的语句,若是把其中的Reduce部分挖出来单独运行,是可以得出结果的,但是若是放在Do里面,Reduce所得的结果就无法输出了……这是为什么?
tt = {10.11,14.31,17.48,25.99};
i = 0;
Do[m1 = (tt[[i]] - d[[1]])*(
tt[[i]] - d[[1]])/c[[1]]/c[[1]];
Reduce[{m1 - 1 ≤ 39.1*k + 16*o + 14*
n + 10.8*b ≤ m1 + 1,k ≥ 0,o ≥ 0,n ≥ 0,b ≥ 0},{k,
o,n,b},Integers],{i,1,4,1}]
这段代码是想分别用tt的五个元素,求解满足不等式的k,n,o,b的整数解,可是一执行语句就报错……为什么?
d[[1]]和c[[1]]是两个参数,在这里无关紧要,换成两个任意常数就行了,我这里是-2.846和3.294。tt这个列表里有四个元素,我的目的是在分别以这四个元素为参量的情况下求解我所给出的不等式的整数解,但是会报错。大家运行一下这段程序就明白了。提示信息的大意是mathematca无法精确求解。
话说这个其实昨天后来我也解决了。我在m1的外面加了个取整函数,又把整个不等式都乘了个10后就不会报出这个错了。(由此看来mathematica在整数规划方面并不是很强啊……)但现在产生了一个新的问题。
修改后的语句,若是把其中的Reduce部分挖出来单独运行,是可以得出结果的,但是若是放在Do里面,Reduce所得的结果就无法输出了……这是为什么?
▼优质解答
答案和解析
请给完整的程序,并且把题目描述清楚,报什么错误也要写出来啊
看了使用mathematica求解...的网友还看了以下:
平安集团在深圳福田CBD修建的平安大厦是深圳新的地标建筑,该地于2007年11月6日以16.568 2020-04-07 …
一下matlab程序错在哪clear allclclamada=0.3;theta=0.7;m=0 2020-05-16 …
风向中的16方位与1到16怎么对应?实验室有个自动气象站,对风向的记录为1到16这16个数字,请问 2020-05-19 …
(a+1)(a^2+1)(a^4+1)(a^8+1)(a^16+1)=(a-1)[(a+1)(a^ 2020-05-22 …
解三元一次方程组,1+1+D-E+F=0,1+16+D+4E+F=0,16+4+4D-2E+F=0 2020-08-03 …
一架飞机于北京时间5月1日8时从北京机场起飞,经16小时后到达伦敦机场,到达时当地时间(区时)为多 2020-08-04 …
命题p:∃x∈R,|x+1|+k<x,命题q:∀x>0,y>0,z>0>且x+y+z=1,有k≤1x 2020-10-31 …
(2011•怀化一模)命题p:∃x∈R,|x+1|+k<x,命题q:∀x>0,y>0,z>0>且x+ 2020-11-01 …
命题p:∃x∈R,|x+1|+k<x,命题q:∀x>0,y>0,z>0>且x+y+z=1,有k≤++ 2020-11-01 …
假若120°E是t日与t+1日的分界线(1)这时t日占全球范围的()A、大于1/2B、小于1/2C、 2020-12-05 …