matlab7.0x课后习题答案
Matlab课后习题解答
0.80
0.95
电阻y
15
18
19
21
22.6
23.8
26
>> x=[0.1,0.3,0.4,0.55,0.7,0.8,0.95];
y=[15,18,19,21,22.6,23.8,26];
p1=polyfit(x,y,1);
p3=polyfit(x,y,3);
p5=polyfit(x,y,5);
if x>=90
disp('优秀');
elseif x>=80
disp('良好');
elseif x>=60
disp('及格');
else
disp('不及格');
end
>> x=85
x =
85
良好
Q3：编写函数，计算
>> sum=0;
>> for i=1:50
a=1;
for j=1:i
a=a*j;
end
ans =
0
(4)
>> syms n
>> limit(sqrt(n+2)-2*sqrt(n+1)+sqrt(n),n,inf)
ans =
0
Q2:用MATLAB软件求下列函数极限：
(1)
>> syms x
>> limit((((1+x)^(1/3)-1)/x),x,0)
ans =
1/3
(4)
>> syms x
(2)
>> syms x
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文件中.8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息.（3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数.注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法？各有什么优点？ 答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷；（2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； （3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； （4）通过数据文件建立，该方法可以调用由其他软件产生数据.3.3在进行算术运算时，数组运算和矩阵运算各有什么要求？答：进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。
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习题及参考答案经典.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 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574（2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []2965318772546， 将(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 的最大值、最小值、均值和标准差。
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=1 ezplot(y^2-(x^3)/(a-x))(5)摆线:()()t b y t t a x cos 1,sin -=-=命令：>> 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)阿基米德螺线>> 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)心形线>> theta=0:0.1:2*pi;>> rho1=1+cos(theta);>> subplot(1,2,1),polar(theta,rho1)练习2.21. 求出下列极限值（1）nn n n 3lim 3+∞→命令：>>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 =（3）x x x 2cot lim 0→命令：syms x ；>> limit(x*cot(2*x),x,0)ans =1/2（4）)(cos lim c m xx ∞→ 命令：syms x m ；limit((cos(m/x))^x,x,inf)ans =1（5）)111(lim 1--→e x x x 命令：syms x>> limit(1/x-1/(exp(x)-1),x,1)ans =(exp(1)-2)/(exp(1)-1)（6）)(2lim x x x x -+∞→命令：syms x>> limit((x^2+x)^(1/2)-x,x,inf)ans =1/2练习2.41. 求下列不定积分，并用diff 验证：（1）⎰+xdx cos 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)⎰+e x dx 1clearsyms 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^2clearsyms 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 ⎰sec 3syms x yy=sec(x)^3;f=int(y,x)f =1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x))clearsyms 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 x ⎰-10clearsyms xy=int(x^(-x),x,0,1)y =int(x^(-x),x = 0 .. 1)vpa(y,10)ans =1.2912859972）xdx e x cos 3202⎰πclearsyms xy=int(exp(2*x)*cos(x)^3,x, clearsyms 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 x e 210221-⎰π>> 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 x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),-inf,inf);>> vpa(y,14)ans =.99999999999999练习2.51判断下列级数的收敛性，若收敛，求出其收敛值。
matlab7.X课后习题
第二章第三题clc;clear;close alla=5-round(100*rand(50,2))/10 a =-4.5000 3.10002.7000 -1.8000-1.1000 2.00000.1000 -0.4000-3.9000 3.5000-2.6000 -2.00000.4000 1.20004.8000 -3.6000-3.2000 -3.50000.6000 -0.9000-1.2000 0-2.9000 -4.0000-4.2000 -3.2000-2.4000 -1.40003.2000 -3.20000.9000 -1.6000-4.4000 1.6000-4.2000 2.10000.9000 1.6000-3.9000 -0.30004.4000 -2.30001.5000 1.9000-3.1000 -3.40004.9000 -0.70003.6000 1.30003.0000 -2.00003.0000 -0.5000-1.0000 0.60002.3000 -1.90003.0000 -1.20004.8000 -2.9000-2.5000 -4.60000.5000 -0.2000-4.3000 -3.80000.3000 3.30000.8000 -4.8000-3.5000 2.3000-0.3000 2.50003.0000 -3.8000-1.7000 -2.4000-3.4000 3.60004.8000 4.9000-1.8000 -3.90001.2000 3.0000-3.3000 2.00000 -1.6000-2.1000 2.20000.7000 0.30002.0000 4.40003.1000 -4.9000第六题clc;clear;close alla=randn(10,10)b=find(-0.5<a<0.5);z=length(b)a =Columns 1 through 9-0.4326 -0.1867 0.2944 -0.3999 -1.6041 -1.0106 0.0000 0.5689 0.6232 -1.6656 0.7258 -1.3362 0.6900 0.2573 0.6145 -0.3179 -0.2556 0.7990 0.1253 -0.5883 0.7143 0.8156 -1.0565 0.5077 1.0950 -0.3775 0.9409 0.2877 2.1832 1.6236 0.7119 1.4151 1.6924 -1.8740 -0.2959 -0.9921 -1.1465 -0.1364 -0.6918 1.2902 -0.8051 0.5913 0.4282 -1.4751 0.2120 1.1909 0.1139 0.8580 0.6686 0.5287 -0.6436 0.8956 -0.2340 0.2379 1.1892 1.0668 1.2540 1.1908 0.2193 0.3803 0.7310 0.1184 -1.0078 -0.0376 0.0593 -1.5937 -1.2025 -0.9219 -1.0091 0.5779 0.3148 -0.7420 0.3273 -0.0956 -1.4410 -0.0198 -2.1707 -0.0195 0.0403 1.4435 1.0823 0.1746 -0.8323 0.5711 -0.1567 -0.0592 -0.0482 0.6771 -0.3510 -0.1315Column 100.38990.0880-0.6355-0.55960.4437-0.94990.78120.5690-0.8217-0.2656z =26第三章第五题clc;clear;close all x=-3:0.1:3;y1=2*x+5;y2=x.^2-3*x+1;figure(1)subplot(2,1,1);plot(x,y1),grid on axis([-3 3 -5 15 ]) subplot(2,1,2);plot(x,y2),grid on axis([-3 3 -5 15 ])第七题clc;clear;close all x=[190 33 45 42 45]; explode=[0 1 0 0 0]; figure(1)pie3(x,explode)title('三维饼图')第八题clc,clear,close all[x,y]=meshgrid(-3:.125:3);z=(x-2).^2+(y-1.2).^2;meshc(x,y,z);axis([-4 4 -4 4 0 30])title('z=(x-2).^2+(y-1.2).^2的网格曲线')第九题clc;clear;close all[x,y]=meshgrid(-5:.1:5);z=x.^2+y.^2+sin(x*y);figure(1)subplot(2,1,1)plot3(x,y,z),grid onsubplot(2,1,2)[C,h]=contour(x,y,z);set(h,'showtext','on','textstep',get(h,'levelstep')*2)第四章第三题clc;clear;close all i=0;for y=[1:1:2006]; while y>=10if rem(y,10)==0 i=i+1;endy=fix(y/10); endenddisp(i)504第五题clc;clear;close all i=1x1=[-3:0.01:3];for x=x1if x>=-3&x<-1y1(i)=(-x^2-4*x-3)/2;elseif x>=-1&x<1y1(i)=-x^2+1;elseif x>=1&x<=3y1(i)=(-x^2+4*x-3)/2;endi=i+1;endfigure(1)plot(x1,y1),grid on第八题clc;clear;close allnumber=input('input the profit\n');if number>40bonus=(number-40)*0.01+20*0.02+10*0.05+10*0.1; elseif number<=40&number>20bonus=(number-20)*0.02+10*0.05+10*0.1;elseif number<=20&number>10bonus=(number-10)*0.05+10*0.1;elsebonus=number*0.1;endData=sprintf('the bonus is %d',bonus); disp(Data);input the profit150000the bonus is 1.501500e+003第九题clc;clear;close allm=1;n=2;for i=1:1:15x(i)=n/m;k=n;n=m+n;m=k;endsum(x)ans =24.5701第五章第一题clc;clear;close alla1=[1 9 8;7 2 5;3 -2 7];b1=inv(a1)det(a1)a2=[1 0 -7 5;0 -26 7 2;7 4 3 5;8 -3 2 15];b2=inv(a2)det(a2)b1 =-0.0543 0.1787 -0.06560.0769 0.0385 -0.11540.0452 -0.0656 0.1380ans =-442b2 =0.1501 0.0654 0.3046 -0.1603-0.0509 -0.0471 -0.0312 0.0336-0.1697 -0.0204 -0.0703 0.0827-0.0676 -0.0416 -0.1593 0.1478ans =-9334第二题clc;clear;close alla1=[1 2 3;2 2 5;3 5 1];b1=[11 12 31]';x1=a1\b1a2=[3 1 0 5;0 6 7 3;0 4 3 0;2 -1 2 6]; b2=[2 4 7 8]';x2=a2\b2x1 =1.80005.2000-0.4000x2 =28.1892-9.459514.9459-14.6216第八题clc;clear;close allp1=[1 -2 -3 4 2];r1=roots(p1)p2=[1 -7 5 31 -30];r2=roots(p2)p3=[1 -1 -25 25];r3=roots(p3)p4=[-2 3 1 5 8 0];r4=roots(p4)r1 =2.41421.4142-1.4142-0.4142r2 =5.0000-2.00003.00001.0000r3 =-5.00005.00001.0000r4 =2.41760.0148 + 1.3215i0.0148 - 1.3215i-0.9473第九题clc;clear;close allp1=[1 -2 -3 4 2];y1=polyval(p1,[-1.5 2.1 3.5]) p2=[1 -7 5 31 -30];y2=polyval(p2,[-1.5 2.1 3.5]) p3=[1 -1 -25 25];y3=polyval(p3,[-1.5 2.1 3.5]) p4=[-2 3 1 5 8 0];y4=polyval(p4,[-1.5 2.1 3.5]) y1 =1.0625 -1.9039 43.5625y2 =-36.5625 11.7711 -10.3125y3 =56.8750 -22.6490 -31.8750y4 =26.2500 24.7733 -468.1250第十三题clc;clear;close alla1=rand(1,50);b1=max(a1)c1=min(a1)m1=mean(a1)s1=std(a1)^2a2=randn(1,50);b2=max(a2)c2=min(a2)m2=mean(a2)s2=std(a2)^2b1 =0.9901c1 =0.0150m1 =0.5246s1 =0.0646b2 =1.4885c2 =-2.2023m2 =-0.2348s2 =0.8153第十六题clc;clear;close allt=[-2:.1:2];x1=sin(50*t)+randn(size(t)); x2=cos(50*t)+randn(size(t)); x3=sin(50*t)+randn(size(t)); r1=corrcoef(x1,x2)r2=corrcoef(x1,x3)X1=randn(size(t));X2=randn(size(t));X3=randn(size(t));R1=corrcoef(X1,X2)R2=corrcoef(X1,X3)r1 =1.0000 -0.0274-0.0274 1.0000r2 =1.0000 0.33740.3374 1.0000R1 =1.0000 -0.0346-0.0346 1.0000R2 =1.0000 0.11180.1118 1.0000第六章第二题clc;clear;close allr=round(100*rand(7,5,10));r1=reshape(r(:,:,1),1,35);a1=mean(r1)r2=reshape(r(:,:,2),1,35);a2=mean(r2)r3=reshape(r(:,:,3),1,35);a3=mean(r3)r4=reshape(r(:,:,4),1,35);a4=mean(r4)r5=reshape(r(:,:,5),1,35);a5=mean(r5)r6=reshape(r(:,:,6),1,35);a6=mean(r6)r7=reshape(r(:,:,7),1,35);a7=mean(r7)r8=reshape(r(:,:,8),1,35);a8=mean(r8)r9=reshape(r(:,:,9),1,35);a9=mean(r9)r10=reshape(r(:,:,10),1,35);a10=mean(r10)A=[a1 a2 a3 a4 a5 a6 a7 a8 a9 a10] [B,A]=sort(A)a1 =54.1714 a2 =48a3 =52.4000 a4 =44.8286 a5 =53.2286 a6 =49.2571 a7 =50.7429 a8 =46.0571 a9 =46.1714 a10 =50.8286A =54.1714 48.0000 52.4000 44.8286 53.2286 49.2571 50.7429 46.0571 46.1714 50.8286B =44.8286 46.0571 46.1714 48.0000 49.2571 50.7429 50.8286 52.4000 53.2286 54.1714A =4 8 9 2 6 7 10 35 1第七章第二题clc;clear;close allx=[-1:.01:1];y=-2*x.^2+4;y1=-2*x.^2+4+sin(2*pi*x);figure(1)plot(x,y,x,y1),grid ontext(-0.3,4,'最大值\rightarrow')text(-1.0,4,'极大值\rightarrow')text(-0.6,2.8,'极小值\rightarrow')text(-0.1,4.8,'极大值\rightarrow')text(0.5,1.7,'极小值\rightarrow')第四题function word(y)i=1;while(any(y))[A(i).Name,y]=strtok(y);A(i).no=i;A(i).length=length(A(i).Name);A(i).value=sum([double(A(i).Name)]'); disp(A(i));i=i+1;endclc;clear;close all; y='I love MATLAB'; word(y);Name: 'I'no: 1length: 1value: 73Name: 'love' no: 2length: 4value: 438Name: 'MATLAB' no: 3length: 6value: 433。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/a428921bf524ccbff12184eb.png)
实验一 MATLAB 运算基础1.先求下列表达式的值，然后显示 MATLAB 工作空间的使用情 况并保存全部变量解:4.完成下列操作：(1) 求[100,999] 之间能被21整除的数的个数。
(1) z i 2sin 85° 1 e 2Z 2 卯(x L)，其中x2 0.45 2i 5Z 3 0.3a 0.3a e e 2sin (a 0.3)3.0, 2.9,L ,2.9, 3.0 Z 4t 2t 2 t 2 1 2t 其中 t=0:0.5:2.5(2) 建立一个字符串向量，删除其中的大写字母。
解: (1)结果：m=100:999;n=fin d(mod(m,21)==0);len gth( n)ans =43(2).建立一个字符串向量例如:ch二'ABC123d4e56Fg9：则要求结果是: ch二'ABC123d4e56Fg9：k=fi nd(ch>='A'&ch<='Z'); ch(k)=[] ch =123d4e56g9 实验二MATLAB 矩阵分析与处理1.设有分块矩阵A E 3 3 °2 3 R 3 2S 2 2 ，其中E 、R 、0、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证A 2 E R RS 0 S 22 3 4 1 1 1 3 4 5 1 1 1 4 5 6 x , 0.95 x 2 0.67 x 3 0.52(1) 求方程的解。
(2) 将方程右边向量元素 b 3改为0.53再求解，并比较 b 3的变化和解的相对变化。
(3)计算系数矩阵A 的条件数并分析结论。
解：M 文件如下： 解：M 文件如下; 5.下面是一个线性方程组:实验三选择结构程序设计1. 求分段函数的值。
x2x 6 x 0且x 3y x2 5x 6 0 x 5且x 2 及x 3x2x 1 其他用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0 时的y值。
matlab习题及答案
matlab习题及答案Matlab习题及答案Matlab是一种强大的数学计算软件，被广泛应用于科学计算、数据分析和工程设计等领域。
在学习和使用Matlab的过程中，习题是一种非常有效的学习方式。
本文将给出一些常见的Matlab习题及其答案，帮助读者更好地掌握Matlab的使用技巧。
一、基础习题1. 计算1到100之间所有奇数的和。
解答：```matlabsum = 0;for i = 1:2:100sum = sum + i;enddisp(sum);```2. 编写一个函数，计算任意两个数的最大公约数。
解答：```matlabfunction gcd = computeGCD(a, b)while b ~= 0temp = b;a = temp;endgcd = a;end```3. 编写一个程序，生成一个5×5的随机矩阵，并计算矩阵的行和列的平均值。
解答：```matlabmatrix = rand(5);row_average = mean(matrix, 2);col_average = mean(matrix);disp(row_average);disp(col_average);```二、进阶习题1. 编写一个程序，实现插入排序算法。
解答：```matlabfunction sorted_array = insertionSort(array)n = length(array);for i = 2:nj = i - 1;while j > 0 && array(j) > keyarray(j+1) = array(j);j = j - 1;endarray(j+1) = key;endsorted_array = array;end```2. 编写一个程序，实现矩阵的转置。
解答：```matlabfunction transposed_matrix = transposeMatrix(matrix) [m, n] = size(matrix);transposed_matrix = zeros(n, m);for i = 1:mfor j = 1:ntransposed_matrix(j, i) = matrix(i, j);endendend```3. 编写一个程序，实现二分查找算法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1、利用基本矩阵产生3*3和15*8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机阵（[-1,1]之间）、正态分布随机阵（均值为1，方差为4）。
解：A1=eye(3); A2=ones(3); A3=zeros(3); A4=2*rand(3)-1; A5=2*randn(3)+1;B1=eye(15,8); B2=ones(15,8); B3=zeros(15,8); B4=2*rand(15,8)-1; B5=2*randn(15,8)+1; 结果：由于数据是随机产生的，所以在没有给出运行结果。
2、利用diag等函数产生下列矩阵：a=[0 0 8;0 -7 5;2 3 0] b=[2 0 4;0 5 0;7 0 8]然后利用reshape函数将它们变换成行向量。
解：产生a的程序：b=diag([8 -7 2]);c=b+diag([5 3],-1);a=fliplr(c)产生b的程序：s=[2 2 8];t=[4 3 7];v=diag(s);p=diag(t)+fliplr(v);b=fliplr(p)运行结果：a =0 0 80 -7 52 3 0b =2 0 40 5 07 0 8利用reshape函数将它们变换成行向量：reshape(a,1,9)ans =0 0 2 0 -7 3 8 5 03、产生一均匀分布在（-5,5）之间的随机阵（50*2），要求精确到小数点后一位。
解：A=5-round(100*rand(50,2))/10部分数据结果：A =2.4000 4.2000-0.1000 2.7000-4.6000 -3.3000-0.5000 -0.40003.50004.20004、编程实现当t∈[-π,π],间隔为1°时求解正弦和余弦值。
解：t=(-1*pi:1/180:pi);y1=sin(t)y2=cos(t)部分数据结果：Columns 10 through 18（y1）-0.0500 -0.0555 -0.0611 -0.0666 -0.0722 -0.0777 -0.0832-0.0888 -0.09435、利用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.35803、在一测量矩阵A（100*3）中，存在有奇异值（假设大于100的置认为是奇异值），编程实现删去奇异值所在的行。
解：A=120*randn(10,3);[i,j]=find(A>100);A(i,:)=[] %删去存在奇异值的行运行结果（每次运行结果是不同的，仅作参考）：A =49.5355 -23.7550 -73.0269-118.4354 39.3214 -88.447291.1482 -28.5962 -209.985521.2336 -74.0239 -9.5871-15.8184 72.1322 22.044471.4429 11.0770 34.89484、在给定的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')subplot(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、绘制饼图。
程序为：x=[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)的三维曲线和网格曲线。