d*********a 发帖数: 255 | 1 【 以下文字转载自 Computation 讨论区 】
发信人: dragoninsea (龙哥), 信区: Computation
标 题: 帮看看一小段matlab程序
发信站: BBS 未名空间站 (Mon Oct 29 23:24:04 2007)
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 | d****y 发帖数: 53 | 2 there are many ways to approximate norminv. See e.g. http://home.online.no/~pjacklam/notes/invnorm/, which should be easy to code up. |
|