大神:
您好,小妹最近一直被这个数学模型困扰,因为是学管理的,对计算机代码不熟悉,写不出代码,如果哪位大神能根据我的模型给出个能实现输出的代码,我将立即兑现承诺,将500金币双手奉上!下面的word版本里是我的数学模型,下面给出一些相关的代码格式,好让大神明白如何实现利润lr和价格w这几个需要输出的结果:
l=10;
v=4;
hd=0.8;
hm=0.8;
n=1;
for i=1:10
j=i;
w1=0.5*a-0.25*b*l*(2*i-1)+0.5*v-0.25*hd*l+0.5*l*hm*(i-(j-1)*n);
w2=a/3-(b*l*(i-1))/3+2*v/3+2*l*hm*(i-(j-1)*n)/3;
w3=a-i*b*l+b*l;
if(w1<=0) w11=0;
else if(w1>=a-i*b*l-hd*l) w11=a-i*b*l-hd*l;
else w11=w1;
end
q1=l*(a/(2*c)+b*l/(4*c)*(2*i-1)-w11/(2*c)-hd*l/(4*c));
lr21=w11*q1-v*q1-l*hm*(i-(j-1)*n)*q1;
lr11=((a-w11+hd*(i-1)*l)*(a-w11+hd*(i-1)*l)*l+(b-hd)*(a-w11+hd*(i-1)*l)*l*l*(2*i-1)+((b-hd)*(b-hd)*l*l*l*(3*i*i-3*i+1))/3)/(4*c);
if(w2<a+i*b*l-b*l)
w22=a+i*b*l-b*l;
else if(w2>=a+i*b*l-hd*l)
w22=a+i*b*l-hd*l;
else w22=w2;
end
q2=((a+b*l*i-hd*l-w22)*(a+b*l*i-hd*l-w22))/(4*c*(b-hd));
lr22=w22*q2-v*q2-l*hm*(i-(j-1)*n)*q2;
lr12=((a-w22+hd*(i-1)*l)*(a-w22+hd*(i-1)*l)*(a-w22+b*i*l-hd*l)/(b-hd)+(b-hd)*(a-w22+hd*(i-1)*l)*(a-w22+b*i*l-hd*l)/(b-hd)*(w22+b*i*l+hd*l-a-2*hd*i*l)/(b-hd)+(b-hd)*(b-hd)*(a-w22+b*i*l-hd*l)/(b-hd)*(((w22-a-hd*l*(i-1))*(w22-a-hd*l*(i-1)))/((b-hd)*(b-hd))+(w22-a-hd*l*(i-1))/(b-hd)*i*l+i*i*l*l)/3)/(4*c);
if(lr21>=lr22)
sprintf('i=j=%3f,w=%3f, q=%3f, lr2=%3f, lr1=3%f',i,w11,q1,lr21,lr11)
else
sprintf('i=j=%3f,w=%3f, q=%3f, lr2=%3f, lr1=3%f',i,w11,q1,lr21,lr11)
end
end
这是一位版内的前辈根据模型原来的C++源代码改编的,因为源代码里公式有很多错误,所以一直实现不了,给出这个代码只是想让大神根据这个格式就可以写出来,不用想的太复杂,如果您是MATLAB高手,一看数学模型就知道该怎么写代码,请您无视上面的代码,如果实现了输出,金币立马全给您。
图片里是模型参数和最后的运行结果。
大神,您在何方,小妹期待您的应助!!!   ![菜鸟求助数模代码]()
模型参数.png
![菜鸟求助数模代码-1]()
运行结果.png |