已经解决了,不过代码有点长,想过用向量来处理,但由于图不是均匀的,我未能成功,大家可以一起想想能否简化
clf
a=1;
b=6;
for x=a:0.1:b;
y=[-3.719,3.719];
plot([x,x],y)
hold on
end
for i=0.0001:0.0001:0.001;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.001:0.001:0.01;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.01:0.01:0.1;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.1:0.01:0.2;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.2:0.02:0.8;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.8:0.01:0.9;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.9:0.01:0.99;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.99:0.001:0.999;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
for i=0.999:0.0001:0.9999;
y1=norminv(i);
x1=[a,b];
plot(x1,[y1,y1])
hold on
end
title('正态概率纸') |