matlab例题

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

五、某公司投资2000万元建成一条生产线。投产后,在时刻t的追加成本和追加收益分别为3/225)(tttG(百万元/年),3/218)(ttH(百万元/年)。试确定该生产线在何时停产可获最大利润?最大利润是多少?提示:利用函数TGHtR0t20-d))t(-)t(()((百万元),由于H(t)-G(t)单调下降,所以H(t)=G(t)时,R(t)取得最大利润。5.解:构造函数f(t)=H(t)-G(t)=13-t-3t2/3=0;令t1/3=x,则f(t)=-t3-3t2+13可得矩阵P=[-1,-3,0,13]求最佳生产时间的源程序如下:p=[-1,-3,0,13];x=roots(p);t=x.^3运行结果如下:t=3.6768+21.4316i3.6768-21.4316i4.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]);gridon;holdon;fplot('5+t+2*t^(2/3)',[0,20],'r');holdoff;%发现t约为4[t,f,h]=fsolve('18-x^(2/3)-5-x-2*x^(2/3)',4)%求得t=4.6465t=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=[165123150123141];y=[187126172125148];P=polyfit(x,y,3)P=1.0e+003*-0.00000.0013-0.17798.4330所以它的线性拟合曲线为:p(x)=1.3x2—177.9x+84335.(1)建立函数文件:functionf=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.06750.1715fmin=1.9920e-010(2)f=inline('-sin(x)-cos(x.^2)');fmax=fminbnd(f,0,pi)fmax=0.73109.(1)矩阵求逆法:A=[235;374;1-71];b=[10;3;5];x=inv(A)*bx=-1.8060-0.53733.0448矩阵除法法:A=[235;374;1-71];b=[10;3;5];x=A\bx=-1.8060-0.53733.0448矩阵分解法:A=[235;374;1-71];b=[10;3;5];[L,U]=lu(A);x=U\(L\b)x=-1.8060-0.53733.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('该员工工作时数是:');ifh120x=(h-120)*84*(1+0.15)+120*84;elseifh60x=h*84-700;elsex=h*84;enddisp([num2str(n),'号员工','的应发工资为',num2str(x)]);实验三第一题%第一题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[];%建立存放所有y值的矩阵forx0=xifx00&x0~=-3y=[y,x0*x0+x0-6];elseifx0=0&x05&x0~=2&x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx%输出所有xy%输出所有y%第一题程序二x=[-5,-3,1,2,2.5,3,5];y=[];fora=1:7ifx(a)0&x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseifx(a)=0&x(a)5&x(a)~=2&x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend

1 / 5
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功