matlab怎么加边界，不是初始条件，是边界条件

lirong-huahua

function dy=xq(theta,y)
V=0.0000138474*(0.17+(1-cos(theta))/2+(1-sqrt(1-0.04*(sin(theta))^2))/0.4);
V1=0.0000138474*(sin(theta)/2+0.1*sin(theta)*cos(theta)/sqrt(1-0.04*(sin(theta))^2));
k=1.4;
w=102.573;
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=(3.2041*(1-y(3))-1.9308*y(1))/2.189255;
dy(3)=(-k/(w*V))*((w*y(3)*V1)-(sqrt(2*k*82369/(k-1)))*0.01099*y(1)*(y(3)^(1/k))*(sqrt(1-y(3)^(1-1/k))));
end

怎样加上边界条件：
若y(1)<0
则y(1)=0

lijingyuak

不清楚你的问题是什么，直接在开头写
if y(1)<0
y(1)=0
end
不就行啦