习题7.10

合集 - python数学建模课后习题(41)

1.习题5.4求解非凸非线性规划2024-10-112.习题5.5求非凸函数的线性规划问题2024-10-123.习题6.42024-10-144.习题6.32024-10-145.习题6.22024-10-146.习题5.7 生产计划问题2024-10-147.习题4.3 求4.3节模型二的解2024-10-158.习题2.12024-10-159.习题2.22024-10-1510.习题2.32024-10-1511.习题2.42024-10-1512.习题2.5（1）2024-10-1513.习题2.5（2）2024-10-1514.习题2.62024-10-1515.习题2.7（1）2024-10-1516.习题2.7（2）2024-10-1517.习题2.82024-10-1518.习题2.92024-10-1519.习题2.102024-10-1520.习题2.112024-10-1521.习题2.122024-10-1522.习题2.132024-10-1523.习题3.22024-10-1524.习题3.32024-10-1525.习题6.52024-10-2626.习题6.62024-10-2627.习题7.12024-11-1228.习题7.32024-11-1229.习题7.42024-11-1230.习题7.72024-11-12

31.习题7.102024-11-12

32.习题8.42024-11-1233.习题8.52024-11-1234.8.8最优捕鱼策略2024-12-0935.8.92024-12-0936.9.2车间生产滚珠2024-12-0937.9.32024-12-0938.9.52024-12-0939.习题10.12024-12-2140.习题10.22024-12-2141.习题10.32024-12-21

收起

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from scipy.interpolate import interp1d, PchipInterpolator, CubicSpline
from scipy.optimize import curve_fit
from scipy.stats import norm

file_path = '7.17.xlsx'
data = pd.read_excel(file_path, header=None)
x_data = data.iloc[:, 0].values
y_data = data.iloc[:, 1].values

x_interp = np.linspace(-2, 4.9, 400)

interp_functions = {
    'Linear': interp1d(x_data, y_data, kind='linear', bounds_error=False, fill_value='extrapolate'),
    'Cubic': interp1d(x_data, y_data, kind='cubic', bounds_error=False, fill_value='extrapolate'),
    'PCHIP': PchipInterpolator(x_data, y_data, extrapolate=True),
    'CubicSpline': CubicSpline(x_data, y_data, bc_type='natural', extrapolate=True)
}

y_interps = {name: func(x_interp) for name, func in interp_functions.items()}

plt.figure(figsize=(10, 6))
for name, y_interp in y_interps.items():
    plt.plot(x_interp, y_interp, label=name)
plt.plot(x_data, y_data, 'o', label='Original Data')
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.title('Interpolation Methods Comparison')
plt.grid(True)
plt.show()

degrees = range(1, 6)
coeffs = {}
residuals = {}
for degree in degrees:
    coefficients, residual = np.polyfit(x_data, y_data, degree, cov=True)
    residual_std = np.sqrt(residual[0, 0])
    coeffs[degree] = coefficients
    residuals[degree] = residual_std

for degree, coeffs_val in coeffs.items():
    print(f"Degree {degree} Polynomial Coefficients: {coeffs_val}")
    print(f"Residual Standard Deviation: {residuals[degree]:.4f}")

best_degree = min(residuals, key=residuals.get)
print(f"Best fitting polynomial degree: {best_degree}")

best_poly = np.poly1d(coeffs[best_degree])
y_poly_fit = best_poly(x_interp)

plt.figure(figsize=(10, 6))
plt.plot(x_interp, y_poly_fit, label=f'{best_degree}-degree Polynomial Fit')
plt.plot(x_data, y_data, 'o', label='Original Data')
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.title('Polynomial Fit')
plt.grid(True)
plt.show()

def normal_dist(x, mu, sigma):
    return norm.pdf(x, mu, sigma)

params, params_covariance = curve_fit(normal_dist, x_data, y_data, p0=[np.mean(x_data), np.std(x_data)])
mu, sigma = params

x_normal_fit = np.linspace(min(x_data), max(x_data), 400)
y_normal_fit = normal_dist(x_normal_fit, mu, sigma)

plt.figure(figsize=(10, 6))
plt.plot(x_normal_fit, y_normal_fit, label=f'Normal Distribution Fit\nmu={mu:.2f}, sigma={sigma:.2f}')
plt.plot(x_data, y_data, 'o', label='Original Data')
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.title('Normal Distribution Fit')
plt.grid(True)
plt.show()

print("学号：3008")

结果如下图所示