matlab简单的程序求改错function Cha2demo5rho=961;G=6.67;L=10;eps=5e-6;deltaP=15e3;k=1.8;n=0.64;K=1.48;n1=n;K1=K*((3*n+1)/(4*n))^n;D1=(32*K1*L*8^(n1-1)*(4*G/pi/rho)^n1)^(1/(3*n1+1));D2=2*D1;delta=1e-4;while abs((D2-D1)/D2)>deltau=(G/rho)/(pi*D1

matlab简单的程序求改错function Cha2demo5rho=961;G=6.67;L=10;eps=5e-6;deltaP=15e3;k=1.8;n=0.64;K=1.48;n1=n;K1=K*((3*n+1)/(4*n))^n;D1=(32*K1*L*8^(n1-1)*(4*G/pi/rho)^n1)^(1/(3*n1+1));D2=2*D1;delta=1e-4;while abs((D2-D1)/D2)>deltau=(G/rho)/(pi*D1
matlab简单的程序求改错
function Cha2demo5
rho=961;G=6.67;L=10;eps=5e-6;deltaP=15e3;k=1.8;n=0.64;K=1.48;n1=n;
K1=K*((3*n+1)/(4*n))^n;
D1=(32*K1*L*8^(n1-1)*(4*G/pi/rho)^n1)^(1/(3*n1+1));
D2=2*D1;
delta=1e-4;
while abs((D2-D1)/D2)>delta
u=(G/rho)/(pi*D1^2/4);
Regen=D1^n1*u(2-n1)*rho/(8^(n1-1)*K1);
if Regen2100
f=FrictFacto(Regen,n1);
end
Le=k*D1/(4*f);
D1=D2;
D2=(2*f*(L+Le)/(rho*deltaP)*(4*G/pi)^2)^0.2;
end
fprint('\t管道直径为D=%.4f%s\n',D2,'m')
fprint('\t摩擦因子为f=%.4f%s',f,'m')
%------------------------------------------------
function f=FtictFactor(Regen,n1)
f0=16/Regen; %4.53
f=fzero(@fFunc,f0,[],Regen,n1)
%------------------------------------------------
function F=fFunc(f)
F=1/sqrt(f)-4.0*log(Regen*f^(1-n1/2))/n1^0.75+0.4/n1^1.2;%4.54

matlab简单的程序求改错function Cha2demo5rho=961;G=6.67;L=10;eps=5e-6;deltaP=15e3;k=1.8;n=0.64;K=1.48;n1=n;K1=K*((3*n+1)/(4*n))^n;D1=(32*K1*L*8^(n1-1)*(4*G/pi/rho)^n1)^(1/(3*n1+1));D2=2*D1;delta=1e-4;while abs((D2-D1)/D2)>deltau=(G/rho)/(pi*D1
clf;
t=6*pi*(0:100)/100;
y=1-exp(0.3*t).*cos(0.7*t);
tt=t(abs(y-1)>0.05);
ts=max(tt);
plot(t,y,'r-','LineWidth',3)
axis([-inf,6*pi,0.6,inf])%此处修改,左边的inf应改为-inf或者按你的要求改小于6*pi.
set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)])
grid on
title('/it y=1-e^{-\alphat}cos{\omegat}')
text(13.5,1.2,'\fontsize{12}{\alpha}=0.3')
text(13.5,1.1,'\fontsize{12}{\omega}=0.7')
hold on;
plot(ts,0.95,'bo','Markersize',10);
hold off;
cell_string{1}='\fontsize{12}\uparrow';
cell_string{2}='\fontsize{16}\fontname{隶书}镇定时间';
cell_string{3}='\fontsize{6}';
cell_string{4}=['\fontsize{14}\rmt_{s}= ',num2str(ts)];%此处修改,num2str(ts)前应有分隔的,.
text(ts,0.85,cell_string)
xlabel('\fontsize{14}\bft\rightarrow')
ylabel('\fontsize{14}\bfy\rightarrow')
求采纳为满意回答.