求stata求解高次方程的命令

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

[tr=rgb(255, 252, 247)][/tr]

正以stata为工具学习统计学， 其中在平均发展速度计算中有一个高次方程，如 x+x2+x3+x4+x5=100 其中X旁边的数字，都是2次项，3次项，4次方，5次方的意思 求解x 在汉米尔顿常见命令和函数中没有找到对应的命令，请教一下

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：Stata tata 高次方 平均发展 发展速度 工具 统计学

Stata不适合用来处理这问题。 建议楼主还是用Mathematica或Matlab…
不过到是可以使用Stata绘图，很粗浅地看一下x在实数上的可能解…

graph tw function x^5+x^4+x^3+x^2+x-100,range(-3 3) yline(0, lstyle(foreground))

函数与水平线交叉的那一点对到的x 即x在实数上的可能解

时间有限，写了个简单的代码，可以初步算出大概的结果：

clear
set obs 1000000
gen d = .
gen x = .
local j = 1
forvalues x = 2(0.000001)3{
  qui replace x = `x' in `j'
  local d = `x' +`x'^2 +`x'^3 +`x'^4 +`x'^5 - 100
  qui replace d = abs(`d') in `j++'
  *dis in y "." _c
}
sort d
br
sort x
line d x   //有点动画的效果

易知f(x)=x^5+x^4+x^3+x^2+x-100在正实数区间上是增函数，且其零解属于(2.239, 2.24)。可使用mata：

mata
for (i=2.239;i<=2.24;i=i+.0000001) {
if (abs(i+i^2+i^3+i^4+i^5-100)<0.00001) i
}
end

sungmoo 发表于 2012-12-1 08:25
易知f(x)=x^5+x^4+x^3+x^2+x-100在正实数区间上是增函数，且其零解属于(2.239, 2.24)。可使用mata：

mat ...

谢谢您的热心回复。

都是大牛

sungmoo 发表于 2012-12-1 08:25
易知f(x)=x^5+x^4+x^3+x^2+x-100在正实数区间上是增函数，且其零解属于(2.239, 2.24)。可使用mata：

mat ...

感谢您对我使用stata解高次方程命令的回答，从您给出的mata命令中，关键是2.239 和2.24两个值的确定。进一步请问这两个值如何求出的呢？

高手啊！