暴雨强度公式推导
用R的nls程序即可:
A=26.457777
B=7.60363
N=0.62901
standard error=0.009808
以下为程序代码及输出:
> baoyudata=data.frame(cbind(
+ T=c(5, 10, 15, 20, 30, 45, 60, 90, 120),
+ M=c(5.38, 4.31, 3.74, 3.34, 2.68, 2.18, 1.86, 1.48, 1.26)))
> baoyudata
T M
1 5 5.38
2 10 4.31
3 15 3.74
4 20 3.34
5 30 2.68
6 45 2.18
7 60 1.86
8 90 1.48
9 120 1.26
> logBY <-function (T, A, B, N) {
+ + log(A) - N*log(T+B)}
> bynls <- nls(log(M) ~ logBY(T, A, B, N), data=baoyudata, start=list(B=12, A=2, N=2))
> summary(bynls)
Formula: log(M) ~ logBY(T, A, B, N)
Parameters:
Estimate Std. Error t value Pr(>|t|)
B 7.60363 0.74279 10.24 5.07e-05 ***
A 26.45777 1.77004 14.95 5.65e-06 ***
N 0.62901 0.01418 44.37 8.77e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.009808 on 6 degrees of freedom
Number of iterations to convergence: 6
Achieved convergence tolerance: 1.234e-08
>
暴雨强度公式有哪些
q=5544*(P0.3-0.42)除以（t+10-7*lgP)0.82+0.07lgP 单位是L/(s.hm2) 其中0.3和0.82+0.07lgP是指数