日志日志比例中最佳拟合线

11次观点(过去30天)
日元yap.
日元yap. 2021年8月24日
评论: 西蒙陈 2021年8月24日
嗨,我想创造一个 直的 在图中第一部分的最佳拟合线,我不想要曲线最佳贴合线,我该怎么办?谢谢你。
rho = 1000;%[kg / m3]
d = 12.6 * 10 ^ -3;%[m]
l = 1.5;%[m]
miu = 0.001;%[PAS]
g = 9.81;
a = pi * d ^ 2/4;
Q0 = [1600,1500,1400,1300,1200,1100,10,900,800,700,600,500,400,300,240,220,......
200,180,160,140,​​120,100,80,70,70,60,50,40,30,20,10];% (L /人力资源)
q = q0 /(1000 * 3600);
%湿湿数字计
P_DPG = [20.1,17.5,15.7,13.1,11.6,9.3,8,6.5,5.3,4.1,2.1,1.3,0.8];%[KPA]
%倒测压仪
h = [6.9,5.9,5,4.1,3.2,2.2.6,1.8,1.3,0.7,0.5];%[厘米]
h_m = h / 100;
p_mtr = rho * g * h_m;
%Capsuhelic仪表
p_cpg = [33,31,28,23,17,12,8];
p = [p_dpg。* 1000,p_mtr,p_cpg];
V = Q / A;
Re = rho * v * d / miu;
f = p * d ./(2 * rho * l * v. ^ 2);
x = re(25:31);
Y = F(25:31);
p = polyfit(x,y,1);
y = polyval(p,x);
图1)
loglog(re,f,'X''行宽'1)
抓住
loglog(x,y,' - '的)
网格
XLIM([10 ^ 2 10 ^ 5])
ylim ([0.001 - 0.1])
Xlabel('雷诺斯号码重新'的)
ylabel('摩擦因子F'的)
标题('f vs'的)

接受答案

西蒙陈
西蒙陈 2021年8月24日
像这样?
p = polyfit(log(x),log(y),1);
y = polyval(p,log(x));
图1)
loglog(re,f,'X''行宽'1)
抓住
loglog(x,exp(y),' - '的)
网格
2评论
西蒙陈
西蒙陈 2021年8月24日
仅仅因为您使用的是loglog刻度,所以您需要方程式:
(logy) = m(logx) + c,使函数polyfit适合于一条直线。
对于从(log x)计算变量y(但不是log(y))的多维亚,这是相同的,因此您需要在图中使用指数转换它。

登录评论。

更多答案(0)

标签

社区宝藏狩猎

找到Matlab Central中的宝藏,并发现社区如何帮助您!

开始狩猎!