刚用上matlab 看看这个错误是什么意思 还希望指正程序如下 想画一个f与u(正态函数x1的均值的函数关系)r=@(x2,x4)(1-2.62*(1-0.36*(x4/x2)^(-0.56))^1.5*(x4/x2)^1.16)p=@(x1,x2,x3,x5)(174.42*(x1/x5)*((x3/(x2-x1))^0.85)
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/06 18:01:57
![刚用上matlab 看看这个错误是什么意思 还希望指正程序如下 想画一个f与u(正态函数x1的均值的函数关系)r=@(x2,x4)(1-2.62*(1-0.36*(x4/x2)^(-0.56))^1.5*(x4/x2)^1.16)p=@(x1,x2,x3,x5)(174.42*(x1/x5)*((x3/(x2-x1))^0.85)](/uploads/image/z/10480544-8-4.jpg?t=%E5%88%9A%E7%94%A8%E4%B8%8Amatlab+%E7%9C%8B%E7%9C%8B%E8%BF%99%E4%B8%AA%E9%94%99%E8%AF%AF%E6%98%AF%E4%BB%80%E4%B9%88%E6%84%8F%E6%80%9D+%E8%BF%98%E5%B8%8C%E6%9C%9B%E6%8C%87%E6%AD%A3%E7%A8%8B%E5%BA%8F%E5%A6%82%E4%B8%8B+%E6%83%B3%E7%94%BB%E4%B8%80%E4%B8%AAf%E4%B8%8Eu%EF%BC%88%E6%AD%A3%E6%80%81%E5%87%BD%E6%95%B0x1%E7%9A%84%E5%9D%87%E5%80%BC%E7%9A%84%E5%87%BD%E6%95%B0%E5%85%B3%E7%B3%BB%EF%BC%89r%3D%40%28x2%2Cx4%29%281-2.62%2A%281-0.36%2A%28x4%2Fx2%29%5E%28-0.56%29%29%5E1.5%2A%28x4%2Fx2%29%5E1.16%29p%3D%40%28x1%2Cx2%2Cx3%2Cx5%29%28174.42%2A%28x1%2Fx5%29%2A%28%28x3%2F%28x2-x1%29%29%5E0.85%29)
刚用上matlab 看看这个错误是什么意思 还希望指正程序如下 想画一个f与u(正态函数x1的均值的函数关系)r=@(x2,x4)(1-2.62*(1-0.36*(x4/x2)^(-0.56))^1.5*(x4/x2)^1.16)p=@(x1,x2,x3,x5)(174.42*(x1/x5)*((x3/(x2-x1))^0.85)
刚用上matlab 看看这个错误是什么意思 还希望指正
程序如下 想画一个f与u(正态函数x1的均值的函数关系)
r=@(x2,x4)(1-2.62*(1-0.36*(x4/x2)^(-0.56))^1.5*(x4/x2)^1.16)
p=@(x1,x2,x3,x5)(174.42*(x1/x5)*((x3/(x2-x1))^0.85))
q=@(x6,x7)(x6*x7)
f=0;
p1=0;p2=0;
u=0.08:0.04/100:0.12;
for i=1:1:10000
x1=normrnd(u,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);
t=p(x1,x2,x3,x5)*(r(x2,x4)/q(x6,x7))^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;
>> plot(u,f)
刚用上matlab 看看这个错误是什么意思 还希望指正程序如下 想画一个f与u(正态函数x1的均值的函数关系)r=@(x2,x4)(1-2.62*(1-0.36*(x4/x2)^(-0.56))^1.5*(x4/x2)^1.16)p=@(x1,x2,x3,x5)(174.42*(x1/x5)*((x3/(x2-x1))^0.85)
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)