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
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/22 04:42:29
![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](/uploads/image/z/12969808-16-8.jpg?t=matlab%E7%AE%80%E5%8D%95%E7%9A%84%E7%A8%8B%E5%BA%8F%E6%B1%82%E6%94%B9%E9%94%99function+Cha2demo5rho%3D961%3BG%3D6.67%3BL%3D10%3Beps%3D5e-6%3BdeltaP%3D15e3%3Bk%3D1.8%3Bn%3D0.64%3BK%3D1.48%3Bn1%3Dn%3BK1%3DK%2A%28%283%2An%2B1%29%2F%284%2An%29%29%5En%3BD1%3D%2832%2AK1%2AL%2A8%5E%28n1-1%29%2A%284%2AG%2Fpi%2Frho%29%5En1%29%5E%281%2F%283%2An1%2B1%29%29%3BD2%3D2%2AD1%3Bdelta%3D1e-4%3Bwhile+abs%28%28D2-D1%29%2FD2%29%3Edeltau%3D%28G%2Frho%29%2F%28pi%2AD1)
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')
求采纳为满意回答.