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')%计算