数学建模matlab作业 (2)

8，函数f(x)定义为如下的行列式： | x x^2 x^3 | | 1 2x 3x^2 | | 0 2 6x |

求f(x)在x=1处的数值导数，并画图比较数值导数与真实导函数之间的差别 syms x; A=[

x x^2 x^3 ; 1 2*x 3*x^2; 0 2 6*x ]; B=det(A); y=diff(B,x) subs(y,1) lims=[-5,5];

fplot('y',lims); hold on;

fplot('6*x^2',lims)

160140120100806040200-20-5-4-3-2-1012345

9，当x取值pi/6,2pi/6,...,5pi/6 时， 求 ( sin(x) )^2 + ( cos(2x) )^2的数值导数和数值二阶导数； syms x;

A=sin(x).^2+cos(2*x).^2; x=pi/6:pi/6:5*pi/6; subs(diff(A),x)

subs((diff(A,2)),x)

10,已知整数1,2,...,10的平方根，请用三次样条插值估算2.3的平方根 N=1:1:10; M=sqrt(N);

Y=interp1(N,M,2.3,'spline')

11, 某一过程测量得:

x：1.2000 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 y：2.9462 3.4094 3.8090 4.1643 4.4794 4.7654 5.0254 5.2668 请用函数y=a+b*log5(x)对此进行曲线拟合；

x=[1.2000,1.4000,1.6000,1.8000,2.0000,2.2000,2.4000,2.6000]; y=[2.9462,3.4094,3.8090,4.1643,4.4794,4.7654,5.0254,5.2668]; fun1=@(c,x) c(2).*log(x)./log(5)+c(1); b=lsqcurvefit(fun1,[1,1],x,y); plot(x,y,'+',x,fun1(b,x))

12，线性规划(linprog)：某工厂生产甲、乙、丙三种产品，单位产品所需工时分别为2、3、1个工时；

单位产品所需原材料分别为3、1、5公斤；单位产品利润分别为2元、3元、5元。 工厂每天可利用的工时为12个，可供应的原材料为15公斤。 试确定使总利润为最大的日生产计划和最大利润 f=[-2,-3,-5];

a=[2,3,1;3,1,5];b=[12;15];

[x,y]=linprog(f,a,b,[],[],[0,0,0]);%linprog是求min的 x=x z=-y

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新医药卫生数学建模matlab作业 (2)全文阅读和word下载服务。