solverp.m

matlab rp 微气泡模型。用于研究非线性科学

clear
clc

global R0 sigema rou0 Pv gamma yita P0 f1 f2 P1 P2
P0=1*10^5;
P1=1*10^5;
P2=1*10^5;
R0=5*1e-6;
f1=3*10^6;
sigema=0.073;
rou0=1000;
Pv=3.2718*1e3;
gamma=1.4;
yita=0.001;

option=odeset('maxstep',1e-8);
figure(1)
[T R]=ode45('rp',[0 38/f1] ,[R0 0]);
v=4/3*pi*R(:,1).^3;
plot(T,v)
title('The solution of classic R-P ==>V')
xlabel ('time sequence')
ylabel ('volume')