Audio

开发平台：

Matlab

pump.asv：源码内容

							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('%dn',Pp1(1));
%plot(Pp1(i));
end
%fprintf('%dn',Pp1(5));
%plot(Pp1);

						