找回密码
 注册
Simdroid-非首页
查看: 166|回复: 6

[符号计算] 威布尔分布的累积分布推导不出来吗?

[复制链接]
发表于 2010-9-15 09:03:50 | 显示全部楼层 |阅读模式 来自 江苏南京
  1. fw=beta*(t-alpha)^(beta-1)/theta^beta*Exp[-((t-alpha)/theta)^beta]
  2. Fw=Integrate[fw,{t,alpha,x},Assumptions->beta>0&&alpha>0&&theta>0&&t>0&&x>0&&t>alpha]
复制代码


这里的Fw怎么推导没结果?
如何得到如图形所示的样子?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
 楼主| 发表于 2010-9-15 09:13:14 | 显示全部楼层 来自 江苏南京
Simdroid开发平台
推导复杂点的公式时经常难以获得自己想要的格式,很难看,不得不自己手工推导再试试。
但是不管是自己推得还是mathematica推的,两个形式不同,还要带入实数验证是否一致,巨麻烦。
回复 不支持

使用道具 举报

 楼主| 发表于 2010-9-15 22:58:38 | 显示全部楼层 来自 江苏南京
他的均值和方差,我都手工推出来了,就是不知道怎么用mathematica弄出来。

我推导对数正态分布也没遇到这些啊,大不了加几个假设就可以了;威布尔分布总是搞不定。
回复 不支持

使用道具 举报

发表于 2010-9-15 23:27:16 | 显示全部楼层 来自 台湾

  1. fw[t_,beta_,alpha_,theta_]:=beta*(t-alpha)^(beta-1)/theta^beta*Exp[-((t-alpha)/theta)^beta];
  2. (*A change in variables is required to evaluate the integral.  Let  y =
  3. (x-alpha)/theta. Then the integral is*)
  4. Integrate[\[Theta]*fw[\[Theta]*y+\[Alpha],\[Beta],\[Alpha],\[Theta]],{y,0,(
  5. x-\[Alpha])/\[Theta]},Assumptions->{\[Beta]>0,x>\[Alpha],\[Theta]>0}]
复制代码
回复 不支持

使用道具 举报

 楼主| 发表于 2010-9-16 00:25:26 | 显示全部楼层 来自 江苏南京
  1. fw[t_,beta_,alpha_,theta_]:=beta*(t-alpha)^(beta-1)/theta^beta*Exp[-((t-alpha)/theta)^beta];
  2. (*A change in variables is required to evaluate the integral.  Let  y =
  3. (x-alpha)/theta. Then the integral is*)
  4. Integrate[\[Theta]*fw[\[Theta]*y+\[Alpha],\[Beta],\[Alpha],\[Theta]],{y,0,(
  5. x-\[Alpha])/\[Theta]},Assumptions->{\[Beta]>0,x>\[Alpha],\[Theta]>0}]
复制代码

chungyuandye 发表于 2010-9-15 23:27


现在好了,万分感谢!!
经过chungyuandye一改,不但结果出来了,还很漂亮。
我用此法继续在mathematica推导了一阶原点矩,二阶原点矩和方差,
发现结果和我推导的完全一致,连gamma函数都能显示出来,这很重要。
再次感谢!

另外的问题,如果有时间麻烦帮助一下:
上边为什么要变换呢?
我在推导过程中怎么判断或决定要变换什么?
回复 不支持

使用道具 举报

发表于 2010-9-18 11:06:46 | 显示全部楼层 来自 上海
Read Documenation PDF and CDF
回复 不支持

使用道具 举报

 楼主| 发表于 2010-9-28 16:28:10 | 显示全部楼层 来自 江苏南京
Read Documenation PDF and CDF
FreddyMusic 发表于 2010-9-18 11:06


看了,不错。
当对威布尔分布变换后就不可行了。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-5-9 05:25 , Processed in 0.039831 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表