用matlab求解一个两重积分方程(未知数在积分下限,含复数积分)这实际上是我求出的一个特征函数,我希望通过这个特征函数求出随机变量的PDF,并由此解决一个给定概率的方程.如图所示:理
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/02 05:34:17
![用matlab求解一个两重积分方程(未知数在积分下限,含复数积分)这实际上是我求出的一个特征函数,我希望通过这个特征函数求出随机变量的PDF,并由此解决一个给定概率的方程.如图所示:理](/uploads/image/z/6838939-19-9.jpg?t=%E7%94%A8matlab%E6%B1%82%E8%A7%A3%E4%B8%80%E4%B8%AA%E4%B8%A4%E9%87%8D%E7%A7%AF%E5%88%86%E6%96%B9%E7%A8%8B%EF%BC%88%E6%9C%AA%E7%9F%A5%E6%95%B0%E5%9C%A8%E7%A7%AF%E5%88%86%E4%B8%8B%E9%99%90%2C%E5%90%AB%E5%A4%8D%E6%95%B0%E7%A7%AF%E5%88%86%EF%BC%89%E8%BF%99%E5%AE%9E%E9%99%85%E4%B8%8A%E6%98%AF%E6%88%91%E6%B1%82%E5%87%BA%E7%9A%84%E4%B8%80%E4%B8%AA%E7%89%B9%E5%BE%81%E5%87%BD%E6%95%B0%2C%E6%88%91%E5%B8%8C%E6%9C%9B%E9%80%9A%E8%BF%87%E8%BF%99%E4%B8%AA%E7%89%B9%E5%BE%81%E5%87%BD%E6%95%B0%E6%B1%82%E5%87%BA%E9%9A%8F%E6%9C%BA%E5%8F%98%E9%87%8F%E7%9A%84PDF%2C%E5%B9%B6%E7%94%B1%E6%AD%A4%E8%A7%A3%E5%86%B3%E4%B8%80%E4%B8%AA%E7%BB%99%E5%AE%9A%E6%A6%82%E7%8E%87%E7%9A%84%E6%96%B9%E7%A8%8B.%E5%A6%82%E5%9B%BE%E6%89%80%E7%A4%BA%EF%BC%9A%E7%90%86)
用matlab求解一个两重积分方程(未知数在积分下限,含复数积分)这实际上是我求出的一个特征函数,我希望通过这个特征函数求出随机变量的PDF,并由此解决一个给定概率的方程.如图所示:理
用matlab求解一个两重积分方程(未知数在积分下限,含复数积分)
这实际上是我求出的一个特征函数,我希望通过这个特征函数求出随机变量的PDF,并由此解决一个给定概率的方程.如图所示:
理想情况的是f(x)能有闭式表达式
没有闭式表达式也没有关系,能求出最后一步r的数值解即可.
感激不尽.
用matlab求解一个两重积分方程(未知数在积分下限,含复数积分)这实际上是我求出的一个特征函数,我希望通过这个特征函数求出随机变量的PDF,并由此解决一个给定概率的方程.如图所示:理
说明
1、f(x)很难求出闭式表达式.
2、使用符号积分计算量很大,而且可能根本就无法计算.
3、可使用integral函数进行数值积分.与quad系列函数相比,integral的优势是可以计算积分限为无穷大的情况.该函数自2012a引入.
4、f(x)的值应为实数,但由于数值计算的误差可能导致结果为复数,所以计算之后对结果取实部.
5、因后续还需要对f(x)进行积分,所以用arrayfun将其写成支持向量输入的形式.
6、PFA的计算,积分上限为无穷大,但取无穷大的上限容易导致出错(我不确定是否可通过算法设置避免),因而从实际计算需要出发,上限取一个有限值,例如1000.
结果
r =
10.3200
参考代码
N=25;M=4;
Fw=@(w)((1-j*w).*(1-j*w/N).*(1+w.^2/N)).^-M;
fx=@(X)real(arrayfun(@(x)1/(2*pi)*integral(@(w)Fw(w).*exp(-j*w*x),-inf,inf),X))
ezplot(fx,[-2 12])
ylabel('f(x)')
P=@(r)integral(fx,r,1000)-0.01;
r=fsolve(P,10.3)