数学建模分配问题图和代码.doc

附录一 图

附录二 程序代码

代码（1）
clc,clear
u=1;
v=1.4;
w=1.3;
a=10000;
xx=u*a/((-log(0.7))^0.5)
yy=v*a/((-log(0.7))^0.5)
zz=w*a/((-log(0.7))^0.5)
x=linspace(0,60000);
y1=1-exp(-(1/xx*x).^2);
y2=1-exp(-(1/yy*x).^2);
y3=1-exp(-(1/zz*x).^2);
plot(x,y1,'b-',x,y2,'b--',x,y3,'r')
legend('y1企业单位','y2机关单位','y3事业单位 ')
title('企业、机关和事业单位满意度函数')
代码（2）
clc,clear
a=[0.520540699 0.636282765 0.629132233 0.693296223 0.740927091 0.84407568 0.980355395 1.103007657 1.280008353 1.458129 0.777175856 0.976488106 0.973183519 1.105332004 1.260764195 1.425830679 1.737382445 1.945204292 2.277792471 2.545222 0.745693558 0.95509605 1.054786 1.187656615 1.382634817 1.6285183 1.935217186 2.48730562 3.377363572 3.416021];
s=10000*a(1,:);
d=10000*a(2,:);
f=10000*a(3,:);
a=10000*a;
g=a';
bar(g,1)
legend('企业单位平均养老金','事业单位平均养老金','机关单位平均养老金')
title('企业、机关与事业单位平均养老金比较图')
colormap([1 0 0;0 1 0;0 0 1])
代码（3）
clc,clear
x=[1 2 3 4 5 6 7 8 9 10];
y=[0.520540699 0.636282765 0.629132233 0.693296223 0.740927091 0.84407568 0.980355395 1.103007657 1.280008353 1.458129];
y=10000*y;
z=[0.777175856 0.976488106 0.973183519 1.105332004 1.260764195 1.425830679 1.737382445 1.945204292 2.277792471 2.545222];
z=10000*z;
w=[0.745693558 0.95509605 1.054786 1.187656615 1.382634817 1.6285183 1.935217186 2.48730562 3.377363572 3.416021];
w=10000*w;
plot(x,y,'bo');
hold on
plot(x,z,'b*')
hold on
plot(x,w,'ro')
hold on
代码（4）
clc,clear
x=[1 2 3 4 5 6 7 8 9 10];
y=[ 0.520540699 0.636282765 0.629132233 0.693296223 0.740927091 0.84407568 0.980355395 1.103007657 1.280008353 1.458129];
y=10000*y;
z=[ 9550 10227 11547 13188 15169 17730 20973 24186 26909 30587];
plot(x,y,'ro');
hold on
plot(x,z,'b*')
hold on
代码（5）
x=[1 2 3 4 5 6 7 8 9 10];
y=[12028 16356 18772 20315 21603 22742 24489 25963 27881 29208];
z=[22439 23919 25522 27301 28445 29239 30605 31588 33834 34671];
w=[20836 22203 23684 25329 26384 27114 28375 29281 30895 32130];
plot(x,y,'ro')
hold on
plot(x,z,'b*')
hold on
plot(x,w,'ro')
hold on
title('修改模型')
代码（6）
clc,clear
x=[1 2 3 4 5 6 7 8 9 10];
y=[0 300 600 1300 2625 4000 4608 5000 5200 5592];
plot(x,y,'ro')
hold on
syms x
x=linspace(1,10);
y=5600./(1+100*exp(-x+0.5));
plot(x,y)
hold on
title('补偿养老金')
代码（7）
clc,clear
x=[1 2 3 4 5 6 7 8 9 10];
y=[0.520540699 0.636282765 0.629132233 0.693296223 0.740927091 0.84407568 0.980355395 1.103007657 1.280008353 1.458129];
y=10000*y;
z=[0.777175856 0.976488106 0.973183519 1.105332004 1.260764195 1.425830679 1.737382445 1.945204292 2.277792471 2.545222];
z=10000*z;
w=[0.745693558 0.95509605 1.054786 1.187656615 1.382634817 1.6285183 1.935217186 2.48730562 3.377363572 3.416021];
w=10000*w;
q=[5664.2 5912.1 6446.6 7408.7 9118.3 12101 16582 21553 26535 28665 ];
syms x
x=linspace(1,10);
y=98.92*x.^2-101.57*x+5636+15000./(1+800*exp(-x+0.5));
w=328.78*x.^2-547.53*x+8523.5;
plot(x,y,x,w)
hold on
title('养老金对比')