数学应用软件试题与答案

..本试卷共6大题,请写出每个问题求解的Matlab代码,并给出求解结果。1、求的最优解及其目标函数值。〔15分f=[-5-12-8-4];A=[-3-5-4023004851];b=[-20;14;200];Aeq=[12-11];beq=[45];vlb=[2100];vub=[];[x,fval]=linprog<f,A,b,Aeq,beq,vlb,vub>x=2.00001.000023.833364.8333fval=-472.0000最优解为472,目标函数为〔2.0000,1.0000,23.8333,64.8333共6页第1页2、求的最优解和极小值,初始值选[0,0]。〔10分functionf=ozk1<x>f=x<1>^4+3*x<1>^2+x<2>^2-2*x<1>-2*x<2>-2*<x<1>^2>*x<2>+6x0=[0,0];x=fminunc<‘ozk1’y=ozk1<x>x=1.00002.0000f=4.0000y=4.0000最优解是〔1,2,极小值为4。共6页第2页3、求的最优解和极小值,初始值选[0,0,0]。〔20分functionf=ozk2<x>f=2*x<1>^2+4*x<2>^2-2*x<1>*x<2>+x<3>^2+5*x<3>function[g,ceq]=mycon<x>g=[2*x<1>*x<3>-9;x<1>+2*x<1>*x<2>-8*x<3>+2];ceq=[];x0=[0;0;0];A=[132];b=[6];Aeq=[2-11];beq=[4];vlb=[0;0;0];vub=[];[x,fval]=fmincon<'ozk2',x0,A,b,Aeq,beq,vlb,vub,'mycon'>x=1.764700.4706fval=8.8028共6页第3页4、用4/5阶龙格-库塔-费尔贝格算法求微分方程在自变量区间[0,15]的数值解,用红色虚线作出函数图像表示。〔15分functiondy=ozk3<t,y>dy=zeros<1,1>;dy<2>=10*sin<y>;[T,Y]=ode45<'ozk3',[0,15],[160]>;plot<T,Y<:,1>,'r--'>共6页第4页5、试用下表所提供的数据求回归方程,其中方程中参数的初值为[0.5,0.5,0.5,0.1]。利用结果预测在=10处的值,并在同一幅图中画出原数据点〔用红色"*"表示和求得的函数关系图。〔20分234579121417212856354247535965687376828699functionf=ozk4<b,x>f=b<1>+b<2>*exp<b<3>.*x>+b<4>*exp<-b<3>.*x>>>x=[234579121417212856];>>y=[354247535965687376828699];>>b0=[0.5;0.5;0.5;0.1];>>[b,r,J]=nlinfit<x,y,'ozk4',b0>b=87.23590.02610.1094-62.8818共6页第5页6、<1>已知表格数据如下,试分别用二次多项式拟合和多元线性回归两种方法求出因变量对自变量的关系式。〔12分12345678921.947.261.970.872.866.450.325.33.2<2>在一天24小时内,从零点开始每间隔2小时测得的某物体变化数据为12.0,10.2,9.0,10.1,18.3,24.5,28.2,27.1,25.6,20.5,18.7,15.2,13.3,推测中午1点时的数据是多少？推测在每0.1小时的变化数据〔由于数据太多,不必计算显示后,描绘其变化曲线。〔8分〔1t=1:1:9;y=[21.947.261.970.872.866.450.325.33.2];A=p