2020年数学建模C题代码

d 源论文 点击%chulisj.m%第一步数据预处理y_jin_jinexlsread(附件1123家有信贷记录企业的相关数据.xlsx,2,E2:E210948);y_jin_shuiexlsread(附件1123家有信贷记录企业的相关数据.xlsx,2,F2:F210948);y_xiao_jinexlsread(附件1123家有信贷记录企业的相关数据.xlsx,3,E2:E162485);y_xiao_shuiexlsread(附件1123家有信贷记录企业的相关数据.xlsx,3,F2:F162485);%读取发票状态jin_zfxlsread(附件1123家有信贷记录企业的相关数据.xlsx,2,H2:H210948);xiao_zfxlsread(附件1123家有信贷记录企业的相关数据.xlsx,3,H2:H210948);y_jin_qyxlsread(附件1123家有信贷记录企业的相关数据.xlsx,2,A2:A210948);y_xiao_qyxlsread(附件1123家有信贷记录企业的相关数据.xlsx,3,A2:A162485);n_jin_jinexlsread(附件2302家无信贷记录企业的相关数据.xlsx,3,E2:E395176);n_jin_shuiexlsread(附件2302家无信贷记录企业的相关数据.xlsx,3,F2:F395176);n_xiao_jinexlsread(附件2302家无信贷记录企业的相关数据.xlsx,2,E2:E330836);n_xiao_shuiexlsread(附件2302家无信贷记录企业的相关数据.xlsx,2,F2:F330836);n_jin_zfxlsread(附件2302家无信贷记录企业的相关数据.xlsx,3,H2:H395176);n_xiao_zfxlsread(附件2302家无信贷记录企业的相关数据.xlsx,2,H2:H330836);n_jin_qyxlsread(附件2302家无信贷记录企业的相关数据.xlsx,3,A2:A395176);n_xiao_qyxlsread(附件2302家无信贷记录企业的相关数据.xlsx,2,A2:A330836);%有信贷——数据处理if jin_zf0 %发票状态为作废时进项金额、税额都作废y_jin_jine10;y_jin_shuie10;elsey_jin_jine1y_jin_jine;y_jin_shuie1y_jin_shuie;if y_jin_jine0 %负数发票时进项金额为0税额不变供以后抵扣y_jin_jine10;endendif xiao_zf0y_xiao_jine10;y_xiao_shuie10;elsey_xiao_jine1y_xiao_jine;y_xiao_shuie1y_xiao_shuie;if y_xiao_jine0 %负数发票时销项金额为0税额不变供以后抵扣y_xiao_jine10;endendxlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_jin_jine1,1, B2:B210948);xlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_jin_shuie1,1, C2:C210948);xlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_xiao_jine1,2, B2:B162485);xlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_xiao_shuie1,2, C2:C162485);%无信贷数据处理if n_jin_zf0n_jin_jine10;n_jin_shuie10;elsen_jin_jine1n_jin_jine;n_jin_shuie1n_jin_shuie;if n_jin_jine0n_jin_jine10;endendif n_xiao_zf0n_xiao_jine10;n_xiao_shuie10;elsen_xiao_jine1n_xiao_jine;n_xiao_shuie1n_xiao_shuie;if n_xiao_jine0n_xiao_jine10;endendxlswrite(处理302家无信贷记录企业的相关数据.xlsx,n_jin_jine1,1, B2:B395176);xlswrite(处理302家无信贷记录企业的相关数据.xlsx,n_jin_shuie1,1, C2:C395176);xlswrite(处理302家无信贷记录企业的相关数据.xlsx,n_xiao_jine1,2, B2:B330836);xlswrite(处理302家无信贷记录企业的相关数据.xlsx,n_xiao_shuie1,2, C2:C330836);%求各企业总收益%有信贷企业y_zongzeros(123,1);y_jin_zongzeros(123,1);y_xiao_zongzeros(123,1);for i1:123for j1:210947if y_jin_qy(j)iy_jin_zong(i)sum(y_jin_jine1(j)y_jin_shuie1(j));endendendfor i1:123for k1:162484if y_xiao_qy(k)iy_xiao_zong(i)sum(y_xiao_jine1(k)y_xiao_shuie1(k));endendendfor i1:123y_zong(i)y_xiao_zong(i)-y_jin_zong(i); %总收益销项总-进项总endxlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_zong,3, B2:B124);%无信贷企业n_zongzeros(302,1);n_jin_zongzeros(302,1);n_xiao_zongzeros(302,1);for i124:425for j1:395175if n_jin_qy(j)in_jin_zong(i-123)sum(n_jin_jine1(j)n_jin_shuie1(j));endendendfor i124:425for k1:330835if n_xiao_qy(k)in_xiao_zong(i-123)sum(n_xiao_jine1(k)n_xiao_shuie1(k));endendendfor i1:302n_zong(i)n_xiao_zong(i)-n_jin_zong(i);endxlswrite(处理302家无信贷记录企业的相关数据.xlsx,n_zong,3, B2:B303);%risksy.m%收益归一化处理求收益信任值y_shouyixlsread(处理123家有信贷记录企业的相关数据1.xlsx,3,B2:B124);n_shouyixlsread(处理302家无信贷记录企业的相关数据.xlsx,3,B2:B303);y_maxmax(y_shouyi);y_minmin(y_shouyi);n_maxmax(n_shouyi);n_minmin(n_shouyi);y_rszeros(123,1);n_rszeros(302,1);for i1:123y_rs(i)(y_shouyi(i)-y_min)/(y_max-y_min);endxlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_rs,3, C2:C124);for k1:302n_rs(k)(n_shouyi(k)-n_min)/(n_max-n_min);endxlswrite(处理302家无信贷记录企业的相关数据.xlsx,n_rs,3, C2:C303);xlswrite(处理302家无信贷记录企业的相关数据.xlsx,1-n_rs,3, D2:D303);%riskxy.m%信誉等级数据预测求信誉信任值[num,str]xlsread(附件1123家有信贷记录企业的相关数据.xlsx,1,C2:C124);y_xyzeros(123,1);%信誉等级ABCD对应的信誉数据分别为[0.75,1],[0.5,0.75],[0.25,0.5],[0,0.25]区间的随机数for i1:123if strcmp(str(i),A)y_xy(i)0.750.25*rand(1);endif strcmp(str(i),B)y_xy(i)0.50.25*rand(1);endif strcmp(str(i),C)y_xy(i)0.250.25*rand(1);endif strcmp(str(i),D)y_xy(i)00.25*rand(1);endendxlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_xy,3, D2:D124);%riskall.m%加权求总信任值及风险y_risksyxlsread(处理123家有信贷记录企业的相关数据1.xlsx,3,C2:C124);y_riskxyxlsread(处理123家有信贷记录企业的相关数据1.xlsx,3,D2:D124);y_riskzeros(123,1);for i1:123y_risk(i)0.6*y_risksy(i)0.4*y_riskxy(i);%收益信任值和信誉信任值的权重为64endxlswrite(处理123家有信贷记录企业的相关数据1.xlsx,y_risk,3, E2:E124);xlswrite(处理123家有信贷记录企业的相关数据1.xlsx,1-y_risk,3, F2:F124);%ls_nihe.m%利率与流失率的拟合函数x xlsread(附件3银行贷款年利率与客户流失率关系的统计数据.xlsx, sheet1, A3:A31);y1xlsread(附件3银行贷款年利率与客户流失率关系的统计数据.xlsx, sheet1, B3:B31);y2xlsread(附件3银行贷款年利率与客户流失率关系的统计数据.xlsx, sheet1, C3:C31);y3xlsread(附件3银行贷款年利率与客户流失率关系的统计数据.xlsx, sheet1, D3:D31);p1polyfit(x,y1,2);p2polyfit(x,y2,2);p3polyfit(x,y3,2);figure(1);draw(p1,x,y1);function [fitresult, gof] createFit1(x, y)%CREATEFIT1(X,Y)% Create a fit.%% Data for untitled fit 1 fit:% X Input : x% Y Output: y% Output:% fitresult : a fit object representing the fit.% gof : structure with goodness-of fit info.%% 另请参阅 FIT, CFIT, SFIT.% 由 MATLAB 于 17-Aug-2023 09:26:20 自动生成%% Fit: untitled fit 1.[xData, yData] prepareCurveData( x, y );% Set up fittype and options.ft fittype( exp2 );opts fitoptions( Method, NonlinearLeastSquares );opts.Display Off;opts.StartPoint [0.903090038453583 0.560082104254179 -2.28174716847795 -23.1065744517619];% Fit model to data.[fitresult, gof] fit( xData, yData, ft, opts );% Plot fit with data.figure;h plot( fitresult, xData, yData );set(h, LineWidth, 1);legend( h, 原始数据, 拟合曲线, Location, southeast,FontSize16);% Label axesset(gca,FontName,Times New Rome,FontSize,16);xlabel(贷款年利率,FontName,宋体,FontSize,18);ylabel(客户流失率,FontName,宋体,FontSize,18);set(gcf,Units,centimeters,Position,[5 5 20 12]);clcclose allclc;clearformat longdata1xlsread(C:\Users\35129\Desktop\暑假集训\2020C6\附件22.xlsx);data2xlsread(C:\Users\35129\Desktop\暑假集训\2020C6\附件1123家有信贷记录企业的相关数据.xlsx,进项发票信息);data3xlsread(C:\Users\35129\Desktop\暑假集训\2020C6\附件1123家有信贷记录企业的相关数据.xlsx,销项发票信息);resultlzeros(123,11);%记录的结果保存在resultl中for i1:123i1;mfind(data2(:,1)i);%索引企业的位置数据 标记redata2(m(1,1):max(m),:);%re是所有同一家企业的数据resortrows(re,[3,4,5]);%按年月日来排序voidfind(re(:,10)0);%标记 报废发票的位置vsize(void,1);%报废发票的个数validsize(m,1)-void;%有效的票的个数%删除完无用发票mmfind((re(:,10)1)(re(:,9)0));%标记出有用的发票reere(mm,:);%筛出报废的发票reesortrows(ree,[3,4,5]);%再次排序%smonth(1,:)ree(1,3:4);%记录起始年月%smonth(2,:)ree(size(ree,1),3:4);%最后一个% length(((smonth(2,1)-2017)*12smonth(2,2))-((smonth(1,1)-2017)*12smonth(1,2)))1;%计算月数updata2(m(1,1):max(m),6);%消费 单位总和snumel(up);%交易次数hsize(ree);unumel(unique(up));%求上游企业的数量 unique表示不相同的个数 numel统计个数resultl(i,1)i;%表示企业代号bzeros(12,1);qqfind(ree(:,3)2017);kkree(qq(1,1):max(qq),:);for j1:12aafind(kk(:,4)j);b(j,1)sum(ree(aa(1,1):max(aa),9));endendclcclose allclc;clearformat longdataxlsread(C:\Users\35129\Desktop\暑假集训\2020C6\附件1123家有信贷记录企业的相关数据.xlsx,进项发票信息);resultlzeros(123,11);%记录的结果保存在resultl中for i1:123i1;mfind(data(:,1)i);%索引企业的位置数据 标记redata(m(1,1):max(m),:);%re是所有同一家企业的数据resortrows(re,[3,4,5]);%按年月日来排序voidfind(re(:,10)0);%标记 报废发票的位置vsize(void,1);%报废发票的个数validsize(m,1)-void;%有效的票的个数%删除完无用发票mmfind((re(:,10)1)(re(:,9)0));%标记出有用的发票reere(mm,:);%筛出报废的发票reesortrows(ree,[3,4,5]);%再次排序%smonth(1,:)ree(1,3:4);%记录起始年月%smonth(2,:)ree(size(ree,1),3:4);%最后一个% length(((smonth(2,1)-2017)*12smonth(2,2))-((smonth(1,1)-2017)*12smonth(1,2)))1;%计算月数updata(m(1,1):max(m),6);%消费 单位总和snumel(up);%交易次数unumel(unique(up));%求上游企业的数量 unique表示不相同的个数 numel统计个数resultl(i,1)i;%表示企业代号nian1find(ree(:,3)2016);Aree(nian1,4);if numel(unique(A))~1count hist(A,unique(A));elsecountnumel(nian1);endnian2find(ree(:,3)2017);Aree(nian2,4);count1 hist(A,unique(A));count[count,count1];nian3find(ree(:,3)2018);Aree(nian3,4);count2 hist(A,unique(A));count[count,count2];nian4find(ree(:,3)2019);Aree(nian4,4);count3 hist(A,unique(A));count[count,count3];nian5find(ree(:,3)2020);Aree(nian5,4);if numel(nian5)~0count4 hist(A,unique(A));count[count,count4];ends10;for j1:numel(count)-1s1s1(count(j1)-count(j))/count(j);endresultl(i,3)12*s/ numel(count);%业务率s1s1/numel(count)*12;resultl(i,4)s1;%年均变化率resultl(i,7)numel(count);enddata1xlsread(C:\Users\35129\Desktop\暑假集训\2020C6\附件1123家有信贷记录企业的相关数据.xlsx,销项发票信息);for i1:123mfind(data1(:,1)i);%索引企业的位置数据 标记redata1(m(1,1):max(m),:);%re是所有同一家企业的数据resortrows(re,[3,4,5]);%按年月日来排序voidfind(re(:,10)0);%标记 报废发票的位置vsize(void,1);%报废发票的个数validsize(m,1)-void;%有效的票的个数%删除完无用发票mmfind((re(:,10)1)(re(:,9)0));%标记出有用的发票reere(mm,:);%筛出报废的发票reesortrows(ree,[3,4,5]);%再次排序%smonth(1,:)ree(1,3:4);%记录起始年月%smonth(2,:)ree(size(ree,1),3:4);%最后一个% length(((smonth(2,1)-2017)*12smonth(2,2))-((smonth(1,1)-2017)*12smonth(1,2)))1;%计算月数updata1(m(1,1):max(m),6);%消费 单位总和snumel(up);%交易次数unumel(unique(up));%求上游企业的数量 unique表示不相同的个数 numel统计个数resultl(i,1)i;%表示企业代号nian1find(ree(:,3)2016);Aree(nian1,4);if numel(unique(A))~1count hist(A,unique(A));elsecountnumel(nian1);endnian2find(ree(:,3)2017);Aree(nian2,4);count1 hist(A,unique(A));count[count,count1];nian3find(ree(:,3)2018);Aree(nian3,4);count2 hist(A,unique(A));count[count,count2];nian4find(ree(:,3)2019);Aree(nian4,4);count3 hist(A,unique(A));count[count,count3];nian5find(ree(:,3)2020);Aree(nian5,4);if numel(nian5)~0count4 hist(A,unique(A));count[count,count4];ends10;for j1:numel(count)-1s1s1(count(j1)-count(j))/count(j);endresultl(i,5)12*s/ numel(count);%业务率s1s1/numel(count)*12;resultl(i,6)s1;%年均变化率resultl(i,8)numel(count);endxlswrite(附件.xlsx,resultl,A1:K123);clcclose allclc;clearformat longdataxlsread(C:\Users\35129\Desktop\暑假集训\2020C6\附件22.xlsx);wwzeros(123,4);aafind(data(:,2)2);qqdata(aa(:,1),3);junzhi1mean(qq);zz1std(qq);createFit1(x,y1);createFit1(x,y2);createFit1(x,y3);