MATLAB上机习题一
MATLAB上机模拟考试题
M A T L A B上机模拟考试题(一)(总8页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--(1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。
（）off(2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。
（） long rat hex short e(3)下列变量名中_____是合法的。
（） *y,a,1 \y,a1234 ,1 bcx ,i,j(4)已知x=0:5,则x有_____个元素。
（）(5)一下运算符中哪个的优先级最高_____。
（） A./ B.^ C.~=D.&(6)使用检测函数isnumeric(10)的结果是_____。
（）(7)三维图形中默认视角是_____度。
（） A.方位角=0 俯仰角=90B.方位角=90 俯仰角=0C.方位角= 仰俯角=30D.方位角=0 仰俯角=180(8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。
（）(9)运行以下命令，则_____描述是正确的。
（）>>syms a b cd >>A=[a b;c d] 占用的内存小于100B B.创建了5个符号变量占用的内存是a b c d的总和 D.不存在（10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。
（11）运行命令bitor（8,7）的结果是_____。
（12）运行以下命令：>>x=0:10; >>y1=sin(x); >>y2=5*sin(x);>>y3=[10*sin(x);20*sin(x)]; >>plot(x,y1,x,y2,x,y3) 则在一个图形窗口中，可以看到_____条曲线。
（13）符号表达式“g=sym(sin(a*z)+cos(w*v))”中的自由符号变量是_____。
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 上机练习题及答案---------------------------------------------------------------------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 上机练习题（一）
1. 计算下式
7
2009/42008/41949222?+-⨯-⨯=
2. 试计算d
3461.1,30,5b d a b c
a b c +=-+===其中
要求只输出d 的值，其他变量的值不输出。
**************************************************************
3. 计算下式：
181(2)sin()cos()6,3
y
x x y π-+-==（，其中 0.1104
2log 805(3)2009
6(4)1,5;
2e
e x x x -⨯+-=其中
4. 完成以上题目后清除工作空间中的所有变量。
**************************************************************
5. 求下列方程的根：
（1）325x x -=
（2）32340x x x -+=
6. 根据课本任意找出4个矩阵在MA TLAB 中输入和输出。
附1：增加的一些常用函数pi：常数
round(x)：四舍五入至最近整数
fix(x) ：无论正负，舍去小数至最近整数min(x): 向量x的元素的最小值
max(x): 向量x的元素的最大值
diff(x): 向量x的相邻元素的差
length(x): 向量x的元素个数
sin(x)：正弦函数
cos(x)：余弦函数
tan(x)：正切函数
asin(x)：反正弦函数
acos(x)：反余弦函数
atan(x)：反正切函数。
MATLAB上机题目
1、熟悉MATLAB 运行环境，找出命令历史窗口、工作空间窗口，记下MATLAB 当前浏览路径。
2、求下列联立方程的解
8102569
583
24754
12743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x
3、MATLAB 环境下，给出下列值的大小
eps=？ Realmax=? Realmin=?
4、生成矩阵
a=102030405060710⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
，取出第三行第二列元素；一、二行，二至五列元素（用一个命令）；
a(17)=?
若a(1,2)=[],则a=？
若a(20)=4.5,则a=？
5、i=inf,j=NaN,则i/j=? i*j=?
6、生成矩阵
b=123435672i i i i +-⎡⎤⎢⎥-⎢⎥⎢⎥+⎣⎦
,求其共轭、转置、共轭转置矩阵 7、检查工作空间中的变量，以及这些变量的详细信息
8、列出2×2阶的单位矩阵I ，4×4阶的魔方矩阵M 和4×2阶的全1矩阵A 、全零矩阵B 。
(2)将这些矩阵拼接为6×6阶的矩阵C ：
⎥⎦
⎤⎢⎣⎡=M B A I C ' (3)取出C 的第2，4，6行，组成3×6阶的矩阵1C ，取出第2，4，6列，组成6×3阶的矩阵2C
9、产生8×6阶的正态分布随机数矩阵1R ,求其各列的平均值和标准差，并求其全体的平均值和标准差。
10、产生4×6阶的均匀分布随机数矩阵R ，要求其元素在1到16之间，并求此矩阵前四列组成的方正的逆阵。
matlab习题(1)
1.输入20个数，求其中最大数和最小数。
要求分别用循环结构和调用MA TLAB的max函数、min函数来实现。
a=rand(1,20);max1=max(a) %用max函数求最大值min1=min(a) %用max函数求最大值%以下用循环方法求解：max2=a(1);min2=a(1);for i=2:20if a(i)>max2max2=a(i);endif a(i)<min2min2=a(i);endendmax2min2调用MATLAB的max函数、min函数来实现:>> a=rand(1,20)zuidazhi=max(a)zuixiaozhi=min(a)a =Columns 1 through 80.7513 0.2551 0.5060 0.6991 0.8909 0.9593 0.5472 0.1386 Columns 9 through 160.1493 0.2575 0.8407 0.2543 0.8143 0.2435 0.9293 0.3500 Columns 17 through 200.1966 0.2511 0.6160 0.4733zuidazhi =0.9593zuixiaozhi =0.13862. 求Fibonacci数列(1)大于4000的最小项。
(2)5000之内的项数。
(1)function t=fibonacc1f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendt=f(n);return(2)function n=fibonacc2f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendn=n-1;return3. 写出下列程序的输出结果:s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0s=s+k(j);endendendSAns=1081、解方程组Ax ＝b ，分别用求逆解法与直接解法求其解。
matlab上机练习（附答案）
1.以下两种说法对吗？（1）MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。
（2）MATLAB指令窗中显示的数据有效位数不超过七位。
2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同？DIARY filename causes a copy of all subsequent command window inputand most of the resulting command window output to be appended to thenamed file. If no file is specified, the file 'diary' is used.DIARY OFF suspends it.DIARY ON turns it back on.DIARY, by itself, toggles the diary state.Use the functional form of DIARY, such as DIARY('file'),when the file name is stored in a string.3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上？“位于搜索路径上的目录”与“当前目录”在MATLAB中的功用相同吗？4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组？用直接键入法？用数组编辑器？用M文件编辑器？5.运用数组算术运算符去掉下面程序里的for/end循环：x=11:15for k=1:length(x)z(k)=x(k)^2+2.3*x(k)^0.5;endx=11:15 z1=x.^2+2.3*x.^0.56.不使用数组算术运算符，重写下面的程序代码：x=[2 1 4]z=1./(1+x.^2)x=2;k=1;while i<=4,z2(k)=1/(1+i^2);i=i+1;x=x+1;end7．某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
matlab上机练习答案
实验一MA TLAB的基本命令与基本函数１已知矩阵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)输入如下矩阵Ａ０π/3Ａ＝π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(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上机练习一参考解答
上机练习一参考解答一、实验目的1、 熟悉Matlab 编程2、 体会数学上恒等，算法上不一定恒等二、实验内容1． Using the Taylor polynomial of degree nine and three-digit rounding arithmetic to find an approximationto 5-e by each of the following methods.(A) ∑=--≈905!)5(n n n e , (B) ∑=-≈=9055!5/11n nn e e An approximate value of 5-e correct to three digits is 31074.6-⨯. Which formula, (A) or (B), gives the most accuracy, and why?1） 算法基础利用x e 的Taylor 公式00!!n nk x n n x x e n n ∞===≈∑∑，x -∞<<+∞ （1）及001/1/1/!!n nk x x n n x x e e n n ∞-====≈∑∑，x -∞<<+∞， （2）其中k 是根据精度要求给定的一个参数。
在本题中将k 取为9， x 取为－5或5即可由公式（1）或（2）得到5-e 的近似计算方法（A ）或（B ）。
2） 程序下述程序用公式（A ）及（B ）分别在Matlab 许可精度下及限定在字长为3的算术运算情况下给出5-e的近似计算结果，其中results_1, results_2为用方法（A ）在上述两种情况下的计算结果，err_1, err_2为相应的绝对误差；类似的，results_3, results_4为用方法（B ）在上述两种情况下的计算结果，err_3, err_4为相应的绝对误差；具体程序如下：% Numerical Experiment 1.1 % by Xu Minghua, May 17, 2008 clc; %Initialize the data x=-5; k=9; m=3; %three-digit rounding arithmetic %------------------------------------ % Compute exp(x) by using Method (A) % with the computer precision results_1=1; power_x=1; for i=1:k factor_x=x/i; power_x=power_x*factor_x; results_1=results_1+power_x; end results_1 err_1=abs(exp(x)-results_1)%------------------------------------% Compute exp(x) by using Method (A) % with the 3-digits precisionresults_2=1;power_x=1;for i=1:kfactor_x=digit(x/i,m);power_x=digit(power_x*factor_x,m); results_2=digit(results_2+power_x,m); endresults_2err_2=abs(exp(x)-results_2)%------------------------------------% Compute exp(x) by using Method (B) % with the computer precisiont=-x;results_3=1;power_x=1;for i=1:kfactor_x=t/i;power_x=power_x*factor_x;results_3=results_3+power_x; endresults_3=1/results_3err_3=abs(exp(x)-results_3)%------------------------------------% Compute exp(x) by using Method (B) % with the 3-digits precisiont=-x; results_4=1;power_x=1;for i=1:kfactor_x=digit(t/i,m);power_x=digit(power_x*factor_x,m);results_4=digit(results_4+power_x,m); endresults_4=digit(1/results_4,m)err_4=abs(exp(x)-results_4)%------------------------------------上述主程序用到一个子程序digit.m, digit(x,m)的作用是将x四舍五入成m位数。
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。
Matlab上机实验答案
0.2000 0.1667 0.1429 0.1250 0.1111
1 1 1 1 1
1 2 3 4 5
1 3 6 10 15
1 4 10 20 35
1 5 15 35 70
>> Hh=det(H) Hh =
3.7493e-12 >> Hp=det(P) Hp =
1
>> Th=cond(H) Th =
0 0 1.0000 0 0
0.8147 0.9058 0.1270 1.0000 0
0.9134 0.6324 0.0975 0 2.0000
>> A^2
ans =
1.0000 0 0 0 0 >> [E,R+R*S;O,S^2]
0 1.0000 0 0 0
0 0 1.0000 0 0
1.6294 1.8116 0.2540 1.0000 0
阵和对角阵，试通过数值计算验证 A >> E=eye(3); >> R=rand(3,2); >> O=zeros(2,3); >> S=diag(1:2); >> A=[E,R;O,S]
E O
R RS 。 S2
A=
1.0000 0 0 0 0
0 1.0000 0 0 0
3. 建立一个 5×5 矩阵，求它的行列式值、迹、秩和范数。 >> A=rand(5) A= 0.2785 0.5469 0.9575 0.9649 0.1576 >> det(A) ans = -0.1322 >> trace(A) ans = 3.4127 >> rank(A) ans = 0.9706 0.9572 0.4854 0.8003 0.1419 0.4218 0.9157 0.7922 0.9595 0.6557 0.0357 0.8491 0.9340 0.6787 0.7577 0.7431 0.3922 0.6555 0.1712 0.7060
Matlab上机题目
Matlab上机题目上机一MATLAB集成环境使用与运算基础1、在计算机上安装MATLAB系统，并熟悉其调试环境，熟悉4个主要的工作窗口（命令窗口、工作区管理窗口、当前目录窗口、历史窗口）。
2、通过Internet访问MATLAB的有关资源，了解其应用。
3、我国人口按2000年第五次全国人口普查的结果为12.9533亿，如果年增长率为1.07%，求公元2010年末的人口数。
提示：计算人口的公式为10(1)nP P r=+，其中，1P为几年后的人口，P为人口初值，r为年增长率，n为年数。
4、当MATLAB启动后，系统会自动建立一个工作空间，这时的工作空间内只包含系统所提供的一些特殊变量，如pi，i等，以后随着用户求解问题过程的深入，会逐渐增加一些用户自己定义的变量。
如不采用clear命令来删除这些内存变量，这些变量会一直存在下去，直到用户关闭MATLAB系统释放工作空间后才会消失。
（1）在工作空间创建2*3 3*3 1*1 2*1四个矩阵A B C D；（2）熟悉变量的查看命令who和whos；（3）熟悉变量的文件保存与获取命令save和lode：比较savesave filenamesave filename A Bsave filename C D -append几种调用格式的区别。
5、输入矩阵A=[1，3，5；5，8，3；6，1，6]、B=[3，6；9，3；4，7]、C=[3，7,9；4，0，7]、D=2：6，体会命令[A,B]、[A;C]、[A,B;D]所产生的结果，总结由小矩阵生成大矩阵的方法。
A=[1,3,5;5,8,3;6,1,6]A =1 3 55 8 36 1 6>> B=[3,6;9,6;4,7]B =3 69 64 7>> C=[3,7,9;4,0,7]C =3 7 94 0 7>> D=2:6D =2 3 4 5 6 >> [A,B]ans =1 3 5 3 65 8 3 9 66 1 6 47 >> [A;C]ans =1 3 5 5 8 3 6 1 6 3 7 9 4 0 7>> [A,B;D] ans =1 3 5 3 6 5 8 3 9 6 6 1 6 4 72345 66、在M 文件编辑器中创建一个名为wl_1的M 函数文件，计算3421( 3.2)/(7.3)6()y x x x x =-+-+在x=5和7时的值，实现在命令窗口中调用M 文件，输入自变量的值时输出函数值。
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上机习题一 解答
9. 判断下面语句的运算结果，并与 MATLAB 计算结果做比较： 1） 4 < 20 2） 4 ≤ 20 3） 4 == 20 4） 4 ~ = 20 5） ' b ' < ' B ' 解答： 1）true 2）true 3）false 4）true 5）false
注意，字符比较实际上比较的是字符的 ASCII 码的大小。
提示，可以使用 isvarname 对上述变量名进行检验。
解答：1、2、5、6、9、10 是合法的，其它不能作为变量名。 注意，第 6 问的函写脚本程序，命名为 barbell.m，完成如下计算： 如图 1 所示的杠铃，
每个球的半径是 10cm，两个求直接链接杆的长度是 15cm，杆的直径为 1cm，计 算杠铃的体积和表面积；使用 disp 函数显示计算结果。 （提示：disp 函数的用法 参考教材的例 2-31， 面积计算时相交部分球冠可近似为圆， 体积计算时可忽略重 叠的球冠体积） 解答： barbell.m（拷贝文件内容到如下方框中） ： %注意，将圆柱与球相交部分面积近似为圆形 clear all clc R=10; %球半径,单位cm l=15; %杆长度,单位cm d=1; %杆直径,单位为cm
n2a + p nRT (V − nb ) = V2
该方程新增了 a 和 b 两个变量，用来表示纯净气体的属性。 1）编写脚本程序，命名为 state_equation.m，利用如下数据： 压强，p 摩尔数，n 体积，V a b 理想气体常数，R 220 mbar 2 mol 1L 5.536 L2bar/mol2 0.03049 L/mol 8.314 J/(mol.K)
1. 运行 MATLAB 软件，拖放、关闭界面上的子窗口，并恢复到原始视图。 解答：略。注意采用下图方式恢复视图
matlab上机题汇总(必做)
上机题汇总1设置matlab 的工作环境，将工作目录设置为d:\work ，添加搜索目录d:\example2在matlab 的命令窗口里完成如下计算，其中t 的值分别取-1,0,1，表达式如下：4/3)2t y e π-=3自行产生一个5行5列的数组，得到最中间的三行三行矩阵。
4用magic 产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置5求方程组的根 x 1+4x 2-3x 3=22x 1+5x 2-x 3=11x 1+6x 2+x 3=126已知：一个多项式的系数向量是p=[1 -6-72 -27]，求这个多项式的根。
7已经两个多项式的系数分别是：[1 2 3 4]和[1 4 9 16]，请求这两个多项式的乘积，及商和余数。
8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式9 A=[2 3 4；1 5 7；6 2 5]用什么函数，保证第一列排序的时候，其他列跟着变化。
最后的结果是：[1 5 7；2 3 4；6 2 5]10补充题：电路分析电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。
下图描述了3个回路电压的方程方程式如下：11121222132423432532()0()()0()0V R I R I I R I I R I R I I R I I R I V -++-=⎧⎪-++-=⎨⎪-++=⎩假设5个电阻值为已知，2个电压值也为已知，求3个电流值。
11自行产生一个5行5列的数组，用两种方法得到最中间的三行三行矩阵。
12根据a=reshape(-4:4,3,3)做一个矩阵，然后（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变)。
(完整版)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.计算⎥⎦⎤⎢⎣⎡=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上机题目一答案
上机一1、（1（2） 查询所有人的成绩。
J=struct(…number ‟,{9901,9902,9903,9904,9905},‟name ‟,{…Xiao li ‟, … …Xiao zhou ‟,‟Xiao xu ‟,‟Xiao liao ‟,‟Li na ‟},… ‟score ‟,{90,83,71,82,64},‟grade ‟,{…A ‟,‟B ‟,‟B ‟,‟A ‟,‟A ‟}); A=J(3),B=[J.score], 运行结果：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 。
A=[27,24,0,-1;23,35,1,14;exp(4),-8,13,20;10,12,19,21;11,18,log(2),0], C=A([4,5,2],[2,4,1]),D=C;D(2,:)=8, G=zeros(size(A)); G(2:4,2:4)=fliplr(C); G(1,1)=1,H=rem(A,3)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的值.a=A./B,b= floor(A./B),c=A./B ， d=B.\A ，e=A.*B, f=A*B ‟, g=2*ones(3,4)-eye(3,4)|A-B>B&B-8,4、要求在闭区间[0，2π ]上产生50个等距采样的一维数组。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB上机习题一
请按以下步骤完成上机实验：
1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内；
2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面；
3）如果习题是问答题，请将答案写在题目后；
4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后；
5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”；
6）上传该文件到FTP的相关目录。
1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。
2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。
3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。
4. 通过使用帮助确定内置变量ispc 的含义。
6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。
5. 创建double类型的变量，并进行计算
1）a=87，b=190，计算a+b、a-b、a*b；
2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致。
6. 计算如下表达式：
1）()sin 60
2）3e
3）3cos 4π⎛⎫ ⎪⎝⎭
4）2
7562323336
+⨯-⨯⨯
7. 设u=2，v=3，计算：
1）4log uv v
2）()
2
2u e v v u +-
3
8. 计算如下表达式：
1）()()3542i i -+
2）()sin 28i -
9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：
1）420<
2）420≤
3）420==
4）4~20=
5）''''b B <
10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b >
2）a b <
3）&&a b b c >>
4）a d ==
5）|a b c >
6）~~d
11. 判断下列变量名哪些是合法的MATLAB 变量：
1）fred
7）fred! 2）book_1
8）book-1 3）2ndplace
9）Second_Place 4）#1
10）No_1 5）vel_5
11）vel.5 6）tan
12）while 提示，可以使用isvarname 对上述变量名进行检验。
12. 编写脚本程序，命名为barbell.m ，完成如下计算：
如图1所示的杠铃，
每个球的半径是10cm ，两个求直接链接杆的长度是15cm ，杆的直径为1cm ，计算杠铃的体积和表面积；使用disp 函数显示计算结果。
（提示：disp 函数的用法参考教材的例2-31，面积计算时相交部分球冠可近似为圆，体积计算时可忽略重叠的球冠体积）
解答：
barbell.m （拷贝文件内容到如下方框中）：
运行屏幕截图：
13. 在高中课本中大家学过描述气体压强（p ）、温度（T ）、体积（V ）和气体摩尔数（n ）之间的理想气体状态方程：
pV nRT =
其中，R 为理想气体常数。
上述状态方程描述的是低压强和高温度时气体的特性。
在1873年，范德瓦尔斯对这个方程进行了修正，使其更好地描述气体在不同压强和温度条件下的状态，该方程被称为范德瓦尔斯方程：
()22n a p V nb nRT V ⎛⎫+-= ⎪⎝
⎭ 该方程新增了a 和b 两个变量，用来表示纯净气体的属性。
1）编写脚本程序，命名为state_equation.m ，利用如下数据：
分别使用理想气体状态方程和范德瓦尔斯方程计算水蒸气的温度，并使用disp 函数显示计算结果。
2）根据范德瓦尔斯方程及第1问的数据，编写getTemperature(p,V,n)函数实现温度计算，其中a、b、R作为常量。
调用该函数测试计算结果。
（提示，注意单位制，1 bar=105 Pa）
解答：
state_equation.m（拷贝到如下方框中）：
运行屏幕截图：
getTemperature.m（拷贝到如下方框中）：
运行屏幕截图：
14. 种群繁殖符合指数规律：
0rt
P P e
其中，P为目前的种群数，P0为原始种群数，r为持续增长率，t为时间。
现有100只兔子，持续增长率是每年90%（r=0.9），编写rabbit(n)函数用来计算n年后兔子的数量，并测试10年结果。
（提示，注意利用2.2.1节中的数学函数）
解答：
rabbit.m（拷贝到如下方框中）：
运行屏幕截图：
15. 一般来说，站得越高会看的越远，那么到底能看多远呢？主要取决于山的高度和地球的半径，如下图所示：
由于月球和地球的半径相差很大，所以能看到的地平线距离相差也很远，现已知地球半径是6378km，月球半径是1737km，编写脚本文件horizon.m计算0到8000m高的山上（珠穆朗玛峰海拔8850m）在地球和月球上看到地平线的距离分别是多少。
解答：
horizon.m（拷贝文件内容到如下方框中）：
运行屏幕截图：。