09.1《Matlab及其的应用》上机实践报告.docVIP

《Matlab及应用》上机实践报告 系 别： 专业班级： 学 号： 姓 名： 一、简易计算题：(第1-4小题每小题5分，第5小题10分，共30分) 1、计算，当和时的值。 方法一： 编写一个函数文件humps.m，即： function y=humps(x) y=x.^3+(x-0.98).^2/(x+1.35).^3-5*(x+1./x); 然后调用，在命令窗口中输入： y1=humps(2) y1 = -4.4723 y2=humps(4) y2 = 42.8096 （或者把y1和y2放在一起，通过下列方式求解） y=feval(humps,[2,4]) y = -4.4423 42.8077 方法二： f=inline( y=x.^3+(x-0.98).^2/(x+1.35).^3-5*(x+1./x)); f(2) ans = -4.4723 f(4) ans = 42.8096 （或者把y1和y2放在一起，通过下列方式求解） y=feval(f,[2,4]) y = -4.4423 42.8077 2、计算。 cos(60*pi/180)-(9-sqrt(2))^(1/3) ans = -1.4649 3、设，，求，。 A=[1 4 8;-3 6 -5;2 -7 -12]; B=[5 4 -2;6 -2 -8;-1 -9 7]; C=A*B C = 21 -76 22 26 21 -77 -20 130 -32 D=A.*B D = 5 16 -16 -18 -12 40 -2 63 -84 所以，。 4、设对称实矩阵，求其特征根。 a=[2 4 9;4 2 4;9 4 18]; p1=poly(a) %求特征多项式的系数 p1 = 1.0000 -22.0000 -37.0000 122.0000 roots(p1) %求特征根 ans = 23.3603 -3.0645 1.7042 所以，原对称实矩阵的特征根为 5、设为符号变量，，，试进行如下运算 （1）；（2）；（3）；（4）； （5）对进行因式分解 syms x fx gx fx=x^4+x^2+1; gx=x^3+4*x^2+5*x+8; fx+gx ans = x^4+5*x^2+9+x^3+5*x %（1）的结果 fx-gx ans = x^4-3*x^2-7-x^3-5*x %（2）的结果 expand(fx*gx) ans = x^7+4*x^6+6*x^5+12*x^4+6*x^3+12*x^2+5*x+8 %（3）的结果 expand(fx/gx) ans = 1/(x^3+4*x^2+5*x+8)*x^4+1/(x^3+4*x^2+5*x+8)*x^2+1/(x^3+4*x^2+5*x+8) %（4）的结果 factor(fx) ans = (x^2+x+1)*(x^2-x+1) %（5）的结果 二、：(每小题分共分) 1、在和处的导数值，的值由下表给出： 1.0 1.1 1.2 1.3 1.4 0.25 0.2268 0.2066 0.1890 0.1736 x=[1.0,1.1,1.2,1.3,1.4]; y=[0.25,0.2268,0.2066,0.1890,0.1736]; p1=polyfit(x,y,3) %3阶公式拟合后的多项式系数 p1 = -0.0667 0.3700 -0.7883 0.7350 x1=linspace(0.8,1.6,100); %在[0.8，1.6]区间取100个点 y1=polyval(p1,x1); %计算拟合多项式的值 plot(x,y,*b,x1,y1,r); %绘制原始数据和拟合曲线 legend(原始数据,3阶多项式); k=polyder(p1) %对拟合后的3阶多项式进行求导 k = -0.2000 0.7400 -0.7883 s1=polyval(k,1.0) %对求导后的多项式计算在的值 s1 = -0.2483 s2=polyval(k,1.2) %对求导后的多项式计算在的值 s2 = -0.1883 2、在某处测得海洋不同深处的水