- 积分
- 314
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-1-25
- 最后登录
- 1970-1-1
![未绑定新浪微博用户 新浪微博达人勋](source/plugin/sina_login/img/gray.png)
|
发表于 2013-3-15 23:23:42
|
显示全部楼层
做过降水集中度分析大大侠帮看一下,这段计算降水集中度c,集中期D 的代码错在哪里,就是算不对,还有那个acrtan会得到负数的,降水期不可能出现负值的,希望懂的人能够指点一下,改怎么处理法,能给个正确的代码就更感谢万分了。本人邮箱:roos_163@163.com
for i=1:72;
SIN(i)=sin(pi/36*i);
end
SIN=SIN';
%求Rx(i)
for i=1:56;
Rx(:,i)=x(:,i).*SIN;
end
Rxi=sum(Rx);
%求72个角度的cos值
for i=1:72;
COS(i)=cos(pi/36*i);
end
COS=COS';
%求Ry(i)
for i=1:56;
Ry(:,i)=x(:,i).*COS;
end
Ryi=sum(Ry);
%把Ri转置成与Rxi具有相同维数,以便实现点除
Ri=Ri';
%求Ci
c=sqrt(Rxi.^2+Ryi.^2)./Ri;
%求Di
D=atan(Rxi./Ryi); |
|