-
Notifications
You must be signed in to change notification settings - Fork 0
/
etaAF.m
12 lines (11 loc) · 923 Bytes
/
etaAF.m
1
2
3
4
5
6
7
8
9
10
11
12
function [ etaA ] = etaAF( x1,zeta1 )
% 顺风向风振加速度的脉动系数 表 J.1.2
etaA_tab = [4.14 3.93 3.75 3.59 3.46 3.35 2.67 3.24 2.12 1.97 1.86 1.76 1.69 1.62 1.56 1.47 1.40 1.34 1.29 1.24 1.20 1.17 1.14 1.11 1.09;...
2.94 2.79 2.66 2.55 2.46 2.38 1.90 1.66 1.51 1.40 1.32 1.25 1.20 1.15 1.11 1.05 0.99 0.95 0.91 0.88 0.85 0.83 0.81 0.79 0.77;...
2.41 2.28 2.18 2.09 2.02 1.95 1.55 1.36 1.23 1.15 1.08 1.03 0.98 0.94 0.91 0.86 0.81 0.78 0.75 0.72 0.70 0.68 0.66 0.65 0.63;...
2.10 1.99 1.90 1.82 1.75 1.69 1.35 1.18 1.07 1.00 0.94 0.89 0.85 0.82 0.79 0.74 0.71 0.68 0.65 0.63 0.61 0.59 0.58 0.56 0.55;...
1.88 1.78 1.70 1.63 1.57 1.52 1.21 1.06 0.96 0.89 0.84 0.80 0.76 0.74 0.71 0.67 0.63 0.61 0.58 0.56 0.55 0.53 0.52 0.50 0.49];
x1Ind = [5:10 20:10:100 120:20:300]';
zetaInd = [0.01 0.02 0.03 0.04 0.05];
etaA = interp2(x1Ind,zetaInd,etaA_tab,x1,zeta1,'linear');
end