基于BP神经网络的个人住房贷款还款信用评估:模型构建与实证分析_第1页
基于BP神经网络的个人住房贷款还款信用评估:模型构建与实证分析_第2页
基于BP神经网络的个人住房贷款还款信用评估:模型构建与实证分析_第3页
基于BP神经网络的个人住房贷款还款信用评估:模型构建与实证分析_第4页
基于BP神经网络的个人住房贷款还款信用评估:模型构建与实证分析_第5页
已阅读5页,还剩17页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

基于BP神经网络的个人住房贷款还款信用评估:模型构建与实证分析一、引言1.1研究背景与意义随着我国经济的快速发展和城市化进程的加速,个人住房贷款业务在金融市场中占据着日益重要的地位。中国央行发布的金融数据显示,2024年10月份个人住房贷款规模企稳,明显好于前9个月月均减少690亿元的水平,也好于前两年同期水平。个人住房贷款业务不仅满足了居民的住房需求,也推动了房地产市场的繁荣,对国民经济的增长起到了积极的促进作用。然而,随着贷款规模的不断扩大,银行面临的信用风险也日益凸显。信用风险是指借款人在贷款期限内无法按时足额偿还贷款本息,从而给银行带来损失的可能性。信用风险的存在不仅会影响银行的资产质量和盈利能力,还可能对整个金融体系的稳定造成威胁。准确的信用评估是银行控制个人住房贷款信用风险的关键。通过科学合理的信用评估,银行可以筛选出信用良好、还款能力强的借款人,降低违约风险,保障贷款资金的安全。信用评估还有助于银行合理确定信贷额度和利率,优化资源配置,提高经营效益。传统的信用评估方法往往依赖于专家经验和简单的统计模型,存在主观性强、准确性低等问题,难以满足现代金融市场的需求。BP神经网络作为一种强大的机器学习算法,具有自学习、自适应和非线性映射能力,能够处理复杂的数据关系,在信用评估领域展现出独特的优势。它可以通过对大量历史数据的学习,自动提取数据中的特征和规律,从而对借款人的信用状况进行准确评估。将BP神经网络应用于个人住房贷款还款信用评估,能够有效提高评估的准确性和效率,为银行的风险管理提供有力支持。1.2国内外研究现状在个人住房贷款还款信用评估领域,国内外学者进行了大量的研究。早期的研究主要采用传统的统计方法,如判别分析、线性回归和Logistic回归等。这些方法基于一定的假设前提,通过对历史数据的统计分析来构建信用评估模型。例如,Durand在1941年正式系统地提出使用数理统计模型辅助消费者授信决策,对贷款的信用风险进行评价,开启了个人信用评估从定性分析向定量分析的转变。然而,传统统计方法在处理复杂的非线性关系时存在局限性,难以准确捕捉数据中的潜在信息。随着信息技术的发展,机器学习和数据挖掘技术逐渐被引入信用评估领域。BP神经网络作为一种重要的机器学习算法,因其强大的非线性映射能力和自学习能力,受到了广泛的关注和应用。国外学者率先将BP神经网络应用于个人信用评估,并取得了较好的效果。研究表明,BP神经网络能够有效处理信用评估中的非线性问题,提高评估的准确性。在国内,相关研究也在不断深入。学者们结合我国个人住房贷款市场的特点,对BP神经网络在信用评估中的应用进行了大量的实证研究。例如,通过对某商业银行个人住房贷款数据的分析,验证了BP神经网络模型在信用风险评估中的有效性。尽管BP神经网络在个人住房贷款还款信用评估中取得了一定的成果,但仍存在一些问题需要解决。例如,BP神经网络容易陷入局部最优解,对初始权重和阈值敏感,训练时间较长等。为了克服这些问题,国内外学者提出了多种改进方法,如结合遗传算法、模拟退火算法等对BP神经网络进行优化,以提高模型的性能和泛化能力。1.3研究方法与创新点本文采用了多种研究方法,以确保研究的科学性和可靠性。通过广泛查阅国内外相关文献,了解个人住房贷款还款信用评估的研究现状和发展趋势,为研究提供理论基础和研究思路。选取了多家银行的个人住房贷款数据作为案例,对实际的信用评估情况进行深入分析,总结存在的问题和挑战。运用实证研究方法,构建基于BP神经网络的个人住房贷款还款信用评估模型,并利用实际数据进行训练和验证,通过对比分析评估模型的性能和效果。本文的创新点主要体现在以下两个方面。将遗传算法、模拟退火算法等多种优化算法与BP神经网络相结合,对BP神经网络的初始权重和阈值进行优化,以提高模型的收敛速度和准确性,克服BP神经网络容易陷入局部最优解的问题。综合考虑多个维度的因素,构建更加全面、合理的个人住房贷款还款信用评估指标体系,使评估结果更加准确地反映借款人的信用状况。二、BP神经网络基本原理2.1BP神经网络结构BP神经网络是一种多层前馈神经网络,其基本结构主要由输入层、隐藏层和输出层组成,各层之间通过带有权重的连接相互关联。输入层作为数据的入口,负责接收外部输入信号,这些信号可以是各种与个人住房贷款还款信用评估相关的特征数据,如借款人的收入水平、信用记录、负债情况等。输入层的神经元数量取决于输入数据的特征维度,每个神经元对应一个输入特征。隐藏层是BP神经网络的核心部分,它位于输入层和输出层之间,可以有一层或多层。隐藏层的主要作用是对输入信号进行非线性变换,通过复杂的内部计算,提取输入数据中的深层次特征和规律,从而学习输入与输出之间的复杂映射关系。隐藏层神经元数量和层数的选择是一个关键问题,它们会影响网络的学习能力和泛化性能。如果隐藏层神经元数量过少,网络可能无法学习到足够的特征,导致欠拟合;而如果神经元数量过多,可能会使网络学习到过多的噪声,出现过拟合现象。通常,隐藏层神经元数量可以根据经验公式进行初步估算,再通过实验进行调整优化。输出层是网络的最终结果输出部分,其神经元数量根据具体的任务需求而定。在个人住房贷款还款信用评估中,输出层的神经元可能只有一个,表示借款人的信用状况(如信用良好或信用不良),也可能有多个,表示不同的信用等级。输出层根据隐藏层传递过来的信息,经过计算生成最终的输出结果,这个结果就是网络对输入数据的预测或评估。在BP神经网络中,各层神经元之间的连接权重决定了信号传递的强度和方向。权重是神经网络的重要参数,通过训练不断调整权重,使得网络能够准确地学习到输入数据与输出结果之间的关系。初始权重通常是随机初始化的,这样可以避免网络在训练过程中陷入局部最优解。除了权重,每个神经元还包含一个偏置项,它类似于数学中的截距,用于调整神经元的激活阈值,使得神经元能够更好地对输入信号做出响应。偏置项和权重一样,在训练过程中也会不断更新,以优化网络的性能。2.2前向传播过程前向传播是BP神经网络计算输出的核心过程,在个人住房贷款还款信用评估模型中,输入层接收借款人的各项特征数据,如年龄、收入、负债比例、信用历史时长等信息。这些数据作为输入信号,被传递到隐藏层。在隐藏层中,每个神经元会接收来自输入层神经元的信号,并对这些信号进行加权求和。假设隐藏层中第i个神经元接收来自输入层第j个神经元的信号,其连接权重为w_{ij},输入层第j个神经元的输出为x_j,则隐藏层第i个神经元的加权输入z_i可以表示为:z_i=\sum_{j=1}^{n}w_{ij}x_j+b_i,其中b_i是隐藏层第i个神经元的偏置项,n是输入层神经元的数量。加权求和后的结果z_i并不会直接作为隐藏层神经元的输出,而是需要经过激活函数f的处理。激活函数的作用是为神经网络引入非线性因素,使得网络能够学习和模拟复杂的非线性关系。如果没有激活函数,无论神经网络有多少层,其整体功能都将等同于一个线性模型,无法处理复杂的现实问题。常见的激活函数有Sigmoid函数、Tanh函数和ReLU函数等。以Sigmoid函数为例,其数学表达式为f(x)=\frac{1}{1+e^{-x}},该函数将输入值映射到(0,1)之间,具有平滑、连续可导的特点。经过激活函数处理后,隐藏层第i个神经元的输出y_i为:y_i=f(z_i)。隐藏层的输出结果会继续作为下一层(如果有多个隐藏层,则为下一个隐藏层;如果是最后一层隐藏层,则为输出层)的输入信号,重复上述加权求和和激活函数处理的过程,直到信号传递到输出层。在输出层,同样进行加权求和和激活函数处理(根据具体问题,输出层的激活函数可能与隐藏层不同),最终得到网络的输出结果\hat{y}。这个输出结果就是BP神经网络对借款人还款信用状况的预测值,例如在二分类问题中,\hat{y}可能表示借款人违约的概率,取值范围在[0,1]之间,当\hat{y}大于某个阈值(如0.5)时,可判断借款人信用不良,存在违约风险;当\hat{y}小于该阈值时,可判断借款人信用良好,还款能力较强。前向传播过程是一个信号不断传递和处理的过程,通过各层神经元的协同工作,将输入数据逐步转化为有意义的输出结果,为后续的误差计算和反向传播提供基础。2.3反向传播与训练算法反向传播是BP神经网络训练过程中的关键步骤,其核心目的是根据网络输出与期望输出之间的误差,调整网络中的权重和偏置,使得误差逐步减小,从而提高网络的预测准确性。在个人住房贷款还款信用评估模型中,当完成前向传播得到网络的输出结果\hat{y}后,需要将其与实际的借款人信用状况(即期望输出y)进行比较,计算两者之间的误差。常用的误差函数是均方误差(MeanSquaredError,MSE),其计算公式为:E=\frac{1}{2}\sum_{k=1}^{m}(\hat{y}_k-y_k)^2,其中m是样本数量,\hat{y}_k是第k个样本的预测输出,y_k是第k个样本的实际输出。通过计算均方误差,可以衡量网络当前的预测性能,误差值越大,说明网络的预测结果与实际情况偏差越大,需要对网络进行调整。在计算出误差后,反向传播开始发挥作用。反向传播利用链式法则,从输出层开始,将误差逐层反向传播到输入层,计算每个权重和偏置对误差的贡献,即计算误差关于各层权重和偏置的梯度。以输出层到隐藏层的权重w_{ij}为例,其梯度\frac{\partialE}{\partialw_{ij}}的计算过程较为复杂,涉及到多个中间变量和链式求导。首先,根据误差函数对输出层神经元的输出求偏导,得到\frac{\partialE}{\partial\hat{y}_i};然后,根据输出层神经元的激活函数对其加权输入求偏导,得到\frac{\partial\hat{y}_i}{\partialz_i};接着,根据隐藏层神经元的输出与输出层神经元的加权输入之间的关系,得到\frac{\partialz_i}{\partialw_{ij}}。最后,通过链式法则\frac{\partialE}{\partialw_{ij}}=\frac{\partialE}{\partial\hat{y}_i}\frac{\partial\hat{y}_i}{\partialz_i}\frac{\partialz_i}{\partialw_{ij}},计算出权重w_{ij}的梯度。同样的方法,可以计算出隐藏层到输入层的权重梯度以及各层的偏置梯度。得到梯度后,使用梯度下降法等优化算法来更新权重和偏置。梯度下降法的基本思想是沿着梯度的反方向更新参数,以减小误差函数的值。对于权重w_{ij},其更新公式为:w_{ij}^{new}=w_{ij}^{old}-\eta\frac{\partialE}{\partialw_{ij}},其中\eta是学习率,它决定了每次权重更新的步长。学习率是一个重要的超参数,如果学习率过大,权重更新的步幅会过大,可能导致网络在训练过程中无法收敛,甚至出现振荡现象;如果学习率过小,权重更新的速度会非常缓慢,训练时间会大大增加。因此,选择合适的学习率对于网络的训练效果至关重要,通常需要通过实验进行调试。偏置项b_i的更新公式与权重类似,b_i^{new}=b_i^{old}-\eta\frac{\partialE}{\partialb_i}。在训练过程中,不断重复前向传播和反向传播的步骤,每次迭代都根据计算得到的梯度更新权重和偏置,使得误差函数逐渐减小。当误差达到预定的阈值或者达到最大迭代次数时,认为网络训练完成,此时得到的权重和偏置就是经过训练优化后的参数,网络能够对新的输入数据进行较为准确的信用评估预测。除了基本的梯度下降法,还有一些改进的训练算法,如带动量的梯度下降法、Adagrad、Adadelta、Adam等,这些算法在不同程度上改进了梯度下降法的性能,能够更快地收敛或者更好地处理不同类型的数据,在实际应用中可以根据具体情况选择合适的训练算法来提高BP神经网络的训练效率和性能。2.4BP神经网络关键要素激活函数在BP神经网络中起着至关重要的作用,它为神经网络引入了非线性因素,使网络能够学习和处理复杂的非线性关系。如果没有激活函数,神经网络将只是一个简单的线性模型,无法解决实际中的复杂问题。常见的激活函数有Sigmoid函数、Tanh函数和ReLU函数等。Sigmoid函数的数学表达式为f(x)=\frac{1}{1+e^{-x}},其输出范围在(0,1)之间,具有平滑、连续可导的特点,常用于二分类问题的输出层,将网络输出映射为概率值,表示样本属于某一类别的可能性。然而,Sigmoid函数存在梯度消失问题,当输入值的绝对值较大时,其导数接近于0,在反向传播过程中,梯度会随着层数的增加而迅速减小,导致网络难以训练深层结构。Tanh函数(双曲正切函数)的数学表达式为f(x)=\frac{e^x-e^{-x}}{e^x+e^{-x}},其输出范围在(-1,1)之间,与Sigmoid函数类似,也具有平滑可导的性质。Tanh函数的输出是以0为中心的,这使得它在一些情况下比Sigmoid函数表现更好,能够加快网络的收敛速度。但同样,Tanh函数也存在梯度消失问题,在深层神经网络中会影响训练效果。ReLU函数(修正线性单元)的数学表达式为f(x)=\max(0,x),即当x大于0时,输出x;当x小于等于0时,输出0。ReLU函数具有计算简单、收敛速度快等优点,能够有效缓解梯度消失问题,因此在隐藏层中得到了广泛应用。然而,ReLU函数也存在一些缺点,例如在x小于0时,神经元会处于完全不激活状态,可能导致部分神经元永远无法被激活,这种现象被称为“死亡ReLU”问题。损失函数用于衡量网络预测结果与真实结果之间的差异,它是网络训练过程中的优化目标。在个人住房贷款还款信用评估中,常用的损失函数有均方误差(MSE)和交叉熵损失函数。均方误差损失函数通过计算预测值与真实值之间差的平方的平均值来衡量误差,其公式为E=\frac{1}{n}\sum_{i=1}^{n}(\hat{y}_i-y_i)^2,其中n是样本数量,\hat{y}_i是第i个样本的预测值,y_i是第i个样本的真实值。均方误差损失函数计算简单,易于理解,适用于回归问题,如预测借款人的违约概率等连续值。交叉熵损失函数主要用于分类问题,它衡量的是两个概率分布之间的差异。在个人住房贷款还款信用评估中,如果将信用状况分为不同的类别(如良好、一般、不良等),则可以使用交叉熵损失函数。对于二分类问题,交叉熵损失函数的公式为E=-\sum_{i=1}^{n}[y_i\log(\hat{y}_i)+(1-y_i)\log(1-\hat{y}_i)],其中y_i是第i个样本的真实类别(0或1),\hat{y}_i是第i个样本预测为正类(1)的概率。交叉熵损失函数能够很好地反映分类问题中预测结果与真实结果的差异,在分类任务中比均方误差损失函数更有效。优化算法的作用是在训练过程中调整网络的权重和偏置,以最小化损失函数。除了前面提到的梯度下降法及其改进算法外,还有许多其他的优化算法可供选择。带动量的梯度下降法在更新权重时,不仅考虑当前的梯度,还考虑之前的权重更新方向,通过引入动量项来加速收敛并减少振荡。Adagrad算法根据每个参数的梯度历史信息自适应地调整学习率,对于频繁更新的参数,学习率会逐渐减小;对于稀疏的参数,学习率会相对较大,从而提高训练效率。Adadelta算法是对Adagrad算法的改进,它通过使用过去梯度平方的指数加权平均数来动态调整学习率,避免了Adagrad算法中学习率单调递减的问题。Adam算法结合了Adagrad和RMSProp算法的优点,它不仅能够自适应地调整学习率,还能处理梯度的稀疏性问题,在许多深度学习任务中表现出色。不同的优化算法具有不同的特点和适用场景,在实际应用中需要根据具体问题和数据特点选择合适的优化算法,以提高BP神经网络的训练效果和效率。三、个人住房贷款还款信用评估指标体系3.1评估指标选取原则全面性原则要求评估指标能够涵盖影响个人住房贷款还款信用的各个方面。个人住房贷款还款信用受到多种因素的综合影响,包括借款人的个人基本信息、财务状况、信用历史以及房产状况等。个人基本信息中的年龄、婚姻状况、教育程度等因素,可能反映出借款人的稳定性和未来收入潜力;财务状况中的收入水平、负债情况等,直接关系到借款人的还款能力;信用历史中的信用卡还款记录、其他贷款还款情况等,体现了借款人的还款意愿和信用意识;房产状况中的房产价值、房产位置等,也会对贷款风险产生影响。只有选取全面的评估指标,才能综合考量这些因素,避免遗漏重要信息,从而准确评估借款人的还款信用。代表性原则意味着所选取的指标应能准确代表其所属的影响因素类别,具有典型性和关键作用。在众多影响还款信用的因素中,有些因素的影响更为直接和显著,这些因素所对应的指标就应作为代表性指标被选取。收入水平是衡量借款人还款能力的关键指标,它直接反映了借款人按时偿还贷款本息的资金来源。较高的收入水平通常意味着较强的还款能力,而收入不稳定或过低则可能增加违约风险。信用历史中的逾期次数也是一个具有代表性的指标,逾期次数直接体现了借款人过去的还款行为和信用状况,逾期次数越多,表明借款人的还款意愿和信用风险越高。选取具有代表性的指标,可以更有效地反映借款人的还款信用状况,提高评估的准确性和效率。可操作性原则要求评估指标的数据易于获取、计算和分析。在实际应用中,评估指标的数据应能够通过合理的渠道收集得到,并且计算方法应简单明了,便于银行或金融机构进行操作。个人基本信息中的年龄、婚姻状况等数据,借款人在申请贷款时可以直接提供;财务状况中的收入水平,可以通过借款人的工资流水、纳税证明等材料获取;信用历史数据可以从征信系统中查询得到。如果指标的数据难以获取或计算复杂,将增加评估的成本和难度,甚至可能导致评估无法进行。确保评估指标的可操作性,有助于提高信用评估的效率和实用性,使其能够在实际业务中得到广泛应用。相关性原则强调评估指标与个人住房贷款还款信用之间应存在紧密的关联。所选取的指标应能够直接或间接地反映借款人的还款能力和还款意愿,与还款信用具有明确的因果关系或逻辑联系。收入水平与还款能力直接相关,收入越高,还款能力越强;信用历史中的逾期记录与还款意愿相关,有逾期记录说明借款人可能存在还款意愿不强的问题。如果选取的指标与还款信用无关或相关性较弱,不仅无法准确评估还款信用,还可能干扰评估结果,导致错误的决策。只有选取与还款信用相关性强的指标,才能构建出科学合理的评估指标体系,为个人住房贷款还款信用评估提供可靠的依据。3.2具体评估指标分析个人基本信息类指标包含多个方面,对还款信用有着重要影响。年龄与还款信用存在一定关联,一般来说,中青年群体收入相对稳定且具有增长潜力,在工作经验和职业发展上处于上升期,收入水平往往较高且稳定,具备较强的还款能力,违约风险相对较低;而年龄过小可能收入不稳定,年龄过大则可能面临收入减少和身体状况不佳等问题,还款能力和稳定性可能受到影响。婚姻状况也不容忽视,已婚借款人通常家庭经济相对稳定,夫妻双方可以共同承担家庭经济责任,在面对贷款还款时,家庭的经济缓冲能力较强,还款稳定性相对较高;未婚借款人可能在经济支持和稳定性方面相对较弱。教育程度同样是重要因素,较高教育程度的借款人往往更容易获得高收入、稳定的工作,其职业发展前景较好,收入水平相对较高且稳定,还款能力更有保障;教育程度较低的借款人可能在就业机会和收入水平上受到一定限制,还款能力相对较弱。财务状况类指标是评估还款信用的关键。收入水平是核心指标之一,它直接决定了借款人的还款能力。稳定且较高的收入意味着借款人有足够的资金按时偿还贷款本息,收入越高,还款能力越强,违约风险越低;相反,收入不稳定或过低,借款人可能难以按时足额还款,增加违约风险。负债比例反映了借款人的债务负担情况,计算公式为总负债除以总资产。负债比例过高,说明借款人的债务负担沉重,可用于偿还住房贷款的资金相对减少,还款能力受到削弱,违约风险相应增加;而负债比例较低,表明借款人的债务压力较小,有更多的资金用于偿还住房贷款,还款能力相对较强。资产状况也不容忽视,丰富的资产如房产、车辆、存款等,一方面可以作为还款的保障,在借款人遇到短期资金困难时,可通过变现资产来偿还贷款;另一方面,资产状况也反映了借款人的经济实力,资产雄厚的借款人还款能力和稳定性相对较高。信用历史类指标体现了借款人过去的信用行为。信用卡还款记录是重要参考,按时足额还款的借款人,表明其具有良好的信用意识和还款习惯,还款意愿较强,在申请个人住房贷款时,违约风险相对较低;而存在逾期还款记录的借款人,说明其还款意愿可能存在问题,未来违约的可能性相对较大。其他贷款还款情况同样重要,如车贷、消费贷等,如果借款人在这些贷款中都能按时还款,说明其信用状况良好,还款能力和意愿得到了一定的验证;反之,若有逾期或违约记录,则会增加个人住房贷款的风险。信用历史长度也对信用评估有影响,较长的信用历史能够提供更多的信用行为数据,有助于更全面、准确地评估借款人的信用状况。如果借款人长期保持良好的信用记录,说明其信用稳定性较高,还款信用更可靠;而信用历史较短的借款人,由于缺乏足够的信用数据,评估的准确性可能会受到一定影响。房产状况类指标对还款信用评估有重要作用。房产价值是关键因素,房产价值越高,在借款人违约时,银行通过处置房产收回贷款的可能性越大,贷款风险相对降低;同时,房产价值也反映了借款人的资产实力,高价值房产的借款人通常经济状况较好,还款能力相对较强。房产位置也不容忽视,位于经济发达地区、交通便利、配套设施完善地段的房产,其市场需求大,价值相对稳定且具有增值潜力,这样的房产作为抵押物,银行面临的风险相对较小;而位于偏远地区、经济发展滞后地段的房产,市场需求小,价值波动较大,在借款人违约时,银行处置房产可能面临困难,贷款风险相对增加。房产类型也会影响还款信用评估,如普通住宅、别墅、商业用房等,不同类型的房产在市场流通性、价值稳定性等方面存在差异,从而对贷款风险产生不同的影响。普通住宅的市场流通性较好,需求相对稳定,作为抵押物的风险相对较低;而商业用房等可能受到市场环境、经营状况等因素的影响较大,风险相对较高。3.3指标数据收集与预处理指标数据的收集渠道丰富多样。银行内部数据是重要来源之一,银行在日常业务开展过程中,积累了大量客户的基本信息、财务信息以及信用信息等。在个人住房贷款申请过程中,借款人需要向银行提供个人身份证明、收入证明、资产证明等资料,这些信息构成了银行内部数据的一部分。银行还会记录借款人在本行的账户交易流水、信用卡使用情况等信息,这些数据对于评估借款人的还款能力和信用状况具有重要价值。征信机构数据也是不可或缺的,我国有专业的征信机构,如中国人民银行征信中心,它们收集了各类金融机构上报的个人信用信息,涵盖了个人的贷款记录、信用卡还款记录、逾期记录等。通过查询征信机构的数据,能够全面了解借款人在整个金融市场的信用状况,为个人住房贷款还款信用评估提供有力支持。公开信息渠道也能获取部分有用数据,例如政府部门的相关网站可能提供个人的社保缴纳记录、纳税记录等信息。社保缴纳记录可以反映借款人的就业稳定性和收入水平,纳税记录则能进一步验证借款人的收入真实性。一些公共事业部门的缴费记录,如水电费、燃气费等,也能在一定程度上反映借款人的生活稳定性和信用状况。互联网大数据平台也逐渐成为数据收集的新渠道,一些互联网金融公司、电商平台等积累了大量用户的消费行为数据、交易数据等。通过分析这些数据,可以挖掘出与借款人还款能力和信用相关的信息,如消费习惯、消费能力等。收集到的数据往往存在各种问题,需要进行预处理以提高数据的可用性。数据清洗是关键步骤,主要是去除重复数据、纠正错误数据以及填补缺失数据。在实际数据收集过程中,可能会出现重复记录的情况,这不仅会占用存储空间,还会影响数据分析的准确性,通过数据清洗可以识别并删除这些重复数据。数据中可能存在错误信息,如数据录入错误、格式错误等,需要对这些错误进行纠正,以确保数据的准确性。数据缺失也是常见问题,对于缺失的数据,可以采用多种方法进行填补,如均值填补法、中位数填补法、回归填补法等。均值填补法是用该变量的均值来填补缺失值;中位数填补法是用中位数来填补;回归填补法是通过建立回归模型,利用其他相关变量来预测缺失值并进行填补。数据归一化也是重要的预处理方法,由于不同指标的数据量纲和取值范围可能不同,这会对模型的训练和评估产生影响。收入水平可能以万元为单位,取值范围较大;而信用历史长度可能以年为单位,取值范围相对较小。如果直接将这些数据输入模型,可能会导致模型更关注取值范围大的指标,而忽视取值范围小的指标。通过数据归一化,可以将不同指标的数据转换到同一尺度上,消除量纲和取值范围的影响。常见的数据归一化方法有最小-最大标准化和Z-score标准化。最小-最大标准化将数据映射到[0,1]区间,公式为x_{new}=\frac{x-x_{min}}{x_{max}-x_{min}},其中x是原始数据,x_{min}和x_{max}分别是该指标数据的最小值和最大值,x_{new}是归一化后的数据。Z-score标准化则是将数据转换为均值为0,标准差为1的标准正态分布,公式为x_{new}=\frac{x-\mu}{\sigma},其中\mu是均值,\sigma是标准差。通过数据归一化处理,可以提高数据的可比性和模型的性能,使基于这些数据训练的BP神经网络模型能够更准确地进行个人住房贷款还款信用评估。四、基于BP神经网络的信用评估模型构建4.1模型设计思路本研究构建的基于BP神经网络的个人住房贷款还款信用评估模型,旨在利用BP神经网络强大的非线性映射和自学习能力,对个人住房贷款申请人的还款信用进行准确评估。模型以第三章所确定的个人住房贷款还款信用评估指标体系为基础,将这些评估指标作为BP神经网络的输入。具体来说,输入层接收借款人的个人基本信息(如年龄、婚姻状况、教育程度等)、财务状况(收入水平、负债比例、资产状况等)、信用历史(信用卡还款记录、其他贷款还款情况、信用历史长度等)以及房产状况(房产价值、房产位置、房产类型等)相关数据。这些数据经过预处理(数据清洗、归一化等)后,能够以统一的标准和格式输入到BP神经网络中,为后续的计算和分析提供准确的数据基础。在网络内部,通过隐藏层对输入数据进行复杂的非线性变换和特征提取。隐藏层神经元利用激活函数(如ReLU、Sigmoid等)对输入信号进行处理,挖掘不同指标之间的潜在关系和模式,从而学习到借款人各项特征与还款信用之间的复杂映射关系。输出层则根据隐藏层传递过来的信息,输出借款人的还款信用评估结果。评估结果可以采用多种形式表示,如信用评分(设定一个具体的分数范围,如0-100分,分数越高表示信用越好)、信用等级(分为优、良、中、差等不同等级)或违约概率(以一个概率值表示借款人违约的可能性,取值范围在0-1之间)。通过这种方式,BP神经网络能够综合考虑多个因素,对个人住房贷款还款信用进行全面、准确的评估,为银行等金融机构的贷款决策提供有力支持。4.2网络结构确定输入层神经元个数的确定直接取决于个人住房贷款还款信用评估指标的数量。在第三章构建的评估指标体系中,涵盖了个人基本信息、财务状况、信用历史和房产状况四个方面的多个指标,共计[X]个指标。因此,输入层神经元个数设定为[X],以确保能够完整地接收和处理这些评估指标数据。每个输入层神经元对应一个评估指标,将相应的指标数据作为输入信号传递给隐藏层。输出层神经元个数依据评估结果的表示方式而定。如果采用信用评分的方式,输出层可以设置1个神经元,该神经元输出的数值即为信用评分;若采用信用等级评估,假设将信用等级分为[Y]个等级(如优、良、中、差四个等级,此时Y=4),则输出层神经元个数为[Y],每个神经元代表一个信用等级,通过神经元的输出值(如经过Softmax函数处理后得到的概率值)来确定借款人所属的信用等级;当以违约概率表示评估结果时,输出层设置1个神经元,输出值在0-1之间,表示借款人违约的概率。隐藏层神经元个数的确定是一个关键且复杂的问题,目前尚无严格的理论方法来精确确定。本研究采用经验法和试凑法相结合的方式。首先,参考经验公式l=\sqrt{m+n}+a(其中l为隐藏层神经元个数,m为输入层神经元个数,n为输出层神经元个数,a为1-10之间的常数),初步估算隐藏层神经元个数的范围。然后,在这个范围内进行试凑,通过多次实验,观察不同隐藏层神经元个数下BP神经网络模型的性能表现,如训练误差、测试误差、收敛速度等指标。在实验过程中,保持其他参数不变,仅改变隐藏层神经元个数,对模型进行训练和测试。经过一系列的实验和分析,最终确定一个能够使模型性能达到最优的隐藏层神经元个数,例如确定为[Z]个神经元。这样确定的隐藏层神经元个数能够在保证模型学习能力的同时,避免过拟合或欠拟合问题,使模型能够准确地学习到输入数据与输出结果之间的复杂关系,从而对个人住房贷款还款信用进行有效的评估。4.3训练算法选择与比较在BP神经网络的训练过程中,训练算法的选择对模型的性能和训练效率有着至关重要的影响。本研究选取了几种常见且具有代表性的训练算法进行对比分析,包括附加动量法、自适应学习速率、拟牛顿法、共轭梯度法和LM算法(Levenberg-Marquardt算法)。附加动量法在传统BP算法的基础上,引入了动量项。其核心思想是在权重更新时,不仅考虑当前误差的梯度,还考虑上一次权重更新的方向和大小,通过动量因子将上一次权重更新的影响传递到本次更新中。公式为\Deltaw(t)=\eta\frac{\partialE}{\partialw}+\alpha\Deltaw(t-1),其中\Deltaw(t)是当前时刻t的权重更新量,\eta是学习率,\frac{\partialE}{\partialw}是误差对权重的梯度,\alpha是动量因子(通常取值在0-1之间),\Deltaw(t-1)是上一时刻的权重更新量。当网络在误差曲面较为平坦的区域时,梯度值较小,普通BP算法的权重更新可能会停滞,而附加动量法由于动量项的存在,可以继承之前的更新趋势,使权重更新不会停止,从而加快收敛速度。同时,当算法陷入局部极小值附近时,动量项可以使权重更新具有一定的惯性,有可能跳过局部极小值区域,继续向全局最小值方向搜索。自适应学习速率算法则是为了解决传统BP算法中固定学习率带来的问题。在传统BP算法中,学习率是一个固定值,在训练开始初期,合适的学习率可能在训练后期变得过大或过小,导致训练不稳定或收敛速度慢。自适应学习速率算法根据训练过程中的误差变化情况自动调整学习率。通常的准则是检查权值更新是否真正降低了误差函数,如果确实如此,则说明当前学习率较小,可以适当增加;若权值更新导致误差增大,即产生了过调现象,那么就应该减小学习率的值。例如,常见的自适应学习速率调整公式为:若E(k)-E(k-1)<0,则\eta(k)=\eta(k-1)\times1.05;若E(k)-E(k-1)>0,则\eta(k)=\eta(k-1)\times0.7,其中E(k)为第k步的误差平方和,\eta(k)为第k步的学习率。这种动态调整学习率的方式能够提高训练效率,使模型更快地收敛到较优解。拟牛顿法是一类优化算法,它通过近似海森矩阵(HessianMatrix)来加速收敛。海森矩阵是误差函数对权重的二阶导数矩阵,计算海森矩阵的逆矩阵在理论上可以更快地找到误差函数的最小值,但直接计算海森矩阵及其逆矩阵的计算量非常大。拟牛顿法通过一些近似方法来避免直接计算海森矩阵,而是利用迭代过程中梯度的信息来近似海森矩阵的逆矩阵。常见的拟牛顿法有DFP(Davidon-Fletcher-Powell)算法和BFGS(Broyden-Fletcher-Goldfarb-Shanno)算法等。拟牛顿法在处理一些复杂的非线性问题时,能够在较少的迭代次数内达到较好的收敛效果,但其计算过程相对复杂,对内存的需求也较大。共轭梯度法也是一种用于求解无约束优化问题的迭代算法,在BP神经网络训练中常用于加速收敛。它的基本思想是通过构造一组共轭方向,使得在这些方向上依次进行搜索,能够更快地找到误差函数的最小值。共轭方向是指对于一个正定矩阵(在BP神经网络训练中,与海森矩阵相关),两个向量的内积在经过该矩阵变换后为零。共轭梯度法在每次迭代时,不仅利用当前点的梯度信息,还结合了之前搜索方向的信息,从而在一定程度上避免了传统梯度下降法中常见的锯齿现象,提高了收敛速度。与拟牛顿法相比,共轭梯度法不需要存储和计算近似海森矩阵,内存需求较小,计算量相对较低,适用于大规模数据集的训练。LM算法是一种专门为解决非线性最小二乘问题而设计的优化算法,在BP神经网络训练中表现出良好的性能。它结合了梯度下降法和高斯-牛顿法的优点,在误差曲面较平缓时,采用高斯-牛顿法,能够快速收敛;在误差曲面较陡峭时,采用梯度下降法,保证算法的稳定性。LM算法通过引入一个阻尼因子\mu来平衡这两种方法,当\mu较小时,算法更接近高斯-牛顿法;当\mu较大时,算法更接近梯度下降法。在训练过程中,根据误差的变化情况动态调整\mu的值。LM算法在收敛速度和精度方面通常优于其他几种算法,尤其适用于误差函数较为复杂的情况,但它对内存的需求较大,计算过程也相对复杂。为了选择最适合个人住房贷款还款信用评估模型的训练算法,进行了详细的实验对比。实验使用相同的训练数据和测试数据,对采用不同训练算法的BP神经网络模型进行训练和测试。在训练过程中,记录每个模型的训练误差、训练次数、训练时间以及收敛速度等指标。通过对这些指标的综合分析,发现LM算法在本次实验中表现最优。LM算法的训练误差最小,能够使模型更准确地拟合训练数据;训练次数相对较少,减少了训练的时间成本;收敛速度最快,能够在较短的时间内达到较好的训练效果。因此,本研究最终选择LM算法作为基于BP神经网络的个人住房贷款还款信用评估模型的训练算法,以确保模型能够高效、准确地进行训练和评估。4.4模型训练与优化在确定了BP神经网络的结构和训练算法后,使用准备好的训练数据对模型进行训练。训练数据包含了大量借款人的个人住房贷款相关信息,这些信息经过预处理后,被输入到BP神经网络中。在训练过程中,数据首先通过输入层进入网络,然后依次经过隐藏层和输出层进行前向传播计算,得到网络的预测输出。将预测输出与实际的还款信用情况(即期望输出)进行比较,计算两者之间的误差。这里采用均方误差(MSE)作为误差函数,公式为E=\frac{1}{2}\sum_{k=1}^{m}(\hat{y}_k-y_k)^2,其中m是样本数量,\hat{y}_k是第k个样本的预测输出,y_k是第k个样本的实际输出。根据计算得到的误差,利用LM算法进行反向传播,计算误差对网络中各权重和偏置的梯度,然后根据梯度来更新权重和偏置,使得误差逐渐减小。在训练过程中,不断重复前向传播和反向传播的步骤,每次迭代都根据新的梯度更新权重和偏置,直到误差达到预定的阈值或者达到最大迭代次数时,认为模型训练完成。在实际训练过程中,为了提高模型的性能和准确性,进行了一系列的优化操作。首先,对模型的参数进行调整。例如,尝试不同的学习率,学习率决定了权重更新的步长,如果学习率过大,权重更新的步幅会过大,可能导致模型在训练过程中无法收敛,甚至出现振荡现象;如果学习率过小,权重更新的速度会非常缓慢,训练时间会大大增加。通过多次实验,找到一个合适的学习率,使得模型在保证收敛速度的同时,能够达到较好的精度。还对动量因子(如果采用附加动量法)、阻尼因子(对于LM算法)等其他参数进行了调整和优化,以找到这些参数的最佳取值组合,使模型性能达到最优。增加训练次数也是一种优化方法。在一定范围内,增加训练次数可以使模型更好地学习到数据中的规律和特征,提高模型的准确性。但训练次数过多也可能导致过拟合现象,即模型对训练数据过度学习,在测试数据上表现不佳。因此,需要在训练过程中密切关注模型在训练集和测试集上的性能表现,当发现模型在测试集上的性能开始下降时,及时停止训练,避免过拟合。还可以采用一些正则化方法来防止过拟合。例如,L1和L2正则化,L1正则化是在损失函数中添加权重的绝对值之和作为惩罚项,公式为E_{L1}=E+\lambda\sum_{i}|w_i|,其中E是原始的损失函数,\lambda是正则化参数,w_i是权重;L2正则化是添加权重的平方和作为惩罚项,公式为E_{L2}=E+\lambda\sum_{i}w_i^2。通过添加正则化项,能够使模型的权重更加稀疏,减少模型对某些特征的过度依赖,从而提高模型的泛化能力。经过一系列的训练和优化操作,得到了性能良好的基于BP神经网络的个人住房贷款还款信用评估模型。该模型能够准确地对新的借款人的还款信用进行评估,为银行等金融机构的贷款决策提供可靠的依据,有效降低个人住房贷款的信用风险。五、实证分析5.1数据来源与样本选取本研究的数据来源于某大型商业银行过去[X]年的个人住房贷款客户信息。该银行在个人住房贷款领域具有丰富的业务经验和广泛的客户基础,其数据涵盖了不同地区、不同收入水平和不同信用状况的借款人信息,具有较高的代表性和可靠性。在数据收集过程中,银行提供了包括借款人的个人基本信息(如年龄、婚姻状况、教育程度等)、财务状况(收入水平、负债比例、资产状况等)、信用历史(信用卡还款记录、其他贷款还款情况、信用历史长度等)以及房产状况(房产价值、房产位置、房产类型等)在内的详细数据。这些数据为构建个人住房贷款还款信用评估模型提供了充足的信息。为了确保数据的质量和有效性,对原始数据进行了严格的筛选和预处理。首先,去除了数据中存在明显错误或缺失值过多的样本,以避免这些异常数据对模型训练产生负面影响。对于存在少量缺失值的样本,采用了均值填充、回归预测等方法进行填补,以保证数据的完整性。对数据进行了归一化处理,将不同指标的数据统一到相同的尺度,消除量纲差异对模型训练的影响。最终,从原始数据中选取了[M]条样本数据作为研究样本。为了充分利用数据进行模型训练和评估,将这些样本数据按照[70:30]的比例划分为训练集和测试集。其中,训练集包含[M1]条样本数据,用于训练BP神经网络模型,使其学习到借款人各项特征与还款信用之间的关系;测试集包含[M2]条样本数据,用于评估模型的性能和泛化能力,检验模型在新数据上的预测准确性。通过合理的数据划分,能够有效避免过拟合现象,确保模型在实际应用中具有良好的表现。5.2模型训练与结果分析使用训练集数据对基于BP神经网络的个人住房贷款还款信用评估模型进行训练。在训练过程中,采用了LM(Levenberg-Marquardt)算法,该算法结合了梯度下降法和高斯-牛顿法的优点,在误差曲面较平缓时,采用高斯-牛顿法,能够快速收敛;在误差曲面较陡峭时,采用梯度下降法,保证算法的稳定性。通过引入阻尼因子\mu来平衡这两种方法,当\mu较小时,算法更接近高斯-牛顿法;当\mu较大时,算法更接近梯度下降法。在训练过程中,根据误差的变化情况动态调整\mu的值,以提高训练效率和模型的准确性。经过多次实验,确定了模型的相关参数。学习率设置为[具体学习率值],该值在保证模型收敛的同时,能够加快训练速度。最大迭代次数设定为[具体最大迭代次数],以确保模型有足够的训练次数来学习数据中的规律。训练目标(均方误差)设定为[具体训练目标值],当模型的训练误差达到该目标值时,认为模型训练达到了较好的效果。在训练过程中,记录了模型的训练误差、训练次数、训练时间以及收敛速度等指标。随着训练的进行,模型的训练误差逐渐减小,表明模型在不断学习和优化,对训练数据的拟合能力逐渐增强。经过[实际训练次数]次迭代后,模型的训练误差达到了[实际训练误差值],满足了预先设定的训练目标,训练过程结束。整个训练过程耗时[具体训练时间]秒,收敛速度较快,表明LM算法在本模型训练中表现出了良好的性能。通过对训练结果的分析,发现模型在训练集上表现出了较高的准确性。这意味着模型能够较好地学习到借款人各项特征与还款信用之间的复杂关系,对训练集中的样本数据能够进行准确的信用评估。但需要注意的是,训练集上的良好表现并不一定意味着模型在未知数据(测试集)上也能有同样出色的表现,因此需要进一步使用测试集对模型进行验证和评估。5.3模型验证与泛化能力测试使用划分好的测试集数据对训练好的BP神经网络模型进行验证和泛化能力测试。将测试集中的样本数据输入到模型中,模型输出对应的还款信用评估结果,然后将这些预测结果与测试集中样本的实际还款信用情况进行对比分析。通过计算预测值与实际值之间的误差,采用均方误差(MSE)和平均绝对误差(MAE)作为评估指标。均方误差的计算公式为MSE=\frac{1}{n}\sum_{i=1}^{n}(\hat{y}_i-y_i)^2,其中n是测试集样本数量,\hat{y}_i是第i个样本的预测值,y_i是第i个样本的实际值;平均绝对误差的计算公式为MAE=\frac{1}{n}\sum_{i=1}^{n}|\hat{y}_i-y_i|。经过计算,模型在测试集上的均方误差为[具体MSE值],平均绝对误差为[具体MAE值]。这两个误差指标相对较小,说明模型的预测值与实际值较为接近,模型在测试集上具有较高的准确性。为了更直观地展示模型的预测效果,绘制了预测值与实际值的对比散点图。在散点图中,横坐标表示实际的还款信用情况,纵坐标表示模型的预测值。如果模型的预测完全准确,所有的散点应该分布在直线y=x上。从实际绘制的散点图来看,大部分散点都集中在直线y=x附近,进一步验证了模型的准确性。通过对测试集的验证,评估模型的泛化能力。泛化能力是指模型对未见过的数据的适应和预测能力。一个具有良好泛化能力的模型,在测试集上应该能够表现出与训练集相似的性能。从测试结果来看,模型在测试集上的误差虽然比训练集上略有增加,但增加幅度较小,表明模型具有较好的泛化能力,能够对新的个人住房贷款借款人的还款信用进行较为准确的评估,具有一定的实际应用价值。5.4与传统评估方法对比将基于BP神经网络的个人住房贷款还款信用评估模型与传统的信用评估方法进行对比,以进一步验证BP神经网络模型的优势。选择了线性判别分析(LDA)和Logistic回归这两种常见的传统信用评估方法作为对比对象。线性判别分析是一种经典的统计分类方法,它通过寻找一个线性变换,将高维数据投影到低维空间,使得同一类别的数据在投影空间中尽可能聚集,不同类别的数据尽可能分开,从而实现分类。在个人住房贷款还款信用评估中,LDA根据借款人的各项特征指标,构建线性判别函数,通过计算判别函数的值来判断借款人的信用状况。Logistic回归是一种广义的线性回归分析模型,常用于解决二分类问题。在个人住房贷款还款信用评估中,Logistic回归假设借款人违约的概率与各项特征指标之间存在一种非线性的关系,通过最大似然估计法来估计模型的参数,从而得到借款人违约概率的预测值。根据设定的阈值,将预测概率与阈值进行比较,判断借款人的信用状况。使用相同的训练集和测试集数据,分别对BP神经网络模型、LDA模型和Logistic回归模型进行训练和测试。在测试过程中,记录并对比三个模型的准确率、召回率、F1值等性能指标。准确率是指分类正确的样本数占总样本数的比例,计算公式为Accuracy=\frac{TP+TN}{TP+TN+FP+FN},其中TP表示真正例,即实际为正类且被正确预测为正类的样本数;TN表示真反例,即实际为反类且被正确预测为反类的样本数;FP表示假正例,即实际为反类但被错误预测为正类的样本数;FN表示假反例,即实际为正类但被错误预测为反类的样本数。召回率是指正确分类为正类的样本数占实际为正类样本数的比例,计算公式为Recall=\frac{TP}{TP+FN}。F1值是准确率和召回率的调和平均值,用于综合衡量模型的性能,计算公式为F1=\frac{2\timesPrecision\timesRecall}{Precision+Recall}。经过实验计算,BP神经网络模型的准确率为[BP神经网络模型准确率值],召回率为[BP神经网络模型召回率值],F1值为[BP神经网络模型F1值];LDA模型的准确率为[LDA模型准确率值],召回率为[LDA模型召回率值],F1值为[LDA模型F1值];Logistic回归模型的准确率为[Logistic回归模型准确率值],召回率为[Logistic回归模型召回率值],F1值为[Logistic回归模型F1值]。从这些指标可以看出,BP神经网络模型在准确率、召回率和F1值等方面均优于LDA模型和Logistic回归模型。这表明BP神经网络模型能够更准确地识别借款人的信用状况,在个人住房贷款还款信用评估中具有更高的可靠性和有效性,能够为银行等金融机构提供更有价值的决策支持。六、案例分析6.1案例背景介绍本案例选取了一位在某商业银行申请个人住房贷款的借款人。借款人李某,32岁,已婚,拥有本科学历,在一家国有企业担任中层管理人员。家庭月收入为20000元,家庭月支出约8000元,包括日常生活开销、子女教育费用等。李某名下有一辆价值15万元的汽车,银行存款50万元,无其他负债。其信用历史良好,信用卡使用5年以来,从未出现过逾期还款记录,也没有其他贷款记录,信用历史长度为5年。李某申请个人住房贷款用于购买一套位于市中心的三居室住宅,房产价值200万元。该房产地理位置优越,周边配套设施完善,交通便利。李某首付60万元,申请贷款金额为140万元,贷款期限为30年,采用等额本息还款方式,年利率为5%。银行在受理李某的贷款申请时,采用传统的信用评估方法,主要依据借款人的收入水平、负债情况、信用记录等因素进行评估。根据银行内部的评估标准,李某的收入稳定,负债比例低,信用记录良好,综合评估后,银行认为李某的还款能力较强,信用风险较低,批准了李某的贷款申请。6.2BP神经网络评估过程与结果将李某的个人住房贷款相关信息作为输入数据,按照之前构建的BP神经网络信用评估模型进行评估。首先,对输入数据进行预处理,包括数据清洗和归一化处理,确保数据的准确性和一致性,并将其转换到适合神经网络输入的尺度。经过数据预处理后,将数据输入到基于LM算法训练好的BP神经网络模型中。模型按照前向传播的过程,依次经过输入层、隐藏层和输出层的计算。在隐藏层中,神经元利用激活函数对输入信号进行非线性变换和特征提取,挖掘不同指标之间的潜在关系和模式。最终,输出层输出评估结果,这里以违约概率的形式表示,输出值为0.05。与银行传统评估方法相比,传统评估主要依赖于简单的指标判断和经验分析,虽然能从直观上判断借款人的还款能力和信用状况,但缺乏对多因素复杂关系的深入挖掘。而BP神经网络评估模型通过对大量历史数据的学习,能够更全面、准确地捕捉借款人各

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论