root/applications/dual/SIDP/bakalarka/SIDPS/matlab/motor2/interpol.m @ 1351

Revision 1351, 0.5 kB (checked in by zimamiro, 13 years ago)
Line 
1function rizeni=interpol(beta,H,C,H_tilde_on)
2i=(1:4)*H_tilde_on*ones(5,1);
3j=ones(1,4)*H_tilde_on*[1;2;3;4;5];
4pom=H-beta(i,j);
5
6if (pom(1)>0)
7    rizeni=C(1);
8else
9    if (pom(end)<0 || pom(end)==0)
10        rizeni=C(end)*2^pom(end);
11    else
12        for i=2:size(H,2)
13            if (pom(i)>0)
14                delka=H(i)-H(i-1);
15                rizeni=(1+pom(i-1)/delka)*C(i-1)+(1-pom(i)/delka)*C(i);
16                break               
17            end           
18        end
19    end
20end
21
22end
Note: See TracBrowser for help on using the browser.