【神经网络】tensorflow实验7--回归问题

1. 实验目的

①掌握一元线性回归模型的实现方法
②掌握多元线性回归模型的实现方法
③掌握三维数据可视化方法

2. 实验内容

①使用TensorFlow建立一元线性回归模型，使用商品房销售数据训练模型，并使用训练好的模型预测房价
②使用TensorFlow建立多元线性回归模型，使用商品房销售数据预测房价，并实现三维数据可视化

3. 实验过程

题目一：

使用9.5小节中的“商品房销售记录表”作为样本数据，训练一元线性回归模型，根据商品房面积预测房价。
提示用户输入商品房面积，并进行输入校验。合理的输入如下：
面积：20-500之间的实数
如果输入正确，根据模型估计房价，并显示。
如果输入数据类型错误，或者输入数据范围不合理，根据错误类型提示，并等待用户重新输入，输错3次，则程序结束。
要求：
(1)编写代码，实现程序功能；
(2)记录实验过程和结果：尝试调试超参数，使模型达到最优的性能，记录实验过程和结果。

import  tensorflow as tf
import numpy as np
import matplotlib.pyplot as plt

#加载面积和房价
x = tf.constant([137.97,104.50,100.00,124.32,79.20,99.00,124.00,114.00,106.69,138.05,53.75,46.91,68.00,63.02,81.26,86.21])
y = tf.constant([145.00,110.00,93.00,116.00,65.32,104.00,118.00,91.00,62.00,133.00,51.00,45.00,78.50,69.65,75.69,95.30])


#求均值
meanX = tf.reduce_mean(x)
meanY = tf.reduce_mean(y)

#求权值
sumXY = tf.reduce_sum((x - meanX) * (y - meanY))
sumXX = tf.reduce_sum((x - meanX) * (x - meanX))

w = sumXY / sumXX

b = meanY - w * meanX

print("权值w为",w.numpy(),"\n偏置值b为",b.numpy())
print("线性模型：y=",w.numpy(),"*x + ",b.numpy())


for i in range(3):
    print("商品房面积和预测房价")
    x1 = input("请输入商品房面积")

    if(i == 3):
        print("错误过多，GameOver")

    else:

        if x1.isdigit():
            x1 = float(x1)
            if(x1 <= 500) & (x1 >= 20):
                y1= w * x1 + b
                print("面积%f的商品房价格为%f"%(x1,y1))
                break
            else:
                print("输入的面积大小错误，请重新输入")
        else:
            print("输入的面积类型错误，请重新输入")

题目二：

使用9.5小节中的“商品房销售记录表”作为样本数据，训练多元线性回归模型，实现一个房价预测系统。
要求：
(1)尝试调试超参数，使模型达到最优的性能，记录实验过程和结果；
(2)创建3D绘图对象来绘制空间点集。x轴表示房屋面积，y轴表示房间数，z轴表示样本的销售价格；
(3)提示用户输入商品房面积和房间数，并进行输入校验。如果输入正确，根据模型预测房价。
合理的输入如下：
面积：20-500之间的实数
房间数：1-10之间的整数
如果输入数据类型错误，或者输入数据范围不合理，根据错误类型给出提示，并等待用户重新输入，输错3次，则程序结束。
（输入时，请注意房间面积对应房间数的合理性）
提示：TensorFlow中矩阵求逆函数tf.linalg.inv()

import tensorflow as tf
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
x1 = np.array([137.97,104.50,100.00,124.32,79.20,99.00,
               124.00,114.00,106.69,138.05,53.75,46.91,68.00,63.02,81.26,86.21])
x2 = np.array([3.00,2.00,2.00,3.00,1.00,2.00,3.00,2.00,2.00,3.00,1.00,1.00,1.00,1.00,2.00,2.00])
y = np.array([145.00,110.00,93.00,116.00,65.32,104.00,
               118.00,91.00,62.00,133.00,51.00,45.00,78.50,69.65,75.69,95.30])
x0 = np.ones(len(x1))
ty_x0 = tf.constant(x0) #创建常数张量
ty_x1 = tf.constant(x1)
ty_x2 = tf.constant(x2)
ty_y = tf.constant(y)

ty_X = tf.stack((ty_x0,ty_x1,ty_x2), axis=1) #将全1数组与x1x2数组堆叠构造16行3列的属性矩阵
ty_y = tf.reshape(ty_y, [16,1]) #十六和一的一维数组

Xt = tf.transpose(ty_X) # X的转置
XtX_1 = tf.linalg.inv(tf.matmul(Xt,ty_X)) #对结果求逆
XtX_1_Xt = tf.matmul(XtX_1,Xt) #求逆后乘以x的转置
W = tf.matmul(XtX_1_Xt,ty_y) #继续乘以y

W = tf.reshape(W,[-1]) #转化为一维数组

# 绘制3D图形
fig = plt.figure()
ax3d = Axes3D(fig)
ax3d.scatter(x1,x2,y)
ax3d.set_xlabel('Area',color = 'r',fontsize = 14)
ax3d.set_ylabel('Room',color = 'r',fontsize = 14)
ax3d.set_zlabel('Price',color = 'r',fontsize = 14)

plt.show()
j = 0
while (True):
                print("请输入房屋面积和房间数，预测房屋销售价格：")
                x1_test = input("商品房面积(20-500)：")
                x2_test = input("房间数(1-10)：")
                if x1_test.isdigit() and x2_test.isdigit():
                                x1_test = float(x1_test)
                                x2_test = int(x2_test)
                                if (x1_test < 20 or x1_test > 500):
                                             print("对不起，您输入的房屋面积超出范围")
                                             j += 1
                                elif(x2_test < 1 or x2_test > 10):
                                             print("对不起你输入的房间数超出范围")
                                             j += 1
                                else:
                                             y_pred = W[1] * x1_test+ W[2] * x2_test + W[0]
                                             y_float = float(y_pred)
                                             print(y_float)
                                             print("预测价格：", round(y_float, 2), "万元")
                else:
                                print("对不起，您的输入无效")
                                j += 1
                if j>= 3:
                                print("对不起，您已经3次输入错误，程序退出")
                                exit()

4.实验小结

① 实验过程中遇到了哪些问题，你是如何解决的？
在实验中用到的方法还不熟悉，在网上进行查找，看ppt
② 分别使用Numpy和TensorFlow进行数据的加载和数组的堆叠
hstac沿水平方向堆叠数组（numpy array）
vstack沿垂直方向堆叠数组（numpy array）
预加载数据：用一个constant常量将数据集加载到计算图中（主要用于小数据集）
placehold feed_dict：从内存中读取数据，占位符填充数据
queue队列：基于队列的输入通道（在计算图计算前从队列中读取数据）
③ 题目一和题目二在代码实现过程中，有重复的代码段，你知道该如何精简自己的代码吗？请进行简要的描述。
重复的代码段可以直接代替
④ 在题目基本要求的基础上，你对每个题目做了那些扩展和提升？或者你觉得在编程实现过程中，还有哪些地方可以进行优化？
对于判断错误的分类更加详细，如先判断面积是否满足，再判断房间数。