步骤如下。
函数定义
功能y=fitd (x，xdata ) )。
全球n r
fori=1:长度(xdata ) ) ) ) ) 0
n=1:N；
x0=1./n.^2. * exp (-n.^2* pi ^2* x * xdata (I )/R^2)；
y(I )=1-6/pi^2*sum(x0 )；
结束
结束
%%%%%%%
主程序
clear clc；
全球n r
data=[
0 0.00985；
0.21 0.298662；
0.42 0.326742；
1 0.369794；
4 0.438064；
7 0.488712；
10 0.518047；
14 0.536574；
18 0.570807；
22 0.593273；
26 0.603525；
30 0.610483；
35 0.627549；
42 0.657646；
49 0.689258；
60 0.741513]；
xdata=data (:2 )；
y data=data (:1 )；
x0=0.001；
N=30；
R=0.01；
x=lsqcurvefit(@fitd，x0，xdata，ydata ) )。
%%%%。 %%%。
x=
3.8190e-004