【2024国赛美赛数学建模常用方法简介，小白也能M奖+！】

程序猿乖乖

于 2024-07-06 11:09:26 发布

阅读量540

点赞数 11

分类专栏：数学建模文章标签：数学建模统一建模语言云计算

本文链接：https://blog.csdn.net/m0_67414667/article/details/140225819

版权

数学建模专栏收录该内容

1 篇文章 0 订阅

订阅专栏

1. 统计建模

线性回归

概述：
线性回归是一种用于探索因变量与一个或多个自变量之间线性关系的统计方法。

示例代码：
使用Python的scikit-learn库进行线性回归建模：

from sklearn.linear_model import LinearRegression
import numpy as np

# 准备数据
X = np.array([[1, 1], [1, 2], [2, 2], [2, 3]])
y = np.dot(X, np.array([1, 2])) + 3

# 创建线性回归模型
model = LinearRegression()

# 拟合模型
model.fit(X, y)

# 打印模型参数
print('Intercept:', model.intercept_)
print('Coefficients:', model.coef_)

应用场景：
预测房价、销售量等因变量与多个自变量线性相关的问题。

逻辑回归

概述：
逻辑回归是一种用于处理二分类问题的统计方法，通过计算概率来预测事件的发生概率。

示例代码：
使用Python的scikit-learn库进行逻辑回归建模：

from sklearn.linear_model import LogisticRegression
import numpy as np

# 准备数据
X = np.array([[1, 2], [2, 3], [3, 4], [4, 5]])
y = np.array([0, 0, 1, 1])

# 创建逻辑回归模型
model = LogisticRegression()

# 拟合模型
model.fit(X, y)

# 打印模型参数
print('Intercept:', model.intercept_)
print('Coefficients:', model.coef_)

应用场景：
预测客户是否流失、疾病风险等二分类问题的建模。

2. 优化建模

线性规划

概述：
线性规划是一种用于解决最大化或最小化线性目标函数在线性约束条件下的优化问题的方法。

示例代码：
使用Python的PuLP库进行线性规划建模：

from pulp import *

# 创建问题
prob = LpProblem("Simple_LP_Problem", LpMaximize)

# 定义决策变量
x = LpVariable("x", lowBound=0)
y = LpVariable("y", lowBound=0)

# 定义目标函数
prob += 3*x + 2*y

# 定义约束条件
prob += 2*x + y <= 8
prob += x + 2*y <= 6

# 求解问题
prob.solve()

# 输出结果
print("Optimal Solution:")
print("x =", value(x))
print("y =", value(y))

应用场景：
资源分配、生产优化等问题的求解。

非线性规划

概述：
非线性规划是一种处理目标函数或约束条件包含非线性项的优化问题的方法。

示例代码：
使用Python的SciPy库进行非线性规划建模：

from scipy.optimize import minimize

# 定义目标函数和约束条件
def objective(x):
    return x[0]**2 + x[1]**2

def constraint1(x):
    return x[0] + x[1] - 3

# 定义初始点
x0 = [1, 1]

# 定义约束条件类型
con1 = {'type': 'eq', 'fun': constraint1}

# 求解非线性规划问题
sol = minimize(objective, x0, constraints=[con1])

# 输出结果
print("Optimal Solution:")
print("x1 =", sol.x[0])
print("x2 =", sol.x[1])

应用场景：
最优化投资组合、工艺参数优化等包含非线性约束或目标函数的问题求解。

3. 概率与统计建模

贝叶斯统计

概述：
贝叶斯统计是一种推断理论，用来处理不确定性，并基于先验知识和数据更新后验概率。

示例代码：
使用Python的PyMC3库进行贝叶斯统计建模：

import pymc3 as pm
import numpy as np

# 准备数据
data = np.random.normal(0, 1, size=100)

# 定义模型
with pm.Model() as model:
    mu = pm.Normal('mu', mu=0, sigma=1)
    sigma = pm.HalfNormal('sigma', sigma=1)
    obs = pm.Normal('obs', mu=mu, sigma=sigma, observed=data)
    
    # 求解模型
    trace = pm.sample(1000, tune=1000)

# 打印结果
print(pm.summary(trace))

应用场景：
风险评估、信用评分、医学诊断等需要结合先验知识和数据进行推断的场景。

马尔可夫链

概述：
马尔可夫链是一种描述状态和状态转移概率的数学工具，用于建模随机过程。

示例代码：
使用Python的numpy库生成简单的马尔可夫链：

import numpy as np

# 定义状态空间和转移矩阵
states = ['Rainy', 'Sunny']
transition_matrix = np.array([[0.7, 0.3],
                              [0.4, 0.6]])

# 模拟马尔可夫链
def generate_weather(days):
    current_state = np.random.choice(states)
    weather = [current_state]
    for _ in range(days - 1):
        current_state = np.random.choice(states, p=transition_matrix[states.index(current_state)])
        weather.append(current_state)
    return weather

# 模拟天气状态变化
weather_sequence = generate_weather(7)
print("Generated Weather Sequence:")
print(weather_sequence)

应用场景：
天气预测、股票价格模拟等依赖当前状态影响未来状态的问题建模。

4. 机器学习与人工智能建模

监督学习 - 支持向量机（SVM）

概述：
支持向量机是一种强大的监督学习算法，用于分类和回归问题，通过找到最佳的分隔超平面来实现分类。

示例代码：
使用Python的scikit-learn库进行支持向量机分类建模：

from sklearn import svm
import numpy as np

# 准备数据
X = np.array([[0, 0], [1, 1]])
y = np.array([0, 1])

# 创建SVM分类器
clf = svm.SVC(kernel='linear')

# 拟合模型
clf.fit(X, y)

# 打印预测结果
print("Predicted Class for [2, 2]:", clf.predict([[2, 2]]))

应用场景：
图像识别、文本分类、生物信息学等分类问题的建模和预测。

无监督学习 - K均值聚类

概述：
K均值聚类是一种常用的无监督学习算法，用于将数据点分成K个簇，使得每个数据点属于距离最近的簇。

示例代码：
使用Python的scikit-learn库进行K均值聚类建模：

from sklearn.cluster import KMeans
import numpy as np

# 准备数据
X = np.array([[1, 2], [1, 4], [1, 0],
              [4, 2], [4, 4], [4, 0]])

# 创建K均值聚类模型
kmeans = KMeans(n_clusters=2, random_state=0)

# 拟合模型
kmeans.fit(X)

# 打印聚类结果
print("Cluster Centers:", kmeans.cluster_centers_)
print("Cluster Labels:", kmeans.labels_)