题目2:拟合
生物学家L.Schroeder Schr1 在研究天蛾幼虫的生长时采用了下面的数据确定w (活幼虫的重量,以克计算)和R (幼虫消耗的氧气,以毫升/小时计算)之间的关系*R=bw^a。
(1) 利用对数最小二乘方程 lnR = lnb + alnw 拟合,确定参数a, b 。
(2) 计算 (1) 中的平方误差。
(3) 修改 (1) 中的对数最小二乘方程lnR = lnb + alnw + c(lnw)^2 ,确定参数a, b, c 。
(4) 计算(3) 中的平方误差
(1)代码(python):
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import math
from numpy.linalg import inv
from numpy import dot
XY = pd.read_excel(r"data.xlsx")
w = np.array(XY['w'])
R = np.array(XY['R'])
x = np.log(w)
y = np.log(R)
#one = np.ones(34)
#X = np.array([x,one])
#Y = np.array([y])
#z = dot(dot