MATLAB上机习题一-2016-解答

合集下载

matlab上机习题详细讲解-试题答案解析

matlab上机习题详细讲解-试题答案解析

学习指导参考P 第一次实验答案1. 设要求以0.01秒为间隔,求出y 的151个点,并求出其导数的值和曲线。

clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3) y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线(a,b,n 自定数据)clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线,要求有图形标注。

clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1,采用不同插值方法进行插值,并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。

matlab上机练习答案

matlab上机练习答案

实验一MA TLAB的基本命令与基本函数1已知矩阵a =11 12 13 1421 22 23 2431 32 33 3441 42 43 44求(1) A(:,1) (2) A(2,:)(3) A(:,2:3) (4) A(2:3,2:3)(5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:)(9) ones(2,2) (10) eye(2)(11) [A,[ones(2,2);eye(2)]](12) diag(A) (13) diag(A,1)(14) diag(A,-1) (15) diag(A,2)2(1)输入如下矩阵A0π/3A=π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵A中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵A中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵A的特征值与特征矢量:称特征矢量为M,而特征值矩阵为L(6) 求Msin(L)M-13已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2)其中μ0=1.785×10-3,a=0.03368,b=0.000221,求水在0,20,40,80℃时的黏度。

程序如下:miu0=1.785e-3;a=0.03368;b=0.000221;t=0:20:80miu=miu0./ (1+a*t+b*t.^2)(2)一个长管,其内表面半径为a,温度为Ta ;外表面半径为Tb;则其径向和切向应力可分别表示为:⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛+----=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛----=r b a b r b a b a a b v T T E r b a b r b a b a a b v T T E b a t b a r ln ln 11)/ln()1(2)(ln ln 1)/ln()1(2)(2222222222ασασ式中r 为管子的径向坐标,E 为管子材料的弹性模量,ɑ为热膨胀系数。

MATLAB 上机 习题及答案

MATLAB 上机 习题及答案

15、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5,要求先求得P(x)=P1(x)+P2(x),然后计算xi=0.2*i各点上的P(xi)(i=0,1,2,…,5)值。

p1=[1.0 0.0 0.0 -2.0 1.0];>> p2=[0.0 0.0 1.0 4.0 -0.5];>> p1x=poly2sym(p1);p2x=poly2sym(p2);>> p=p1x+p2xp =x^4+2*x+1/2+x^2>> x=0:5;>> x.^4+2*x+1/2+x.^2ans =0.5000 4.5000 24.5000 96.5000 280.5000 660.50001、试个MATLAB的工作空间中建立以下2个矩阵:A=[1 2]1234B⎡⎤=⎢⎥⎣⎦,求出矩阵A和B的乘积,并将结果赋给变量C。

>> A=[1 2]A =1 2>> B=[1 23 4]B =1 23 4>> C=A*BC =7 102、利用MATLAB提供的帮助信息,了解inv命令的调用格式,并作简要说明。

help invINV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.See also SLASH, PINV, COND, CONDEST, LSQNONNEG, LSCOV. Overloaded methodshelp gf/inv.mhelp zpk/inv.mhelp tf/inv.mhelp ss/inv.mhelp lti/inv.mhelp frd/inv.mhelp sym/inv.mhelp idmodel/inv.m3、使用help命令查询函数plot的功能以及调用方法,然后利用plot命令绘制函数y=sin(x)的图形,其中0xπ≤≤。

matlab上机习详解_试题答案

matlab上机习详解_试题答案

P 第一次实验答案1. 设要求以0.01秒为间隔,求出y 的151个点,并求出其导数的值和曲线。

clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线(a,b,n 自定数据)clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点,画出以x 为横坐标,y为纵坐标的曲线,要求有图形标注。

clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1,采用不同插值方法进行插值,并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。

Matlab上机练习题及答案

Matlab上机练习题及答案

Matlab 上机练习题及答案---------------------------------------------------------------------1、矩阵Y=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3472123100451150425,给出元素1的全下标和单下标,并用函数练习全下标和单下标的转换,求出元素100的存储位置。

取出子矩阵⎥⎦⎤⎢⎣⎡21301,并求该矩阵的维数。

解:命令为:Y=[5,2,4;0,15,1;45,100,23;21,47,3]Y(2,3)Y(10)sub2ind([43],2,3)[i,j]=ind2sub([43],10)find(Y==100)sub2ind([43],3,2)B=Y(2:2:4,3:-2:1)或B=Y([24],[31])[m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[10-1;241;-205],B=[0-10;213;112]求2A+B 、A 2-3B 、A*B 、B*A 、A.*B ,A/B 、A\B 解:命令为:A=[10-1;241;-205]B=[0-10;213;112]E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B---------------------------------------------------------------------3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵,并计算两者的乘积。

解:命令为:A=eye(3,4)B=8*ones(4)C=A*B---------------------------------------------------------------------4、创建矩阵a=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601,取出其前两列构成的矩阵b ,取出前两行构成矩阵c ,转置矩阵b构成矩阵d ,计算a*b 、c<d ,c&d,c|d ,~c|~d 解:命令为:a=[-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7]b=a(:,[12])c=a([12],:)d=b’e=a*b f=c<d g=c&d h=c|d i=~c|~d---------------------------------------------------------------------5、求!201∑=n n 解:命令文件为sum=0;s=1;for n=1:20s=n*s;sum=sum+s;end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值,其中a 是一个数字,由键盘输入,表达式中位数最多项a 的个数,也由键盘输入。

MATLAB期末上机试题带答案

MATLAB期末上机试题带答案

MATLAB 期末上机考试试题带答案版姓名: 学号: 成绩:1.请实现下图:xyy=sin(x)x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图:x=linspace(0,2*pi,100); y1=sin(x);subplot(2,2,1)plot(x,y1,'k--')grid onxlabel('x')ylabel('y')title('sin(x)')legend('y=sin(x)')y2=cos(x);subplot(2,2,2)plot(x,y2,'r--')grid onxlabel('x')ylabel('y')title('cos(x)')legend('y=cos(x)')y3=tan(x);subplot(2,2,3)plot(x,y3,'k-')grid onxlabel('x')ylabel('y')title('tan(x)')legend('y=tan(x)')y4=cot(x);subplot(2,2,4)plot(x,y4)grid onxlabel('x')ylabel('y')title('cot(x)')legend('y=cot(x)')3.解方程组:a=[3 2 1;1 -1 3;2 4 -4];b=[7;6;-2] ;x=a\b4.请实现下图:yxx=linspace(0,4*pi,1000);y1=sin(x);y2=sin(2*x);plot(x,y1,'--',x,y2,'b*')grid onxlabel('x');ylabel('y');title('耿蒙蒙')legend('sin(x)','sin(2*x)')5.请在x,y在(-2,2)内的z=xexp (-x2-y2) 绘制网格图[x,y]=meshgrid(-2:0.1:2);z=x.*exp (-x.^2-y.^2);mesh(x,y,z)6.请实现peaks函数:-55xPeaksy[x,y]=meshgrid(-3:1/8:3); z=peaks(x,y); mesh(x,y,z) surf(x,y,z) shading flataxis([-3 3 -3 3 -8 8])xlabel('x');ylabel('y');title('Peaks')7.请在x=[0,2],y=[-0.5*pi,7.5*pi],绘制光栅的振幅为0.4的三维正弦光栅。

matlab试题及答案

matlab试题及答案

matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是:A. 矩阵B. 向量C. 标量D. 数组答案:A2. 下列哪个命令可以用来绘制函数图形?A. `plot`B. `graph`C. `draw`D. `chart`答案:A3. MATLAB中,以下哪个是正确的矩阵转置操作?A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案:B二、简答题1. 简述MATLAB中矩阵的基本操作。

答案:在MATLAB中,矩阵是最基本的数据结构,可以进行加、减、乘、除等基本运算。

矩阵的创建可以使用方括号`[]`,例如`A = [1 2;3 4]`。

矩阵的转置使用单引号`'`,例如`A'`。

矩阵的求逆使用`inv`函数,例如`inv(A)`。

2. MATLAB中如何实现循环结构?答案:MATLAB中实现循环结构主要有两种方式:`for`循环和`while`循环。

`for`循环用于已知迭代次数的情况,例如:```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况,例如:```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B,请计算它们的乘积C,并求C的行列式。

A = [1 2; 3 4]B = [5 6; 7 8]答案:首先计算矩阵乘积C:```matlabC = A * B;```然后计算C的行列式:```matlabdetC = det(C);```结果为:```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数,计算并返回一个向量的范数。

答案:```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本,实现以下功能:- 随机生成一个3x3的矩阵。

MATLAB作业1参考答案_2016_

MATLAB作业1参考答案_2016_

MATLAB 作业1参考答案1 、试用符号元素工具箱支持的方式表达多项式5432()34236f x x x x x x =+++++,并令11s x s -=+,将f(x) 替换成s 的函数。

【求解】可以先定义出f 函数,则由subs() 函数将x 替换成s 的函数>> syms s x;f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;F=subs(f,x,(s-1)/(s+1)) F =(s-1)^5/(s+1)^5+3*(s-1)^4/(s+1)^4+4*(s-1)^3/(s+1)^3+2*(s-1)^2/(s+1)^2+3*(s-1)/(s+1)+62、用数值方法可以求出63626302124822i i S ===++++++∑,试不采用循环的形式求出和式的数值解。

由于数值方法采用double 形式进行计算的,难以保证有效位数字,所以结果不一定精确。

试采用符号运算的方法求该和式的精确值。

【求解】用符号运算的方式可以采用下面语句>> sum(2.^[1:63])ans =1.844674407370955e+019>> sum(sym(2).^[1:63])ans =18446744073709551614由于结果有19 位数值,所以用double 型不能精确表示结果,该数据类型最多表示16 位有效数字。

其实用符号运算方式可以任意保留有效数字,例如可以求200 项的和或1000 项的和可以由下面语句立即得出。

>> sum(sym(2).^[1:200])ans = 3213876088517980551083924184682325205044405987565585670602750>> sum(sym(2).^[1:1000])ans = 2143017214372534641896850098120003621122809623411067214887500776740 70 2102249872244986396757631391716255189345835106293650374290571384628 08 7196915514939714960786913554964846197084214921012474228375590836430 60 9294996716388253479753511833108789215412582914239295537308433532085 96633052487736744113361387503、编写一个矩阵相加函数mat_add() ,使其具体的调用格式为A=mat_add(A1 ,A2 ,A3 ,…) ,要求该函数能接受任意多个矩阵进行加法运算。

(完整版)matlab上机练习题答案

(完整版)matlab上机练习题答案

(完整版)matlab 上机练习题答案1.计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积>〉 a=[6 9 3;2 7 5]; 〉〉 b=[2 4 1;4 6 8]; 〉〉 a 。

*b ans =12 36 3 8 42 402。

对于B AX =,如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ,求解X 。

〉〉 A=[4 9 2;7 6 4;3 5 7]; 〉> B=[37 26 28]’; >〉 X=A\B X = -0.5118 4.0427 1.33183。

⎥⎦⎤⎢⎣⎡-=463521a ,⎥⎦⎤⎢⎣⎡-=263478b ,观察a 与b 之间的六种关系运算的结果 >〉 a=[1 2 3;4 5 6]; 〉> b=[8 –7 4;3 6 2]; >〉 a 〉b ans =0 1 0 1 0 1 〉> a 〉=b ans =0 1 0 1 0 1 >> a 〈b ans =1 0 1 0 1 0 〉> a<=b ans =1 0 1 0 1 0 >〉 a==b ans =0 0 0 0 0 0 〉〉 a~=b ans =1 1 1 1 1 14计算多项式乘法(x 2+2x +2)(x 2+5x +4)>> c=conv ([1 2 2],[1 5 4]) c =1 7 16 18 8 5计算多项式除法(3x 3+13x 2+6x +8)/(x +4) 〉〉 d=deconv ([3 13 6 8],[1 4]) d =3 1 26求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解〉〉 a=[2 4 7 4;9 3 5 6]; 〉> b=[8 5]’; 〉> x=pinv (a )*b x =—0.2151 0.4459 0。

Matlab上机练习一答案

Matlab上机练习一答案

Matlab上机练习⼀答案Matlab 上机练习⼀班级学号姓名按要求完成题⽬,并写下指令和运⾏结果。

(不需要画图)1、求下列联⽴⽅程的解=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59942、设 ??++=)1(sin 35.0cos 2x x x y ,把x=0-2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线。

>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')3、产⽣8×6阶的正态分布随机数矩阵R1, 求其各列的平均值、和、乘积。

并求该矩阵全体数的平均值。

(mean var )a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)或者u=reshape(a,1,48);p1=mean(u)p2=var(u)4、绘制)(222y x e x z +-=在定义域x=[-2,2],y=[-2,2]内的曲⾯。

(利⽤meshgrid )x=-2:2;y=x;[X,Y]= meshgrid(x,y);Z=X^2*exp(-(X^2+Y^2));mesh(X,Y,Z)5、求代数⽅程3x 5+4x 4+7x 3+2x 2+9x+12=0的所有根。

(利⽤roots 函数)p=[3 4 7 2 9 12];roots(p)6、把1开五次⽅,并求其全部五个根。

Matlab上机题及答案

Matlab上机题及答案

1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。

输出全部水仙花数。

for m=100:999m1=fix(m/100); %求m的百位数字m2=rem(fix(m/10),10); %求m的十位数字m3=rem(m,10); %求m的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend2.从键盘输入若干个数,当输入0时结束输入,求这些数的平均值和它们之和。

sum=0;n=0;val=input('Enter a number (end in 0):');while (val~=0)sum=sum+val;n=n+1;val=input('Enter a number (end in 0):');endif (n > 0)summean=sum/nend3. 若一个数等于它的各个真因子之和,则称该数为完数,如6=1+2+3,所以6是完数。

求[1,500]之间的全部完数。

for m=1:500s=0;for k=1:m/2if rem(m,k)==0s=s+k;endendif m==sdisp(m);endend4. 从键盘上输入数字星期,在屏幕上显示对应英文星期的单词。

function weekn=input('input the number:');if isempty(n)errror('please input !!')endif n>7|n<1error('n between 1 and 7')endswitch ncase 1disp('Monday')case 2disp('Tuesday')case 3disp('Wednesday')case 4disp('Thursday')case 5disp('Friday')case 6disp('Saturday')case 7disp('Sunday')end5. 某公司销售电脑打印机的价格方案如下:()如果顾客只买一台打印机,则一台的基本价格为$150。

matlab2016年课后习题及答案详解【精选】

matlab2016年课后习题及答案详解【精选】

2012年12月7日星期五第2章 MATLAB 矩阵运算基础2.1 在MATLAB 中如何建立矩阵,并将其赋予变量a ?⎥⎦⎤⎢⎣⎡194375>> a=[5 7 3;4 9 1]2.2 有几种建立矩阵的方法?各有什么优点?可以用四种方法建立矩阵:①直接输入法,如a=[2 5 7 3],优点是输入方法方便简捷;②通过M 文件建立矩阵,该方法适用于建立尺寸较大的矩阵,并且易于修改;③由函数建立,如y=sin(x),可以由MATLAB 的内部函数建立一些特殊矩阵;④通过数据文件建立,该方法可以调用由其他软件产生数据。

2.3 在进行算术运算时,数组运算和矩阵运算各有什么要求?进行数组运算的两个数组必须有相同的尺寸。

进行矩阵运算的两个矩阵必须满足矩阵运算规则,如矩阵a 与b 相乘(a*b )时必须满足a 的列数等于b 的行数。

2.4 数组运算和矩阵运算的运算符有什么区别?在加、减运算时数组运算与矩阵运算的运算符相同,乘、除和乘方运算时,在矩阵运算的运算符前加一个点即为数组运算,如a*b 为矩阵乘,a.*b 为数组乘。

2.5 计算矩阵与之和。

⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6];>> a+bans =777914131512142.6 求的共轭转置。

⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x >> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i 2.7 计算与的数组乘积。

Matlab上机题库及详细答案

Matlab上机题库及详细答案
例17:输入一个字符,若为大写字母,则输出其对应的小写字母;若为小写字母,则输出其对应的大写字母;若为数字字符则输出其对应的数值,若为其他字符则原样输出。
解:c=input('请输入一个字符','s');
if c>='A' & c<='Z'
disp(setstr(abs(c)+abs('a')-abs('A')));
elseif c>='a'& c<='z'
disp(setstr(abs(c)- abs('a')+abs('A')));
(2)均值为0.6方差为0.1的5阶正态分布随机矩阵
解:>>x=20+(50-20)*rand(5);
>>y=0.6+sqrt(0.1)*randn(5)
例13:将101~125等25个数填入一个5行5列的表格中,使其每行每列及对角线的和均为565。
解:M=100+magic(5)
M =117 124 101 108 115
-0.4606
0.3848
例2、用简短命令计算并绘制在0x6范围内的sin(2x)、sinx2、sin2x。
解:x=linspace(0,6)
y1=sin(2*x),y2=sin(x.^2),y3=(sin(x)).^2;
plot(x,y1,x, y2,x, y3)
例3:画出指数衰减曲线y1=exp(-t/3)*sin(3*t)和它的包络y2=exp(-t/3),t的取值范围是(0,4pi)。

matlab上机考试题及答案

matlab上机考试题及答案

matlab上机考试题及答案1. 题目:编写一个MATLAB函数,计算并返回一个向量中所有元素的平方和。

答案:函数定义如下:```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目:使用MATLAB的内置函数,找出一个矩阵中的最大元素及其位置。

答案:可以使用`max`函数来找出矩阵中的最大元素,同时使用`find`函数来获取其位置。

示例代码如下:```matlabA = [1, 2, 3; 4, 5, 6; 7, 8, 9];[maxValue, linearIndex] = max(A(:));[row, col] = ind2sub(size(A), linearIndex);```3. 题目:给定一个向量,使用MATLAB编写代码,实现向量元素的逆序排列。

答案:可以使用`flip`函数来实现向量的逆序排列。

示例代码如下:```matlabvector = [1, 2, 3, 4, 5];reversedVector = flip(vector);```4. 题目:编写一个MATLAB脚本,计算并绘制一个正弦波的图像。

答案:可以使用`sin`函数生成正弦波数据,并使用`plot`函数绘制图像。

示例代码如下:```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);xlabel('x');ylabel('sin(x)');title('Sine Wave');```5. 题目:给定一个3x3的矩阵,使用MATLAB编写代码,计算其行列式。

答案:可以使用`det`函数来计算矩阵的行列式。

示例代码如下:```matlabmatrix = [1, 2, 3; 4, 5, 6; 7, 8, 9];determinant = det(matrix);```结束语:以上是MATLAB上机考试的题目及答案,希望能够帮助大家更好地掌握MATLAB的编程技巧和函数使用。

Matlab上机题库及详细答案

Matlab上机题库及详细答案
rate=8/100;
case num2cell(25:49) %价格大于等于2500但小于5000
Байду номын сангаасrate=10/100;
otherwise %价格大于等于5000
rate=14/100;
end
price=price*(1-rate) %输出商品实际销售价格
例19已知,当n=100时,求的值。
解:程序如下:
y=0;n=100;for i=1:n;y=y+1/(2*i-1);End
例20:一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数
解:for m=100:999
m1=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);
if m==m1*m1*m1+m2*m2*m2+m3*m3*m3;disp(m);end
D=diag(1:5);D*A %用D左乘A,对A的每行乘以一个指定常数
例15:输入x,y的值,并将它们的值互换后输出。
程序如下:
x=input('Input x please.');y=input('Input y please.');z=x;x=y;y=z;
disp(x);disp(y);
例16::求一元二次方程ax2 +bx+c=0的根。
switch fix(price/100)
case {0,1} %价格小于200
rate=0;
case {2,3,4} %价格大于等于200但小于500
rate=3/100;
case num2cell(5:9) %价格大于等于500但小于1000

MatLab练习册1(答案)(1)

MatLab练习册1(答案)(1)

MatLab 练习题(1) 参考答案编题: 周德强 (长江大学数学学院)1, 基本运算符输入练习>> (5*2+1.3-0.8)*10^2/25 答案:422,计算y =23log )10sin(4)43.0(2+⋅+-e π 答案:2.2618>>sin(10*pi)*exp(-0.3+4^2)+log(23)/log(4) %将23log 4化成自然对数的形式 3,换行输入练习: >>z = 10*sin(pi/3)* ...cos(pi/3) 答案:4.33014,函数求值运算 >>x= sin(223/3); >>y = x^2; >>z = y*10;>>x+2y-5z 答案:-37.60255,ie z i z i z ⋅==+=π23215,2,72,计算3221z z z z z +=答案:-2.1379 + 1.6552i>>z1=2+7*i >>z2=2*i>>z3=5*exp(2*pi*i) >>z=(z1*z2)/(z2+z3)6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x 答案:(填写程序语句) >>x=(3:5.5:44) 或 >>x=[3:5.5:44] 或 >>x=3:5.5:44 7, 建立等差一维数组x :首项为0,末项为π,项数为15 答案:(填写程序语句) >>x=linspace(0,pi,15)8, 计算行列式的值187624323=A 答案:-8>>A = [3,2,3; 4,2,6; 7,8,1] %注意n 行n 列的矩阵输入方式>>det(A)9, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=187624323A ,矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积.答案:>>A = [3,2,3; 4,2,6; 7,8,1] >>B = [1,1,1; 2,2,2; 3,3,3] >> A*B16 16 16 26 26 2626 26 26>> A.*B %体会点运算 3 2 3 8 4 12 21 24 310,写出计算 Sin(30o )的程序语句.答案:(填写程序语句) >> sin(pi*30/180) 或 sin(pi/6)11,,计算f(2), f(1)f(2)+f 2(3) 答案:>> x=2>> f2=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14) f2 =-4.4792 >> x=1>> f1=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14) >> x=3>> f3=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14) >> f1*f2+f3^2-12.602312,,计算f 在x=[2 3 4 6]各处的函数值; 答案:-4.47920.9484 2.2741 4.3335方法1: >> x=[2 3 4 6];>> f=(x.^3 - 2*x.^2 + x - 6.3)./(x.^2 + 0.05*x - 3.14) %注意用点运算方法2: >> syms x %定义x为符号变量>> f=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14);%写出f的符号表达式>> f1=subs(f,x,[2 3 4 6]) %计算f在x=[2 3 4 6]各处的函数值MATLAB提供基本的算术运算有:加 (+)、减 (-)、乘 (*)、除 (/)、幂次方 (^),范例为:5+3, 5-3, 5*3, 5/3, 5^3MATLAB常用数学函数●三角函数和双曲函数●指数函数●复数函数●圆整函数和求余函数● 矩阵变换函数● 其他函数例:>>y = sin(10)*exp(-0.3*4^2) y = -0.0045【例】复数ie z i z i z 63212,21,43π=+=+=表达,及计算321z z z z =。

(完整版)matlab上机题目一

(完整版)matlab上机题目一

上机一1、(1(2) 查询所有人的成绩。

2、已知:A=[27,24,0,-1;23,35,1,14;e 4,-8,13,20;10,12,19,21;11,18,ln2,0], 求:(1)求A 中第4、5、2行与第2、4、1列组成的C 阵,并将C 中第2行元素全改成8,得到D 阵;(2) 创建与A 同维的G 阵,使C 阵左右翻转后位于G 中第2行至第4行、第2列至第4列位置,且G 中第1行第1列元素为1,其余元素为0。

(3) 求A 阵除以3的余数阵H 。

3、已知:A=[5,17,-2,3;-1,0,21,11;14,9,-10,8],B=[2,5,8,-1;12,-2,7,8;11,13,16,4] 求:(1)将A 除以B 阵的对应元素,并将结果朝负无穷大方向取整。

(2)求A./B , B.\A ,A.*B 和A*B ’的值。

(3)求2*ones(3,4)-eye(3,4)|A-B>B&B-8的值.4、要求在闭区间[0,2π ]上产生50个等距采样的一维数组。

试用两种不同的指令实现。

5、对矩阵,⎥⎦⎤⎢⎣⎡=4321A 问如下MATLAB 指令:)(),(),5.0(^),5.0(.^A sqrtm A sqrt A A 所得的结果相同吗?他们中的哪个结果为复数矩阵? 6、说出MATLAB 指令A (3,1,2,: )=1:4所产生数组的维数、大小和长度;然后对A 进行降维处理;最后指出所有非零元素的“全下标”位置。

7、现有五个数组: (2×3)的实数数组A ; (2×3)的复数数组B ; (2×3)的字符串数组C ; (2×3)的异质(元胞)数组D ; (2×3)的结构(构架)数组E 。

问,它们各占多少字节?8、写出使以下这段文字成为字符串的MATLAB 程序。

注意保持这段文字的格式。

在英式用法中,引号通常是单引号,如‘Fire!’。

MATLAB机试题及答案

MATLAB机试题及答案

上机题汇总1设置matlab 的工作环境,将工作目录设置为d:\work ,添加搜索目录d:\example设置当前目录:在Matlab 工具栏Current Directory 中输入或者浏览 设置搜索目录命令在资源管理器中创建work 文件夹 addpath('d:\work'); savepath;【也可以在file/ Set Path 路径设置窗口中完成】2在matlab 的命令窗口里完成如下计算,其中t 的值分别取-1,0,1,表达式如下:4/3)2ty eπ-=y=sqrt(2)/2*exp(-4*t).*sin(4*sqrt(3*t)+pi/3)3自行产生一个5行5列的数组,得到最中间的三行三列矩阵。

答:>>A=magic(5) I=[2 3 4];J=[2 3 4]; M=A(I,J)4用magic 产生一个5*5的矩阵,将这个矩阵的第二行与第三行互换位置答:>>A=magic(5) I=[1 3 2 4 5];J=[1 2 3 4 5]; M=A(I,J)5求方程组的根x 1+4x 2-3x 3=2 2x 1+5x 2-x 3=11 x 1+6x 2+x 3=12答:>>A=[1 4 -3;2 5 -1;-1 3 4];>>b=[2;11;12];>>x=A\b 或x=inv(A)*b6已知:一个多项式的系数向量是p=[1 -6-72 -27],求这个多项式的根。

答:>> p=[1 -6-72 -27] >> r=roots(p)7已经两个多项式的系数分别是:[1 2 3 4]和[1 4 9 16],请求这两个多项式的乘积,及商和余数。

答:>>p1=[1 2 3 4];p2=[1 4 9 16]; >>C=conv(p1,p2) >>[q,r]=deconv(p1,p2)8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式答:>>r=[-5 -3+4i -3-4i] >>p=poly(r) >>disp(poly2sym(p))9 A=[2 3 4;1 5 7;6 2 5]用什么函数,保证第一列排序的时候,其他列跟着变化。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

4. 在命令窗口输入 ispc,按键盘的 F1 键会弹出帮助,阅读并写出该内置变量的 含义。 (也可通过主菜单中[Help]->[Product Help]菜单打开帮助,并搜索 ispc 来完 成) 解答:
表示 MATLAB 的版本是否是 Windows 系统平台下的版本
5. 观察 MATLAB 命令窗口中关键字、字符串、注释、错误提示的显示颜色。提 示:输入相关表达式来观察,例如字符串可以这样输: ss=’abc’ 解答: 关键字:蓝色 字符串:紫色 注释:绿色需要 MATLAB7.0 以上版本,请确认系统装安装有。第 1、2 机房需要在开机时 选择【APPLY】系统。 截图方法:按键盘上右上角的[Print Screen]键,在作业文档中按 Ctrl+V 键粘贴即可。
1. 打开 MATLAB,观察界面中包含哪些子窗口,写出窗口名字 解答: (不需截图) 命令窗口( Command Window ) 、工作区窗口( Workspace ) 、历史命令窗口 (Command History) 、当前文件夹窗口(Current Folder) 、...
u
(e 2)
+ v)
2
v2 − u
3) 解答:
u = 2; v = 3;
u − 3v uv
%1) 4*u*v/log(v) %2) (exp(u)+v)^2/(v^2-u) %3) sqrt(u-3*v)/(u*v) 运行屏幕截图:
9. 计算如下表达式: 1) ( 3 − 5i )( 4 + 2i ) 2) sin ( 2 − 8i ) 解答:
2. 使用 path 命令将计算机的 C 盘添加到系统路径,使用 cd 命令将 C 盘设置为 当前工作路径 解答:(不需截图,仅需要写出相应的指令表达式即可) path('C:\', path) cd('C:\') %或者 path(path, 'C:\')
3. 在命令窗口中输入 iskeyword,观察 MATLAB 内置的关键字有哪些,尝试用其 中的某个作为变量名并赋值,观察 matlab 的输出。 (这些关键字是不允许作 为变量名的,MATLAB 会提示错误) 解答:
n2a + p nRT (V − nb ) = V2
该方程新增了 a 和 b 两个变量,用来表示纯净气体的属性。 请利用如下数据: 压强,p 摩尔数,n 体积,V a b 理想气体常数,R 220 mbar 2 mol 1L 5.536 L2bar/mol2 0.03049 L/mol 8.314 J/(mol.K)
提示,可以使用 isvarname 对上述变量名进行检验。 解答: 合法的变量名为:1、2、5、6、9、10 注意变量名命名规则:只能以字母开头,包含字母、数字或下划线,长度不能超 过 63 个,内置关键字不能作为变量名。
7. 创建 double 类型的变量,并进行计算 1)a=87,b=190,计算 a+b、a-b、a*b; 2)创建 uint8 类型的变量,数值与(1)中相同,进行相同的计算。 观察二者的结果有何不同,并解释为什么
7 5 + 6 × − 22 3 4) 2 3 × 3 3× 6 解答:
%1) sind(60) %或者 sin(60*pi/180) %2) exp(3) %3) cos(3/4*pi) %4) (5+6*7/3-2^2)/(2/3*3/(3*6)) 运行屏幕截图:
9. 设 u=2,v=3,计算: 1) 4 uv log v
MATLAB 上机习题一
请按以下步骤完成上机实验: 1)在 FTP 上下载“MATLAB 上机习题一.doc”文件,所有习题列在该文件内; 2)在 MATLAB 中完成所有习题,并将屏幕截图粘贴到相应习题后面; 3)如果习题是问答题,请将答案写在题目后; 4)如果有的习题要求提供脚本文件,请将脚本文件内容拷贝到相应习题后; ; 5)将文件保存并重命名为自己的学号,例如“20110771.doc” 6)上传该文件到 FTP 的相关目录。
%1) (3-5i)*(4+2i) %2) sin(2-8i) 运行屏幕截图:
10. 在高中课本中大家学过描述气体压强(p) 、温度(T) 、体积(V)和气体摩 尔数(n)之间的理想气体状态方程:
pV = nRT
其中, R 为理想气体常数。 上述状态方程描述的是低压强和高温度时气体的特性。 在 1873 年,范德瓦尔斯对这个方程进行了修正,使其更好地描述气体在不 同压强和温度条件下的状态,该方程被称为范德瓦尔斯方程:
分别使用理想气体状态方程和范德瓦尔斯方程计算水蒸气的温度,并使用 disp 函数显示计算结果。 (提示,注意单位制, 1 bar=105 Pa,1mbar=10-3bar,公式 中的各物理量单位应使用国际单位制) 解答:
%定义参数 p = 220 * 1E-3 * 1E5; %转换为Pa V = 1*1E-3; %转换为m^3 h = 15; n = 2; R = 8.314; a = 5.536 * 1E-6 * 1E5; %转换为m^2Pa/mol^2 b = 0.030491 * 1E-3; %转换为m^3/mol %1) 理想气体状态方程 T1 = p*V/(n*R); %2) 范德瓦尔斯方程 T2 = (p+n^2*a/V^2)*(V-n*b)/(n*R); %显示结果 disp(T1) disp(T2) 运行屏幕截图:
解答:
%1) a = 87; b = 190; a+b a-b a*b %2) a = uint8(87); b = uint8(190); a+b a-b a*b
屏幕运行截图:
8. 计算如下个式,写出 MATLAB 中的相应表达式: 1) sin ( 60 )
2) e3
3 3) cos π 4
6. 判断下列变量名哪些是合法的 MATLAB 变量: 1)fred 2)book_1 3)2ndplace 4)#1 5)vel_5 6)tan 7)fred! 8)book-1 9)Second_Place 10)No_1 11)vel.5 12)while
13)a123456789b123456789c123456789d123456789e123456789f1234567891234
相关文档
最新文档