lirong-huahua 发表于 2013-1-1 12:59:21

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

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 发表于 2017-2-21 20:37:58

不清楚你的问题是什么,直接在开头写
if y(1)<0
y(1)=0
end
不就行啦
页: [1]
查看完整版本: matlab怎么加边界,不是初始条件,是边界条件