返回信息流有个方程是认知无线电里能量检测时求门限用的:
Pfa=[gammainc(u,th/2,'upper')]/gamma(u)
这里u的是样值的数目,是已知的假设为30吧,pfa也是已知的0.1,要求的是th
可是用slove和fslove都解不出来,这个用matlab该怎么弄才能得出th的值?
高手帮帮忙了!
[em18]
这是一条镜像帖。来源:北邮人论坛 / matlab / #6058同步于 2009/7/27
该镜像源已超过 30 天没有更新,可能在源站已被删除。
Matlab机器人发帖
求助:不完全的gamma函数的逆函数
wjnwxn
2009/7/27镜像同步10 回复
订阅后,新回复会通过你的通知中心匿名送达。
9 条回复
u太大了,gammainc(u, th/2, 'uppper')在gamma(u)面前太渺小了,趋于0,但是Pfa=0.1,很明显,在fsolve默认的搜索范围内的th无法满足Pfa=gammainc(u, th/2, 'upper')/gamma(u)的要求。
【 在 wjnwxn 的大作中提到: 】
: 有个方程是认知无线电里能量检测时求门限用的:
: Pfa=[gammainc(u,th/2,'upper')]/gamma(u)
: 这里u的是样值的数目,是已知的假设为30吧,pfa也是已知的0.1,要求的是th
: ...................
非常感谢!!
【 在 ReLive 的大作中提到: 】
: lamda = gaminv(1-Pf, u, 1) * 2;
: 用help里gammainc和gaminv说明里的函数原型可推出
[em18]
用Matlab验证发现这种算法不对。p=gamminc(x,a)=gamcdf(x,a,1),但是这个函数gaminv(p,a,b)是求gamcdf(x,a,b)中的参数x的逆函数。
即x=gaminv(p,a,b),还是没法求第二个参数a。
求高人继续指点!!!!!!
想请教各位技术大牛,如果电脑上已经安装了matlab 2009a,想升级成2012a(里面有新函数gammaincinv)的话,是必须重新安装呢,还是有什么简单的办法可以升级?或者有没有办法让2009a多认识一个新函数?谢谢!~
可以把新函数对应的文件直接复制到MATLAB 09a对应的路径下就好了。
【 在 joylioe 的大作中提到: 】
: 想请教各位技术大牛,如果电脑上已经安装了matlab 2009a,想升级成2012a(里面有新函数gammaincinv)的话,是必须重新安装呢,还是有什么简单的办法可以升级?或者有没有办法让2009a多认识一个新函数?谢谢!~