试题三(非线性方程迭代法收殓性影子价格lag函数).docVIP

考试课程 数学实验 2001.01.05 班级 姓名 学号 说明：（1）第1，2，3题每题10分，直接将答案填在试题纸上； （2）第4题20分，将简要解题过程和结果写在试题纸上。 1． 给定方程。用迭代公式求该方程在区间[0，2]内的解，则迭代是 1 阶收敛的。若取迭代初值为0, 则 1.8971943063 （精确到10位小数）。 该方程的Newton迭代公式是： x(k+1)=x(k)-(x(k)-sin2(x(k))-1)/(1-sin(2*x(k))) 迭代是 2 阶收敛的。 判敛法：(先大致估计x*取值) 1)对于任意迭代公式x(k+1)=f(x(k))，其不动点为x* 收敛的充要条件(局部收敛性定理)：f(x)在x*的一个邻域内连续、可微，且，则对于该邻域内任意初值x0，序列{xk}均收敛于x* 2)若已判断其收敛，则将f(x)在x*作泰勒展开，若f’(x*)=f’’(x*)=……=f(p-1)(x*); f(p) (x*)≠0，则此迭代公式为p阶收敛 x0=0; x(1)=x0; for j=1:49 x(j+1)=sin(x(j))^2+1; end; j x(j+1) NEWTON法(2阶收敛): 对于方程f(x)=0，其牛顿法迭代公式为： x(k+1)=x(k)-(x(k)-sin2(x(k))-1)/(1-sin(2*x(k))) 2． 设 用数值解法算出 y(1)= 0.2714 （精确到4位小数）, 你用的方法是 龙格-库塔方法 ，调用的 Matlab命令是 ode45(@f, ts,y0) ，算法精度为 4阶 。 %待解常微分方程组函数M文件源程序： function dy=ff (x,y) dy=[y(2);y(2)*sin(x)-y(1)*exp(x)]; %应用欧拉方法和龙格-库塔方法求解该常微分方程： ts=:0.1:1; y0=[1,0]; [x,y]=ode45(@ff, ts,y0); %龙格-库塔方法求数值解 [, y(:,1)] 输出结果： 1.000000000000000 0.271423374882902 3. 假定取95%置信水平。从甲种汽油抽取5桶的含硫量为（%）：1.2, 1.5, 0.9, 0.8 1.3, 甲种汽油含硫量的置信区间为（精确到4位小数） [0.7823，1.4977] ；从乙种汽油抽取4桶的含硫量为（%）：1.2, 0.8, 1.5, 1.7, 要检验：甲的含硫量小于乙，作的零假设是 H0：μ甲≥μ乙Matlab命令是 ttest2(x,y) ，检验结果是 接受H0 ，作这个检验的前提是 数据来自正态总体，相互独立。 1) x=[1.2 1.5 0.9 0.8 1.3] [mux sigmax mucix sigmacix]=normfit(x) 输出结果： mucix = 0.782280071552989 1.497719928447011 2) x=[1.2 1.5 0.9 0.8 1.3] y=[1.2 0.8 1.5 1.7] [h,sig,ci]=ttest2(x,y,[],-1) 输出结果： h=0 4. 某厂用甲、乙两种原料添加填充剂制成一种新型材料，经实验知道，新型材料的强度y主要取决于单位体积内原料甲的含量x1(千克) 和原料乙的含量x2(千克)，并得到以下数据： y 105 150 128 125 131 129 146 129 161 138 x1 13.2 21.9 20.0 14.1 18.4 20.7 16.0 18.5 20.8 15.0 x2 15.0 20.9 13.7 19.3 17.8 14.1 24.0 18.2 25.0 25.0 在投入正式生产时得知，每千克原料甲含有毒物质5克，每千克原料乙含有毒物质1克，而单位体积新型材料中规定有毒物质不得超过100克；每千克原料甲售价300元，每千克原料乙售价400元，而生产单位体积新型材料用于购买甲、乙两种原料的成本预算限额为12,000元；由于技术原因，单位体积新型材料中原料乙的含量不得超过25千克。（精确到4位小数） （1）确定新型材料的强度y与单位体积内原料甲的含量x1和原料乙的含