程序如下：

function q=Ex0001()
% Ex0001 收缩喷管的流量压力关系图
k=1.3;
zeta1=0:0.01:1;
if (zeta1>0.546)&(zeta1<1)
q=sqrt((2/(k-1))*((k+1)/2).^((k+1)/(k-1))*(zeta1.^(2/k)-zeta1.^((k+1)/k)));
elseif (zeta1>0)&(zeta1<0.546)
q=1;
else zeta1==1
q=0;
end
axis([0 2 0 2]);
plot(zeta1,q)

[此贴子已经被作者于2006-5-26 21:02:17编辑过]

希望的图形是要与x,y轴相接的啊

function q=fact(zeta1)
% Ex0001 收缩喷管的流量压力关系图
k=1.3;
zeta1=0:0.01:1;
if (zeta1>0.546)&(zeta1<=1)
q=sqrt((2/(k-1))*((k+1)/2).^((k+1)/(k-1))*(zeta1.^(2/k)-zeta1.^((k+1)/k)));
elseif (zeta1>=0)&(zeta1<=0.546)
q=1;
end
axis([0 2 0 2]);
plot(zeta1,q)

这样改过以后，运行不出来曲线图，只是一个空的框框
而把这些语句直接输入命令窗口就可以运行出曲线，
但是却不是与x轴相交的。
不懂为什么会这样啊！

function q=Ex0001()
% Ex0001 收缩喷管的流量压力关系图
k=1.3;
zeta1=0:0.01:1;
for m=1:1:length(zetal);
if (zeta1(m)>0.546)&(zeta1(m)>1)
q(m)=sqrt((2/(k-1))*((k+1)/2)^((k+1)/(k-1))*(zeta1(m)^(2/k)- zeta1.^((k+1)/k)));
elseif (zeta1(m)>>0)&(zeta1(m)><0.546)
q=1;
else zeta1(m)>==1
q=0;
end
end
axis([0 2 0 2]);
plot(zeta1,q)

试试上面的程序！！！

谢谢哦

已经解决了

可以呢