MATLAB课后习题
matlab课后习题答案（附图）
matlab课后习题答案（附图）习题2.1画出下列常见曲线的图形y （1）⽴⽅抛物线3x命令：syms x y;ezplot('x.^(1/3)')(2)⾼斯曲线y=e^(-X^2);命令：clearsyms x y;ezplot('exp(-x*x)')(3)笛卡尔曲线命令：>> clear>> syms x y;>> a=1;>> ezplot(x^3+y^3-3*a*x*y)(4)蔓叶线命令：>> clear>> syms x y;>> a=1ezplot(y^2-(x^3)/(a-x))(5)摆线:()()tsin-=,=-by1命令：>> clear>> t=0:0.1:2*pi;>> x=t-sin(t);>>y=2*(1-cos(t)); >> plot(x,y)7螺旋线命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t;>>plot3(x,y,z)(8)阿基⽶德螺线命令：clear>> theta=0:0.1:2*pi;>> rho1=(theta);>> subplot(1,2,1),polar(theta,rho1)(9) 对数螺线命令：cleartheta=0:0.1:2*pi;rho1=exp(theta);subplot(1,2,1),polar(theta,rho1)(12)⼼形线命令：>> clear >> theta=0:0.1:2*pi; >> rho1=1+cos(theta); >> subplot(1,2,1),polar(theta,rho1)练习2.21. 求出下列极限值（1）nnn n3→命令：>>syms n>>limit((n^3+3^n)^(1/n)) ans = 3（2）)121(lim n n n n ++-+∞→命令：>>syms n>>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf) ans = 0（3）x x x 2cot lim 0→命令：syms x ；>> limit(x*cot(2*x),x,0) ans = 1/2 （4）)(coslimcm xx ∞→命令：syms x m ； limit((cos(m/x))^x,x,inf) ans = 1（5）)111(lim 1--→exx x命令：syms x>> limit(1/x-1/(exp(x)-1),x,1) ans =(exp(1)-2)/(exp(1)-1) （6）)(2lim x x xx -+∞>> limit((x^2+x)^(1/2)-x,x,inf)ans = 1/2练习2.41. 求下列不定积分，并⽤diff 验证：（1）+x dxcos 1>>Clear >> syms x y >> y=1/(1+cos(x)); >> f=int(y,x) f =tan(1/2*x) >> y=tan(1/2*x); >> yx=diff(y ,x); >> y1=simple(yx) y1 =1/2+1/2*tan(1/2*x)^2 (2)+exdx1clear syms x yy=1/(1+exp(x));f=int(y,x) f =-log(1+exp(x))+log(exp(x)) syms x yy=-log(1+exp(x))+log(exp(x)); yx=diff(y,x); y1=simple(yx) y1 = 1/(1+exp(x)) (3)dx x x ?sin 2syms x yy=x*sin(x)^2; >> f=int(y,x) f =x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2 clearsyms x y y=x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2; yx=diff(y,x); >> y1=simple(yx) y1 = x*sin(x)^2 (4)xdx ?sec3syms x y y=sec(x)^3;f=int(y,x) f =1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)) clear syms x yy=1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x)); yx=diff(y,x); y1=simple(yx) y1 =1/cos(x)^32. 求下列积分的数值解 1）dx x-10clearsyms xy=int(x^(-x),x,0,1) y =int(x^(-x),x = 0 .. 1) vpa(y,10) ans =1.291285997 2）xdx e x cos3202?πclearsyms xy=int(exp(2*x)*cos(x)^3,x, clear syms xy=int((1/(2*pi)^(1/2))*exp(-x^2/2),x,0,1) y =7186705221432913/36028797018963968*erf(1/2*2^(1/2))*2^(1/2)*pi^(1/0,2*pi) y =22/65*exp(pi)^4-22/65vpa(ans,10)(3)dx xe21221-π>> clear >> syms x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),0,1); >> vpa(y,14) ans =.341344746068552（4）>> clear >> syms x>> y=int(x*log(x^4)*asin(1/x^2),1,3); Warning: Explicit integral could not be found. > In sym.int at 58 >> vpa(y,14) ans = 2.45977212823752（5） >> clear >> syms x1判断下列级数的收敛性，若收敛，求出其收敛值。
MATLAB课后习题
MATLAB课后习题第⼀部分 MATLAB 运算基础1. 先求下列表达式的值，然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。
(1) 0122sin 851z e=+ (2) 221ln(1)2z x x =++，其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ?≤=-≤，其中t =0:0.5:2.52. 已知：1234413134787,2033657327A B --==-求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B-???(1) 求它们的乘积C 。
(2) 将矩阵C 的右下⾓3×2⼦矩阵赋给D 。
(3) 查看MATLAB ⼯作空间的使⽤情况。
4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建⽴⼀个字符串向量，删除其中的⼤写字母。
第⼆部分 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证22E R RS A OS +??=?。
2. 产⽣5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其⾏列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
为什么？ 3. 建⽴⼀个5×5矩阵，求它的⾏列式值、迹、秩和范数。
MATLAB课后习题集
第2章 MATLAB 概论1、 与其他计算机语言相比较，MA TLAB 语言突出的特点是什么？2、 MATLAB 系统由那些部分组成？3、 安装MATLAB 时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？4、 MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？5、 如何启动M 文件编辑/调试器？6、 存储在工作空间中的数组能编辑吗？如何操作？7、 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？8、 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？9、 在MA TLAB 中有几种获得帮助的途径？第3章 MATLAB 数值运算3.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？3.2 有几种建立矩阵的方法？各有什么优点？3.3 在进行算术运算时，数组运算和矩阵运算各有什么要求？3.4 数组运算和矩阵运算的运算符有什么区别？3.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
3.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
3.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
3.8 “左除”与“右除”有什么区别？3.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
3.11 ⎥⎦⎤⎢⎣⎡-=463521a ，⎥⎦⎤⎢⎣⎡-=263478b ，观察a 与b 之间的六种关系运算的结果。
3.12 []7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。
matlab课后习题及答案
第一章 5题已知a=4.96，b=8.11，计算)ln(b a eba +-的值。
解：clear clc a=4.96; b=8.11;exp(a-b)/log(a+b) ans =0.0167 6题已知三角形的三边a=9.6,b=13.7, c=19.4,求三角形的面积。
提示：利用海伦公式area =))()((c s b s a s s ---计算，其中S=(A+B+C)/2. 解：clear clc a=9.6; b=13.7; c=19.4; s=(a+b+c)/2area=sqrt(s*(s-a)*(s-b)*(s-c)) s =21.3500 第二章 8题已知S=1+2+2^2+2^3+……+2^63,求S 的值 解：clear clc S=0;for i=0:1:63 S=S+2^i; end S S =1.8447e+019 9题分别用for 和while 循环结构编写程序，计算∑=-1001n 1n 2）（的值。
解：clear clc s=0;for n=1:100 s=s+(2*n-1); end s s =10000 clear clc n=1; s=0;while n<=100 s=s+(2*n-1); n=n+1; end s s =10000 第三章 2题在同一坐标下绘制函数x ,,2x-,2x xsin(x)在()∏∈,0x 的曲线。
解：clear clcx=0:0.2:pi; y1=x; y2=x.^2; y3=-(x.^2); y4=x.*sin(x);plot(x,y1,'-' ,x,y2,'-' ,x,y3,'-' ,x,y4,'-')0.511.522.53-10-8-6-4-202468109题用不同的线型和颜色在同一坐标内绘制曲线y1=2ex5.0 、y2=sin(2∏x )的图形。
matlab课后习题答案
习题二1.如何理解“矩阵就是MATＬAB最基本得数据对象”？答:因为向量可以瞧成就是仅有一行或一列得矩阵,单个数据(标量)可以瞧成就是仅含一个元素得矩阵,故向量与单个数据都可以作为矩阵得特例来处理。
ﻩ因此,矩阵就是MATLＡB最基本、最重要得数据对象。
2.设A与Ｂ就是两个同维同大小得矩阵,问:(1)A＊B与A。
＊B得值就是否相等?答:不相等。
(2)A。
/B与Ｂ。
＼A得值就是否相等?答:相等。
(3)A／B与Ｂ\A得值就是否相等？答:不相等、(4)A／Ｂ与B＼A所代表得数学含义就是什么？答:A/Ｂ等效于B得逆右乘A矩阵,即A＊inv(B),而B＼Ａ等效于Ｂ矩阵得逆左乘A矩阵,即inv(Ｂ)＊A。
3.写出完成下列操作得命令。
(1)将矩阵A第2～５行中第１,３, 5列元素赋给矩阵B。
答:B=A(2:５,1:２:５); 或B=A(２:5,［１ 3 5］)(2)删除矩阵A得第７号元素、答:A(7)=[］(3)将矩阵A得每个元素值加3０、答:Ａ＝Ａ+30;(4)求矩阵A得大小与维数、答:size(A);ｎdiｍｓ(A);(5)将向量ｔ得0元素用机器零来代替。
答:t(finｄ(t==0))＝eps;(6)将含有12个元素得向量x转换成矩阵、答:reｓhapｅ(x,３,4);(7)求一个字符串得ASCIＩ码。
答:aｂｓ(‘１23’); 或ｄouble(‘12３');(8)求一个ASCＩI码所对应得字符、答:ｃｈａr(49);4.下列命令执行后,L1、L２、L３、L4得值分别就是多少?A=1:9;B=1０-A;。
.、L1=Ａ==B;Ｌ2=A<=5;L3=A〉3&Ａ〈７;L4＝ｆｉnd(Ａ>3&A<7);答:Ｌ１得值为[0, 0, 0, 0, １, 0, 0,0, 0］L２得值为［1, 1, １, １, 1,０, 0, 0, ０]Ｌ3得值为［０, 0, ０, １,１, １, 0, 0, 0］L4得值为［4,5, 6]5.已知完成下列操作:(1)取出A得前3行构成矩阵Ｂ,前两列构成矩阵C,右下角子矩阵构成矩阵Ｄ,B与C得乘积构成矩阵E。
MATLAB课后习题
MATLAB课后习题第1章MATLAB操作基础1 MA TLAB的主要功能包括数值计算和符号计算功能，绘图功能，编程语⾔以及应⽤⼯具箱。
2 MA TLAB的优点在于可以在窗⼝直接输⼊运算命令，并的到计算结果；有许多现成的函数可以直接调⽤；可以直接绘制图形对结果更直接的显⽰。
3 MA TLAB中help命令在窗⼝中输⼊后将会显⽰当前帮助系统中的所包含的所有项⽬，即搜索路径中所有的⽬录项⽬，所以可以搜索到。
4 主要有控制系统⼯具箱，信号处理⼯具箱，神经⽹络⼯具箱，最优化⼯具箱，⾦融⼯具箱，统计学⼯具箱，功能性⼯具箱。
5 inv 为求⽅阵的逆矩阵，plot为对所给函数作图，max为取最⼤值，round为四舍五⼊到最临近的整数。
第2章MATLAB矩阵及其运算2-1 合法的变量有global My_exp AbCd AB_C_2-2 直接输⼊法；利⽤M⽂件建⽴；利⽤冒号表达式建⽴⼀个向量；建⽴⼤矩阵。
2-3 A*B要求A阵的列数必须与B阵的⾏数相等，A.*B表⽰对应元素相乘，矩阵运算左除与右除值⼀般不同，点除也有类似的规律，但对于标量则都相同。
实验指导1.(1) w=sqrt(2)*(1+0.34245*10^(-6))w = 1.4142(2) a=3.5;b=5;c=-9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)x =0.9829(3) a=3.32;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)y = -128.4271(4) t=[2,1-3i;5,-0.65];z=0.5*exp(2*t)*log(t+sqrt(1+t.*t))z =1.0e+004 * 0.0048 + 0.0002i 0.0048 - 0.0034i1.58992.0090 - 1.3580i2.(1) A=[-1 5 -4;-0 7 8;3 61 7]A =-1 5 -40 7 8B=[8 3 -1;2 5 3;-3 2 0]B =8 3 -12 5 3-3 2 0A+6*Bans =47 23 -1012 37 26-15 73 7I=[1 0 0;0 1 0;0 0 1]I =1 0 00 1 00 0 1A*A-B+Ians =-18 -217 1722 533 10921 867 526(2)A*Bans =14 14 16-10 51 21125 328 180A.*Bans =-8 15 40 35 24-9 122 0B*Aans =-11 0 -157 228 53(3)[A,B]ans =-1 5 -4 8 3 -10 7 8 2 5 33 61 7 -3 2 0[A([1,3],:);B^2]ans =-1 5 -43 61 773 37 117 37 13-20 1 93 (1)A=[23 10 -0.778 0;41 -45 65 5;32 5 0 32;6 -9.54 54 3.14]A =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.00006.0000 -9.5400 54.0000 3.1400B=A(1:3,:)B =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000 C=A(:,1:2)C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D=A(2:4,3:4)D =65.0000 5.00000 32.000054.0000 3.1400E =1.0e+003 *0.9141 -0.22391.20802.71231.1330 -0.2103(2)Eans =0 10 00 1E&Dans =1 10 11 1E|Dans =1 11 11 1~Dans =0 01 00 0~Eans =0 00 00 04， H=hilb(5); P=pascal(5); Hh=det(H)Hh = 3.7493e-012Hp=det(P)Hp = 1 Th=cond(H)Th = 4.7661e+005Tp=cond(P)Tp = 8.5175e+003条件数越趋近于1，矩阵的性能越好，所以帕斯卡矩阵性能更好。
MATLAB课后习题
5、利用rand函数产生（0,1）间的均匀分布的10*10随机矩阵A，然后统计A中大于等于的元素的个数。
解：A=rand(10);B=A >= ;C=sum(B);count=sum(C)运行结果（每次运行结果是不同的，仅作参考）：count=326、利用randn函数产生均值为0，方差为1的10*10随机矩阵A，然后统计A中大于且小于的元素的个数。
%解：A=randn(10);B=(A<&(A>;C=sum(sum(B))运行结果（每次运行结果是不同的，仅作参考）：C=481、解：if and(a<1,b<=语句1;elseif and(a<1,b>`语句2;elseif and(a>=1,b<=语句3;else语句4;2、有一矩阵A，找出矩阵中值等于1的元素，并将它们重新排列成列向量B。
解：A=2*rand(4);k=find(A<=1);￥A(k)=[];%删除下标为k的元素B=A'运行结果（每次运行结果是不同的，仅作参考）B =)3、在一测量矩阵A（100*3）中，存在有奇异值（假设大于100的置认为是奇异值），编程实现删去奇异值所在的行。
解：A=120*randn(10,3);[i,j]=find(A>100);A(i,:)=[] %删去存在奇异值的行【4、在给定的100*100矩阵中，删去整行为0的行，删去整列为0的列。
解：A=diag([1 2 3 4],1)B=any(A)[i,j]=find(B==0)A(:,i)=[] %删除全为0的列B=any(A')[i,j]=find(B==0)A(j,:)=[] %删除全为0的行—运行结果：初始值：A =0 1 0 0 00 0 2 0 00 0 0 3 00 0 0 0 40 0 0 0 0操作后：A =1 0 0 0；0 2 0 00 0 3 00 0 0 41、将窗口分割成四格，分别绘制正弦、余弦、正切和余切函数曲线，并加上适当的标注。
完整word版,Matlab课后习题
习题 11. 执行下列指令，观察其运算结果, 理解其意义： (1) [1 2;3 4]+10-2i(2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4])(8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi)(10) [1 2;3 4]>=[4,3;2 1](11)find([10 20;30 40]>=[40,30;20 10])(12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示：a 为行号，b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2)2. 执行下列指令，观察其运算结果、变量类型和字节数，理解其意义： (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun)3. 本金K 以每年n 次，每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加，当增加到rK 时所花费的时间为)01.01ln(ln p n rT +=(单位：年)用MA TLAB 表达式写出该公式并用下列数据计算：r =2, p =0.5, n =12.4．已知函数f (x )=x 4-2x 在(-2, 2)内有两个根。
取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。
matlab课后习题及答案详解
matlab课后习题及答案详解第1章练习题1.安装matlab时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？在安装matlab时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即matlab选项）必须安装。
第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。
2.matlab操作方式桌面存有几个窗口？如何并使某个窗口瓦解桌面沦为单一制窗口？又如何将瓦解过来的窗口再次置放至桌面上？与其他计算机语言相比较，matlab语言注重的特点就是什么？matlab系统由那些部分共同组成？在matlab操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的close按钮，一个是可以使窗口成为独立窗口的undock按钮，点击undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择dock……菜单项就可以将独立的窗口重新防止的桌面上。
matlab具备功能强大、使用方便、输出简便、库函数多样、开放性弱等特点。
matlab系统主要由开发环境、matlab数学函数库、matlab语言、图形功能和应用程序接口五个部分组成。
3.如何设置当前目录和搜寻路径，在当前目录上的文件和在搜寻路径上的文件存有什么区别？命令历史窗口除了可以观测前面键入的命令外，除了什么用途？当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的setpath菜单项来完成。
在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被matlab运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。
命令历史窗口除了用作查阅以前键入的命令外，还可以轻易执行命令历史窗口中选取的内容、将选取的内容拷贝到剪贴板中、将选取内容轻易拷贝到m文件中。
matlab课后习题答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答：A(7)=[](3)将矩阵A的每个元素值加30。
答：A=A+30;(4)求矩阵A的大小和维数。
答：size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答：char(49);4.下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5.已知完成下列操作：(1)取出A的前3行构成矩阵B，前两列构成矩阵C，右下角32⨯子矩阵构成矩阵D，B与C 的乘积构成矩阵E。
matlab课后习题答案
习题二之阿布丰王创作2.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵,单个数据（标量）可以看成是仅含一个元素的矩阵,故向量和单个数据都可以作为矩阵的特例来处置.因此,矩阵是MATLAB最基本、最重要的数据对象.3.设A和B是两个同维同年夜小的矩阵,问：(1)A*B和A.*B的值是否相等？答：不相等.(2)A./B和B.\A的值是否相等？答：相等.(3)A/B和B\A的值是否相等？答：不相等.(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵,即A*inv(B),而B\A等效于B矩阵的逆左乘A矩阵,即inv(B)*A.4.写出完成下列把持的命令.(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B.答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素.答：A(7)=[](3)将矩阵A的每个元素值加30.答：A=A+30;(4)求矩阵A的年夜小和维数.答：size(A);ndims(A);(5)将向量t的0元素用机器零来取代.答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34 矩阵.答：reshape(x,3,4);(7)求一个字符串的ASCII码.答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符.答：char(49);5.下列命令执行后,L1、L2、L3、L4的值分别是几多？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]6.已知完成下列把持：(1)取出A的前3行构成矩阵B,前两列构成矩阵C,右下角32⨯子矩阵构成矩阵D,B与C的乘积构成矩阵E.答：B=A(1:3,:);C=A(:,1:2);D=A(2:4,3:4);E=B*C;(2)分别求E<D、E&D、E|D、~E|~D和find(A>=10&A<25).答：E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,E|D=111111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,~E|~D=001000⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦find(A>=10&A<25)=[1; 5].7.当A=[34, NaN, Inf, -Inf, -pi, eps, 0]时,分析下列函数的执行结果：all(A),any(A),isnan(A),isinf(A),isfinite(A).答：all(A)的值为0any(A)的值为1isnan(A)的值为[0, 1, 0, 0, 0, 0, 0]isinf(A)的值为[0, 0, 1, 1, 0, 0, 0]isfinite(A)的值为[1, 0, 0,0, 1, 1, 1]8.用结构体矩阵来存储5名学生的基本情况数据,每名学生的数据包括学号、姓名、专业和6门课程的成果.答：student(1).id='0001';student(1).name='Tom';student(1).major='computer';student(1).grade=[89,78,67,90,86,85];9.建立单位矩阵B并回答有关问题.B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3);B{2,2}={12,34,2;54,21,3;4,23,67};(1)size(B)和ndims(B)的值分别是几多？答：size(B)的值为2,2.ndims(B)的值为2.(2)B(2)和B(4)的值分别是几多？答：B(2)=147258369⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,B(4)=[12][34][2][54][21][3][4][23][67]⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(3)B(3)=[]和B{3}=[]执行后,B的值分别是几多？答：当执行B(3)=[]后,B={1, [1, 4, 7; 2, 5, 8; 3, 6, 9], {12, 34, 2; 54, 21, 3; 4, 23, 67}}当执行B{3}=[]后,B={1,[]; [1, 4, 7; 2, 5, 8; 3, 6, 9],{12, 34, 2; 54, 21, 3; 4, 23, 67}}习题三1.写出完成下列把持的命令.(1)建立3阶单位矩阵A.答：A=eye(3);(2)建立5×6随机矩阵A,其元素为[100,200]范围内的随机整数.答：round(100+(200-100)*rand(5,6));(3)发生均值为1,方差为0.2的500个正态分布的随机数.答：1+sqrt(0.2)*randn(5,100);(4)发生和A同样年夜小的幺矩阵.答：ones(size(A));(5)将矩阵A对角线的元素加30.答：A+eye(size(A))*30;(6)从矩阵A提取主对角线元素,并以这些元素构成对角阵B.答：B=diag(diag(A));2.使用函数,实现方阵左旋90o或右旋90o的功能.例如,原矩阵为A,A左旋后获得B,右旋后获得C.答：B=rot90(A);C=rot90(A,-1);3.建立一个方阵A,求A的逆矩阵和A的行列式的值,并验证A与A-1是互逆的.答：A=rand(3)*10;B=inv(A);C=det(A);先计算B*A,再计算A*B,由计算可知B*A=A*B,即A·A-1=A-1·A 是互逆.4.求下面线性方程组的解.答：A=[4,2,-1;3,-1,2;12,3,0];b=[2;10;8];x=inv(A)*b方程组的解为x=6.000026.666727.3333 -⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦5.求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹.(1)112351423052111509A-⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥⎣⎦(2)0.434328.9421B⎡⎤=⎢⎥-⎣⎦答：(1)取主对角线元素：diag(A);上三角阵：triu(A)；下三角阵：tril(A);秩：rank(A);范数：norm(A,1); 或 norm(A);或 norm(A,inf);条件数：cond(A,1); 或 cond(A,2); 或 cond(A,inf)迹：trace(A);(2)【请参考(1)】.6.求矩阵A的特征值和相应的特征向量.答：[V,D]=eig(A);习题四1.从键盘输入一个4位整数,按如下规则加密后输出.加密规则：每位数字都加上7,然后用和除以10的余数取代该数字；再把第一位与第三位交换,第二位与第四位交换.答：a=input('请输入4位整数：');A=[a/1000,a/100,a/10,a];A=fix(rem(A,10));A=rem(A+7,10);b=A(3)*1000+A(4)*100+A(1)*10+A(2);disp(['加密后的值为：',num2str(b)]);2.分别用if语句和switch语句实现以下计算,其中a、b、c的值从键盘输入.答：(1) 用if语句实现计算：a=input('请输入a的值：');b=input('请输入b的值：');c=input('请输入c的值：');x=input('请输入x的值：');if x>=0.5 & x<1.5y=a*x^2+b*x+c;endif x>=1.5 & x<3.5y=a*((sin(b))^c)+x;endif x>=3.5 & x<5.5y=log(abs(b+c/x));enddisp(['y=',num2str(y)]);(2) 用switch语句实现计算：a=input('请输入a的值：');b=input('请输入b的值：');c=input('请输入c的值：');x=input('请输入x的值：');switch fix(x/0.5)case {1,2}y=a*x^2+b*x+c;case num2cell(3:6)y=a*((sin(b))^c)+x;case num2cell(7:10)y=log(abs(b+c/x));enddisp(['y=',num2str(y)]);3.发生20个两位随机整数,输出其中小于平均值的偶数.答：A=fix(10+89*rand(1,20));sum=0;for i=1:20sum=sum+A(i);endB=A(find(A<(sum/20)));C=B(find(rem(B,2)==0));disp(C);4. 输入20个数,求其中最年夜数和最小数.要求分别用循环结构和调用MATLAB 的max 函数、min 函数来实现.答：(1) 用循环结构实现：v_max=0;v_min=0;for i=1:20x=input(['请输入第', num2str(i), '数：']);if x>v_maxv_max=x;end;if x<v_minv_min=x;end;enddisp(['最年夜数为：', num2str(v_max)]);disp(['最小数为：', num2str(v_min)]);(2) 用max 函数、min 函数实现：for i=1:5A(i)=input(['请输入第', num2str(i), '数：']);enddisp(['最年夜数为：', num2str(max(A))]);disp(['最小数为：', num2str(min(A))]);5. 已知：236312222s =+++++,分别用循环结构和调用MATLAB 的sum 函数求s 的值.答：s=0;for i=0:63s=s+2^i;ends(2) 调用sum函数实现：s=0:63;s=2.^s;sum(s)6.当n分别取100、1000、10000时,求下列各式的值.(1)111111(1)(ln2) 234nn+-+-++-+=(2)1111() 3574π-+-+=(3) 11111() 4166443n+++++=(4)224466(2)(2)133557(21)(21)2n nn nπ⎛⎫⨯⨯⨯⨯⎛⎫⎛⎫⎛⎫⎛⎫=⎪⎪⎪⎪ ⎪⨯⨯⨯-⨯+⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭要求分别用循环结构和向量运算（使用sum或prod函数）来实现.答：(1) 用循环结构实现：sum=0;for k=1:100sum=sum+(-1)^(k+1)/k;endsum使用sum函数：x=[];for k=1:10000x=[x, (-1)^(k+1)/k];endsum(x)sum=0;for k=1:100sum=sum+(-1)^(k+1)/(2*k-1);endsum使用sum函数：x=[];for k=1:100x=[x, (-1)^(k+1)/(2*k-1)];endsum(x)(3)用循环结构实现：sum=0;for k=1:100sum=sum+1/(4^k);endsum使用sum函数实现：x=[];for k=1:100x=[x, 1/(4^k)];endsum(x)(4)用循环结构实现：t=1;for k=1:100t=t*(((2*k)*(2*k))/((2*k-1)*(2*k+1)));endt使用prod函数实现：x=[];for k=1:100x=[x, ((2*k)*(2*k))/((2*k-1)*(2*k+1))];endprod(x)7. 编写一个函数文件,求小于任意自然数n的斐波那契（Fibnacci）数列各项.斐波那契数列界说如下：答：function x=fibnacci(n)for i=1:nif i<=2x(i)=1;elsex(i)=x(i-1)+x(i-2);endend8.编写一个函数文件,用于求两个矩阵的乘积和点乘,然后在命令文件中调用该函数.答：函数文件myfnc.m：function [x, y]=myfnc(A, B)tryx=A*B;catchx=[];endy=A.*B;命令文件myexe.m：A=input('请输入矩阵A：');B=input('请输入矩阵B：');[x, y]=myfnc(A, B);if length(x)==0display('两矩阵的维数不匹配,无法进行乘积运算！');elsedisp('矩阵A 和矩阵B 的乘积为：');xenddisp('矩阵A 和矩阵B 的点乘为：');y9. 先用函数的递归调用界说一个函数文件求1n m i i =∑,然后调用该函数文件求100501021111k k k k k k ===++∑∑∑. 答：函数文件myfnc.m ：function sum=myfnc(n, m)if n<=1sum=1;elsesum= myfnc (n-1, m)+n^m;end在命令窗口中调用myfnc.m 文件,计算100501021111k k k k k k ===++∑∑∑： sum=myfnc(100, 1)+myfnc(50, 2)+myfnc(10,-1) 10.写出下列法式的输出结果.① 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);endendends答：执行结果为s=108②命令文件exe.m执行后的结果为：x =41220y=。
matlab课后习题答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答：A(7)=[](3)将矩阵A的每个元素值加30。
答：A=A+30;(4)求矩阵A的大小和维数。
答：size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8) 求一个ASCII 码所对应的字符。
答：char(49);4. 下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5. 已知23100.7780414565532503269.5454 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦ 完成下列操作：(1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32⨯子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
matlab仿真课后习题
第一章习题3.请指出以下的变量名（函数名、M文件名）中，哪些是合法的Abc 2004x lil-1 wu_2004 a&b _xyz 解：合法的变量名有：Abc wu_20044．指令窗操作（1）求[12+2×（7-4）]÷32的运算结果解：>> [12+2*(7-4)]/3^2ans =2（2）输入矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。
解：>> A=[1,2,3;4,5,6;7,8,9]A =1 2 34 5 67 8 9（3）输入以下指令，观察运算结果；clear;x=-8::8;y=x';X=ones(size(y))*x;Y=y*ones(size(x));R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel('x'),ylabel('y'),zlabel('z')解：7．指令行编辑（1）依次键入以下字符并运行：y1=2*sin*pi)/(1+sqrt(5))解：>>y1=2*sin*pi)/(1+sqrt(5))y1 =(2)通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算；y2=2*cos*pi)/(1+sqrt(5))解：>>y2=2*cos*pi)/(1+sqrt(5))y2 =11.编写题4中（3）的M脚本文件，并运行之。
解：第二章习题1.在指令窗中键入x=1::2和y=2::1,观察所生成的数组。
解:>> x=1::2 x =>> y=2::1 y =Empty matrix: 1-by-02．要求在[0，2π]上产生50个等距采样数据的一维数组，试用两种不同的指令实现。
解： y1=0:2*pi/49:2*pi y2=linspace(0,2*pi,50)3.计算e -2t sint ，其中t 为[0，2π]上生成的10个等距采样的数组。
matlab课后习题答案
第2章 MATLAB 矩阵运算基础2.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ >> a=[5 7 3;4 9 1]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 142.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.0000i2.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 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
MATLAB课后习题集附标准答案
第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB 语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB 系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB 时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装. 第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可. 矚慫润厲钐瘗睞枥庑赖。
4、MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？聞創沟燴鐺險爱氇谴净。
答：在MATLAB 操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view 菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。
5、如何启动M 文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动.在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器.酽锕极額閉镇桧猪訣锥。
6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。
7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中. 謀荞抟箧飆鐸怼类蒋薔。
MATLAB课后习题资料
第一部分 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z e=+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
第二部分 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
5、利用rand函数产生（0,1）间的均匀分布的10*10随机矩阵A，然后统计A中大于等于0.6的元素的个数。
解：A=rand(10);B=A >= 0.6;C=sum(B);count=sum(C)运行结果（每次运行结果是不同的，仅作参考）：count=326、利用randn函数产生均值为0，方差为1的10*10随机矩阵A，然后统计A中大于-0.5且小于0.5的元素的个数。
解：A=randn(10);B=(A<0.5)&(A>-0.5);C=sum(sum(B))运行结果（每次运行结果是不同的，仅作参考）：C=481、解：if and(a<1,b<=0.5)语句1;elseif and(a<1,b>0.5)语句2;elseif and(a>=1,b<=0.5)语句3;else语句4;2、有一矩阵A，找出矩阵中值等于1的元素，并将它们重新排列成列向量B。
解：A=2*rand(4);k=find(A<=1);A(k)=[];%删除下标为k的元素B=A'运行结果（每次运行结果是不同的，仅作参考）B =1.47691.83481.53101.15241.36671.09321.28891.29521.3580 3、在一测量矩阵A（100*3）中，存在有奇异值（假设大于100的置认为是奇异值），编程实现删去奇异值所在的行。
解：A=120*randn(10,3);[i,j]=find(A>100);A(i,:)=[] %删去存在奇异值的行4、在给定的100*100矩阵中，删去整行为0的行，删去整列为0的列。
解：A=diag([1 2 3 4],1)B=any(A)[i,j]=find(B==0)A(:,i)=[] %删除全为0的列B=any(A')[i,j]=find(B==0)A(j,:)=[] %删除全为0的行运行结果：初始值：A =0 1 0 0 00 0 2 0 00 0 0 3 00 0 0 0 40 0 0 0 0操作后：A =1 0 0 00 2 0 00 0 3 00 0 0 41、将窗口分割成四格，分别绘制正弦、余弦、正切和余切函数曲线，并加上适当的标注。
程序为：x=0:pi/50:2*pi;k=[1 26 51 76 101];x(k)=[];%删除正切和余切的奇异点figure(1)subplot(2,2,1)plot(x,sin(x),'k--'),grid onlegend('\ity=sin(x)')title('y=sin(x)')xlabel('x'), ylabel('y')1subplot(2,2,2)plot(x,cos(x),'r--'),grid onlegend('\ity=cos(x)')title('y=con(x)')xlabel('x'), ylabel('y')subplot(2,2,3)plot(x,tan(x),'k'),grid onlegend('\ity=tan(x)')title('y=tan(x)')xlabel('x'), ylabel('y')subplot(2,2,4)plot(x,cot(x),'b-'),grid onlegend('\ity=cot(x)')title('y=cot(x)')xlabel('x'), ylabel('y')运行如下：2、绘制多峰函数peaks和三角函数多条曲线。
多峰函数peaks：[x,y]=meshgrid(-3:0.15:3);z=peaks(x,y);x1=x(1,:);figure(1)plot(x1,z),grid ontitle('二维多峰函数')图形为：[x,y]=meshgrid(-3:0.15:3);z=peaks(x,y);figure(1)plot3(x,y,z),grid ontitle('三维多峰函数')三角函数多条曲线：程序为：t=-pi:pi/20:pi;y1=sinh(t); %双曲正弦y2=cosh(t); %双曲余弦figure(1)subplot(2,1,1)plot(t,y1,'r--',t,y2,'k-'),grid onlegend('\ity1=sinh(t)','\ity2=cosh(t)')title('三角函数1')xlabel('t'), ylabel('y')subplot(2,1,2) plot(t,sin(t),'k-'),grid onhold on %保持原有图像函数plot(t,cos(t),'r--')legend('\ity2=cos(t)','\ity1=sin(t)')title('三角函数2')xlabel('t'), ylabel('y')运行图形为：3、将图形窗口分成两个，分别绘制以下函数在[-3,3]区间上的曲线，并利用axis调整轴刻度，使他们具有相同缩放尺度。
y1=2x+5；y2=x2-3x+1。
程序为：x=-3:0.1:3;y1=2*x+5;y2=x.^2-3*x+1;figure(1)subplot(2,2,1)plot(x,y1,'r-'),grid onlegend('\ity1=2*x+5')title('y1=2x+5')xlabel('x'), ylabel('y1')subplot(2,2,2)plot(x,y2,'k-'),grid onlegend('\ity2=x.^2-3*x+1')title('y2=x^2-3x+1')xlabel('x'), ylabel('y2')subplot(2,2,3)plot(x,y1,'r-'),grid onlegend('\ity1=2*x+5')title('调整后的y1=2x+5')axis([-3 3 -10 10])xlabel('x'), ylabel('y1')subplot(2,2,4)plot(x,y2,'k-'),grid onlegend('\ity2=x.^2-3*x+1')title('调整后的y2=x^2-3x+1')axis([-3 3 -10 10]) %调整坐标轴xlabel('x'), ylabel('y2')运行后的图形：4、绘制饼图。
程序为：2x=[190 33 45 42 45];explode=[0 1 0 0 0];figure(1)subplot(2,1,1)colormap hsvpie(x,explode)gtext('生活费')gtext('资料费')gtext('电话费')gtext('衣服')gtext('其它')title('二维饼图')subplot(2,1,2)colormap hsvpie3(x,explode)title('三维饼图')图形为：5、画出函数z=（x-2）2+（y-1.2）2+sin(xy)的三维曲线和网格曲线。
程序为：[x,y]=meshgrid(0:0.5:10); %为三维绘图产生x，y数据矩阵z=(x-2).^2+(y-1.2).^2;figure(1)subplot(2,1,1)mesh(x,y,z),grid on %绘制网格曲线title('网格曲线')subplot(2,1,2)plot3(x,y,z),grid ontitle('三维曲线')6、画出下列函数的曲面及等高线图z=x2+y2+sin(xy)。
程序为：[x,y]=meshgrid(0:pi/10:2*pi);z=x.^2+y.^2+sin(x*y);figure(1)subplot(2,1,1)surfc(x,y,z), grid ontitle('曲面和等高线')subplot(2,1,2)[c,h]=contour(x,y,z);set(h,'showtext','on','textstep',get(h,'levelstep')*2);title('等高线')1、将图形窗口分成两个，分别绘制正割和余割曲线，并加上标注。
程序为：x1=0:pi\10:2*pi;figure(1)subplot(2,1,1)plot(x,sec(x),'k-'),grid onlegend('\ity=sec(x)')title('y=sec(x)')xlabel('x'), ylabel('y')subplot(2,1,2)plot(x,csc(x),'k-'),grid onlegend('\ity=csc(x)')title('y=csc(x)')xlabel('x'), ylabel('y')运行后图形为：2、画出对数和指数曲线并加上标注。
x=0.01:0.1:10;y1=log10(x);y2=exp(x);figure(1)subplot(2,1,1)plot(x,y1,'k-'),grid onlegend('\ity1=log-{10}(x)')title('y1=log-{10}(x)')xlabel('x'), ylabel('y1')subplot(2,1,2)plot(x,y2,'k-'),grid onlegend('\ity2=exp(x)')title('y2=exp(x)')xlabel('x'), ylabel('y2')3、设有函数y=exp(x+5)+x.^3，在半对数坐标系中绘制曲线。