如何用MATLAB求取积分曲线，有奇点

debuger

如何求得上述表达式的数值积分曲线，这里w*为一临界量当w=w*出现奇点，谢谢各位，这是喷管中的一个分析式

happy

那些是已知量？那些是变量？那些是自变量？把问题说清楚

i8i9i89

剔出奇点，分段求积分就可以了

debuger

x 为自变量
p,w为未知量 ，w*为临界速度已知，在由以下方程和上述方程耦合，在其奇点处即dA/dx=0,w=w* 时，dw/dx，dp/dx均未定；如何求取其积分曲线

happy

试试三楼的建议
将方程变成分段微分方程
然后求解，分段微分方程怎么求解前面有，自己搜索一下

御用左前锋

>> i0=5.83;
r=0.367;
ig1=6.09;
ig2=3.09;
ig3=1.71;
ig4=1;
y=0.85;
m=3880;
n=600:50:4000;
ua1=0.377.*r.*n./ig1./i0;
ua2=0.377.*r.*n./ig2./i0;
ua3=0.377.*r.*n./ig3./i0;
ua4=0.377.*r.*n./ig4./i0;
Ttq=-19.313+295.27.*(n/1000)-165.4.*(n/1000).*(n/1000)+40.874.*(n/1000).*(n/1000).*(n/1000)-3.8445.*(n/1000).*(n/1000).*(n/1000).*(n/1000);
Ft1=Ttq.*ig1.*i0.* y./r;
Ft2=Ttq.*ig2.*i0.* y./r;
Ft3=Ttq.*ig3.*i0.* y./r;
Ft4=Ttq.*ig4.*i0.* y./r;
f=0.013;
W=38800;
Ff=f.*W;
CDA=2.77;
Fw=CDA.*ua4.*ua4./21.15;
If=0.218;
Iw=5.396;
q1=1+Iw./r./r./m+If.*ig1.*ig1.*i0.*i0.*y./r./r; q2=1+Iw./r./r./m+If.*ig2.*ig2.*i0.*i0.*y./r./r; q3=1+Iw./r./r./m+If.*ig3.*ig3.*i0.*i0.*y./r./r; q4=1+Iw./r./r./m+If.*ig4.*ig4.*i0.*i0.*y./r./r;
a1=(Ft1-(Fw+Ff))./ig1./m;
a2=(Ft2-(Fw+Ff))./ig2./m;
a3=(Ft3-(Fw+Ff))./ig3./m;
a4=(Ft4-(Fw+Ff))./ig4./m;
plot(ua1,1./a1,ua2,1./a2,ua3,1./a3,ua4,1./a4);
xlabel('ua');ylabel('1/a');
我画出了这几条曲线，然后想求出从绿线的起点往下作一条直线，到ua=70一条直线，以及X轴，曲线，所围成的图形之间的面积，怎么办？