powerProfileFit: Fit a velocity profile to a power law.

In tunelipt/wutils: Useful functions for data processing in a wind tunnel

Description

Uses the function powerFit to find the parameters of the power law velocity profile:

u/uref = (z/zref)^p

Usage

powerProfileFit(z, u, zref = 1)

Arguments

z	Heights were the velocity was measured.
u	Velocity measured.
zref	Reference height.

Value

Vector c(p, uref).

Examples

z <- seq(10, 400, by=10)
zref <- 300
u <- 10*(z/zref)^0.2 + rnorm(10, sd=0.2)

fit <- powerProfileFit(z,u,zref=300)
plot(u, z, xlab='Velocity (m/s)', ylab='Height (mm)')
lines(fit['u']*(z/zref)^fit['p'], z)
print(fit)

tunelipt/wutils documentation built on Nov. 5, 2019, 11:01 a.m.