lz 看你再等答案，你需要的图是不是这个意思？ 
clear
r=@(x)(1-2.62*(1-0.36*(x(2)/x(1))^(-0.56))^1.5*(x(2)/x(1))^1.16);
p=@(x)(174.42*(x(1)/x(4))*((x(3)/(x(2)-x(1)))^0.85));
q=@(x)(x(1)*x(2)); p1=0;p2=0;
u=0.08:0.04/100:0.12;
F=zeros(1,length(u));
for i=1:1:101
f=0;%每次计算f清0
for j=1:1:10000
x1=normrnd(u(i),0.00125);
x2=normrnd(0.3,0.005);
x3=normrnd(0.1,0.00125);
x4=normrnd(0.1,0.005);
x5=normrnd(1.5,0.075);
x6=normrnd(16,0.2667);
x7=normrnd(0.75,0.0125);
xr=[x2,x4];
xp=[x1,x2,x3,x5];
xq=[x6,x7];
t=p(xp)*(r(xr)/q(xq))^0.5;
if(t>=1.6&&t<1.8)
f=f+0.1;
p1=p1+0.0001;
elseif(t>=1.8)
f=f+0.9;
p2=p2+0.0001;
end
end
F(i)=f;
end
plot(u,F)
希望对你有帮助 - 提问者评价
谢谢!
|
