五、某公司投资2000万元建成一条生产线。
投产后,在时刻t 的追加成本和追加收益分别为3/225)(t t t G ++=(百万元/年),3/218)(t t H -=(百万元/年)。
试确定该生产线在何时停产可获最大利润?最大利润是多少?
提示:利用函数⎰=T
G H t R 0t 20-d ))t (-)t (()((百万元),由于H (t )-G (t )单调
下降,所以H (t )=G (t )时,R (t )取得最大利润。
5.解:构造函数f(t)=H(t)-G(t)=13-t-3t 2/3=0 ;
令t 1/3=x,则f(t)=-t 3-3t 2+13
可得矩阵P=[-1,-3,0,13]
求最佳生产时间的源程序如下:
p=[-1,-3,0,13];
x=roots(p);
t=x.^3
运行结果如下:
t =
3.6768 +21.4316i
3.6768 -21.4316i
4.6465
再分别将t 的三个值带入函数f(t),比较大小后,得到最大利润与最佳生产时间。
求最大利润的程序代码如下:
① t=3.6768 +21.4316i;
x=0:0.01:t;
y=13-x-3*x.^(2/3);
trapz(x,y)
运行结果: ans =
25.2583
② t=3.6768 -21.4316i;
x=0:0.01:t;
y=13-x-3*x.^(2/3);
trapz(x,y)
运行结果: ans =
25.2583
③ t=4.6465;
x=0:0.01:t;
y=13-x-3*x.^(2/3);
trapz(x,y)
运行结果: ans =
26.3208
比较以上三组数据,可知最佳生产时间t=4.6465年,可获得的最大利润 26.3208(百万元/年)。
clear; close;
fplot('18-t^(2/3)',[0,20]);grid on;hold on;
fplot('5+t+2*t^(2/3)',[0,20],'r');hold off;
%发现t 约为4
[t,f,h]=fsolve('18-x^(2/3)-5-x-2*x^(2/3)',4)
%求得t=4.6465
t=linspace(0,t,100); y=18-t.^(2/3)-5-t-2*t.^(2/3); trapz(t,y)-20
%最大利润6.3232(百万元)
习题六
1. (科学计算)
A=randn(10,5)
(1)mean(A) ;均值std(A) ;标准方差
(2)max(max(A)) ;最大元素min(min(A)) ;最小元素(3)B=sum(A,2) ;A每行元素的和sum(B) ;A全部元素之和(4)sort(A) ;A的每列元素按升序排列
sort(A,2,’descend’) ;A的每行元素按将序排列
3.
x=[165 123 150 123 141];
y=[187 126 172 125 148];
P=polyfit(x,y,3)
P =
1.0e+003 *
-0.0000 0.0013 -0.1779 8.4330
所以它的线性拟合曲线为:p(x)=1.3x2—177.9x+8433
5.
(1)
建立函数文件:
function f=fxy(u)
x=u(1);y=u(2);
f=3.*x.^2+2*x.*y+y.^2
在命令窗口中输入以下命令:
[U,fmin]=fminsearch('fxy',[1,1])
结果:
U =
1.0e-004 *
-0.0675 0.1715
fmin =
1.9920e-010
(2)
f=inline('-sin(x)-cos(x.^2)'); fmax=fminbnd(f,0,pi) fmax =
0.7310
9.
(1)
矩阵求逆法:
A=[2 3 5;3 7 4;1 -7 1];
b=[10;3;5];
x=inv(A)*b
x =
-1.8060
-0.5373
3.0448
矩阵除法法:
A=[2 3 5;3 7 4;1 -7 1];
b=[10;3;5];
x=A\b
x =
-1.8060
-0.5373
3.0448
矩阵分解法:
A=[2 3 5;3 7 4;1 -7 1];
b=[10;3;5];
[L,U]=lu(A);
x=U\(L\b)
x =
-1.8060
-0.5373
3.0448
(2)方法同(1)
11.
(1)
f=inline('x-sin(x)./x');
x=fzero(f,0.5)
x =
0.8767
(2)
f=inline('(sin(x).^2).*exp(-0.1.*x)-0.5.*abs(x)');
x=fzero(f,1.5)
x =
1.6738
实验八第三题
%第三题
h=6:2:18;
x=6.5:2:17.5;
t1=[18,20,22,25,30,28,24];
t2=[15,19,24,28,34,32,30];
T1=spline(h,t1,x)
T2=spline(h,t2,x)
实验九第三题
绘图
硅谷公司
%第三题
n=input('请输入员工工号:');
h=input('该员工工作时数是:');
if h>120
x=(h-120)*84*(1+0.15)+120*84;
elseif h<60
x=h*84-700;
else
x=h*84;
end
disp([num2str(n),'号员工','的应发工资为',num2str(x)]); 实验三第一题
%第一题程序一
x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];
y=[]; %建立存放所有y值的矩阵for x0=x
if x0<0&x0~=-3
y=[y,x0*x0+x0-6];
elseif x0>=0&x0<5&x0~=2&x0~=3
y=[y,x0*x0-5*x0+6];
else
y=[y,x0*x0-x0-1];
end
end
x %输出所有x
y %输出所有y
%第一题程序二
x=[-5,-3,1,2,2.5,3,5];
y=[];
for a=1:7
if x(a)<0&x(a)~=-3
y=[y,(x(a))^2+x(a)-6];
elseif x(a)>=0&x(a)<5&x(a)~=2&x(a)~=3
y=[y,(x(a))^2-5*x(a)+6];
else
y=[y,x(a)*x(a)-x(a)-1];
end
end。