09.1《Matlab及应用》上机实践报告.doc

蒆《Matlab及应用》上机实践报告肃系别:专业班级:螁学号:姓名:肈一、简易计算题:(第1-4小题每小题5分,第5小题10分,共30分)蒆1、计算,当和时的值。蒄方法一:,即:袇functiony=humps(x)薆y=x.^3+(x-).^2/(x+).^3-5*(x+1./x);薁然后调用,在命令窗口中输入:羁>>y1=humps(2)蚆y1=蚆->>y2=humps(4)蒈y2=(或者把y1和y2放在一起,通过下列方式求解)莃>>y=feval('humps',[2,4])膀y=蒇-:薈>>f=inline('y=x.^3+(x-).^2/(x+).^3-5*(x+1./x)');膆>>f(2)羆ans=羀->>f(4)羅ans=(或者把y1和y2放在一起,通过下列方式求解)螈>>y=feval(f,[2,4])羈y=肅-、计算。螇>>cos(60*pi/180)-(9-sqrt(2))^(1/3)膅ans=膃-、设,,求,。薄>>A=[148;-36-5;2-7-12];蚀>>B=[54-2;6-2-8;-1-97];蕿>>C=A*B莅C=蚁21-7622莂2621-77莈-20130-32蒅>>D=A.*B肂D=衿516-16肆-18-1240薅-263-84蒂所以,。薁袅蚅4、设对称实矩阵,求其特征根。袃>>a=[249;424;9418];罿>>p1=poly(a)%求特征多项式的系数袈p1=-->>roots(p1)%求特征根螁ans=-,原对称实矩阵的特征根为蒆袄螂5、设为符号变量,,,试进行如下运算袁(1);(2);(3);(4);葿(5)对进行因式分解羄>>symsxfxgx膃>>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)的结果蚈薇二、数据拟合与插值:(每小题6分,共12分)莄1、用三阶公式计算在和处的导数值,的值由下表给出:>>x=[,,,,];节>>y=[,,,,];羇>>p1=polyfit(x,y,3)%3阶公式拟合后的多项式系数螃p1=芃-->>x1=linspace(,,100);%在[,]区间取100个点蚆>>y1=polyval(p1,x1);%计算拟合多项式的值螃>>plot(x,y,'*b',x1,y1,'r');%绘制原始数据和拟合曲线蚄>>legend('原始数据','3阶多项式');蒁螈>>k=polyder(p1)%对拟合后的3阶多项式进行求导袃k=袀-->>s1=polyval(k,)%对求导后的多项式计算在的值蒇s1=羃->>s2=polyval(k,)%对求导后的多项式计算在的值蚁s2=芆-、在某处测得海洋不同深处的水温如下:蒆深度(m),求在深度为500m、1000m、1500m处的水温。羂>>x=[446,714,950,1422,1634];蒈>>y=[,,,,];虿>>y1=interp1(x,y,500,'linear')%计算在500处的插值结果螆y1=>>y2=interp1(x,y,1000,'linear')%计算