Matlab 二重积分 求解

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请大家帮我看看这个公式（求二重积分）对不对，我试了好多次都没整明白，Ris和Rif都是289x1的收益率序列，按照下面的编写能得到F（289x1）的矩阵吗？请大侠指点，谢谢！       for j=1:289
            F(j) = quad2d(@(Ris(j),Rif(j)) copulapdf('Clayton',[u(:) v(:)],theta).*tpdf(Ris(j),omegais).*tpdf(Rif(j),omegaif),-3,3,-4,@(Ris(j)) (Ris(j)-Rp(j))/h);
      end


其中：u=tpdf(Ris(j),omegais)； v=tpdf(Rif(j),omegaif);积分区间是[-3,3][-4,(Ris-Rp)/h];

Matlab总是报错说：         ??? Error: File: rpquad2d.m Line: 22 Column: 22
                                         Unbalanced or unexpected parenthesis or bracket.

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB matla atlab 二重积分 Lab 收益率 Error

符号计算
syms x y;                    %定义两个符号变量
a=int(int(x^y,x,0,1),y,1,2)  %积分x,0,1   ,y,1,2
b=simple(a)                  %化简
c=vpa(b,4)                   %得到4位近似解，也可以任意N位解
数值计算
%%二重积分f= @(x,y)exp(sin(x))*ln(y)，y从5*x积分到x^2，x从10积分到20
1 (7.X后版本才有此函数quad2d)
   y1=quad2d(@(x,y) exp(sin(x)).*log(y),10,20,@(x)5*x,@(x)x.^2)

2  y2 =quadl(@(x) arrayfun(@(x) quadl(@(y)exp(sin(x)).*log(y),5*x,x.^2),x),10,20)

3  y3 = dblquad(@(x,y)exp(sin(x)).*log(y).*(y>=5*x & y<=x.^2),10,20,50,400)