d*********a
发帖数: 255 | 1 M=36;
theta=zeros(1,M);
pho=[1 0.5;0.5 1];
syms u;
syms v;
s=[norminv(u);norminv(v)];
for i=1:M
theta(i)=int(int(10*(exp(-0.5*s'*(inv(pho)-eye(2))*s)/sqrt(det(pho))),u,0,1)
,v,0,1);
end
theta
===============================================
??? Error using ==> <
Function '<' not defined for variables of class 'sym'.
Error in ==> C:\MATLAB6p1\toolbox\stats\norminv.m
On line 36 ==> k = find(sigma <= 0 | p < 0 | p > 1 | isnan(p));
如何修改才可以正确运行？或者有没有什么简单函数可以近似norminv函数？谢谢 | r****y
发帖数: 1437 | 2 你这个程序差的也太远了吧，错到天上了去自己还看不出来？
just using icdf should be enough.
1)
【在 d*********a 的大作中提到】
: M=36;
: theta=zeros(1,M);
: pho=[1 0.5;0.5 1];
: syms u;
: syms v;
: s=[norminv(u);norminv(v)];
: for i=1:M
: theta(i)=int(int(10*(exp(-0.5*s'*(inv(pho)-eye(2))*s)/sqrt(det(pho))),u,0,1)
: ,v,0,1);
: end
|
|
