求强人指导---为什么print出来的相关系数绝对值会大于1呢？？？？？

dzxconan

!---------子程序：计算相关系数-----
subroutine xiangguan(a,b,r,sa,sb,avea,aveb)
        parameter (nx=180,ny=89,m=54,undef=32767)
real a(nx,ny,m),b(m),sa(nx,ny),sb,sab(nx,ny),r(nx,ny),avea(nx,ny),aveb,t(nx,ny)

do k=1,m
do j=1,ny
do i=1,nx
if(a(i,j,k)/=undef)then
avea(i,j)=avea(i,j)+a(i,j,k)/m
else
avea(i,j)=undef
endif
enddo
enddo
aveb=aveb+b(k)/m
enddo

do k=1,m
do j=1,ny
do i=1,nx
if(a(i,j,k)/=undef)then
sa(i,j)=sa(i,j)+(a(i,j,k)-avea(i,j))**2/m
sab(i,j)=sab(i,j)+(a(i,j,k)-avea(i,j))*(b(k)-aveb)/m
else
sa(i,j)=undef
sab(i,j)=undef
endif
enddo
enddo
sb=sb+(b(k)-aveb)**2/m
enddo

do j=1,ny
do i=1,nx
if(sa(i,j)/=undef)then
r(i,j)=sab(i,j)/(sqrt(sa(i,j))*sqrt(sb))
else
r(i,j)=undef
endif
if((abs(r(i,j))>=1).and.(r(i,j)/=undef))then
write(*,*)r(i,j),i,j,sa(i,j),sb,sab(i,j)
endif
enddo
enddo


end

dzxconan

真悲剧 发帖还扣金币。。。。。。。。。。。。。。

做个霸气的木头

相关系数大于1肯定是你的算法哪里出了问题，这种楼主自己细心检查程序与算法的对应关系。

dzxconan

做个霸气的木头 发表于 2013-5-3 13:11

                               
登录/注册后可看大图

相关系数大于1肯定是你的算法哪里出了问题，这种楼主自己细心检查程序与算法的对应关系。

能帮忙看看吗？？对应了好多遍 我都是用协方差除以两个的标准差的乘积。。。。O(∩_∩)O谢谢啦

dzxconan

做个霸气的木头 发表于 2013-5-3 13:11

                               
登录/注册后可看大图

相关系数大于1肯定是你的算法哪里出了问题，这种楼主自己细心检查程序与算法的对应关系。

而且我第一次调用没有问题 第二次就有问题啦。。。。。。

topmad

dzxconan 发表于 2013-5-3 13:14

                               
登录/注册后可看大图

而且我第一次调用没有问题 第二次就有问题啦。。。。。。

很可能是某个值赋值的时候出了问题，你检查一下赋值的顺序，看看逻辑上有没有错误

dzxconan

topmad 发表于 2013-5-3 13:54

                               
登录/注册后可看大图

很可能是某个值赋值的时候出了问题，你检查一下赋值的顺序，看看逻辑上有没有错误

O(∩_∩)O谢谢啦~~我想应该是第二次忘了给初值 它的第一次的值还在吧？？？

river

应该没有赋初值的问题吧，楼主这个程序一个初值都没赋啊

topmad

river 发表于 2013-5-3 15:41

                               
登录/注册后可看大图

应该没有赋初值的问题吧，楼主这个程序一个初值都没赋啊

程序我还没有看   你再看看是不是算法上面的问题 过程中间涉及到的一些单位还是什么东西  细心一点

dzxconan

topmad 发表于 2013-5-3 16:14

                               
登录/注册后可看大图

程序我还没有看   你再看看是不是算法上面的问题 过程中间涉及到的一些单位还是什么东西  细心一点

嗯嗯 初值问题下回要细心 O(∩_∩)O谢谢~~