SCM Repository
Annotation of /SplineFitBATMAN_MATLABcode/intersectSpecSpline.m
Parent Directory
|
Revision Log
Revision 346 - (view) (download)
1 : | jie | 346 | function ppmx = intersectSpecSpline(Xn, offset, x, y, ppm, minflag, ppmerror) |
2 : | % written 120213 Dr. Jie Hao, Imperial College London | ||
3 : | |||
4 : | [n,m]=size(Xn); | ||
5 : | offy=[0:n-1].*offset; | ||
6 : | off = repmat(offy',1,m); | ||
7 : | D = Xn + off; | ||
8 : | [y iy]=sort(y); | ||
9 : | x=x(iy); | ||
10 : | |||
11 : | % ppmx = zeros(2,size(D,1)); | ||
12 : | ppmx = intersectSpecSplineRatio(Xn, offset, x, y, ppm); | ||
13 : | for i = 1:size(ppmx,2) | ||
14 : | ind = find(ppm>=ppmx(1,i)-ppmerror & ppm<= ppmx(1,i)+ppmerror); | ||
15 : | if (~isempty(ind)) | ||
16 : | if minflag | ||
17 : | [ny, nyi] = min(D(i,ind)); | ||
18 : | else | ||
19 : | [ny, nyi] = max(D(i,ind)); | ||
20 : | end | ||
21 : | |||
22 : | ppmx(1,i) = ppm(ind(nyi)); | ||
23 : | ppmx(2,i) = ny; | ||
24 : | end | ||
25 : | end |
R-Forge@R-project.org | ViewVC Help |
Powered by ViewVC 1.0.0 |