中南大学材料院matlab操作题集答案

的491操作题集操作题1《Matlab与科学计算》P28例2.1要求计算水在温度为0℃，20℃，40℃，60℃，80℃时的黏度，已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2)其中μ0为0℃水的黏度，值为1.758*10-3。要求分别以长短格式显示计算结果。操作题2《Matlab与科学计算》P31例2.2数值数组【1：5】和字符串转换示例。操作题3《Matlab与科学计算》P44例2.9比较用左除和右除法分别求解恰定方程的解。操作题4《Matlab与科学计算》P48例2.14计算矩阵magic(3)的指数，并比较不同函数的结果。指数运算expm(a)操作题5《Matlab与科学计算》P50例2.18计算下面矩阵[-149-50-154；537180546；-27-9-25]的特征值条件数。特征值函数condeig(A)操作题6《Matlab与科学计算》P62例2.29对矩阵a=1111123413610141020进行三角抽取使之分别变成000011101000及123413001361014100141020操作题7《Matlab与科学计算》P62例2.30用3种方法建立向量v为[12620]多项式的伴随矩阵。操作题8《Matlab与科学计算》P64例2.31已知a=[21-3-1；3107；-124-2；10-15]；求该矩阵的3次方及数组的3次方。操作题9《Matlab与科学计算》P66例2.32数组逻辑运算演示。已知a=[1:3；4：6；7：9]，b=[010；101；001]，x=5,y=ones（3）*5,求x=a，ab=a&b，～b操作题10《Matlab与科学计算》P71例2.41对多项式p=[2-56-19]求根，计算其在x=3的值，进行微分计算并写出其表达式。操作题11《Matlab与科学计算》P71例2.41用5阶多项式对[0,pi/2]上的正弦函数进行最小二乘拟合。操作题12《Matlab与科学计算》P79定义a=[1/x,1/(x+1);1/(x+2),1/(x+3)]与b=[x,1;x+2,0]为符合矩阵，并求b-a,错误!未找到引用源。,错误!未找到引用源。,a2,det(a),错误!未找到引用源。,错误!未找到引用源。等操作题13《Matlab与科学计算》P88例3.7已知x∈（-2,2）,y∈（-2,2）,错误!未找到引用源。利用函数gradient绘制一个矢量图。操作题14《Matlab与科学计算》P101例4.1绘制如图所示带有显示属性设置的二维图形。操作题15《Matlab与科学计算》P104例4.5绘制如图所示的条形图bar和矢量图feather。操作题16《Matlab与科学计算》P104例4.65-绘制以下函数的图形。函数值范围：Lim=[0,2*pi,-1,1];函数为：y=sin(x);0123456-1-0.8-0.6-0.4-0.200.20.40.60.81fplot('[sin(x),cos(x)]',lim)操作题17《Matlab与科学计算》P105例4.7绘制如图所示饼状图(x=[2，4，6，8])注意{}pie操作题18《Matlab与科学计算》P107例4.9绘制如图所示三维螺旋线,x=[0,10*pi],y=sin(x),z=cos(x)。操作题19《Matlab与科学计算》P107例4.10绘制参数为矩阵的三维图.已知x∈（-2,2）,y∈（-2,2）,错误!未找到引用源。。Plot3操作题20《Matlab与科学计算》P109例4.11作函数在x[-8,8],y[-8,8]上的三维网格图形。Meshmeshgrid操作题21《Matlab与科学计算》P110例4.13已知作函数错误!未找到引用源。,x[-4,4],y[-4,4],使用meshc函数绘制三维面图。操作题22《Matlab与科学计算》P113例4.19绘制如图所示的柱面图。Mesh(a,b,c)操作题23《Matlab与科学计算》P113例4.20绘制如图所示地球表面的气温分布示意图。Surf(a,b,c)操作题24《Matlab与科学计算》P118例4.24坐标标注函数应用示意图，如图所示。Legendtitlexlable操作题25《Matlab与科学计算》P123例4.30在同一张图上绘制几个三角函数图(给出图例）；x范围x=0:0.1*pi:2*pi;函数为：y=sin(x);z=cos(x);需在同一张图上绘制出：sin(x);cos(x)；sin(x)+cos(x)。操作题26《Matlab与科学计算》P124例4.31-在4个子图中绘制不同的三角函数图。函数范围：x=0:0.1*pi:2*pi;函数为：sin(x);cos(x)；sin(x)+cos(x)；sin(x).*cos(x)subplot(m,n,p);操作题27《Matlab与科学计算》P222例7.3已知x=0:0.1:10,y=sin(x),对上述数据点作插值，并加密数据点绘制插值点及插值函数。Interp1pchip（H插值）spline（三次样方插值）cubic(三次插值）操作题28《Matlab与科学计算》P227例7.7设y=span{1，x，x2}，用最小二乘法拟合如表所示的数据。（polyfit功能函数进行拟合），并绘出数据点及拟合曲线polyfit(x,y,2)x0.51.01.52.02.53.0y1.752.453.814.808.008.60操作题29《Matlab与科学计算》P228例7.8用最小二乘法求一个形如y=a+bx2的经验公式，使它与表所示的数据相拟合。并绘出数据点及拟合曲线xi1925313844yi19.032.349.073.398.8A=[zeros(5,1),x']B=y'ab=a\b操作题30《Matlab与科学计算》P233例7.10分别采用矩形及梯形公式求积分∫3π0e-0.5tsin(t+π/6)dtFun.m......D=pi/1000;t=0:D:3*pi;y=fun(t);Z=trapz(y)*Dnt=lengtg(t);sc=cumsum(y)*D;scf=sc(nt)操作题31《Matlab与科学计算》P237例7.11采用自适应Simpson求积公式求下列积分值。∫10(x/（x2+4)dxquad操作题32《Matlab与科学计算》P237例7.12采用quadl求积公式如下积分。∫31e-x/2dxQuadl('fun',1,3,1e-10)操作题33《Matlab与科学计算》P246例7.17求解下列方程组。0.4096x1+0.1234x2+0.3678x3+0.2943x4=0.40430.2246x1+0.3872x2+0.4015x3+0.1129x4=0.15500.3645x1+0.1920x2+0.3781x3+0.0643x4=0.42400.1784x1+0.4002x2+0.2786x3+0.3927x4=-0.2557A=[]b=[]Ab=a\b操作题34《Matlab与科学计算》P246例7.18对下列矩阵进行LU分解。A=[123241467][l,u]=lu[a]操作题35《Matlab与科学计算》P265例7.39求下列方程组的符合解：x1-0.7sinx1-0.2cosx2=0x2-0.7cosx1+0.2sinx2=0Functionf=fc(x)Y(1)=Y(2)=Y=[y(1).y(2)];Fslove('fc',x0)操作题36《Matlab与科学计算》P273例7.42求微分方程y′=-y+x+1数值解。其中（0≤x≤1），y(0)=1Ode23('......',[0,1],1)操作题37《Matlab与科学计算》P273例7.45采用R-K方法求微分方程y′=-2y+2x2+2x数值解。其中（0≤x≤0.5），y(0)=1Ode23('......',[0,0.5],1)操作题38《Matlab与科学计算》P274例7.46解如下刚性方程：错误!未找到引用源。=错误!未找到引用源。+错误!未找到引用源。初始条件为：错误!未找到引用源。操作题39《Matlab与科学计算》P275例7.47求解常微分方程的符号解：y’=y+sin(t),y(pi/2)=0。Dsolve('......','......')操作题40《Matlab与科学计算》P307例9.21设有三台机器，用来生产规格相同的铝合金薄板。取样、测量薄板的厚度精确至千分之一厘米。得结果如下：机器10.2360.2380.2480.2450.243机器20.2570.2530.2550.2540.261机器30.2580.2640.2590.2670.262检验各台机器所生产的薄板的厚度有无显著的差异？Anova1(x')操作题41《Matlab与科学计算》P308例9.22一火箭使用了4种燃料，3种推进器做射程试验，每种燃料与每种推进器的组合各发射火箭两次，得到结果如下。燃料（A）推进器（B）B1B2B3A158.200056.200065.300052.600041.200060.8000A249.100054.100051.600042.800050.500048.4000A360.100070.900039.200058.300073.200040.7000A475.800058.200048.700071.500051.000041.4000考察推进器和燃料这两个因素对射程是否有显著的影响？Anova2(x,2)操作题42《Matlab与科学计算》P309例9.23为了研究某一化学反应过程中，温度X对产品得率Y的影响，测得数据如下。温度x100110120130140150160170180190得率y45515461667074788589试做y=a+bx型的回归。polyfit操作题43分别用符合求解和数值求解symsxquad('log(1+x)',0,1)ans=0.3863int('log(1+x)',x,0,1)ans=log(4)-1操作题44用ode23,ode45,ode113求解多阶常微分方程错误!未找到引用源。,错误!未找到引用源。,错误!未找到引用源。,x∈[1,10]。建立求解函数文件myfun03functiondy=myfun03(x,y)dy=zeros(3,1);%初始化变量dy，改行可以没有dy(1)=y(2);%dy(1)表示y的一阶导数,其等于y的第二列值dy(2)=y(3);%dy(2)表示y的二阶导数dy(3)=2*y(3)/x^3+3*y(2)/x^3+3*exp(2*x)/x^3%dy(3)表示y的三阶导数求解过程：[x45,y45]=ode45('myfun03',[1,10],[11030]);操作题45求，在x=0.5附近的根。操作题46求；的梯度函数δx,δy和法向量操作题47任意构造一个m文件，并存于d盘某文件，把此文件夹路径设置为搜索路径。操作题48计算多项式p（x）=x3+21x2+20x的根，并由多项式构建原多项式。操作题49求解下面的极限。(1)(2)(3)（4）Limit(f,x,x0)操作题50求函数错误!未找到引用源。的1阶导数。symsxf=(5*x^3+3*x^2-2*x+7)/(-4*x^3+8*x+3);diff(f)ans=(15*x^2+6*x-2)/(-4*x^3+8*x+3)+((12*x^2-8)*(5*x^3+3*x^2-2*x+7))/(-4*x^3+8*x+3)^2操作题51何用符合求积的方法计算下面的积分。Int(f,x,a,b)（1）int((x+sin(x))/(1+cos(x)),x)ans=x*tan(x/2)(2)int(sqrt(log(1/x)),x,0,1)ans=pi^(1/2)/2(3)a=[cos(x)x^2;2^xlog(2+x)];int