中山大学南方学院数值计与Matlab附答案.docxVIP

中山大学南方学院期末考试试卷（2013级电软系各专业2015～2016学年度第1学期）课程名称数值计算与Matlab A 卷考试形式开卷考核类型考查本试卷共一大题，卷面满分100分。题号一二三四五六七八总分复核人得分得分评卷人一、上机题目（请完成下面的上机题目，按照要求完成。共100 分）题目一、给出某地一月内温度的测定值（内含平均温度、最低、最高温度），并求其平均温度、最低、最高温度的平均值、标准差等。（参考下表) 周口四月份每天的温度报表温度日期 平均温度（℃） 最高温度（℃）最低温度（℃）110.21982140199415.8191259.715768.0164715.2208814.820791011971114327121322.629171420.727161522.331191613.015121716.925121813711620781582210621122417.825132518.624142622.729172723.528182824.529182918.627143018.92715源代码：A=[10.2 19 8;14.3 19 10;13.0 19 9;15.8 19 12;9.7 15 7;8.0 16 4;15.2 20 8;14.8 20 7;10.3 17 7;12.1 19 7;14.7 22 9;15.3 27 12;22.6 29 17;20.7 27 16;22.3 31 19;13.0 15 12;16.9 25 12;13.5 16 10;8.7 11 6;7.3 11 4;11.8 15 8;10.9 14 8;15.6 21 12;17.8 25 13;18.6 24 14;22.7 29 17;23.5 28 18;24.5 29 18;18.6 27 14;18.9 27 15]mean(A)std(A)题目二、用MATBAB软件，用二分法求方程在区间[1,2]内根的近似值，为使误差不超过10^-5时所需要的二分次数。代码：f=inline(x^3+4*x^2-10);a=1;b=2;p=b-a;i=0;d=0.00001while pd c=(a+b)/2;if f(a)*f(c)0 b=c;elseif f(a)*f(c)0 a=celse a=c;b=cend p=p/2; i=i+1;endx=(a+b)/2disp(i)结果：题目三、取，用牛顿迭代法求的跟的近似值。代码：function [i,A]=new(f,fd,x0,eps)x1=x0-f(x0)/fd(x0);i=1;A=[];A=[A,x0,x1];while abs(x1-x0)epsx0=x1;x1=x0-f(x0)/fd(x0); i=i+1; A=[A,x1];endf=inline(x0^3+4*x0^2-10);fd=inline(3*x0^2+8*x0);x0=1.5;eps=0.0005;[i,A]=new(f,fd,x0,eps)结果：题目四、弦割法求方程在区间[-2，-1]内的一个实根近似值，使.function [i,A]=new1(f,x0,x1,eps,NMAX)sn=(f(x1)-f(x0))/(x1-x0);x2=x1-f(x1)/sn;i=1;A=[];A=[A,x0,x1,x2];while (abs(x2-x1)eps iNMAX)x0=x1;x1=x2;sn=(f(x1)-f(x0))/(x1-x0);x2=x1-f(x1)/sn; i=i+1; A=[A,x2];endf=inline(x^3-3*x^2-x-9);x0=-2;x1=-1;NMAX=100;eps=0.00001;[i,A]=new1(f,x0,x1,eps,NMAX)结果：题目五、求的根。代码：clc; clearall;t=-5:1:5;x1=(sin(t)).^2;x2=-0.1.*t;y=x1.*exp(x2)-0.5.*abs(t);roots(y)结果：题目六、用高斯列主元算法编程求解下列线性方程组的解。代码：function X=Q6_Gauss(A,b) %A是一个N阶矩阵，b是n维向量，X是解[N N]=size(A);X=zeros(1,N+1);Aug=[A b]; %增广矩阵for p=1:N-1 [Y,j]=max(abs(Aug(p:N,p))); C=Aug(p,:);Aug(j+p