www.pudn.com > SpeakerRecognitionsystem.rar > pump.m

n=50; 
b2=1; 
b4=10e-3; 
h=6.626176e-34; 
a=0.22; 
r=11; 
c=1/12; 
Pp1(1)=10e-6; 
Pp2(1)=4e-6; 
Ps(1)=2e-6; 
Pi(1)=3e-6; 
theta(1)=0; 
for i=2:n 
    %for i=1:n 
    P1(i-1)=-a.*Pp1(i-1)-4.*r.*sqrt(Pp1(i-1).*Pp2(i-1).*Ps(i-1).*Pi(i-1)).*sin(theta(i-1)); 
    P2(i-1)=-a.*Pp2(i-1)-4.*r.*sqrt(Pp1(i-1).*Pp2(i-1).*Ps(i-1).*Pi(i-1)).*sin(theta(i-1)); 
    P3(i-1)=-a.*Ps(i-1)+4.*r.*sqrt(Pp1(i-1).*Pp2(i-1).*Ps(i-1).*Pi(i-1))*sin(theta(i-1)); 
    P4(i-1)=-a.*Pi(i-1)+4.*r.*sqrt(Pp1(i-1).*Pp2(i-1).*Ps(i-1).*Pi(i-1))*sin(theta(i-1)); 
 %Pp1(i)=P1(i)+Pp1(i-1);  
 %Pp2(i)=P2(i)+Pp2(i-1); 
 %Ps(i)=P3(i)+Ps(i-1); 
 %Pi(i)=P4(i)+Pi(i-1); 
 Ws=Ps(i-1)./h; 
 Wi=Pi(i-1)./h; 
 Wc=0.5*(Ws+Wi); 
 Wp1=Pp1(i-1)./h; 
 Wp2=Pp2(i-1)./h; 
 B=b2.*((Ws-Wc)^2-(0.5*(Wp1-Wp2))^2)+b4.*c*((Ws-Wc)^4-(0.5*(Wp1-Wp2))^4); 
qtheta(i-1)=B+2.*r.*(sqrt((Pp1(i-1).*Pp2(i-1).*Pi(i-1))/Ps(i-1))+... 
sqrt((Pp1(i-1).*Pp2(i-1).*Ps(i-1))/Pi(i-1))-sqrt((Ps(i-1).*Pp2(i-1).*Pi(i-1))/Pp1(i-1))-sqrt((Ps(i-1).*Pp1(i-1).*Pi(i-1))/Pp2(i-1))).*cos(theta(i-1)) 
%+r.*(Pp1(i-1)+Pp2(i-1)-Ps(i-1)-Pi(i-1)); 
Pp1(i)=P1(i-1)+Pp1(i-1); 
Pp2(i)=P2(i-1)+Pp2(i-1); 
Ps(i)=P3(i-1)+Ps(i-1); 
Pi(i)=P4(i-1)+Pi(i-1); 
theta(i)=qtheta(i-1)+theta(i-1); 
%fprintf('%d\n',Pp1(1)); 
%plot(Pp1(i)); 
end 
%fprintf('%d\n',Pp1(5)); 
%plot(Pp1);