|
Revision 168
(checked in by cbc, 16 years ago)
|
Adding diwasp customizations.
|
| Line | |
|---|
| 1 |
function S=smoothspec(Sin,kernel) |
|---|
| 2 |
% smooths a directional spectrum using the |
|---|
| 3 |
% first dimension is frequency |
|---|
| 4 |
% kernel is 2*3 matrix with smoothing parameters |
|---|
| 5 |
|
|---|
| 6 |
f1=kernel(1,3); |
|---|
| 7 |
f2=kernel(1,2); |
|---|
| 8 |
f3=kernel(1,1); |
|---|
| 9 |
d1=kernel(2,3); |
|---|
| 10 |
d2=kernel(2,2); |
|---|
| 11 |
d3=kernel(2,1); |
|---|
| 12 |
tot=2*f1+2*f2+f3+2*d1+2*d2+d3; |
|---|
| 13 |
|
|---|
| 14 |
nf=size(Sin,1); |
|---|
| 15 |
nd=size(Sin,2); |
|---|
| 16 |
|
|---|
| 17 |
S=Sin; |
|---|
| 18 |
|
|---|
| 19 |
S(3:nf-2,3:nd-2)=(f1*Sin(1:nf-4,3:nd-2)+f2*Sin(2:nf-3,3:nd-2)+f3*Sin(3:nf-2,3:nd-2)+f2*Sin(4:nf-1,3:nd-2)+f1*Sin(5:nf,3:nd-2)+... |
|---|
| 20 |
d1*Sin(3:nf-2,1:nd-4)+d2*Sin(3:nf-2,2:nd-3)+d3*Sin(3:nf-2,3:nd-2)+d2*Sin(3:nf-2,4:nd-1)+d1*Sin(3:nf-2,5:nd))./tot; |
|---|
| 21 |
|
|---|
| 22 |
|
|---|
