matlab方程求解，HELP~~~

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在做Copula的时候，需要解一个二元正态的方程，如下图：
其中，相关系数rho为已知，rho=0.8  ,需要求y的解
由于积分部分没有原函数，求高手指点如何求解，甚为感激，盼君赐教！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB atlab matla Help Atl 如何 matlab 原函数

找找书吧，肯定有方法的，积分函数传递参数的做法应该可以实现。

感谢楼上，我觉得重点不是参数传递吧，而是算法的问题，我用二分法计算的，效率很低！

期待epoh老师解答~~

yishions 发表于 2011-12-3 13:47
期待epoh老师解答~~

刚试了:
  Warning: Explicit integral could not be found

表示符号方法没结果
请尝试数值方法.
dblquad.m
  Numerically evaluate double integral over rectangle
q = dblquad(fun,xmin,xmax,ymin,ymax,tol,method)

因为被积函数的原函数是不可以用初等函数表示的，所以符号方法应该是不好算。我用的就是数值方法，就是运行速度出奇的慢~~~

修改后程序约4分钟跑完.

修改部分:

  数据:原rho(251)=-0.0307,改为rho(251)=0.0307

         因为p=-0.0307,f1 = -0.0025,f2 = -9.7692e-005

         f1*f2>0,会造成程序中止

  程序:改用dblquad,及-Inf改为-50,a & b维持原来.

         2*1-p^2  --> 2*(1-p^2)

%%%%%in command window
datarho=xlsread('allrho.xls');
rho=datarho(:,1);
T = length(rho);
for i=1:T
i
p=rho(i);
yy(i)=HalfInterval(-3.3,-1.6,p,0.00001);  
end
%%%%
HalfInterval.m, FindRoots.m (need pw)
   HalfInterval.rar (777 Bytes)

2011-12-4 08:44:24 上传

epoh哥真是个天才，像用-50代替-Inf打死我也想不出来，跑下试试

程序完美运行，精度也很高，你让我怎么感谢你

yishions 发表于 2011-12-4 12:35
程序完美运行，精度也很高，你让我怎么感谢你

哈哈!问题解决就好
主要还是你的两个function写得很好
我只是小修而已