24小时热门版块排行榜    

查看: 849  |  回复: 1

dax071124

新虫 (初入文坛)

[求助] 菜鸟求助matlab循环叠加问题,大神看过来,听说标题要长@@@@@@@@@@@@@@@@@@@

a(1)=0.4;a(2)=0.5;
S(1)=100;S(2)=70;
q(1)=0.05;q(2)=0.03;
v(1)=0.1;v(2)=0.2;
r=0.25;rol=-0.5;T=1;n=2;
g=zeros(n,1);b=g;c=g;d=g;
for i=1:n
    b=0;c=0;d=0;
    b=b+a(i).*S(i).*exp(-q(i).*T);
    c=c+a(i).*S(i);
    for j=1:n
      if i==j
          rol(i,j)=1;
      else
          rol(i,j)=rol;
      end
    d=d+a(i).*a(j).*S(i).*S(j).*exp((-q(i)-q(j)+rol(i,j).*v(i).*v(j)).*T);
    end
    qcall=-1/T.*log(b./c)
    vcall=sqrt(1./T.*log(d./(b.^2)))         
end

qcall =

    0.0500


vcall =

    0.8079

??? Subscripted assignment dimension mismatch.

>> -log((40*exp(-0.05)+35*exp(-0.03))/75)

ans =

    0.0406
不知为何程序运行出来的和我根据公式算出来的精确值不一样,还有出错提示是什么意思?如何修改

P_(ZE3TG]H5E`_()8LA}(NU.jpg



P_(ZE3TG]H5E`_()8LA}(NU.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dax071124

新虫 (初入文坛)

没人知道吗?都是给定的常数啊,就是一个循环求和,我找不到错。
2楼2012-10-03 21:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dax071124 的主题更新
信息提示
请填处理意见