以房子的面积X作为输入特征，对房子的总价Y建立线性模型(写出模型表达式) ，并给出损失函数和梯度函数（写出梯度函数），并且在训练集上通过梯度下降求解参数，画出损失函数和迭代次数的图，其中面积和总价已知，给我python代码

在Python中，我们可以使用scikit-learn库来创建一个简单的线性回归模型，假设我们的数据是一个二维数组`X`表示房屋面积，一维数组`y`表示房屋总价。线性模型的基本形式可以表示为：

\[ \hat{y} = w_0 + w_1 * x_i \]

其中，`w_0`是截距（偏置项），`w_1`是斜率（权重），`x_i`是每个样本的面积。

损失函数通常选择均方误差(Mean Squared Error, MSE)，对于一组预测值`yp`和真实值`y`，MSE计算公式为：

\[ Loss = \frac{1}{n} \sum_{i=1}^{n} (y_i - \hat{y}_i)^2 \]

梯度下降算法用于最小化损失函数。对于线性模型，梯度函数（损失关于权重的导数）为：

\[ \frac{\partial Loss}{\partial w_0} = \frac{-2}{n} \sum_{i=1}^{n} (y_i - (\hat{y}_i)) \]
\[ \frac{\partial Loss}{\partial w_1} = \frac{-2}{n} \sum_{i=1}^{n} (y_i - (\hat{y}_i)) * x_i \]

以下是使用Python和scikit-learn实现上述过程的示例代码：

import numpy as np
from sklearn.linear_model import LinearRegression
import matplotlib.pyplot as plt

# 假设我们有数据
X = # 房屋面积数据
y = # 房屋总价数据

# 初始化参数
n_samples, n_features = X.shape
w0 = 0
w1 = 0
learning_rate = 0.01
epochs = 1000
losses = []

# 创建模型
model = LinearRegression()

for epoch in range(epochs):
    # 前向传播
    y_pred = w0 + w1 * X
    loss = ((y - y_pred) ** 2).mean()
    
    # 梯度计算
    dw0 = -(2 / n_samples) * sum(y_pred - y)
    dw1 = -(2 / n_samples) * sum((y_pred - y) * X)

    # 更新参数
    w0 -= learning_rate * dw0
    w1 -= learning_rate * dw1
    
    losses.append(loss)

    # 显示损失随迭代次数的变化（这里仅用于可视化）
    if epoch % 100 == 0:
        print(f"Epoch {epoch+1}: Loss = {loss:.4f}")

# 可视化
plt.plot(range(epochs), losses, label='Loss Function')
plt.xlabel('Iterations')
plt.ylabel('Loss')
plt.title('Loss vs. Iterations')
plt.legend()
plt.show()

注意：在这个代码片段中，实际的数据`X`和`y`需要你自己提供，并且`n_samples`和`n_features`应该对应数据的实际维度。此外，这个例子仅展示了一次性遍历所有数据进行梯度下降的情况，实际上在训练大型数据集时，我们会使用批量梯度下降或随机梯度下降。