晶片成本和晶片面积函数(matlab)

长度单位：cm，价格单位：美元

s_chip_min=0.5;
s_chip_max=4;
r=15;
c=4000;
b=2;
a=4;
s_chip=[s_chip_min:0.1:s_chip_min];
chip_num=(pi*r*r./s_chip) - (2*pi*r./sqrt(2*s_chip));
chip_ratio=(1+b*s_chip./a).^(-4);
cost=c./(chip_num.*chip_ratio);

a1=polyfit(s_chip,cost,2);
b1= polyval(a1,s_chip);

r1= sum((cost-b1).^2);       
plot(s_chip,cost,'*', 'LineWidth',1)          
hold on
plot(s_chip,b1, 'r', 'LineWidth',1)
disp(a1);