关于matlab的代码分析

关于matlab的代码分析

y=0.14085;
u=0;
c=0.3144;
dc=0.354;
hmin=-0.25679;
hmax=0.24321;
s=0.4150;
deltau0=0.5;
a0=0.404;

这段代码定义一些参数，包括常数y、u、c、dc、hmin、hmax、s、deltau0和a0，并初始化它们的值。

syms z b;

这句代码使用MATLAB中的符号变量函数syms定义符号变量z和b，以便在后续计算中使用。

f1=  (1/(b-z)^7 - 2*(c^6)/(b-z)^13)*(exp((-(z - u)^2)/(2*(y^2))));
f2=  (1/(dc^7) - 2*(c^6)/(dc^13))*(exp((-(z - u)^2)/(2*(y^2))));

这两行代码计算了法向载荷与界面间距Fn的分子部分，其中f1表示在界面上方的载荷部分，f2表示在界面下方的载荷部分。

b1=int(f1,z,hmin,b-dc);
b2=int(f2,z,b-dc,hmax);

这两行代码使用MATLAB中的积分函数int来计算Fn的分母部分的积分值，其中b1计算界面上方的Fn分母部分的积分值，b2计算界面下方的Fn分母部分的积分值。

Fn=((9.577*s*(c^6))/y)*(b1+b2);

这句代码通过将Fn的分子部分与分母部分相乘并除以y得到Fn的表达式。

figure;
fplot(b,Fn,[-0.08,5]);
xlabel('x');
ylabel('Fn');
legend('Fn');

这段代码先使用MATLAB中的figure函数创建一个新的图形窗口，然后使用fplot函数画出平均摩擦力图，其中b表示x轴，Fn表示y轴，[-0.08,5]指定了x轴的范围。接着使用xlabel和ylabel函数设置x轴和y轴的标签为’x’和’Fn’，最后使用legend函数在图例中显示’Fn’。