您需要添加另外两个循环来打印另外两组二维网格:
for rr=1:1:length(x1a_dis)
x1a_new = repmat(x1a_dis(rr),size(z2));
hold on
plot3(x1a_new,y1,z2,'r');
end
for rr=1:1:length(x1a_dis)
x1a_new = repmat(x1a_dis(rr),size(z2));
hold on
plot3(z2,x1a_new,y1,'r');
end
for rr=1:1:length(x1a_dis)
x1a_new = repmat(x1a_dis(rr),size(z2));
hold on
plot3(y1,z2,x1a_new,'r');
end
注意:您可以通过以下几种方式简化代码:
一
台词:
y1(ii,jj) = exp(log(1-x(ii)^(S(jj)*T))/(S(jj)*T));
可以简化为:
y1(ii,jj) = (1-x(ii)^(S(jj)*T)).^(1/(S(jj)*T));
使用日志和经验规则。
二
循环(更改编号1后):
y1 = zeros(length(x),length(S));
for ii = 1:length(x)
for jj = 1:len