matlab题库

操作题集操作题1《Matlab与科学计算》P28例2.1要求计算水在温度为0℃，20℃，40℃，60℃，80℃时的黏度，已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2)其中μ0为0℃水的黏度，值为1.758*10-3。要求分别以长短格式显示计算结果。formatshortemuw0=1.785e-3;a=0.03368;b=0.000221;t=0:20:80;muw=muw0./(1+a*t+b*t.^2)muw=1.7850e-0031.0131e-0036.6092e-0044.6772e-0043.4940e-004操作题2《Matlab与科学计算》P31例2.2数值数组【1：5】和字符串转换示例。a=[1:5];b=num2str(a)操作题3《Matlab与科学计算》P44例2.9比较用左除和右除法分别求解恰定方程的解。操作题4《Matlab与科学计算》P48例2.14计算矩阵magic(3)的指数，并比较不同函数的结果。b=magic(3);expm(b)expmdemo2(b)操作题5《Matlab与科学计算》P50例2.18计算下面矩阵[-149-50-154；537180546；-27-9-25]的特征值条件数。a=[-149-50-154;537180546;-27-9-25][V,D,s]=condeig(a)操作题6《Matlab与科学计算》P62例2.29对矩阵a=1111123413610141020进行三角抽取使之分别变成000011101000及123413001361014100141020a=pascal(4)tril(a,-1)tril(a,2)tril是下三角抽取，triu是上三角抽取操作题7《Matlab与科学计算》P62例2.30用3种方法建立向量v为[12620]多项式的伴随矩阵。v=[12620]a1=compan(v)a2=[-v(2:4);eye(2),zeros(2,1)]a3=-v(2:4);a3(2:3,1:2)=eye(2)操作题8《Matlab与科学计算》P64例2.31已知a=[21-3-1；3107；-124-2；10-15]；求该矩阵的3次方及数组的3次方。a=[21-3-1;3107;-124-2;10-15]a^3a.^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，～ba=[1:3;4:6;7:9];b=[010;101;001];x=5;y=ones(3)*5xa=x=aab=a&bB=~b操作题10《Matlab与科学计算》P71例2.41对多项式p=[2-56-19]求根，计算其在x=3的值，进行微分计算并写出其表达式。poly2sym数值2符号；polyval求多项式的值；roots求多项式的根；conv多项式的乘法(向量之卷积)conv(p,d);polyder多项式微分；polyfit多项式拟合。操作题11《Matlab与科学计算》P71例2.41用5阶多项式对[0,pi/2]上的正弦函数进行最小二乘拟合。x=0:pi/20:pi/2;y=sin(x);a=polyfit(x,y,5);x1=0:pi/30:pi*2;y1=sin(x1);y2=a(1)*x1.^5+a(2)*x1.^4+a(3)*x1.^2+a(4)*x1.^2+a(5)*x1+a(6);plot(x1,y1,'b-',x1,y2,'r*')legend('原曲线','拟合曲线')axis([0,7,-1.2,4])操作题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),等symsxa=[1/x1/(x+1);1/(x+2)1/(x+3)]b=[x,1;x+2,0]b-aa'a/ba^2a.^2det(a)%符号矩阵的行列式运算inv（b）符号矩阵的逆rank（a）符号矩阵的秩,操作题13《Matlab与科学计算》P88例3.7已知x∈（-2,2）,y∈（-2,2）,利用函数gradient绘制一个矢量图。[x,y]=meshgrid(-2:.2:2,-2:.2:2);z=x.*exp(-x.^2-y.^2);[px,py]=gradient(z,.2,.2);contour(z),%等高线绘制函数holdonquiver(px,py)%矢量图绘制函数操作题14《Matlab与科学计算》P101例4.1绘制如图所示带有显示属性设置的二维图形。操作题15《Matlab与科学计算》P104例4.5绘制如图所示的条形图和矢量图。操作题16《Matlab与科学计算》P104例4.65-绘制以下函数的图形。函数值范围：Lim=[0,2*pi,-1,1];函数为：y=sin(x);10.80.60.40.20-0.2-0.4-0.6-0.8-10123456操作题17《Matlab与科学计算》P105例4.7绘制如图所示饼状图(x=[2，4，6，8])操作题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）,。操作题20《Matlab与科学计算》P109例4.11作函数[-8,8],y[-8,8]上的三维网格图形。[x,y]=meshgrid(-8:0.5:8,-8:0.5:8)z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2)mesh(z)操作题21《Matlab与科学计算》P110例4.13已知作函数[-4,4],y[-4,4],使用meshc函数绘制三维面图。操作题22《Matlab与科学计算》P113例4.19绘制如图所示的柱面图。x=0:pi/20:pi*3;r=5+cos(x);[a,b,c]=cylinder(r,30);mesh(a,b,c)操作题23《Matlab与科学计算》P113例4.20绘制如图所示地球表面的气温分布示意图。操作题24《Matlab与科学计算》P118例4.24坐标标注函数应用示意图，如图所示。在x,x操作题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)。x=0:0.1*pi:2*pi;y=sin(x);z=cos(x);plot(x,y,'-*')holdonplot(x,z,'-o')plot(x,y+z,'-+')legend('sin(x)','cos(x)','sin(x)+cos(x)',0)操作题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)x=0:0.1*pi:2*pi;subplot(2,2,1);plot(x,sin(x),'-*')title('sin(x)')subplot(2,2,2)plot(x,cos(x),'-o')title('cos(x)')subplot(2,2,3)plot(x,sin(x).*cos(x),'-x')title('sin(x)*cos(x)')subplot(2,2,4)plot(x,sin(x)+cos(x),'-h')title('sin(x)+cos(x)')操作题27《Matlab与科学计算》P222例7.3已知x=0:0.1:10,y=sin(x),对上述数据点作插值，并加密数据点绘制插值点及插值函数。x=0:0.05:10;y=sin(x);x1=0:0.25:10;y1=interp1(x,y,x1);plot(x,y,'*',x1,y1)操作题28《Matlab与科学计算》P227例7.7设y=span{1，x，x2}，用最小二乘法拟合如表所示的数据。（polyfit功能函数进行拟合），并绘出数据点及拟合曲线x0.51.01.52.02.53.0y1.752.453.814.808.008.60x=[0.5:0.5:3.0];y=[1.752.453.814.808.008.60];a=polyfit(x,y,2)x1=0.5:0.05:3.0;y1=a(1)*x1.^2+a(2)*x1+a(3);plot(x1,y1,'-or',x,y,'-+b')操作题29《Matlab与科学计算》P228例7.8用最小二乘法求一个形如y=a+bx2的经验公式，使它与表所示的数据相拟合。并绘出数据点及拟合曲线xi1925313844yi19.032.349.073.398.8xi=19:6:44xi=1925313743yi=[19.032.349.073.398.8];formatshorta=polyfit(xi,yi,2)a=0.0635-0.59327.2811x1=19:.1:44;y1=a(1)*x1.^2+a(2)*x1+a(2);%此处发生了一点错误，将a（3）误输为a(2),后面会有改正plot(x1,y1,'-x')plot(x1,y1,'-')x2=xi.^2;%采用求解矩阵的方法来求解此拟合问题x2=[ones(5,1),x2']x2=1361162519611136911849ab=x2\yi'ab=-1.35220.0540y3=ab(1)+ab(2)*x1.^2;holdon;plot(x1,y3,'--g')y1=a(1)*x1.^2+a(2)*x1+a(3);plot(x1,y1,'-')%从图中可以看出，两种拟合方法较为吻合操作题30《Matlab与科学计算》P233例7.10分别采用矩形及梯形公式求积分∫3π0e-0.5tsin(t+π/6)dty=inline('exp(-0.5*t).*sin(t+pi/6)')quad(y,0,3*pi)操作题31《Matlab与科学计算》P237例7.11采用自适应Simpson求积公式求下列积分值。∫10(x/（x2+4)dx操作题32《Matlab与科学计算》P237例7.12采用quadl求积公式如下积分。∫31e-x/2dxquad操作题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操作题34《Matlab与科学计算》P246例7.18对下列矩阵进行LU分解。A=[123241467]U=UPshangdanjiaoL=low操作题35《Matlab与科学计算》P265例7.39求下列方程组的符合解：x1-0.7sinx1-0.2cosx2=0x2-0.7cosx1+0.2sinx2=0[x1,x2]=solve('x1-0.7*sin(x1)-0.2*cos(x2)=0','x2-0.7*cos(x1)