matlab二重循环遇到问题,%离散等轴断面孔型高度H2=zeros(200,1);B2=zeros(200,1);H4=zeros(200,1);B4=zeros(200,1);H2=140:0.2:180;B2=140:0.2:180;H4=120:0.2:160;B4=120:0.2:160;%求中间非等轴断面的高和宽%利用乌萨托夫斯基

matlab二重循环遇到问题,%离散等轴断面孔型高度H2=zeros(200,1);B2=zeros(200,1);H4=zeros(200,1);B4=zeros(200,1);H2=140:0.2:180;B2=140:0.2:180;H4=120:0.2:160;B4=120:0.2:160;%求中间非等轴断面的高和宽%利用乌萨托夫斯基
matlab二重循环遇到问题,
%离散等轴断面孔型高度
H2=zeros(200,1);
B2=zeros(200,1);
H4=zeros(200,1);
B4=zeros(200,1);
H2=140:0.2:180;
B2=140:0.2:180;
H4=120:0.2:160;
B4=120:0.2:160;
%求中间非等轴断面的高和宽
%利用乌萨托夫斯基宽展公式
%第一道平箱等效宽和高
m0=0.96;H0=200;B0=200;D1=735;D2=720;
H1=145*ones(200,1);
B1=zeros(200,1);
for i=1:200
for j=0:0.01:30
Hc0=200*m0;Hc1(i)=H1(i)*m0;dta1=B0/Hc0;
ke1(i)=Hc0/(D1-Hc1(i));
yita1(i)=Hc1(i)/Hc0;
if yita1(i)>0.5
w1(i)=-10^(-1.269*dta1*ke1(i)^0.556);
else
w1(i)=-10^(-3.457*dta1*ke1(i)^0.968);
end
beta1(i)=yita1(i)^w1(i);
B1(i)=beta1(i)*H0;
%验算轧件在第二孔型中的充满情况 %宽变高 高变宽?
dta2(i)=H1(i)/(B1(i)*m0); %此处需要注意 dta轧件断面形状系数
Hc2(i)=H2(i)*m0;
Bc1(i)=B1(i)*m0;
ke2(i)=Bc1(i)/(D2-Hc2(i));
yita2(i)=Hc2(i)/Bc1(i);
if yita2(i)>0.5
w2(i)=-10^(-1.269*dta2(i)*ke2(i)^0.556);
else
w2(i)=-10^(-3.457*dta2(i)*ke2(i)^0.968);
end
beta2(i)=yita2(i)^w2(i);
b2(i)=beta2(i)*H1(i);
if abs(b2(i)-H2(i))>0.5
H1(i)=H1(i)+j;
else
H1(i)=H1(i);
end
end
end
H1;
B1;

matlab二重循环遇到问题,%离散等轴断面孔型高度H2=zeros(200,1);B2=zeros(200,1);H4=zeros(200,1);B4=zeros(200,1);H2=140:0.2:180;B2=140:0.2:180;H4=120:0.2:160;B4=120:0.2:160;%求中间非等轴断面的高和宽%利用乌萨托夫斯基
%离散等轴断面孔型高度
H2=linspace(140,180,200);
B2=linspace(140,180,200);
H4=linspace(120,160,200);
B4=linspace(120,160,200);
%求中间非等轴断面的高和宽
%利用乌萨托夫斯基宽展公式
%第一道平箱等效宽和高
m0=0.96;H0=200;B0=200;D1=735;D2=720;
H1=145*ones(200,1);
Hc1=ones(200,1);
B1=zeros(200,1);
for i=1:200
for j=linspace(0,30,200);
Hc0=200*m0;Hc1(i)=H1(i)*m0;dta1=B0/Hc0;
ke1(i)=Hc0/(D1-Hc1(i));
yita1(i)=Hc1(i)/Hc0;
if yita1(i)>0.5
w1(i)=-10^(-1.269*dta1*ke1(i)^0.556);
else
w1(i)=-10^(-3.457*dta1*ke1(i)^0.968);
end
beta1(i)=yita1(i)^w1(i);
B1(i)=beta1(i)*H0;
%验算轧件在第二孔型中的充满情况 %宽变高 高变宽?
dta2(i)=H1(i)/(B1(i)*m0); %此处需要注意 dta轧件断面形状系数
Hc2(i)=H2(i)*m0;
Bc1(i)=B1(i)*m0;
ke2(i)=Bc1(i)/(D2-Hc2(i));
yita2(i)=Hc2(i)/Bc1(i);
if yita2(i)>0.5
w2(i)=-10^(-1.269*dta2(i)*ke2(i)^0.556);
else
w2(i)=-10^(-3.457*dta2(i)*ke2(i)^0.968);
end
beta2(i)=yita2(i)^w2(i);
b2(i)=beta2(i)*H1(i);
if abs(b2(i)-H2(i))>0.5
H1(i)=H1(i)+j;
else
H1(i)=H1(i);
end
end
end
plot(B1,H1)
按上面语句试试,经过运算,有复数产生.