0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

如何利用Matlab的lsqnonlin函数处理复杂函数的拟合问题

冬至子 来源:FEM and FEA 作者:追逐繁星的Mono 2023-07-10 11:42 次阅读

1. 简介

在工程问题的计算中,我们经常需要处理一些离散数据的拟合问题,而最小二乘法是处理曲线拟合问题的常用方法。目前,许多软件都提供有基于最小二乘法进行曲线拟合的功能,例如在Origin和Excel中均可直接利用离散数据进行曲线拟合。然而,这些软件只能处理一些简单函数的拟合问题,当需要拟合的函数较为复杂时,或者无法用简单的表达式来表述时,则往往无法直接进行拟合。为此,本文将对最小二乘法的基本原理做简单介绍,随后介绍如何利用Matlablsqnonlin函数处理复杂函数的拟合问题。

1.1 曲线拟合的最小二乘法原理

利用最小二乘法进行曲线拟合的本质为寻找某个近似函数 φ ( x ),使得该函数与离散点之间尽可能地逼近。若将偏差定义为近似函数的近似值 φ ( xi )与离散点 yi *之间的差值:

1.jpg

求解上述线性方程组即可得到拟合多项式的系数。

2. 利用Matlab处理曲线拟合问题

基于上述计算原理,Matlab提供了polyfit函数用于处理多项式曲线的拟合问题,对于一些较为复杂但仍可通过简单表达式进行表述的函数,也可以利用Matlab的拟合工具箱(Curve fitting)进行拟合。但在某些情况,当拟合函数非常复杂,以致于无法用简单表达式进行表述时(例如分段函数以及涉及到条件语句),则无法使用拟合工具箱进行拟合。对于此类问题,可以使用Matlab优化工具箱中的lsqnonlin函数进行解决。

2.1 lsqnonlin函数

lsqnonlin函数用于求解以下述形式表示的非线性最小二乘法拟合问题:

在使用该函数进行最小二乘法拟合时,lsqnonlin函数并不需要用户提供min || f ( x )||(平方和),而是需要用户提供自定义函数fun,用于计算矢量形式表示的 f ( x ):

lsqnonlin函数常用语法为:

x = lsqnonlin(fun,x0)
x = lsqnonlin(fun,x0,lb,ub)

其中fun为用户自定义函数,x0为计算采用的初始值,lsqnonlin函数首先利用x0通过自定义函数fun计算 fi (x)的取值并计算平方和,随后通过优化算法调整x的取值直至得到平方和的最小值。此外,lb和ub还可以用于定义x的取值范围,使得x满足lb≤x≤ub。

例如,对于节1.1中所述的多项式,根据最小二乘法的定义,则自定义函数 f ( x )应表示为:

注意此时 f ( x )中的x为以向量形式表示的多项式 P ( x )的系数:

1.jpg

在计算时,用户需要指定多项式系数的初始值,则lsqnonlin函数将利用最小二乘法计算多项式系数。

下面,本文将以笔者所在领域常用的NASGRO方程为例,介绍如何利用lsqnonlin函数处理此类复杂函数的曲线拟合问题。

2.2 NASGRO方程简介

在进行基于断裂力学的损伤容限分析时,应力强度因子和裂纹扩展速率模型是最为重要的输入。一般来说,应力强度因子可以通过经验公式或数值方法进行计算,而裂纹扩展速率模型则需要通过裂纹扩展速率试验获得的试验数据拟合得到。例如,大量的试验结果表明,在裂纹扩展的中速率区域,应力强度因子幅值ΔK和裂纹扩展速率d a /dN满足良好的对数线性关系,可以通过Paris公式进行描述:

1.jpg

其中Cm为材料常数。

尽管Paris公式已经得到广泛的应用,但是Paris公式仅仅描述了裂纹在中速率区域的扩展行为,没有描述近门槛区域和接近断裂的高速率区域的扩展行为,也没有考虑应力比R和裂纹闭合效应对裂纹扩展速率的影响,因此给出的计算结果将过于保守。另一个常用的裂纹扩展速率模型为Newman提出的NASGRO模型,该模型基于Forman模型改进了裂纹扩展速率模型,同时比Paris和Walker模型更加全面,不仅考虑了应力强度因子门槛值和断裂韧度,还体现了应力比以及裂纹闭合效应对裂纹扩展速率d a /dN的影响,如图2.1所示,其表达式如下:

其中R为应力比,ΔK为应力强度因子幅值,ΔKth为应力强度因子幅值门槛值,Kmax为最大应力强度因子,可表示为:

1.jpg

图片

图2.1 NASGRO方程

NASGRO方程中的应力强度因子门槛值ΔKth可采用下面的经验公式进行估算:

其中A0为裂纹张开函数中的多项式系数,ΔK1R =1时的应力强度因子门槛值,Cth是对于正应力(上标为p=positive)和负应力比(上标为m=minus, negative)取不同值的材料常数,*a*~0~是内在小裂纹尺寸(典型值为0.0381mm)。在基于NASGRO方程开发的疲劳裂纹扩展分析软件NASGRO中,正应力比下*C*~th~^p^和Δ*K*~1~是保存在数据库里的值,负应力比下*C*~th~^m^的默认值为0.1。

2.3 NASGRO方程拟合

图2.2为疲劳裂纹扩展分析软件NASGRO材料库中某铝合金材料的裂纹扩展速率数据,已知试验时采用的试样为中心平板试样(M(T)),σmaxσF的比值为0.3,塑性约束因子α为2.0,材料断裂韧度Kc为65.7,应力比R为1时的门槛值ΔK1为1.23,C th ^p^为1.06,C th ^m^为0.1,下面需要通过拟合试验数据获得NASGRO方程的参数 C , m , p , q

图片

图2.2 疲劳裂纹扩展数据

拟合NASGRO方程的难点主要有以下几点:

(1)裂纹扩展速率d a /dN不仅与应力强度因子幅值ΔK有关,还与使用的应力比R有关,因此实际上为多变量的拟合问题;

(2)裂纹张开函数f为分段函数,并且使用了计算最大值的max函数,该函数在拟合时无法用简单函数进行表述。

针对以上问题,NASGRO软件给出的拟合方法为首先给参数pq确定一个初始值,并利用最小二乘法确定参数Cm ,随后根据工程经验来获得可接受的结果,如果对拟合效果不满意,可以调整任意参数,直至获得满意的结果。

显然,这样的拟合策略具有很大的随意性,如果参数pq选取不当,很可能对拟合效果有很大的影响。下面,本文将介绍如何利用lsqnonlin函数在不提前定义参数pq的情况下对NASGRO方程进行拟合。

根据lsqnonlin函数的介绍,首先需要构造自定义函数 f ( x )使其满足最小二乘法计算的基本原理,由于Paris公式具有对数线性的关系,因此尝试将NASGRO方程两边取对数,可得:

上式可以用如下所示的通式表示:

系数bj为与 C , n , pq有关(b 0 =log( C ), b 1 = n , b 2 = p , b 3 =- q )的系数,而gj为与Δ KR和NASGRO中所有剩余参数有关的函数。

根据最小二乘法的定义,应选取参数bj使得:

参考lsqnonlin函数对目标函数的定义,则自定义函数 f ( x )应表示为:

1.jpg

y= R .^2 * (R >0) + R * (R <= 0)

而裂纹张开函数f中涉及到求取最大值的计算以及分段函数的处理,也可以通过上述语法实现,具体的计算过程可参见程序代码(参见附录)。

此外,由于自定义函数 f ( x )为关于系数 bj ( j =1,2,3,4)的函数,为了将试验数据(不同应力比R下的应力强度因子幅值ΔK和裂纹扩展速率d a /d N )传递到函数 f ( x )中进行计算,可以将试验数据定义为全局变量,以便被 f ( x )调用。

通过编写程序,可以计算得到NASGRO方程的系数如表1所示。

图片

拟合曲线与试验数据如图2.3所示。

图片

图2.3 试验数据及拟合曲线对比

附录1 NASGRO方程曲线拟合程序

NASGRO_LSQ.m

NASGRO_LSQ用于定义采用最小二乘法拟合NASGRO方程时的自定义函数 f ( x ),输入参数Coeff为NASGRO方程系数 bj ,输出参数为拟合函数与试验数据误差的平方和。

function F=NASGRO_LSQ(Coeff)
%程序用于计算最小二乘法拟合NASGRO方程的目标函数
%程序返回一个N×1的数值,其中N为数据对的个数
%Coeff为拟合时待求的系数(共4个系数)
%4个系数分别为log(C)、n、p和-q
%DataX(:,1)为应力强度因子幅值
%DataX(:,2)为应力比
%DataY为裂纹扩展速率


%*********全局变量传递**************
global S_max_flow alpha DKth1 a0 Cth_p Cth_m a Kc
global DataX DataY
%S_max_flow为施加的最大应力与流动应力的比值
%alpha为塑性约束因子
%DKth1为应力比为1时对应的门槛值
%a0为与门槛值有关的常数
%Cth_p为正应力比下与门槛值有关的常数
%Cth_m为负应力比下与门槛值有关的常数
%a为计算门槛值时使用的裂纹长度,建议取为远大于a0的值
%Kc为材料断裂韧度
%DataX为应力强度因子幅值
%DataY为裂纹扩展速率
%***********************************


%********Newman裂纹张开函数计算**********
R=DataX(:,2);    %应力比
DK=DataX(:,1);   %应力强度因子幅值
%计算系数A0(与应力比和应力强度因子幅值无关)
A0=(0.825-0.34*alpha+0.05*alpha^2)*...
    (cos(pi/2*S_max_flow))^(1/alpha);
A1=(0.415-0.071*alpha)*S_max_flow;
A3=2*A0+A1-1;
A2=1-A0-A1-A3;
%计算向量形式的裂纹张开函数
f1=max(A0+A1*R+A2*R.^2+A3*R.^3,R);
f2=A0-2*A1;
f3=A0+A1*R;
f=f1.*(R >=0)+f2.*(R< -2)+...
    f3.*(R >=-2&R< 0);  %裂纹张开函数
%****************************************


%********应力强度因子门槛值计算**********
DKth_p1=DKth1*sqrt(a/(a+a0))*((1-R)./(1-f)).^(1+R*Cth_p)./...
        (1-A0).^((1-R)*Cth_p);        %正应力比下的门槛值
DKth_p2=DKth1*sqrt(a/(a+a0))*((1-R)./(1-f)).^(1+R*Cth_m)./...
        (1-A0).^(Cth_p-R*Cth_m);      %负应力比下的门槛值
DKth=DKth_p1.*(R >=0)+...
    DKth_p2.*(R< 0);  %应力强度因子门槛值
%****************************************


%******根据NASGRO方程计算函数F***********
F1=log10((1-f)./(1-R).*DK);   %DataX(1,:)为应力强度因子幅值
F2=log10(1-DKth./DK);
F3=log10(1-1./(1-R).*(DK./Kc));
%****************************************


%******根据NASGRO方程计算裂纹扩展速率***********
y=Coeff(1)+Coeff(2)*F1+Coeff(3)*F2+Coeff(4)*F3;
%***********************************************


%*****构造基于最小二乘法的目标函数F**************
%最小二乘法应保证目标函数F中所有原始之和达到最小
F=(y-log10(DataY)).^2;
%***********************************************
end

NASGRO_nonfit.m

NASGRO_nonfit用于读取处理试验数据、定义材料参数以及调用lsqnonlin函数进行曲线拟合,并通过拟合得到的参数绘制NASGRO方程曲线。

%程序用于非线性拟合NASGRO方程的系数
clear,clc
%***************计算参数输入*****************
global S_max_flow alpha DKth1 a0 Cth_p Cth_m a Kc  %定义全局变量
global DataX DataY 
S_max_flow=0.3;         %最大应力与流动应力的比值
alpha=2;                %塑性约束因子
Kc=65.7;                %材料断裂韧度
DKth1=1.23;            %应力比R=1对应的门槛值
Cth_p=1.06;             %应力比为正时的门槛值计算参数
Cth_m=0.1;              %应力比为负时的门槛值计算参数
a0=0.0381;              %门槛值计算参数
a=a0*1000;              %计算门槛值时采用的裂纹长度
%************************************************


%**************试验数据导入*****************
%数据格式:应力强度因子,裂纹扩展速率,应力比
FCG_data=importdata('FCG_data.txt');
%数据处理
data_num=size(FCG_data,1);         %数据对总个数
FCG_data=sortrows(FCG_data,3);     %将数据按照应力比顺排  
tbl=tabulate(FCG_data(:,3));       %数据统计
R_eval=tbl(:,1);                        %试验出现的应力比(统计数据)
R_num=size(R_eval,1);                   %试验使用的应力比数量
DataX(:,1)=FCG_data(:,1);               %应力强度因子幅值
DataX(:,2)=FCG_data(:,3);               %应力比
DataY=FCG_data(:,2);                    %裂纹扩展速率
clear tbl
%********************************************




%*************最小二乘法拟合*****************
Coeff0=[-9 3 1 -1];      %拟合方程系数的迭代初始值
[Coeff,resnorm,residual,exitflag,output]=lsqnonlin(@NASGRO_LSQ,Coeff0);
%NASGRO方程系数输出
fprintf('NASGRO方程系数C为%s\\n',10^Coeff(1));
fprintf('NASGRO方程系数m为%s\\n',Coeff(2));
fprintf('NASGRO方程系数p为%s\\n',Coeff(3));
fprintf('NASGRO方程系数q为%s\\n',-1*Coeff(4));
%********************************************




%*************绘制试验数据及拟合曲线***************
Output_DK=cell(R_num,1);       %用于输出数据的元胞数组初始化
Output_DaDn=cell(R_num,1);
%绘制试验数据
figure
for i=1:R_num     %试样采用的应力比数量
    data_index=find(FCG_data(:,3)==R_eval(i));
    loglog(FCG_data(data_index,1),FCG_data(data_index,2),'o');
    hold on
end
%绘制拟合曲线
for i=1:R_num     
%计算裂纹张开函数
R=R_eval(i);     %应力比
[f,A0,~]=Crack_opening_f(R,S_max_flow,alpha);
%计算应力强度因子门槛值
[DKth]=...
    DK1th_calc(DKth1,R,f,A0,a0,Cth_p,Cth_m,a);
%计算用于绘制NASGRO方程的应力强度因子幅值
%应力强度因子幅值应大于门槛值并且最大应力强度因子应小于断裂韧度
DK=linspace(DKth,(1-R)*Kc,500);   %应力强度因子幅值
%计算裂纹扩展速率
C=10^Coeff(1);
m=Coeff(2);
p=Coeff(3);
q=-1*Coeff(4);
Crack_rate=C*(((1-f)/(1-R))*DK).^m.*...
    (1-(DKth./DK)).^p./...
    (1-(1/(1-R)*(DK/Kc))).^q;
loglog(DK,Crack_rate)
%拟合曲线数据输出
Output_DK{i}=DK';
Output_DaDn{i}=Crack_rate';
end
legend_str=cell(2*R_num,1);   %定义图例的元胞数组预分配
for i=1:R_num     
    legend_str{i}=['R=' num2str(R_eval(i))];    %获取应力比用于图例显示
    legend_str{i+R_num}=legend_str{i};
end
legend(legend_str);
title('NASGRO方程拟合');
xlabel('应力强度因子幅值ΔK');ylabel('裂纹扩展速率da/dN');
hold off
%**************************************************

附录2 NASGRO方程绘制程序

程序仅用于计算绘制NASGRO方程时需要的参数,未参与自定义函数的计算。

Crack_opening_f.m

function [f,A0,A1,A2,A3] =...
    Crack_opening_f(R,S_max_flow,alpha)
%程序用于计算Newman裂纹张开函数
%R为应力比
%S_max_flow为最大应力与流动应力的比值S_max/S_flow
%alpha为塑性约束因子


A0=(0.825-0.34*alpha+0.05*alpha^2)*...
    (cos(pi/2*S_max_flow))^(1/alpha);
A1=(0.415-0.071*alpha)*S_max_flow;
A3=2*A0+A1-1;
A2=1-A0-A1-A3;


if R >=0
    f=A0+A1*R+A2*R^2+A3*R^3;
    f=max([f;R]);
elseif R< -2
    f=A0-2*A1;
else
    f=A0+A1*R;
end
end

DK1th_calc.m

function [DKth]=...
    DK1th_calc(DKth1,R,f,A0,a0,Cth_p,Cth_m,a)
%程序通过应力比为1的门槛值计算不同应力比下的门槛值
%DKth1为应力比为1时对应的门槛值
%R为应力比
%f为裂纹张开函数,A0为裂纹张开函数中的多项式系数
%a0为材料常数,Cth_p和Cth_m分别为正负应力比下的材料常数
%a为裂纹长度


if R >=0
    DKth=DKth1*sqrt(a/(a+a0))*((1-R)/(1-f))^(1+R*Cth_p)/...
        (1-A0)^((1-R)*Cth_p);
else
    DKth=DKth1*sqrt(a/(a+a0))*((1-R)/(1-f))^(1+R*Cth_m)/...
        (1-A0)^(Cth_p-R*Cth_m);
end
end

试验数据

FCG_data.txt

应力强度因子幅值 裂纹扩展速率 应力比

11.9896  2.43E-05  0.7
11.6632  2.51E-05  0.7
10.5603  1.71E-05  0.7
10.916  1.43E-05  0.7
10.5604  1.43E-05  0.7
10.1043  1.10E-05  0.7
9.25035  8.18E-06  0.7
8.89987  5.91E-06  0.7
8.61001  5.10E-06  0.7
8.46855  6.09E-06  0.7
8.05821  4.95E-06  0.7
7.83889  4.40E-06  0.7
7.41798  3.91E-06  0.7
7.54187  3.48E-06  0.7
6.98103  3.38E-06  0.7
6.86652  1.99E-06  0.7
6.25152  1.99E-06  0.7
6.01454  2.44E-06  0.7
6.11514  1.17E-06  0.7
5.66032  1.62E-06  0.7
5.35634  1.82E-06  0.7
5.47599  1.24E-06  0.7
5.18191  1.39E-06  0.7
5.32693  1.24E-06  0.7
4.93077  1.44E-06  0.7
5.09686  1.07E-06  0.7
4.48922  9.23E-07  0.7
4.61486  8.20E-07  0.7
4.43997  7.73E-07  0.7
4.36708  6.48E-07  0.7
4.06469  6.87E-07  0.7
4.64051  4.41E-07  0.7
4.82333  4.04E-07  0.7
3.59998  4.04E-07  0.7
3.78335  3.49E-07  0.7
3.8678  3.70E-07  0.7
3.91075  3.39E-07  0.7
3.86786  2.38E-07  0.7
3.62002  1.72E-07  0.7
3.31402  1.93E-07  0.7
3.2238  2.05E-07  0.7
3.20604  2.38E-07  0.7
3.10161  2.24E-07  0.7
3.05068  2.05E-07  0.7
2.83944  2.31E-07  0.7
2.6722  1.82E-07  0.7
2.59947  1.57E-07  0.7
2.61388  1.24E-07  0.7
2.68703  1.17E-07  0.7
2.7929  9.53E-08  0.7
2.31499  1.28E-07  0.7
2.315  1.10E-07  0.7
2.34071  9.81E-08  0.7
2.26445  1.10E-07  0.7
2.20281  1.10E-07  0.7
2.00554  8.22E-08  0.7
1.96175  7.75E-08  0.7
1.92954  7.10E-08  0.7
1.88741  6.50E-08  0.7
1.86669  6.31E-08  0.7
1.82592  6.89E-08  0.7
1.74707  4.98E-08  0.7
1.70894  3.93E-08  0.7
1.68088  3.60E-08  0.7
1.64418  3.50E-08  0.7
1.63523  6.32E-09  0.7
1.58198  5.00E-09  0.7
43.3791  6.61E-04  0
39.9322  6.61E-04  0
35.9575  3.36E-04  0
33.2835  3.36E-04  0
31.4968  1.97E-04  0
31.3239  1.43E-04  0
28.8351  1.23E-04  0
28.6768  8.15E-05  0
26.6911  9.17E-05  0
26.5444  6.83E-05  0
24.9811  4.94E-05  0
23.7706  4.02E-05  0
23.2517  3.37E-05  0
22.6188  3.37E-05  0
21.7614  3.90E-05  0
20.8215  3.57E-05  0
20.8217  2.51E-05  0
21.4044  2.36E-05  0
21.7617  2.66E-05  0
19.9226  1.81E-05  0
19.062  1.71E-05  0
18.6457  1.98E-05  0
19.2735  2.16E-05  0
19.5951  2.51E-05  0
17.7421  1.76E-05  0
17.9392  1.52E-05  0
17.0698  1.61E-05  0
16.6051  1.61E-05  0
16.3325  1.47E-05  0
16.5139  1.27E-05  0
15.888  1.24E-05  0
15.7137  1.10E-05  0
15.6275  6.65E-06  0
14.465  1.16E-05  0
14.4652  9.20E-06  0
14.3065  6.65E-06  0
13.8403  9.76E-06  0
13.8402  1.13E-05  0
13.2425  8.18E-06  0
13.3159  7.48E-06  0
13.2427  6.09E-06  0
12.5315  7.27E-06  0
12.3939  7.71E-06  0
12.1233  6.85E-06  0
11.9902  6.27E-06  0
12.5316  5.10E-06  0
11.7285  4.95E-06  0
10.8563  5.91E-06  0
10.7371  5.41E-06  0
10.7966  4.95E-06  0
10.6193  4.67E-06  0
10.3875  4.15E-06  0
9.99386  3.28E-06  0
9.40513  3.28E-06  0
9.30192  2.91E-06  0
9.25075  2.67E-06  0
8.99896  2.51E-06  0
8.90013  2.75E-06  0
8.19305  1.93E-06  0
7.97006  1.71E-06  0
7.79604  1.62E-06  0
7.796  1.87E-06  0
7.54204  1.87E-06  0
6.9813  1.24E-06  0
7.05883  1.07E-06  0
6.53398  7.96E-07  0
6.829  7.08E-07  0
6.67993  5.76E-07  0
6.49812  5.27E-07  0
6.32121  5.93E-07  0
6.08171  4.29E-07  0
6.32129  4.16E-07  0
6.32131  3.81E-07  0
6.35634  3.29E-07  0
5.85128  3.10E-07  0
5.66068  3.01E-07  0
5.75522  2.67E-07  0
5.26876  2.67E-07  0
5.26878  2.38E-07  0
5.41621  2.38E-07  0
5.35678  2.05E-07  0
5.35682  1.72E-07  0
5.29803  1.57E-07  0
5.21103  1.53E-07  0
5.15379  1.72E-07  0
5.21097  1.99E-07  0
5.09716  2.24E-07  0
5.06913  1.99E-07  0
4.95844  1.99E-07  0
4.82348  1.82E-07  0
4.71814  1.93E-07  0
4.74422  2.31E-07  0
4.58968  2.38E-07  0
4.61513  1.77E-07  0
4.9312  1.48E-07  0
4.6407  1.57E-07  0
4.66641  1.36E-07  0
4.48956  1.28E-07  0
4.39151  1.40E-07  0
4.29561  1.40E-07  0
4.31936  1.72E-07  0
4.29556  1.99E-07  0
4.20179  1.62E-07  0
4.20181  1.40E-07  0
4.24848  1.21E-07  0
4.15571  1.21E-07  0
4.11007  1.36E-07  0
4.06492  1.57E-07  0
4.04256  1.40E-07  0
3.97618  1.40E-07  0
3.9325  1.67E-07  0
3.84666  1.32E-07  0
3.88939  1.04E-07  0
3.88942  8.47E-08  0
4.04264  8.47E-08  0
3.70087  1.32E-07  0
3.7009  1.10E-07  0
3.60016  1.10E-07  0
3.62006  1.32E-07  0
3.52153  1.24E-07  0
3.52155  1.04E-07  0
3.38806  1.21E-07  0
3.50218  8.98E-08  0
3.6002  8.22E-08  0
3.5022  7.75E-08  0
3.42573  7.53E-08  0
3.36947  7.53E-08  0
3.31411  9.25E-08  0
3.2239  9.25E-08  0
3.18854  7.31E-08  0
3.25972  6.89E-08  0
3.31417  5.77E-08  0
3.36953  4.84E-08  0
3.08471  5.13E-08  0
3.08473  4.43E-08  0
3.1536  4.17E-08  0
3.31424  3.40E-08  0
3.20635  1.83E-08  0
3.01752  1.14E-08  0
3.01755  9.28E-09  0
3.17125  8.49E-09  0
2.96804  6.14E-09  0
3.08497  5.79E-09  0
3.10206  4.85E-09  0
3.05116  3.40E-09  0
2.87152  1.37E-09  0
2.79336  1.37E-09  0
21.6401  2.50E-04  0.5
21.8805  1.97E-04  0.5
20.8201  1.97E-04  0.5
19.9211  1.23E-04  0.5
20.7061  9.17E-05  0.5
17.9382  6.83E-05  0.5
17.9383  5.56E-05  0.5
17.6438  4.94E-05  0.5
16.6044  5.09E-05  0.5
16.6965  3.90E-05  0.5
15.8873  4.02E-05  0.5
15.6265  3.37E-05  0.5
15.0343  3.27E-05  0.5
14.5447  2.51E-05  0.5
14.3059  2.29E-05  0.5
13.0967  2.10E-05  0.5
12.0561  1.71E-05  0.5
12.3253  1.47E-05  0.5
12.4623  1.13E-05  0.5
11.7281  1.13E-05  0.5
11.1597  1.24E-05  0.5
11.0983  1.13E-05  0.5
10.5023  1.24E-05  0.5
10.5024  9.76E-06  0.5
10.33  9.20E-06  0.5
10.4447  8.42E-06  0.5
10.2732  7.71E-06  0.5
9.93858  6.65E-06  0.5
9.50927  6.65E-06  0.5
8.89981  7.06E-06  0.5
8.99867  5.74E-06  0.5
8.8022  5.25E-06  0.5
8.28374  4.27E-06  0.5
8.32966  3.38E-06  0.5
7.92594  3.91E-06  0.5
7.71025  3.28E-06  0.5
6.90432  4.40E-06  0.5
7.21618  2.51E-06  0.5
6.75371  2.37E-06  0.5
6.67961  1.99E-06  0.5
6.21708  2.23E-06  0.5
6.04791  1.77E-06  0.5
6.04792  1.66E-06  0.5
5.38609  1.10E-06  0.5
5.29767  8.96E-07  0.5
5.23953  8.20E-07  0.5
5.12513  7.73E-07  0.5
5.01323  7.08E-07  0.5
5.26863  4.97E-07  0.5
4.9038  5.93E-07  0.5
4.66611  7.29E-07  0.5
4.41561  4.68E-07  0.5
4.31924  3.59E-07  0.5
4.20166  3.59E-07  0.5
3.99806  3.19E-07  0.5
3.66015  2.45E-07  0.5
3.50207  2.18E-07  0.5
3.48282  1.77E-07  0.5
3.18843  1.77E-07  0.5
3.03396  1.17E-07  0.5
3.01723  1.53E-07  0.5
2.83949  1.44E-07  0.5
2.80833  1.32E-07  0.5
2.65752  1.36E-07  0.5
2.45989  1.48E-07  0.5
2.54274  1.14E-07  0.5
2.59952  9.53E-08  0.5
2.44638  1.04E-07  0.5
2.39295  1.14E-07  0.5
2.3798  9.25E-08  0.5
2.2896  8.98E-08  0.5
2.22728  9.53E-08  0.5
2.21503  8.47E-08  0.5
2.16666  8.47E-08  0.5
2.13107  9.25E-08  0.5
2.11936  7.53E-08  0.5
2.09609  7.10E-08  0.5
2.05032  7.10E-08  0.5
2.07309  6.31E-08  0.5
2.0731  5.77E-08  0.5
2.05034  5.77E-08  0.5
2.01667  5.77E-08  0.5
1.99453  5.28E-08  0.5
1.99454  4.84E-08  0.5
1.99454  4.56E-08  0.5
1.98357  4.17E-08  0.5
1.92958  3.82E-08  0.5
1.98358  3.71E-08  0.5
1.98361  2.68E-08  0.5
1.9727  2.38E-08  0.5
2.0279  2.18E-08  0.5
2.02791  2.06E-08  0.5
1.98367  1.14E-08  0.5
1.86688  4.71E-09  0.5
1.84641  3.21E-09  0.5
1.83625  3.12E-09  0.5
声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • MATLAB仿真
    +关注

    关注

    4

    文章

    176

    浏览量

    19900
  • 最小二乘法
    +关注

    关注

    0

    文章

    22

    浏览量

    8433
收藏 人收藏

    评论

    相关推荐

    关于labview拟合函数利用问题

    关于labview拟合函数的问题,具体描述如下:比如说有关于X,Y两组对应的数据,现在把X作为x坐标轴,Y作为y坐标轴拟合函数,怎么做到根据这个
    发表于 04-19 15:06

    labview调用matlab中的高斯拟合函数cftool

    本帖最后由 像风一样飞翔! 于 2019-11-28 16:32 编辑 matlab有强大的拟合工具包cftool,能否用labview调用matlab拟合
    发表于 11-28 16:28

    matlab自定义函数调用的方法

    matlab自定义函数调用的方法 命令文件/函数文件+ 函数文件 - 多
    发表于 11-29 13:14 88次下载

    matlab主要函数与命令

    matlab主要函数与命令
    发表于 06-19 14:18 29次下载
    <b class='flag-5'>matlab</b>主要<b class='flag-5'>函数</b>与命令

    matlab基本函数

    matlab基本函数 abs :绝对值   acos :反余弦   acosh :反双曲余弦函数   acot :反余切   acoth :反双曲
    发表于 06-18 14:49 1661次阅读

    函数的可视化与Matlab

    函数的可视化与Matlab作2.1 实验与观察:函数的可视化2.1.1 Matlab二维绘图命令1.周期函数与线性p-周期
    发表于 10-17 00:30 2188次阅读
    <b class='flag-5'>函数</b>的可视化与<b class='flag-5'>Matlab</b>作

    matlab取整函数命令代码

    本内容介绍了matlab取整函数。在matlab中有四个取整函数具体:Matlab取整函数有:
    发表于 12-20 14:41 4932次阅读

    Matlab数据拟合基础函数的使用

    Matlab数据拟合基础函数的使用,教你如何拟合两个未知量
    发表于 11-10 16:45 0次下载

    matlab牛逼函数总结

    matlab牛逼函数总结,想大家会因为MATLAB函数不懂吧!自己看看吧!
    发表于 11-06 17:18 0次下载

    MATLAB常用函数总结(表格)

    MATLAB常用函数总结,MATLAB函数速查手册,方便应用MATLAB函数
    发表于 01-21 14:31 0次下载

    matlab函数资料

    matlab 函数 matlab的一些基本的界面设计的函数
    发表于 03-23 17:13 0次下载

    MATLAB 常用函数参考

    matlab的常用函数,常用的函数都可以查。
    发表于 05-19 14:31 0次下载

    Sigmoid函数拟合法分析及其高效处理

    的影响。完成算法在FPGA上的硬件实现,给出一种使用三阶多项式处理Sigmoid函数拟合结果及流水线架构,处理精度达到10-5数量级,最大频率达到127.327 MHz,满足了高速、
    发表于 11-15 13:15 1.2w次阅读
    Sigmoid<b class='flag-5'>函数</b>的<b class='flag-5'>拟合</b>法分析及其高效<b class='flag-5'>处理</b>

    matlab中subplot函数的功能

    subplot是MATLAB中的函数,是将多个图画到一个平面上的工具。在matlab的命令窗口中输入doc subplot或者help subplot即可获得该函数的帮助信息。在
    发表于 11-17 17:44 2.8w次阅读
    <b class='flag-5'>matlab</b>中subplot<b class='flag-5'>函数</b>的功能

    MATLAB信号处理常用函数详解

    MATLAB是一款功能强大的数学软件,尤其在信号处理领域,它提供了众多的函数和工具箱,使得信号的分析、处理、仿真变得简单而高效。本文将详细介绍MAT
    的头像 发表于 05-17 14:31 1697次阅读