(完整版)matlab经典习题及解答
第1章 MATLAB 概论
1.1 与其他计算机语言相比较，MATLAB 语言突出的特点是什么？
MATLAB 具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB 系统由那些部分组成？
MATLAB 系统主要由开发环境、MATLAB 数学函数库、MATLAB 语言、图形功能和应用程序接口五个部分组成。
1.4 MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？
在MATLAB 操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view 菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。
1.5 如何启动M 文件编辑/调试器？
在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动。在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器。
1.6 存储在工作空间中的数组能编辑吗？如何操作？
存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。
1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？
命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。
1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？
当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。
1.9 在MATLAB 中有几种获得帮助的途径？ 在MATLAB 中有多种获得帮助的途径：
（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器； （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息；
（3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。
第2章 MATLAB 矩阵运算基础
2.1 在MATLAB 中如何建立矩阵??
?
???194375，并将其赋予变量a ？
>> 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 的行数。
在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如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+b ans =
7 7 7 9 14 13 15 12 14
2.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 ’ ans =
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 计算????
??=572396a 与??
?
???=864142b 的数组乘积。
>> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8]; >> a.*b ans =
12 36 3 8 42 40
2.8 “左除”与“右除”有什么区别？
在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。
2.9 对于B AX =，如果??????????=753467294A ，???
?
?
?????=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7]; >> B=[37 26 28]’; >> X=A\B X =
4.0427 1.3318
2.10 已知：????
??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >> a=[1 2 3;4 5 6;7 8 9]; >> a.^2 ans =
1 4 9 16 25 36 49 64 81 >> a^
2 ans =
30 36 42 66 81 96 102 126 150
2.11 ????
??-=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
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
2.12 []7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。
相当于a=[1 1 0 1 1]。
2.13 在sin(x )运算中，x 是角度还是弧度？
在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。
2.14 角度[]604530=x ，求x 的正弦、余弦、正切和余切。 >> x=[30 45 60]; >> x1=x/180*pi; >> sin(x1) ans =
0.5000 0.7071 0.8660 >> cos(x1) ans =
0.8660 0.7071 0.5000 >> tan(x1) ans =
0.5774 1.0000 1.7321 >> cot(x1) ans =
1.7321 1.0000 0.5774
2.15 用四舍五入的方法将数组[2.4568 6.3982
3.9375 8.5042]取整。 >> b=[2.4568 6.3982 3.9375 8.5042]; >> round(b) ans =
2 6 4 9
2.16 矩阵????
??????=728365219a ，分别对a 进行特征值分解、奇异值分解、LU 分解、QR 分解及Chollesky 分解。 >> [v,d]=eig(a,b) v =
-0.4330 -0.2543 -0.1744 -0.5657 0.9660 -0.6091 -0.7018 0.0472 0.7736 d =
13.5482 0 0 0 4.8303 0 0 0 3.6216 >> a=[9 1 2;5 6 3;8 2 7]; >> [u,s,v]=svd(a)
-0.5601 0.5320 -0.6350 -0.4762 -0.8340 -0.2788 -0.6779 0.1462 0.7204 s =
15.5234 0 0 0 4.5648 0 0 0 3.3446 v =
-0.8275 0.3917 -0.4023 -0.3075 -0.9156 -0.2592 -0.4699 -0.0907 0.8781 >> [l,u]=lu(a) l =
1.0000 0 0 0.5556 1.0000 0 0.8889 0.2041 1.0000 u =
9.0000 1.0000 2.0000 0 5.4444 1.8889 0 0 4.8367 >> [q,r]=qr(a) q =
-0.6903 0.3969 -0.6050 -0.3835 -0.9097 -0.1592 -0.6136 0.1221 0.7801 r =
-13.0384 -4.2183 -6.8260 0 -4.8172 -1.0807 0 0 3.7733 >> c=chol(a) c =
3.0000 0.3333 0.6667 0 2.4267 1.1447 0 0 2.2903
2.17 将矩阵4257a
??=??
??
、??????=3817b 和???
???=2695c 组合成两个新矩阵： （1）组合成一个4
3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为
按列顺序排列的c 矩阵元素，即
??
???
????
???237912685
574
（2）按照a 、b 、c 的列顺序组合成一个行矢量，即
>> a=[4 2;5 7]; >> b=[7 1;8 3]; >> c=[5 9;6 2]; % (1)
>> d=[a(:) b(:) c(:)] d =
4 7
5 5 8
6 2 1 9
7 3 2 % (2)
>> e=[a(:);b(:);c(:)]' e =
4 5 2 7 7 8 1 3 5 6 9 2 或利用（1）中产生的d >> e=reshape(d,1,12) ans =
4 5 2 7 7 8 1 3 5 6 9 2
第3章 数值计算基础
3.2 求解多项式x
3
-7x 2+2x +40=0的根。
>> r=[1 -7 2 40]; >> p=roots(r); -0.2151 0.4459 0.7949 0.2707
3.4 计算多项式乘法(x 2
+2x +2)(x 2
+5x +4)。 >> c=conv([1 2 2],[1 5 4]) c =
1 7 16 18 8
3.5 计算多项式除法(3x 3
+13x 2
+6x +8)/(x +4)。 >> d=deconv([3 13 6 8],[1 4]) d =
3 1 2
3.7 计算多项式95141242
3
4
++--x x x x 的微分和积分。 >> p=[4 –12 –14 5]; >> pder=polyder(p); >> pders=poly2sym(pder) >> pint=polyint(p);
pders =
12*x^2-24*x-14 pints =
x^4-4*x^3-7*x^2+5*x
3.8 解方程组????
??????=??????????66136221143092x 。 >> a=[2 9 0;3 4 11;2 2 6]; >> b=[13 6 6]'; >> x=a\b x = 7.4000 -0.2000 -1.4000
3.11 矩阵????
??????-=943457624a ，计算a 的行列式和逆矩阵。 >> a=[4 2 -6;7 5 4 ;3 4 9]; >> ad=det(a) >> ai=inv(a) ad = -64 ai =
-0.4531 0.6562 -0.5937 0.7969 -0.8437 0.9062 -0.2031 0.1562 -0.0937
3.12 y =sin(x )，x 从0到2，
x =0.02，求y 的最大值、最小值、均值和标准差。
>> x=0:0.02*pi:2*pi; >> y=sin(x); >> ymax=max(y) >> ymin=min(y) >> ymean=mean(y) >> ystd=std(y) ymax = 1 ymin = -1 ymean = 2.2995e-017 ystd = 0.7071
3.13 []54321=x ，[]108642=y ，计算x 的协方差、y 的协方差、x 与y 的互协方差。
>> x=[1 2 3 4 5]; >> y=[2 4 6 8 10]; >> cx=cov(x) >> cy=cov(y) >> cxy=cov(x,y) cx = 2.5000 cy = 10 cxy =
2.5000 5.0000 5.0000 10.0000
3.15 有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi ，用三次样条法进行插值。 >> x0=0:pi/5:4*pi; >> y0=sin(x0).*exp(-x0/10); >> x=0:pi/20:4*pi; >> y=spline(x0,y0,x); >> plot(x0,y0,'or',x,y,'b')
第4章 符号数学基础
4.1 创建符号变量有几种方法？
MATLAB 提供了两种创建符号变量和表达式的函数：sym 和syms 。
sym 用于创建一个符号变量或表达式，用法如x=sym(‘x’) 及 f=sym(‘x+y+z’)，syms 用于创建多个符号变量，用法如syms x y z 。
f=sym(‘x+y+z’) 相当于 syms x y z f= x+y+z
4.2 下面三种表示方法有什么不同的含义？ （1）f=3*x^2+5*x+2 （2）f='3*x^2+5*x+2' （3）x=sym('x')
（1）f=3*x^2+5*x+2
表示在给定x 时，将3*x^2+5*x+2的数值运算结果赋值给变量f ，如果没有给定x 则指示错误信息。 （2）f='3*x^2+5*x+2'
表示将字符串'3*x^2+5*x+2'赋值给字符变量f ，没有任何计算含义，因此也不对字符串中的内容做任何分析。 （3）x=sym('x') f=3*x^2+5*x+2
表示x 是一个符号变量，因此算式f=3*x^2+5*x+2就具有了符号函数的意义，f 也自然成为符号变量了。
4.3 用符号函数法求解方程a t 2
+b*t +c=0。 >> r=solve('a*t^2+b*t+c=0','t') r =
[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] 4.5 求矩阵??
?
?
??=22211211a a a a A 的行列式值、逆和特征根。 >> syms a11 a12 a21 a22; >> A=[a11,a12;a21,a22]
>> AD=det(A) % 行列式 >> AI=inv(A) % 逆 >> AE=eig(A) % 特征值 A = [ a11, a12] [ a21, a22] AD =
a11*a22-a12*a21 AI =
[ -a22/(-a11*a22+a12*a21), a12/(-a11*a22+a12*a21)] [ a21/(-a11*a22+a12*a21), -a11/(-a11*a22+a12*a21)] AE =
[ 1/2*a11+1/2*a22+1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)] [ 1/2*a11+1/2*a22-1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)]
4.6 因式分解：65552
3
4
-++-x x x x >> syms x;
>> f=x^4-5*x^3+5*x^2+5*x-6; >> factor(f) ans =
(x-1)*(x-2)*(x-3)*(x+1)
4.7 ?
???
??
??
=)sin()log(12
x x e x x a
f ax ，用符号微分求df/dx 。
>> syms a x;
>> f=[a, x^2, 1/x; exp(a*x), log(x), sin(x)];
df =
[ 0, 2*x, -1/x^2] [ a*exp(a*x), 1/x, cos(x)]
第5章 基本图形处理功能
5.1 绘制曲线13++=x x y ，x 的取值范围为[-5,5]。 >> x=-5:0.2:5; >> y=x.^3+x+1; >> plot(x,y)
5.2 有一组测量数据满足-at
e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况
下的曲线。
>> t=0:0.5:10; >> y1=exp(-0.1*t); >> y2=exp(-0.2*t); >> y3=exp(-0.5*t);
>> plot(t,y1,'-ob',t,y2,':*r',t,y3,'-.^g')
5.3 在5.1题结果图中添加标题-at
e =y ，并用箭头线标识出各曲线a 的取值。
>> title('\ity\rm=e^{-\itat}','FontSize',12)
>> text(t(6),y1(6),'\leftarrow\ita\rm=0.1','FontSize',11) >> text(t(6),y2(6),'\leftarrow\ita\rm=0.2','FontSize',11) >> text(t(6),y3(6),'\leftarrow\ita\rm=0.5','FontSize',11)
5.4 在5.1题结果图中添加标题-at
e y 和图例框。
>> title('\ity\rm=e^{-\itat}','FontSize',12) >> legend('a=0.1','a=0.2','a=0.5')
5.5表中列出了4个观测点的6次测量数据，将数据绘制成为分组形式和堆叠形式的条形图。
第1次 第2次 第3次 第4次 第5次 第6次 观测点1 3 6 7 4 2 8 观测点2 6 7 3 2 4 7 观测点3 9 7 2 5 8 4 观测点4 6
4
3
2
7
4
>> y=[3 6 9 6;6 7 7 4;7 3 2 3;4 2 5 2;2 4 8 7;8 7 4 4];
>> bar(y)
>> bar(y ,’stack’)
5.6 x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。 >> x=[66 49 71 56 38]; >> L=[0 0 0 0 1]; >> pie(x,L)
5.7 2
2
y x
xe
z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网
线图、表面图和带渲染效果的表面图。
>> [x,y]=meshgrid([-2:.2:2]); >> z=x.*exp(-x.^2-y.^2); >> mesh(x,y,z)
>> subplot(2,2,1), plot3(x,y,z) >> title('plot3 (x,y,z)')
>> title('mesh (x,y,z)')
>> subplot(2,2,3), surf(x,y,z)
>> title('surf (x,y,z)')
>> subplot(2,2,4), surf(x,y,z), shading interp
>> title('surf (x,y,z), shading interp')
5.8 绘制peaks函数的表面图，用colormap函数改变预置的色图，观察色彩的分布情况。 >> surf(peaks(30));
>> colormap(hot)
>> colormap(cool)
>> colormap(lines)
5.9 用sphere函数产生球表面坐标，绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。 >> [x,y,z]=sphere(30);
>> mesh(x,y,z)
>> mesh(x,y,z),hidden off
>> surf(x,y,z)
>> z(18:30,1:5)=NaN*ones(13,5); >> surf(x,y,z)
显示。
axis square
axis off
第6章 M文件编程
6.1 命令文件与函数文件的主要区别是什么？
（1）命令文件是一系列命令的组合，函数文件的第一行必须用function说明；
（2）命令文件没有输入参数，也不用返回参数，函数文件可以接受输入参数，也可以返回参数；
（3）命令文件处理的变量为工作空间变量，函数文件处理的变量为函数内部的局部变量，也可以处理全局变量。
6.2 如何定义全局变量？
用关键字global可以把一个变量定义为全局变量，在M文件中定义全局变量时，如果在当前工作空间已经存在了相同的变量，系统将会给出警告，说明由于将该变量定义为全局变量，可能会使变量的值发生改变，为避免发生这种情况，应该在使用变量前先将其定义为全局变量。
6.4 if语句有几种表现形式？
（1）if 表达式
语句体
end
（2）if 表达式
语句体1
else
语句体2
end
（3）if 表达式1
语句体1
elseif 表达式2
语句体2
end
（4）if 表达式1
语句体1
elseif 表达式2
语句体2
else
end
6.5 说明break语句和return语句的用法。
break语句用于终止for循环和while循环的执行。如果遇到break语句，则退出循环体，继续执行循环体外的下一行语句。在嵌套循环中，break只存在于最内层的循环中。
return语句用于终止当前的命令序列，并返回到调用的函数或键盘，也用于终止keyboad方式，在MATLAB中，被调用的函数运行结束后会自动返回调用函数，使用return语句时将return插入被调用函数的某一位置，根据某种条件迫使被调用函数提前结束并返回调用函数。
6.7 编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。
function c=collatz(n)
% collatz
% Classic “3n+1” Ploblem from number theory
c=n;
while n>1
if rem(n,2)==0
n=n/2;
else
n=3*n+1;
end
c=[c n];
end
6.8 有一组学生的考试成绩（见表），根据规定，成绩在100分时为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时为良好，成绩在60~79分之间为及格，成绩在60分以下时为不及格，编制一个根据成绩划分等级的程序。
Name=['王','张','刘','李','陈','杨','于','黄','郭','赵'];
Marks=[72,83,56,94,100,88,96,68,54,65];
% 划分区域：满分(100)，优秀(90-99)，良好(80-89)，及格(60-79)，不及格(<60)。
n=length(Marks);
for i=1:n
a{i}=89+i;
b{i}=79+i;
c{i}=69+i;
d{i}=59+i;
end;
c=[d,c];
% 根据学生的分数，求出相应的等级。
for i=1:n
switch Marks(i)
case 100 %得分为100时
Rank(i,:)=' 满分';
case a %得分在90~99之间
Rank(i,:)=' 优秀';
case b %得分在80~89之间
case c %得分在60~79之间
Rank(i,:)=' 及格';
otherwise %得分低于60。
Rank(i,:)='不及格';
end
end
% 将学生姓名，得分，级等信息打印出来。
disp(' ')
disp(['学生姓名 ',' 得分 ',' 等级']);
disp('--------------------------')
for i=1:10;
disp([' ',Name(i),' ',num2str(Marks(i)),' ',Rank(i,:)]);
end
学生姓名得分等级
-----------------------------------------
王 72 及格
张 83 良好
刘 56 不及格
李 94 优秀
陈 100 满分
杨 88 良好
于 96 优秀
黄 68 及格
郭 54 不及格
赵 65 及格
6.10 keyboard命令的作用是什么？当在函数中使用keyboard命令时，是否可以在工作空间浏览器中观察和修改函数中的局部变量？退出keyboard状态的命令是什么？编一个函数文件，在其中设置keyboard，观察变量情况。
当运行中的程序遇到keyboard命令时会暂时停止运行，将控制权交给键盘，这时命令窗口的提示符变为K>>，可以在命令窗口或工作空间浏览器中观察和修改变量内容，如果keyboard命令设置在函数中，则可以观察和修改函数中的局部变量，检查或修改完成后，键入return命令，退出keyboard状态，控制权交还给系统，程序将继续运行。
书上第四章习题
4.1
function [X1,X2]=matlab401(a,b,c)
d=b^2-4*a*c;
if d>0
X1=(-b+sqrt(d))/(2*a);
X2=(-b-sqrt(d))/(2*a);
elseif d==0
X1=(-b)/(2*a);
X2=X1;
else
r=(-b)/(2*a);
X1=r+k*i;
X2=r-k*i;
end
disp([X1,X2])
测试数据
(1)
>>matlab401(1,2,1) %两个相等的实根
(2)
>>matlab401(1,-1,-2) %两个不等的实根
(3)
>>matlab401(1,4,6.25) %两个不等的虚根
4.2
function fenS=matlab402(n)
fenZ=zeros(1,n);
fenM=ones(1,n);
fenZ(1)=2;
fenZ(2)=3;
fenM(1)=1;
fenM(2)=2;
fenSum=fenZ(1)/fenM(1)+fenZ(2)/fenM(2); for i=3:n
fenZ(i)=fenZ(i-1)+fenZ(i-2);
fenM(i)=fenM(i-1)+fenM(i-2);
fenS(i)=fenZ(i)/fenM(i);
fenSum=fenSum+fenS(i);
end
disp(fenSum)
测试数据
(1)
>>matlab402(3) %前3项的和
(2)
>>matlab402(15) %前15项的和
4.3
clear;clc;close all;
funy=zeros(1,601);
i=1;
for x=-3:0.01:3
if (x>-1)&(x<=-3)
y=(-x.^2-4*x-3)/2;
funy(i)=y;
i=i+1;
elseif (x>=-1)&(x<1)
y=-x.^2+1;
i=i+1;
else
y=(-x.^2+4*x-3)/2;
funy(i)=y;
i=i+1;
end
end
x=-3:0.01:3;
plot(x,funy,'r*-')
运行
>>matlab403
4.4
clear;clc;close all;
A=[-4 -2 0 2 4;-3 -1 1 3 5];
B=abs(A);
C=find(B>3);%?ò3?A?D?????μ′óóú3μ??a???ù?úμ?????D=[];%′?·?A?D?????μ′óóú3μ??a??
for i=1:length(C)
D(length(D)+1)=A(C(i));
end
disp(C)
disp(D)
运行
>>matlab404
结果C=[1 9 10]’
D=[-4 4 5]
4.5
clear;clc;close all;
jc=1;i=1;
while jc<1000
i=i+1;
jc=jc*i;
end
n=i-1
运行
>>matlab405
结果n=6
4.6
function matlab406(m,n)
% ?ò3???m,n??·??§?ú2??ü±?3??3yμ?êy2￠ê?3?
%num=zeros(1,n-m);
matlab考试题及答案
%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (i
matlab入门经典范例
num1=[13]; den1=conv([1,1],[1,0]); G1=tf(num1,den1); num2=[5.096,13]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.098,1],[1]); den2=conv( ssys1,ssys2); G2=tf(num2,den2); figure(1) margin(G1); hold on margin(G2); num1=[13]; den1=conv([1,1],[1,0]); G1=tf(num1,den1); num2=[5.096,13]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.098,1],[1]); den2=conv( ssys1,ssys2); G2=tf(num2,den2); figure(1) margin(G1); hold on margin(G2); num=[4.56,10]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.114,1],[1]); den=conv( ssys1,ssys2); G=tf(num,den); figure(1) bode(G) num=[4.56,10]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.114,1],[1]); den=conv( ssys1,ssys2); G=tf(num,den); figure(1) msrgin(G); num1=[13]; den1=conv([1,1],[1,0]);
MATLAB经典数学建模教程
第 1 节Matlab 基本知识 一、Matlab 的主要功能 Matlab是一种功能非常强大的工程语言，诞生于20世纪70年代，1984年正式推向市场。2002年8月，Matlab6.5开始发布。是进行科学研究和产品开发必不可少的工具。 ●数值和符号计算 矩阵(数组)的四则运算（Matrix＋Laboratory）、数值差分、导数、积分、求解微分方程、微分方程的优化等 ●数字图像、数字信号处理 ●工程和科学绘图 ●控制系统设计 ●财务工程 ●建模、仿真功能 二、Matlab 的界面 1.命令窗口（Command Window）： Matlab各种操作命令都是由命令窗口开始，用户可以在命令窗口中输入Matlab命令，实现其相应的功能。此命令窗口主要包括文本的编辑区域和菜单栏（如：四则运算；“;”禁止显示变量的值；↑↓遍历以前的命令）。在命令窗口空白区域单击鼠标右键，打开快捷菜单，各项命令功能如下： Evaluate Selection :打开所选文本对应的表达式的值。 Open Selection :打开文本所对应的MatLab文件。 Cut :剪切编辑命令。 Paste :粘贴编辑命令。 2. M-文件编辑/调试（Editor/Debugger）窗口 Matlab Editor/Debugger窗口是一个集编辑与调试两种功能于一体的工具环境。 M-文件（函数文件） ●什么是M-文件：它是一种和Dos环境中的批处理文件相似的脚本文件，对于简单问题， 直接输入命令即可，但对于复杂的问题和需要反复使用的则需做成M－文件（Script File）。 ●创建M-文件的方法： Matlab命令窗的File/New/M-file。 在Matlab命令窗口运行edit。 ●M－文件的扩展名：*.m ●执行M-文件：F5 ●M文件的调试 选择Debug菜单，其各项命令功能如下： Step :逐步执行程序。 Step in :进入子程序中逐步执行调试程序。
matlab基础练习题带答案
Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。
利用MATLAB进行时域分析
自动控制原理与系统课程实验报告 实验题目：利用MATLAB进行时域分析 班级：机电1131班姓名：刘润学号：38号 一、实验目的及内容 时域分析法是一种直接在时间域中对系统进行分析的方法，具有直观、准确的优点，并且可以提供系统时间响应的全部信息。在此实验中，主要介绍时域法进行系统分析，包括一阶系统、二阶系统以及高阶系统,以及系统的性能指标。通过实验，能够快速掌握、并利用MATLAB及控制系统箱对各种复杂控制系统进行时域分析。 二、实验设备 三、实验原理 典型的二阶系统在不同的阻尼比的情况下，它们的阶跃响应输出特性的差异是很大的。若阻尼比过小，则系统的振荡加剧，超调量大幅度增加；若阻尼比过大，则系统的响应过慢，又大大增加了调整时间，下面通过此实验课题分析输出响应变化规律： 已知二阶振荡环节的传递函数为：G（s）=ωn*ωn/（s*s+2*ζ*ωn*s+ωn*ωn）， 其中ωn=0.4，ζ从0变化到2，求此系统的单位阶跃响应曲线，并分析当ζ发生变化时，二阶系统的响应有什么样的变化规律。
四、实验步骤编出程序如下图： 五、实验结果画出图表如下图：
六、结果分析 （1）当ξ=0（无阻尼）（零阻尼）时： 无阻尼时的阶跃响应为等幅振荡曲线。如图ξ=0曲线。 （2）当0<ξ<1（欠阻尼）时： 对应不同的ξ，可画出一系列阻尼振荡曲线，且ξ越小，振荡的最大振幅愈大。如图ξ=0.4曲线。 (3)当ξ=1（临界阻尼）时： 临界阻尼时的阶跃响应为单调上升曲线。如图ξ=1曲线。 （4）当ξ>1（过阻尼）时： 过阻尼时的阶跃响应也为单调上升曲线。不过其上升的斜率较临界阻尼更慢。如图ξ=1.6曲线 七、教师评语
matlab经典习题及解答
第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 页脚内容1
命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 在MATLAB中有多种获得帮助的途径： （1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器； （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。 第2章MATLAB矩阵运算基础 页脚内容2
matlab课程设计题目
课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图
形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。
matlab经典编辑习题集及解答
第1章 MATLAB 概论 1.1 与其他计算机语言相比较，MATLAB 语言突出的特点是什么？ MATLAB 具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB 系统由那些部分组成？ MATLAB 系统主要由开发环境、MATLAB 数学函数库、MATLAB 语言、图形功能和应用程序接口五个部分组成。 1.4 MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB 操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view 菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M 文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动。在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB 中有几种获得帮助的途径？ 在MATLAB 中有多种获得帮助的途径： （1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器； （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按T ab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 2.1 在MATLAB 中如何建立矩阵?? ? ???194375，并将其赋予变量a ？ >> 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 ）
MATLAB习题参考答案
第1章 MATLAB简介 1、MA TLAB的主要特点有：①语言简洁，编程效率高。②人机界面友善，交互性好。③绘图功能强大，便于数据可视化。④学科众多、领域广泛的MATLAB工具箱。⑤源程序的开放性。 MATLAB的典型应用领域有：①自动控制②汽车③电子④仪器仪表⑤生物医学⑥信号处理⑦通信等。 2、填空题 ⑴命令窗口、命令历史窗口、当前目录窗口 ⑵查阅、保存、编辑 ⑶清除图形窗、清除命令窗口中显示内容、清除MATLAB工作空间中保存的变量。 3、如果想查看某一变量具体内容或者对其修改操作，可以在工作空间中双击该变量名称，可以打开数组编辑器，在数组编辑器中可以查看变量的具体内容，也可以对其修改。如果想删除MATLAB内存中的变量，可以在工作空间中选中该变量，然后利用工作空间窗口的菜单命令或工具条中的快捷图标进行删除。 4、1+2+3+4+5+... （+ 后面可以直接跟...，也可以在+和...中加一个空格。） 6+7+8+9 1+2+3+4+5 ... （5后面必须跟一个空格，不能直接跟...，否则报错，这在预置一个+6+7+8+9 大数组时很重要。） 第2章矩阵与数值数组 1、填空题： ⑴非数、无穷大、机器零阈值，浮点数相对精度,eps= 2.2204e-016。 ⑵全下标、单下标。 2、阅读程序题： （本题主要考察数组的寻访、赋值和简单运算，提示：带；的语句不显示结果） ⑴ans = 2 3 7 Sa = 10 20 30 A = 1 20 5 30 9 10 4 6 8 10 ⑵ A = 1 3 5 7 9 11 13 15
2 4 6 8 10 12 14 16 ans = 1 5 9 13 2 6 10 14 3 7 11 15 4 8 12 16 A = 0 0 5 7 0 0 13 15 2 4 0 0 10 12 0 0 ⑶ ans = -1 -4 6 4 ans = 3 0 5 -2 ans = 3 6 9 12 ans = 3 6 9 12 3、A=magic(4)； L=A<10 L = 0 1 1 0 1 0 0 1 1 1 1 0 1 0 0 1 4、省略。 第3章字符串、元胞和构架数组 1. 直接创建法：S1=['Where there is life,' 'there is hope. '] %注意第2行要加入6个空格使其与第一 行字符数相等，否则报错。 S2=str2mat('Where there is life,','there is hope.') S3=strvcat('Where there is life,','there is hope.') %后两种方法则不用考虑两行 字符的数目 2. 填空题：A（2，3）；A{2，3} 3. 阅读程序题： ⑴subch = ABc12
matlab入门教程
MATLAB入门教程 1．MATLAB的基本知识 1-1、基本运算与函数 在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter键即可。例如： >> (5*2+1.3-0.8)*10/25 ans =4.2000 MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer）并显示其数值於萤幕上。 小提示： ">>"是MATLAB的提示符号（Prompt），但在PC中文视窗系统下，由於编码方式不同，此提示符号常会消失不见，但这并不会影响到MATLAB的运算结果。 我们也可将上述运算式的结果设定给另一个变数x： x = (5*2+1.3-0.8)*10^2/25 x = 42 此时MATLAB会直接显示x的值。由上例可知，MATLAB认识所有一般常用到的加（+）、减（-）、乘（*）、除（/）的数学运算符号，以及幂次运算（^）。 小提示：MATLAB将所有变数均存成double的形式，所以不需经过变数宣告（Variable declaration）。MATLAB同时也会自动进行记忆体的使用和回收，而不必像C语言,必须由使用者一一指定.这些功能使的MATLAB易学易用，使用者可专心致力於撰写程式，而不必被软体枝节问题所干扰。 若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可，如下例： y = sin(10)*exp(-0.3*4^2); 若要显示变数y的值，直接键入y即可： >>y y =-0.0045 在上例中，sin是正弦函数，exp是指数函数，这些都是MATLAB常用到的数学函数。 下表即为MATLAB常用的基本数学函数及三角函数： 小整理：MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle)
matlab2007教程
第 1 章基础准备及入门本章有三个目的：一是讲述MATLAB正常运行所必须具备的基础条件；二是简明地介 绍MATLAB及其操作桌面Desktop的基本使用方法；三是全面介绍MATLAB的帮助系统。 本章的前两节讲述：MATLAB的正确安装方法和MATLAB 环境的启动。因为指令窗是MATLAB最重要的操作界面，所以本章用第 1.3、1.4 两节以最简单通俗的叙述、算例讲述指令窗的基本操作方法和规则。这部分内容几乎对MATLAB各种版本都适用。第1.5到第1.8节专门介绍MATLAB最常用的另五个交互界面：历史指令窗、当前目录浏览器、工作空间浏览器、数组编辑器、M文件编辑器。鉴于实际应用中，帮助信息和求助技能的重要性。本章专设第1.9节专门叙述MATLAB的帮助体系和求助方法。 作者建议：不管读者此前是否使用过MATLAB，都不要忽略本章。 1.1 MATLAB的安装和工具包选择 MATLAB只有在适当的外部环境中才能正常运行。因此，恰当地配置外部系统是保证MATLAB运行良好的先决条件。MATLAB本身可适应于许多机种和系统，如 PC机和 Unix 工作站等。但本节只针对我国使用最广的PC机系统给予介绍。 对PC机用户来说，常常需要自己安装MATLAB。MATLAB R2007a（即旧编号MATLAB7.4）版要求Win2000或WinXP平台。下面介绍从光盘上安装MATLAB的方法。一般说来，当MATLAB光盘插入光驱后，会自启动“安装向导”。假如自启动没有实现，那么可以在<我的电脑>或<资源管理器>中双击setup.exe应用程序，使“安装向导”启动。安装过程中出现的所有界面都是标准的，用户只要按照屏幕提示操作，如输入用户名、单位名、口令等就行。 在安装MATLAB.R2007a时，会出现一个界面，该界面上有两个选项：Typical和Custom。假如你不熟悉MATLAB，或假如你机器的硬盘的自由空间远大于3G，或假如你需要用到光盘上MATLAB的所有功能及工具包，那么你就点选“Typical”。否则，点选“Custom”。 在点选“Custom”后，会引出如图1.1-1的界面。你可以根据需要，在“Select products to install”栏中勾选相应的组件。注意：MATLAB软件光盘总包含很多工具包，它们有的是通用的，有的则专业性很强。对一般用户来说，完全不必采取全部安装，而应根据需要有所选择。否则将占据很多硬盘空间。表1.1-1对各组件的描述供用户选择时参考。
MATLAB习题及参考答案经典.doc
习题： 1， 计算?? ????=572396a 与??? ???=864142b 的数组乘积。 2， 对于B AX =，如果???? ? ?????=753467294A ，??????????=282637B ，求解X 。 3， 已知：?? ?? ? ?????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。(应用sin,cos,tan.cot) 5， 将矩阵?? ?? ??=7524a 、??????=3817b 和??? ???=2695c 组合成两个新矩阵： （1）组合成一个4?3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ?? ??? ???? ???237 912685 574 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []296531877254 6， 将(x -6)(x -3)(x -8)展开为系数多项式的形式。(应用poly,polyvalm) 7， 求解多项式x 3-7x 2+2x +40的根。(应用roots) 8， 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。(应用poly,polyvalm) 9， 计算多项式9514124234++--x x x x 的微分和积分。(应用polyder,polyint ，poly2sym)
10， 解方程组???? ? ?????=??????????66136221143092x 。(应用x=a\b) 11， 求欠定方程组?? ? ???=???? ??5865394742x 的最小范数解。(应用pinv) 12， 矩阵???? ? ?????-=943457624a ，计算a 的行列式和逆矩阵。(应用det,inv) 13， y =sin(x )，x 从0到2π，?x =0.02π，求y 的最大值、最小值、均值和标准差。(应用max,min,mean,std) 14， 参照课件中例题的方法，计算表达式() 2 2 e 1053y x y x z ---=的梯度并绘图。(应用meshgrid, gradient, contour, hold on, quiver) 15， 用符号函数法求解方程a t 2+b*t +c=0。(应用solve) 16， 用符号计算验证三角等式：(应用syms,simple) 17， 求矩阵?? ? ? ??=2221 1211a a a a A 的行列式值、逆和特征根。(应用syms,det,inv,eig) 18， 因式分解：6555234-++-x x x x (应用syms, factor) 19， ? ??? ?? ?? =)sin()log(12 x x e x x a f ax ，用符号微分求df/dx 。(应用syms,diff) 20， 符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t 的变化范围为[0,2π]。(应用syms,ezplot) 21， 绘制曲线13++=x x y ，x 的取值范围为[-5,5]。(应用plot) 22， 有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线,在图中添加标题-at e =y ，并用箭头线标识出各曲线a 的取值，并添加标题-at e =y 和图例框。(应用plot,title,text,legend) 23 24， x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。
Matlab基础教程
1-1、基本运算与函数 在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter键即可。例如： >> (5*2+1.3-0.8)*10/25 ans =4.2000 MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer）并显示其数值於萤幕上。 小提示： ">>"是MATLAB的提示符号（Prompt），但在PC中文视窗系统下，由於编码方式不同，此提示符号常会消失不见，但这并不会影响到MATLAB的运算结果。 我们也可将上述运算式的结果设定给另一个变数x： x = (5*2+1.3-0.8)*10^2/25 x = 42 此时MATLAB会直接显示x的值。由上例可知，MATLAB认识所有一般常用到的加（+）、减（-）、乘（*）、除（/）的数学运算符号，以及幂次运算（^）。 小提示： MATLAB将所有变数均存成double的形式，所以不需经过变数宣告（Variable declaration）。MATLAB同时也会自动进行记忆体的使用和回收，而不必像C语言,必须由使用者一一指定.这些功能使的MATLAB易学易用，使用者可专心致力於撰写程式，而不必被软体枝节问题所干扰。 若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可，如下例： y = sin(10)*exp(-0.3*4^2);
若要显示变数y的值，直接键入y即可： >>y y =-0.0045 在上例中，sin是正弦函数，exp是指数函数，这些都是MATLAB常用到的数学函数。 下表即为MATLAB常用的基本数学函数及三角函数： 小整理：MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数 rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开
Matlab考题题整理带答案
MATLAB 考试试题(1) 产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏） 答案： a=10*rand(1,10)-5; b=sort(a,'descend') 1.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5] 2. 已知变量：A=’ilovematlab’；B=’matlab’, 请找出： （A）B在A中的位置。 （B）把B放在A后面，形成C=‘ilovematlabmatlab’ 3. 请修改下面的程序，让他们没有for循环语句！ A=[1 2 3; 4 5 6; 7 8 9]; [r c]=size(A); for i=1:1:r for j=1:1:c if (A(i,j)>8 | A(i,j)<2) A(i,j)=0; end end end 4. 请把变量A=[1 2 3; 4 5 6; 7 8 9]写到文件里(output.xls)，写完后文件看起来是这样的 1 2 3 4 5 6 7 8 9 5.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。 6.编写M文件，从Yahoo网站批量读取60000.SH至600005.SH在2008年9月份的每日收盘价（提示：使用字符串函数）。 7. 将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。 8.已知资产每日回报率为0.0025，标准差为0.0208，资产现在价值为0.8亿，求5％水平下资产的10天在险价值(Var)。 9.a=[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB中最简单的方法计算b，注意最简单哦。
matlab例题
五、某公司投资2000万元建成一条生产线。投产后，在时刻t 的追加成本和追加收益分别为3/225)(t t t G ++=（百万元/年），3/218)(t t H -=（百万元/年）。试确定该生产线在何时停产可获最大利润？最大利润是多少？ 提示：利用函数?=T G H t R 0t 20-d ))t (-)t (()(（百万元），由于H （t ）-G （t ）单调 下降，所以H （t ）=G （t ）时，R （t ）取得最大利润。 5.解：构造函数f(t)=H(t)-G(t)=13-t-3t 2/3=0 ; 令t 1/3=x,则f(t)=-t 3-3t 2+13 可得矩阵P=[-1,-3,0,13] 求最佳生产时间的源程序如下： p=[-1,-3,0,13]; x=roots(p); t=x.^3 运行结果如下： t = 3.6768 +21.4316i 3.6768 -21.4316i 4.6465 再分别将t 的三个值带入函数f(t)，比较大小后，得到最大利润与最佳生产时间。 求最大利润的程序代码如下： ① t=3.6768 +21.4316i; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 25.2583 ② t=3.6768 -21.4316i; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 25.2583 ③ t=4.6465; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 26.3208 比较以上三组数据，可知最佳生产时间t=4.6465年，可获得的最大利润 26.3208（百万元/年）。 clear; close; fplot('18-t^(2/3)',[0,20]);grid on;hold on; fplot('5+t+2*t^(2/3)',[0,20],'r');hold off; %发现t 约为4
matlab经典编程例题
以下各题均要求编程实现，并将程序贴在题目下方。 1．从键盘输入任意个正整数，以0结束，输出那些正整数中的素数。 clc;clear; zzs(1)=input('请输入正整数：');k=1; n=0;%素数个数 while zzs(k)~=0 flag=0;%是否是素数，是则为1 for yz=2:sqrt(zzs(k))%因子从2至此数平方根 if mod(zzs(k),yz)==0 flag=1;break;%非素数跳出循环 end end if flag==0&zzs(k)>1%忽略0和1的素数 n=n+1;sus(n)=zzs(k); end k=k+1; zzs(k)=input('请输入正整数：'); end disp(['你共输入了' num2str(k-1) '个正整数。它们是：']) disp(zzs(1:k-1))%不显示最后一个数0 if n==0 disp('这些数中没有素数！')%无素数时显示 else disp('其中的素数是：') disp(sus) end 2．若某数等于其所有因子（不含这个数本身）的和，则称其为完全数。编程求10000以内所有的完全数。 clc;clear;
wq=[];%完全数赋空数组 for ii=2:10000 yz=[];%ii的因子赋空数组 for jj=2:ii/2 %从2到ii/2考察是否为ii的因子 if mod(ii,jj)==0 yz=[yz jj];%因子数组扩展，加上jj end end if ii==sum(yz)+1 wq=[wq ii];%完全数数组扩展，加上ii end end disp(['10000以内的完全数为：' num2str(wq)])%输出 3．下列这组数据是美国1900—2000年人口的近似值（单位：百万）。 （1）若. 2c + = y+ 与试编写程序计算出上式中的a、b、c; 的经验公式为 t at bt y （2）若.bt 的经验公式为 y= 与试编写程序计算出上式中的a、b; y ae t （3）在一个坐标系下，画出数表中的散点图（红色五角星），c + =2中 ax bx y+拟合曲线图(蓝色实心线),以及.bt y=(黑色点划线)。 ae （4）图形标注要求:无网格线,横标注“时间t”，纵标注“人口数(百万)”，图形标题“美国1900—2000年的人口数据”。 （5）程序中要有注释，将你的程序和作好的图粘贴到这里。 clf;clc;clear %清除图形窗、屏幕、工作空间 t=1900:10:2000; y=[76 92 106 123 132 151 179 203 227 250 281]; p1=polyfit(t,y,2);%二次多项式拟合
matlab教程详解 (5)
第五章符号计算 符号计算的特点：一，运算以推理解析的方式进行，因此不受计算误差积累问题困扰；二，符号计算，或给出完全正确的封闭解，或给出任意精度的数值解（当封闭解不存在时）；三，符号计算指令的调用比较简单，经典教科书公式相近；四，计算所需时间较长，有时难以忍受。 在MATLAB中，符号计算虽以数值计算的补充身份出现，但涉及符号计算的指令使用、运算符操作、计算结果可视化、程序编制以及在线帮助系统都是十分完整、便捷的。 MATLAB的升级和符号计算内核Maple的升级，决定着符号计算工具包的升级。但从用户使用角度看，这些升级所引起的变化相当细微。即使这样，本章还是及时作了相应的更新和说明。如MATLAB 6.5+ 版开始启用Maple VIII的计算引擎，从而克服了Maple V计算“广义Fourier变换”时的错误（详见第5.4.1节）。 5.1符号对象和符号表达式 5.1.1符号对象的生成和使用 【例5.1.1-1】符号常数形成中的差异 a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] % <1> a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) % <2> a3=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)],'e') % <3> a4=sym('[1/3,pi/7,sqrt(5),pi+sqrt(5)]') % <4> a24=a2-a4 a1 = 0.3333 0.4488 2.2361 5.3777 a2 = [ 1/3, pi/7, sqrt(5), 6054707603575008*2^(-50)] a3 = [ 1/3-eps/12, pi/7-13*eps/165, sqrt(5)+137*eps/280, 6054707603575008*2^(-50)] a4 = [ 1/3, pi/7, sqrt(5), pi+sqrt(5)] a24 = [ 0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] 【例5.1.1-2】演示：几种输入下产生矩阵的异同。 a1=sym([1/3,0.2+sqrt(2),pi]) % <1> a2=sym('[1/3,0.2+sqrt(2),pi]') % <2> a3=sym('[1/3 0.2+sqrt(2) pi]') % <3> a1_a2=a1-a2 % a1 = [ 1/3, 7269771597999872*2^(-52), pi] a2 = [ 1/3, 0.2+sqrt(2), pi] a3 = [ 1/3, 0.2+sqrt(2), pi] a1_a2 =