怎么用matlab来估计广义极值分布的参数吗

Matlab在统计（极值分布）中的问题，在线等，急~

代码如下

X=random('ev',1,1.5,[1000 1]);
>> histfit(X,20,'ev')

matlab里有利用mle来估计gp的参数的函数：gpfit()
y = x(x>u)-u;paramEsts = gpfit(y); %得到的第一个参数是形状参数，第二个是尺度参数可参考下面链接
http://cn.mathworks.com/help/stats/examples/modelling-tail-data-with-the-generalized-pareto-distribution.html

第一种方法：直接调用random函数，如下：datan1=random('norm',3.5,sqrt(5.8))datag1=random('gam',0.3,4.0)datap1=random('poiss',4.2)random函数的第一个字符串参数表明了分布类型，后面的参数依次是该分布按习惯顺序的对应参数，例如正态分布就是μ（均值）和σ（标准差），伽马分布就是α（形状参数）和β（尺度参数），泊松分布就是单参数λ（均值或方差）。第二种方法：用各自分布对应的函数，如下：datan2=normrnd(3.5,sqrt(5.8))datag2=gamrnd(0.3,4.0)datap2=poissrnd(4.2)参数依次是该分布按习惯顺序的对应参数，用法同上，可以认为是random函数结合不同分布的字符串退化为各自分布的独立函数。

---

芝罘区19588068627： 怎么用matlab来估计广义极值分布的参数吗 - ？
邵波痱子： 第一种方法:直接调用random函数,如下:datan1=random('norm',3.5,sqrt(5.8))datag1=random('gam',0.3,4.0)datap1=random('poiss',4.2)random函数的第一个字符串参数表明了分布类型,后面的参数依次是该分布按习惯顺序的对应参数,例如...

芝罘区19588068627： 怎么用matlab算广义积分? - ？
邵波痱子： 用matlab算广义积分方法:%符号积分 int(f,v) int(f,v,a,b)%数值积分 trapz(x,y)%梯形法沿列方向求函数Y关于自变量X的积分 cumtrapz(x,y)%梯形法沿列方向求函数Y关于自变量X的累计积分 quad(fun,a,b,tol)%采用递推自适应Simpson法计算积分 quad1(fun,a,b,tol)%采用递推自适应Lobatto法求数值积分 dbquad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%二重(闭型)数值积分指令 triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%三重(闭型)数值积分指令 .

芝罘区19588068627： 如何用matlab求极值???加50分 跪求! - ？
邵波痱子： 在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_heyoungguang(xx) %%x,k,l分别用xx(1),xx(2),xx(3)代替x...

芝罘区19588068627： 如何用matlab求函数的极值点 - ？
邵波痱子： Rosenbrock函数简介 Rosenbrock函数是数学优化中的常用函数,他可以用来测试优化算法的性能.又叫“Rosenbrock山谷”、“Rosenbrock香蕉函数”、“香蕉函数”.它的定义式如下图所示:匿名函数表示该函数,matlab的匿名函数代码为:ff=@(x)(100*(x(2)-x(1).^2)^2+(1-x(1))^2.单纯形法求极值 检查目标函数 画rosenbrock函数三维图形,如下图所示是画图代码:下面是画出的rosenbrock函数三维图形

芝罘区19588068627： 如何用Matlab求极值 - ？
邵波痱子： 极值: data数据 find(diff(sign(diff(data)))==-2)+1 找极值位置 find(diff(sign(diff(data)))==2)+1 找极值位置 data(find(diff(sign(diff(data)))==-2)+1) data(find(diff(sign(diff(data)))==2)+1) 返极值极值

芝罘区19588068627： 用MATLAB求极值与积分. - ？
邵波痱子： 1) fun=inline('x(1)^2+5*x(2)^2-6*x(1)+10*x(2)+6'); [x,g]=fminsearch(fun,[0,0]); x =3 -1 g =-82) syms x s=int((1+log(x))/x,1,exp(1)); vpa(s) ans = 1.5

芝罘区19588068627： 如何用matlab求多元函数的极值 - ？
邵波痱子： on来求解.写成向量形式,约束A*x<=b,初始值x0附近可以搜出出局部极小值.输入: z=@(x)7.2*sqrt(25+(15-x(1))^2)+(7.2+(107/(20-x(1))))*sqrt((8-x(2))^2+(20-x(1))^2)+7.2*x(2); A=[1,0;0,1]; b=[15;8]; x0=[1;1]; [x,feval]=fmincon(z,x0,A,b) 得到局部极小值的坐标和函数值: x =15.00006.6994 feval =231.9944

芝罘区19588068627： matlab中 函数极值怎么求 - ？
邵波痱子： 极值: data是你的数据, find(diff(sign(diff(data)))==-2)+1 找到极大值的位置 find(diff(sign(diff(data)))==2)+1 找到极小值的位置 data(find(diff(sign(diff(data)))==-2)+1)和 data(find(diff(sign(diff(data)))==2)+1) 返回的是极大值和极小值

芝罘区19588068627： 用matlab求所有极值 - ？
邵波痱子： 随便编了一个,肯定还有更好的方法. x=-6:0.01:6; y=sinc(x); plot(x,y); hold on m=0;n=0; for i=2:length(x)-1if y(i-1)<y(i) && y(i)>y(i+1)m=m+1;x1(m)=x(i);y1(m)=y(i);endif y(i-1)>y(i) && y(i)<y(i+1)n=n+1;x2(n)=x(i);y2(n)=y(i);end end plot(x1,y1,'ro',x2,y2,'bo') hold off

芝罘区19588068627： matlab广义积分 - ？
邵波痱子： 我觉得,从数学上讲,此积分是无穷的,函数只在正负1时有最小值,在别处都在x轴上方,且趋近于1,所以x取的越多,积分越大,取无穷,积分肯定无穷.这是画的图,你看一下 x=-10:10; plot(x,exp((-x.^2)./(1+x.^4))) 希望别人给出一些好的办法 还有我的运行这样digits(32) syms x y=vpa(int(exp((-x.^2)./(1+x.^4)),-inf,+inf)) Warning: Explicit integral could not be found.> In sym.int at 58 y =(Inf)