matlab画基尼系数,matlab 拟合洛伦兹曲线求基尼系数

匿名用户

1级

2018-11-23 回答

%x1=0:0.1:1;

x1=linspace(0,1,8);

x2=0:0.001:1;

t=[0,4.12222,5.57722,13.91623,18.05866,22.62142,14.21962,21.48463]*0.01;%收入数百分比

xlen=length(x1);

y=zeros(1,xlen);

for i=1:xlen

y(i)=sum(t(1:i));

end

% c=polyfit(x1,y,7);

c=fit(x1',y','smoothingspline');

%d=polyval(c,x2);

d=c(x2);

%plot([0,1],[0,1],x1,y,['-'])

plot(x1,x1,'b-.',x1,y,'*',x2,d,'-')

title('洛伦兹曲线')

xlabel('人口百分比累计'),ylabel('收入百分比累计')

axis equal

axis([0,1,0,1])

grid on

%计算基尼系数

area1=trapz(x2,d);

area2=trapz(x1,x1);

JN=area1/area2;

disp(['收入数百分比基尼系数=',num2str(JN)])

%x1=0:0.1:1;

x1=linspace(0,1,8);

x2=0:0.001:1;

t=[0,4.12222,5.57722,13.91623,18.05866,22.62142,14.21962,21.48463]*0.01;%累计收入百分比

for i=1:length(x1)

y(i)=sum(t(1:i));

end

% c=polyfit(x1,y,7);

c=fit(x1',y','smoothingspline');

%d=polyval(c,x2);

d=c(x2);

figure

plot(x1,x1,'b-.',x1,y,'*',x2,d,'-')

% plot([0,1],[0,1],x1,y,['-'])

title('洛伦兹曲线')

xlabel('人口百分比累计')

ylabel('收入百分比累计')

axis equal

axis([0,1,0,1])

grid on

%计算基尼系数

area1=trapz(x2,d);

area2=trapz(x1,x1);

JN=area1/area2;

disp(['累计收入百分比基尼系数=',num2str(JN)])

追问:

f4d42d3eaec138a6bc29be36323a7e43.png

这个怎么说?程序不对啊

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值