fortran 求趋势系数

believefuture

program main
integer i,j
integer a(662,50),p(662),q(662)                      ******** 662个站点，50年的数据，所以时间的平均是25.5
real m(662)
real d(662)
open (1,file='8.31.txt')
do i=1,662
read(1,*) (a(i,j),j=1,50)
end do
do i=1,662
do j=1,50
m(i)=m(i)+a(i,j)
end do
m(i)=m(i)/50
!print *,m(i)
end do
print *,1111
do i=1,662
do j=1,50
p(i)=p(i)+(a(i,j)-m(i))*(j-25.5)
q(i)=q(i)+(a(i,j)-m(i))*(a(i,j)-m(i))*(j-25.5)*(j-25.5)
end do
end do
do i=1,662
d(i)=p(i)/sqrt(q(i)*1.0)
end do
open (2,file='1.txt')
do i=1,662
write(2,'(662f7.2)') d(i)
end do
end

river

求和竟然都不赋初值0的，那可不就是大于1了么

believefuture

river 发表于 2014-8-31 19:32
求和竟然都不赋初值0的，那可不就是大于1了么

我以为默认为0的 那也就是要将m（i），q(i),p(i),这些数组都先赋值为0是吧？

believefuture

believefuture 发表于 2014-8-31 19:35
我以为默认为0的 那也就是要将m（i），q(i),p(i),这些数组都先赋值为0是吧？

你好，我把这些数组赋值为0后，值还是没变还是>1，不知道哪错了program main
integer i,j
integer a(662,50)
real p(662),q(662),m(662),d(662)
open (1,file='8.31.txt')
do i=1,662
read(1,*) (a(i,j),j=1,50)
end do
do i=1,662
m(i)=0.0
p(i)=0.0
q(i)=0.0
d(i)=0.0
end do
do i=1,662
do j=1,50
m(i)=m(i)+a(i,j)
end do
m(i)=m(i)/50
!print *,m(i)
end do
print *,1111
do i=1,662
do j=1,50
p(i)=p(i)+(a(i,j)-m(i))*(j-25.5)
q(i)=q(i)+(a(i,j)-m(i))*(a(i,j)-m(i))*(j-25.5)*(j-25.5)
end do
end do
do i=1,662
d(i)=p(i)/sqrt(q(i)*1.0)
end do
open (2,file='2.txt')
do i=1,662
write(2,'(662f7.2)') d(i)
end do
end

river

赋初值是你这么赋的啊·····举个例子，比如你前面求和的那个赋初值
do i=1,662
m(i)=0.0
do j=1,50
m(i)=m(i)+a(i,j)
end do
m(i)=m(i)/50
!print *,m(i)
end do

believefuture

river 发表于 2014-8-31 19:55
赋初值是你这么赋的啊·····举个例子，比如你前面求和的那个赋初值
do i=1,662
m(i)=0.0

大哥不要烦我，我还是不对

program main
integer i,j
integer a(662,50)
real p(662),q(662),m(662),d(662)
open (1,file='8.31.txt')
do i=1,662
read(1,*) (a(i,j),j=1,50)
end do
do i=1,662
m(i)=0.0
do j=1,50
m(i)=m(i)+a(i,j)
end do
m(i)=m(i)/50
!print *,m(i)
end do
print *,1111
do i=1,662
p(i)=0.0
q(i)=0.0
d(i)=0.0
do j=1,50
p(i)=p(i)+(a(i,j)-m(i))*(j-25.5)
q(i)=q(i)+(a(i,j)-m(i))*(a(i,j)-m(i))*(j-25.5)*(j-25.5)
end do
end do
do i=1,662
d(i)=p(i)/sqrt(q(i)*1.0)
end do
open (2,file='3.txt')
do i=1,662
write(2,'(662f7.2)') d(i)
end do
end

言深深

p(i)=p(i)+(a(i,j)-m(i))*(j-25.5)
q(i)=q(i)+(a(i,j)-m(i))*(a(i,j)-m(i))*(j-25.5)*(j-25.5)

请问这个公式是哪里找来的呢?

believefuture

言深深 发表于 2014-8-31 20:44
p(i)=p(i)+(a(i,j)-m(i))*(j-25.5)
q(i)=q(i)+(a(i,j)-m(i))*(a(i,j)-m(i))*(j-25.5)*(j-25.5)

因为我是50年的数据，所以时间的平均值是（50+1）/2=25.5，公式是看的气象家园的帖子里的

言深深

believefuture 发表于 2014-8-31 21:10
因为我是50年的数据，所以时间的平均值是（50+1）/2=25.5，公式是看的气象家园的帖子里的

趋势，不就是最小二乘法计算斜率？貌似公式不对吧，建议找统计书籍阅读。

river

believefuture 发表于 2014-8-31 20:18
大哥不要烦我，我还是不对

program main

首先还有一个赋初值的问题，d(i)=0.0为什么要放在别的循环里啊，直接放在
do i=1,662
d(i)=p(i)/sqrt(q(i)*1.0)
end do
前面就可以了。剩下的就是你的趋势系数公式很奇怪感觉，你还是看一下现代气象统计预报里面的相关公式吧