国科大《Python科学计算与数据处理》作业3及参考答案(SciPy).docxVIP

1.求解非线性方程组，cosa=1-d22r2

importnumpyasnp

fromscipy.optimizeimportroot

d=140

L=156

defequations(x):

a,r=x

f1=np.cos(a)-1+d**2/(2*r**2)

f2=a*r-L

return[f1,f2]

defjacobian(x):

a,r=x

return[[-np.sin(a),-d**2/r**3],[r,a]]

x0=[1.6,97.5]

#无雅可比

sol1=root(equations,x0)

print(无雅可比:,sol1.x,|函数调用次数:,sol1.nfev)

#带雅可比

sol2=root(equations,x0,jac=jacobian)

print(带雅可比:,sol2.x,|函数调用次数:,sol2.nfev)

运行结果：

无雅可比:[1.594063897|函数调用次数:12

带雅可比:[1.594063897|函数调用次数:10

2.用curve_fit()函数对高斯分布进行拟合，x∈[0,